From 5e9d7ec2e4fb2b0903260b5eaa2746ef72cd47ae Mon Sep 17 00:00:00 2001 From: Sylvain Pion Date: Thu, 11 Mar 1999 15:19:05 +0000 Subject: [PATCH] - Comments. --- .../include/CGAL/Interval_arithmetic.h | 8 ++++++-- 1 file changed, 6 insertions(+), 2 deletions(-) diff --git a/Packages/Interval_arithmetic/include/CGAL/Interval_arithmetic.h b/Packages/Interval_arithmetic/include/CGAL/Interval_arithmetic.h index 7d3b7d250ba..289d794886f 100644 --- a/Packages/Interval_arithmetic/include/CGAL/Interval_arithmetic.h +++ b/Packages/Interval_arithmetic/include/CGAL/Interval_arithmetic.h @@ -34,6 +34,10 @@ // - The first one supposes the rounding mode is set -> +infinity before // nearly all operations, and might set it -> +infinity when leaving, whereas // the second leaves the rounding -> nearest. +// +// Note: When rounding is towards +infinity, to make an operation rounded +// towards -infinity, it's enough to take the opposite of some of the operand, +// and the opposite of the result (see operator+, operator*,...). #include #include @@ -210,7 +214,7 @@ Interval_nt_advanced::operator* (const Interval_nt_advanced & d) const // d>=0 [inf*d.inf; sup*d.sup] // d<=0 [sup*d.inf; inf*d.sup] // d~=0 [sup*d.inf; sup*d.sup] - /* register */ double a = inf, b = sup; + double a = inf, b = sup; if (d.inf < 0) { a=b; @@ -224,7 +228,7 @@ Interval_nt_advanced::operator* (const Interval_nt_advanced & d) const // d>=0 [inf*d.sup; sup*d.inf] // d<=0 [sup*d.sup; inf*d.inf] // d~=0 [inf*d.sup; inf*d.inf] - /* register */ double a = sup, b = inf; + double a = sup, b = inf; if (d.inf < 0) { a=b;