// Copyright (c) 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 #ifndef CGAL_SQRT_EXTENSION_CHINESE_REMAINDER_TRAITS #define CGAL_SQRT_EXTENSION_CHINESE_REMAINDER_TRAITS #include #include #include namespace CGAL { template class Chinese_remainder_traits >{ private: typedef CGAL::Sqrt_extension EXT; typedef Chinese_remainder_traits CRT_NT; typedef Chinese_remainder_traits CRT_ROOT; public : typedef EXT Type; typedef typename CRT_NT::Scalar_type Scalar_type; struct Chinese_remainder{ void operator() ( const Scalar_type& m1, const Scalar_type& m2, const Scalar_type& m, const Scalar_type& s, const Scalar_type& t, const Type& u1, const Type& u2, Type& u) const { NT a0,a1; ROOT root; typename CRT_NT::Chinese_remainder chinese_remainder_nt; chinese_remainder_nt(m1,m2,m,s,t,u1.a0(),u2.a0(),a0); CGAL_postcondition_msg(!((u1.is_extended() && !u2.is_extended()) || (!u1.is_extended() && u2.is_extended())), " chinese remainder not possible for these numbers"); if(u1.is_extended() || u2.is_extended()){ chinese_remainder_nt(m1,m2,m,s,t,u1.a1(),u2.a1(),a1); typename CRT_ROOT::Chinese_remainder chinese_remainder_root; chinese_remainder_root(m1,m2,m,s,t,u1.root(),u2.root(),root); u=Type(a0,a1,root); }else{ u=Type(a0); } } void operator() ( const Scalar_type& m1, const Type& u1, const Scalar_type& m2, const Type& u2, Scalar_type& m, Type& u) const { Scalar_type s,t; CGAL::extended_euclidean_algorithm(m1,m2,s,t); m = m1 * m2; this->operator()(m1,m2,m,s,t,u1,u2,u); } }; }; } // namespace CGAL #endif // CGAL_SQRT_EXTENSION_CHINESE_REMAINDER_TRAITS