int main(){
// Enforce IEEE double precision and to nearest before using modular arithmetic
CGAL::Protect_FPU_rounding<true> pfr(CGAL_FE_TONEAREST);
#ifdef CGAL_USE_LEDA
test_modular_gcd_utils<CGAL::LEDA_arithmetic_kernel>();
#endif // CGAL_USE_LEDA
#ifdef CGAL_USE_CORE
test_modular_gcd_utils<CGAL::CORE_arithmetic_kernel>();
#endif // Lis_HAVE_CORE
return 0;
}
int main()
{
// Enforce IEEE double precision and rounding mode to nearest
CGAL::Protect_FPU_rounding<true> pfr(CGAL_FE_TONEAREST);
test_modular_traits<int>();
#ifdef CGAL_USE_LEDA
test_modular_traits<leda::integer>();
test_modular_traits<CGAL::Polynomial< leda::integer > >();
test_modular_traits<CGAL::Lazy_exact_nt< leda::integer > >();
test_modular_traits<CGAL::Sqrt_extension< leda::integer , leda::integer > >();
#endif
#ifdef CGAL_USE_CORE
test_modular_traits<CORE::BigInt>();
test_modular_traits<CGAL::Polynomial< CORE::BigInt > >();
test_modular_traits<CGAL::Lazy_exact_nt< CORE::BigInt > >();
test_modular_traits<CGAL::Sqrt_extension< CORE::BigInt , CORE::BigInt > >();
#endif
#ifdef CGAL_USE_GMP
test_modular_traits<CGAL::Gmpz>();
#endif
#ifdef CGAL_USE_GMPXX
test_modular_traits< mpz_class >();
#endif
// test Sqrt_extension
test_modular_traits<CGAL::Sqrt_extension< int , int > >();
assert(
(!CGAL::Modular_traits<CGAL::Sqrt_extension<double,double> >
::Is_modularizable::value));
// test Polynomial
test_modular_traits<CGAL::Polynomial< int > >();
assert(
!CGAL::Modular_traits<CGAL::Polynomial<double> >
::Is_modularizable::value);
// test_modular_traits<CGAL::MP_Float >();
test_modular_traits< CGAL::Lazy_exact_nt<int> >();
assert(
!CGAL::Modular_traits<CGAL::Lazy_exact_nt< double > >
::Is_modularizable::value);
}
int main(){
// Set wrong rounding mode to test modular arithmetic
CGAL::Protect_FPU_rounding<true> pfr(CGAL_FE_UPWARD);
#ifdef CGAL_HAS_DEFAULT_ARITHMETIC_KERNEL
typedef CGAL::Arithmetic_kernel AK;
test_AK_1<AK>();
test_AK_2<AK>();
test_AK_4<AK>();
test_AK_5<AK>();
test_AK_6<AK>();
#else
std::cout << "No default arithmetic kernel has been found.\nNothing was tested" << std::endl;
#endif // CGAL_HAS_DEFAULT_ARITHMETIC_KERNEL
return 0;
}
virtual void on_draw()
{
pixfmt pf(rbuf_window());
typedef agg::pixfmt_alpha_blend_gray<gray_blender, agg::rendering_buffer, 3, 2> pixfmt_r;
typedef agg::pixfmt_alpha_blend_gray<gray_blender, agg::rendering_buffer, 3, 1> pixfmt_g;
typedef agg::pixfmt_alpha_blend_gray<gray_blender, agg::rendering_buffer, 3, 0> pixfmt_b;
pixfmt_r pfr(rbuf_window());
pixfmt_g pfg(rbuf_window());
pixfmt_b pfb(rbuf_window());
agg::renderer_base<pixfmt> rbase(pf);
agg::renderer_base<pixfmt_r> rbr(pfr);
agg::renderer_base<pixfmt_g> rbg(pfg);
agg::renderer_base<pixfmt_b> rbb(pfb);
agg::rasterizer_scanline_aa<> ras;
agg::scanline_p8 sl;
rbase.clear(agg::rgba(1,1,1));
agg::ellipse er(width() / 2 - 0.87*50, height() / 2 - 0.5*50, 100, 100, 100);
ras.add_path(er);
agg::render_scanlines_aa_solid(ras, sl, rbr,
agg::gray8(0, unsigned(m_alpha.value())));
agg::ellipse eg(width() / 2 + 0.87*50, height() / 2 - 0.5*50, 100, 100, 100);
ras.add_path(eg);
agg::render_scanlines_aa_solid(ras, sl, rbg,
agg::gray8(0, unsigned(m_alpha.value())));
agg::ellipse eb(width() / 2, height() / 2 + 50, 100, 100, 100);
ras.add_path(eb);
agg::render_scanlines_aa_solid(ras, sl, rbb,
agg::gray8(0, unsigned(m_alpha.value())));
agg::render_ctrl(ras, sl, rbase, m_alpha);
}
// This file needs Open MP.
// Use ,e.g. , gcc-4.3.1 with -fopenmp and -lgomp
int main () {
int tid;
//Beginning of parallel section. Fork a team of threads.
//Specify variable scoping
#pragma omp parallel private(tid)
{
// Enforce IEEE double precision and rounding mode to nearest
CGAL::Protect_FPU_rounding<true> pfr(CGAL_FE_TONEAREST);
tid = omp_get_thread_num();
int old_prime = CGAL::internal::primes[0];
int new_prime = CGAL::internal::primes[tid+1];
assert(CGAL::Residue::get_current_prime() == old_prime);
CGAL::Residue::set_current_prime(new_prime);
assert(CGAL::Residue::get_current_prime() == new_prime);
typedef CGAL::Residue NT;
typedef CGAL::Field_tag Tag;
typedef CGAL::Tag_true Is_exact;
CGAL::test_algebraic_structure<NT,Tag, Is_exact>();
}
}
void test_AK_6(){
CGAL::set_pretty_mode(std::cout);
CGAL::set_pretty_mode(std::cerr);
typedef typename AK::Integer Integer;
typedef typename AK::Rational Rational;
CGAL_USE_TYPE(Integer);
CGAL_USE_TYPE(Rational);
// Enforce IEEE double precision and to nearest before
// using modular arithmetic
CGAL::Protect_FPU_rounding<true> pfr(CGAL_FE_TONEAREST);
typedef CGAL::Polynomial< CGAL::Residue > Poly;
typedef CGAL::Polynomial_traits_d<Poly> PT;
std::cerr << std::endl;
std::cerr <<
"Test for coefficient type CGAL::Residue"
<< std::endl;
std::cerr <<
"----------------------------------------------------------------------"
<< std::endl;
CGAL::Test_Pol::test_multiple_dimensions(PT());
}
int main() {
{
typedef CGAL::Interval_nt<true> NT;
typedef CGAL::Field_with_sqrt_tag Tag;
typedef CGAL::Tag_false Is_exact;
CGAL::test_algebraic_structure<NT,Tag, Is_exact>();
CGAL::test_algebraic_structure<NT,Tag, Is_exact>(NT(4),NT(6),NT(15));
CGAL::test_algebraic_structure<NT,Tag, Is_exact>(NT(-4),NT(6),NT(15));
CGAL::test_algebraic_structure<NT,Tag, Is_exact>(NT(4),NT(-6),NT(15));
CGAL::test_algebraic_structure<NT,Tag, Is_exact>(NT(-4),NT(-6),NT(15));
CGAL::test_algebraic_structure<NT,Tag, Is_exact>(NT(4),NT(6),NT(-15));
CGAL::test_algebraic_structure<NT,Tag, Is_exact>(NT(-4),NT(6), NT(15));
CGAL::test_algebraic_structure<NT,Tag, Is_exact>(NT(4),NT(-6),NT(-15));
CGAL::test_algebraic_structure<NT,Tag, Is_exact>(NT(-4),NT(-6),NT(-15));
CGAL::test_real_embeddable<NT>();
}
{
CGAL::Protect_FPU_rounding<true> pfr(CGAL_FE_UPWARD);
typedef CGAL::Interval_nt<false> NT;
typedef CGAL::Field_with_sqrt_tag Tag;
typedef CGAL::Tag_false Is_exact;
CGAL::test_algebraic_structure<NT,Tag, Is_exact>();
CGAL::test_algebraic_structure<NT,Tag, Is_exact>(NT(4),NT(6),NT(15));
CGAL::test_algebraic_structure<NT,Tag, Is_exact>(NT(-4),NT(6),NT(15));
CGAL::test_algebraic_structure<NT,Tag, Is_exact>(NT(4),NT(-6),NT(15));
CGAL::test_algebraic_structure<NT,Tag, Is_exact>(NT(-4),NT(-6),NT(15));
CGAL::test_algebraic_structure<NT,Tag, Is_exact>(NT(4),NT(6),NT(-15));
CGAL::test_algebraic_structure<NT,Tag, Is_exact>(NT(-4),NT(6), NT(15));
CGAL::test_algebraic_structure<NT,Tag, Is_exact>(NT(4),NT(-6),NT(-15));
CGAL::test_algebraic_structure<NT,Tag, Is_exact>(NT(-4),NT(-6),NT(-15));
CGAL::test_real_embeddable<NT>();
}
typedef CGAL::Interval_nt<true> Interval;
{ // DEFAULT CONSTRUCTOR I(0)
Interval I; // TODO: Standard constructor does not initialize I with 0.
/* assert(I==.0);
assert(I.inf()==.0);
assert(I.sup()==.0);*/
}
{ // CONSTRUCTOR FROM INT
Interval I(1);
assert(I.inf()==1.0);
assert(I.sup()==1.0);
I = Interval(1,2);
assert(I.inf()==1.0);
assert(I.sup()==2.0);
}
{ // CONSTRUCTOR FROM double
Interval I(1.0);
assert(I.inf()==1.0);
assert(I.sup()==1.0);
I= Interval(1.0,2.0);
assert(I.inf()==1.0);
assert(I.sup()==2.0);
}
{ // assign
// TODO: No assign available in Interval_nt
/* Interval I;
I.assign(2.0,3.0);
assert(I.inf()==2.0);
assert(I.sup()==3.0);
I.assign(-2.0,-1.0);
assert(I.inf()==-2.0);
assert(I.sup()==-1.0); */
}
{ //comparison
Interval I,J;
I=Interval(2);
J=Interval(2);
assert( (I==J));
assert(!(I!=J));
I=Interval(2);
J=Interval(3);
assert(!(I==J));
assert( (I!=J));
// I < J
I=Interval(1,2);
J=Interval(3,4);
assert( (I<J));
assert(!(J<I));
assert( (I<=J));
assert(!(J<=I));
assert(!(I>J));
assert( (J>I));
assert(!(I>=J));
assert( (J>=I));
// OVERLAP
I=Interval(1,3);
J=Interval(2,4); // TODO: must explicitly convert to bool to get exception
CGAL_catch_error((bool)(I==J),CGAL::Uncertain_conversion_exception&);
CGAL_catch_error((bool)(I!=J),CGAL::Uncertain_conversion_exception&);
CGAL_catch_error((bool)(I< J),CGAL::Uncertain_conversion_exception&);
CGAL_catch_error((bool)(I> J),CGAL::Uncertain_conversion_exception&);
CGAL_catch_error((bool)(I<=J),CGAL::Uncertain_conversion_exception&);
CGAL_catch_error((bool)(I>=J),CGAL::Uncertain_conversion_exception&);
CGAL_catch_error((bool)(J> I),CGAL::Uncertain_conversion_exception&);
CGAL_catch_error((bool)(J> I),CGAL::Uncertain_conversion_exception&);
CGAL_catch_error((bool)(J>=I),CGAL::Uncertain_conversion_exception&);
CGAL_catch_error((bool)(J<=I),CGAL::Uncertain_conversion_exception&);
// I<=J
I=Interval(1,2);
J=Interval(2,3);
assert( (I<=J));
assert( (J>=I));
assert(!(I> J));
assert(!(J< I));
CGAL_catch_error((bool)(I==J),CGAL::Uncertain_conversion_exception&);
CGAL_catch_error((bool)(I!=J),CGAL::Uncertain_conversion_exception&);
CGAL_catch_error((bool)(I< J),CGAL::Uncertain_conversion_exception&);
CGAL_catch_error((bool)(J> I),CGAL::Uncertain_conversion_exception&);
CGAL_catch_error((bool)(I>=J),CGAL::Uncertain_conversion_exception&);
CGAL_catch_error((bool)(J<=I),CGAL::Uncertain_conversion_exception&);
// degenerated I
I=Interval(1,1);
J=Interval(1,1);
assert( (I==J));
assert(!(I!=J));
assert(!(I> J));
assert(!(I< J));
assert( (I>=J));
assert( (I<=J));
// "I==J"
I=Interval(1,2);
J=Interval(1,2);
CGAL_catch_error((bool)(I==J),CGAL::Uncertain_conversion_exception&);
CGAL_catch_error((bool)(I!=J),CGAL::Uncertain_conversion_exception&);
CGAL_catch_error((bool)(I< J),CGAL::Uncertain_conversion_exception&);
CGAL_catch_error((bool)(I> J),CGAL::Uncertain_conversion_exception&);
CGAL_catch_error((bool)(I>=J),CGAL::Uncertain_conversion_exception&);
CGAL_catch_error((bool)(I<=J),CGAL::Uncertain_conversion_exception&);
}
{// external functions on Intervals
// functions (abs, square, sqrt, pow)
{//abs
{
Interval I=CGAL_NTS abs(Interval(2.0,3.0));
assert(I.inf()==2.0);
assert(I.sup()==3.0);
}{
Interval I=CGAL_NTS abs(Interval(-4.0,3.0));
assert(I.inf()==0.0);
assert(I.sup()==4.0);
}{
Interval I=CGAL_NTS abs(Interval(-4.0,-2.0));
assert(I.inf()==2.0);
assert(I.sup()==4.0);
}
}{//square
{
Interval I=CGAL_NTS square(Interval(2.0,3.0));
assert(I.inf()==4.0);
assert(I.sup()==9.0);
}{
Interval I=CGAL_NTS square(Interval(-4.0,3.0));
assert(I.inf()==0.0);
assert(I.sup()==16.0);
}{
Interval I=CGAL_NTS square(Interval(-4.0,-2.0));
assert(I.inf()==4.0);
assert(I.sup()==16.0);
}
}{//sqrt
{
Interval I=CGAL_NTS sqrt(Interval(2.0,3.0));
assert(1.4<I.inf()&& I.inf()<1.5);
assert(1.7<I.sup()&& I.sup()<1.8);
assert(CGAL::in(CGAL_NTS sqrt(2.0),I));
assert(CGAL::in(CGAL_NTS sqrt(3.0),I));
}{
Interval I=CGAL_NTS sqrt(Interval(-4.0,3.0));
assert(I.inf()==0.0);
assert(1.7<I.sup()&& I.sup()<1.8);
assert(CGAL::in(sqrt(3.0),I));
}{
// TODO: Throws no exception
// CGAL_catch_error(sqrt(Interval(-4.0,-2.0)),
// ...);
}
}{//pow TODO: Not available for Interval_nt
/* {
Interval I=CGAL::pow(Interval(2.0,3.0),3);
assert(I.inf()==8.0);
assert(I.sup()==27.0);
}{
Interval I=CGAL::pow(Interval(-2.0,3.0),3);
assert(I.inf()==-8.0);
assert(I.sup()==27.0);
}{
Interval I=CGAL::pow(Interval(-4.0,-2.0),3);
assert(I.inf()==-64.0);
assert(I.sup()==-8.0);
}*/
}
}{// functions min max
{//min
{
Interval I=(CGAL::min)(Interval(-2.0,-1.0),Interval(1.0,4.0));
assert(I.inf()==-2.0);
assert(I.sup()==-1.0);
}{
Interval I=(CGAL::min)(Interval(2.0,3.0),Interval(1.0,4.0));
assert(I.inf()==1.0);
assert(I.sup()==3.0);
}
}{//max
{
Interval I=(CGAL::max)(Interval(-2.0,-1.0),Interval(1.0,4.0));
assert(I.inf()==1.0);
assert(I.sup()==4.0);
}{
Interval I=(CGAL::max)(Interval(2.0,3.0),Interval(1.0,4.0));
assert(I.inf()==2.0);
assert(I.sup()==4.0);
}
}
}{// functions width, median, singleton
{//width TODO: Not available for Interval_nt
/* assert(CGAL::width(Interval(2.0,2.0))==0.0);
assert(CGAL::width(Interval(2.0,3.0))==1.0);
assert(CGAL::width(Interval(-2.0,3.0))==5.0);*/
}{//median TODO: Is simulated by to_double
assert(CGAL_NTS to_double(Interval(2.0,2.0))==2.0);
assert(CGAL_NTS to_double(Interval(2.0,3.0))==2.5);
assert(CGAL_NTS to_double(Interval(-2.0,3.0))==0.5);
}{ // TODO: Is simulated by is_point
assert((Interval(3.0,3.0)).is_point() ==true);
assert((Interval(2.0,3.0)).is_point() ==false);
}
}{// functions in, in_zero
{//in
assert(CGAL::in(2.0,Interval( 2.0,2.0))==true);
assert(CGAL::in(2.5,Interval( 2.0,3.0))==true);
assert(CGAL::in(4.0,Interval(-2.0,3.0))==false);
}{//in_zero TODO: Not available for Interval_nt
/* assert(CGAL::in_zero(Interval( 2.0,2.0))==false);
assert(CGAL::in_zero(Interval( 2.0,3.0))==false);
assert(CGAL::in_zero(Interval(-2.0,3.0))==true);
assert(CGAL::in_zero(Interval(-2.0,0.0))==true);
assert(CGAL::in_zero(Interval( 0.0,3.0))==true);*/
}
}{ // functions equal, subset, proper_subset, overlap, intersect, hull
{//equal TODO: Is simulated by is_same
assert(Interval(2.0,2.0).is_same(Interval(2.0,2.0))==true);
assert(Interval(1.0,2.0).is_same(Interval(1.0,2.0))==true);
assert(Interval(0.0,2.0).is_same(Interval(1.0,2.0))==false);
}{//subset TODO: Not available for Interval_nt
/* assert(CGAL::subset(Interval(2.0,2.0),Interval(2.0,2.0))==true);
assert(CGAL::subset(Interval(1.0,2.0),Interval(0.0,3.0))==true);
assert(CGAL::subset(Interval(0.0,2.0),Interval(1.0,2.0))==false);*/
}{//proper_subset TODO: Not available for Interval_nt
/* assert(CGAL::proper_subset(Interval(2.0,2.0),Interval(2.0,2.0))
==false);
assert(CGAL::proper_subset(Interval(1.0,2.0),Interval(0.0,3.0))
==true);
assert(CGAL::proper_subset(Interval(0.0,2.0),Interval(1.0,2.0))
==false); */
}{//overlap TODO: Is simulated by do_overlap
assert(Interval( 2,2).do_overlap(Interval(2,2))==true);
assert(Interval(-1,3).do_overlap(Interval(1,2))==true);
assert(Interval( 0,2).do_overlap(Interval(2,3))==true);
assert(Interval( 2,3).do_overlap(Interval(3,4))==true);
assert(Interval(-2,1).do_overlap(Interval(3,4))==false);
assert(Interval( 2,3).do_overlap(Interval(5,6))==false);
}{//intersect TODO: Not available for Interval_nt
/* {
Interval I=CGAL::intersect(Interval(2,2),Interval(2,2));
assert(CGAL::equal(I,Interval(2.0,2.0)));
}{
Interval I=CGAL::intersect(Interval(0,2),Interval(1,3));
assert(CGAL::equal(I,Interval(1,2)));
}*/
}{//hull TODO: Not available for Interval_nt
/* {
Interval I=CGAL::hull(Interval(2,2),Interval(2,2));
assert(CGAL::equal(I,Interval(2,2)));
}{
Interval I=CGAL::hull(Interval(0,2),Interval(1,3));
assert(CGAL::equal(I,Interval(0,3)));
}{
Interval I=CGAL::hull(Interval(-3,-1),Interval(1,3));
assert(CGAL::equal(I,Interval(-3,3)));
}*/
}
}
// IO ------------------------------------------------------------------------
/* {
// input/output
Interval tmp1,tmp2;
// tmp IS_GENERAL = sqrt 2
tmp1 = Interval(1.5,2.1);
std::ostringstream os;
os << LiS::oformat(tmp1);
std::istringstream is(os.str());
is >> LiS::iformat(tmp2);
assert(tmp1.inf()==tmp2.inf());
assert(tmp1.sup()==tmp2.sup());
}*/
return 0;
}
