add assertion to make sure that the rounding mode is good prior to th…

…e call
CGAL · Nov 20, 2024 · 42d435c · 42d435c
1 parent 4ceec7d
commit 42d435c
Showing 1 changed file with 1 addition and 0 deletions.
diff --git a/Number_types/include/CGAL/FPU.h b/Number_types/include/CGAL/FPU.h
@@ -618,6 +618,7 @@ inline double IA_sqrt_toward_zero(double d) {
 #ifdef CGAL_ALWAYS_ROUND_TO_NEAREST
   return (d > 0.0) ? nextafter(std::sqrt(d), 0.) : 0.0;
 #else
+  CGAL_assertion(FPU_get_cw()==CGAL_FE_UPWARD);
   FPU_set_cw(CGAL_FE_DOWNWARD);
   double i = (d > 0.0) ? CGAL_IA_FORCE_TO_DOUBLE(CGAL_BUG_SQRT(CGAL_IA_STOP_CPROP(d))) : 0.0;
   FPU_set_cw(CGAL_FE_UPWARD);