// Copyright (c) 2006-2008 Max-Planck-Institute Saarbruecken (Germany). // All rights reserved. // // This file is part of CGAL (www.cgal.org); you can redistribute it and/or // modify it under the terms of the GNU Lesser General Public License as // published by the Free Software Foundation; either version 3 of the License, // or (at your option) any later version. // // Licensees holding a valid commercial license may use this file in // accordance with the commercial license agreement provided with the software. // // This file is provided AS IS with NO WARRANTY OF ANY KIND, INCLUDING THE // WARRANTY OF DESIGN, MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE. // // $URL$ // $Id$ // // // 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 // 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 std::unary_function { BFI operator()(const Input& pair){ return Sqrt()(Cast()(pair.second)); } }; public: typedef Cache Cache_type; static Cache_type cache; }; template typename Sqrt_extension_bfi_cache::Cache_type Sqrt_extension_bfi_cache::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