// Copyright (c) 2006-2008 Max-Planck-Institute Saarbruecken (Germany). // All rights reserved. // // This file is part of CGAL (www.cgal.org) // // $URL: https://github.com/CGAL/cgal/blob/v5.2/Number_types/include/CGAL/Sqrt_extension/convert_to_bfi.h $ // $Id: convert_to_bfi.h 0779373 2020-03-26T13:31:46+01:00 Sébastien Loriot // SPDX-License-Identifier: LGPL-3.0-or-later OR LicenseRef-Commercial // // // Author(s) : Michael Hemmer // This files adds an optional static cache to convert_to_bfi for Sqrt_extension #ifndef CGAL_SQRT_EXTENSION_CONVERT_TO_BFI_H #define CGAL_SQRT_EXTENSION_CONVERT_TO_BFI_H #include #include #include #include #include #include // Disbale SQRT_EXTENSION_TO_BFI_CACHE by default #ifndef CGAL_USE_SQRT_EXTENSION_TO_BFI_CACHE #define CGAL_USE_SQRT_EXTENSION_TO_BFI_CACHE 0 #endif #if CGAL_USE_SQRT_EXTENSION_TO_BFI_CACHE namespace CGAL { namespace INTERN_SQRT_EXTENSION { template class Sqrt_extension_bfi_cache { typedef std::pair Input; typedef BFI Output; typedef typename Coercion_traits::Cast Cast; typedef typename Algebraic_structure_traits::Sqrt Sqrt; struct Creator : public CGAL::cpp98::unary_function { BFI operator()(const Input& pair){ return Sqrt()(Cast()(pair.second)); } }; public: typedef Cache Cache_type; static Cache_type& cache() { CGAL_STATIC_THREAD_LOCAL_VARIABLE_0(Cache_type, cache_); return cache_; } }; } // namespace INTERN_SQRT_EXTENSION template typename Get_arithmetic_kernel::Arithmetic_kernel::Bigfloat_interval convert_to_bfi(const CGAL::Sqrt_extension& x) { typedef typename Get_arithmetic_kernel::Arithmetic_kernel AK; typedef typename AK::Bigfloat_interval BFI; typedef Bigfloat_interval_traits BFIT; long precision = typename BFIT::Get_precision()(); BFI result; if(x.is_extended()){ typedef INTERN_SQRT_EXTENSION::Sqrt_extension_bfi_cache Get_cache; BFI a0(convert_to_bfi(x.a0())); BFI a1(convert_to_bfi(x.a1())); BFI root(Get_cache::cache()(std::make_pair(precision,x.root()))); result = a0+a1*root; }else{ result = convert_to_bfi(x.a0()); } #ifndef NDEBUG BFI result_; typedef typename Get_arithmetic_kernel::Arithmetic_kernel AT; typedef typename AT::Bigfloat_interval BFI; if(x.is_extended()){ BFI a0(convert_to_bfi(x.a0())); BFI a1(convert_to_bfi(x.a1())); BFI root(CGAL::sqrt(convert_to_bfi(x.root()))); result_ = a0+a1*root; }else{ result_ = convert_to_bfi(x.a0()); } CGAL_assertion(lower(result) == lower(result_)); CGAL_assertion(upper(result) == upper(result_)); #endif return result; } } // namespace CGAL #endif #endif // CGAL_SQRT_EXTENSION_CONVERT_TO_BFI_H