songsenand
/
cgal
mirror of https://github.com/CGAL/cgal
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

wip

This commit is contained in:
Brian Spilsbury 2021-03-03 00:32:09 +09:00
parent 6d048b92e1
commit 7f8f7e2f23
2 changed files with 8 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
Number_types/include/CGAL/Interval_nt.h
View File

@ -1142,16 +1142,20 @@ namespace INTERN_INTERVAL_NT {
// it helps significantly, it might even hurt by introducing a
// dependency.
}
#else
#else // no __AVX512F__
// TODO: Alternative for computing CGAL_IA_SQRT_DOWN(d.inf()) exactly
// without changing the rounding mode:
// - compute x = CGAL_IA_SQRT(d.inf())
// - compute y = CGAL_IA_SQUARE(x)
// - if y==d.inf() use x, else use -CGAL_IA_SUB(CGAL_IA_MIN_DOUBLE,x)
FPU_set_cw(CGAL_FE_DOWNWARD);
#ifdef CGAL_ALWAYS_ROUND_TO_NEAREST
double i = (d.inf() > 0.0) ? nextafter(std::sqrt(d.inf()), 0.) : 0.0;
#else // no CGAL_ALWAYS_ROUND_TO_NEAREST
FPU_set_cw(CGAL_FE_DOWNWARD);
double i = (d.inf() > 0.0) ? CGAL_IA_SQRT(d.inf()) : 0.0;
FPU_set_cw(CGAL_FE_UPWARD);
#endif
#endif // no CGAL_ALWAYS_ROUND_TO_NEAREST
#endif // no __AVX512F__
return Interval_nt<Protected>(i, CGAL_IA_SQRT(d.sup()));
}

2
Triangulation_3/include/CGAL/Delaunay_triangulation_cell_base_with_circumcenter_3.h
View File

@ -130,7 +130,7 @@ public:
circumcenter_ = new Point(this->Cb::circumcenter(gt));
} else {
CGAL_expensive_assertion(
this->Cb::circumcenter(gt) == *circumcenter);
this->Cb::circumcenter(gt) == *circumcenter_);
}
return *circumcenter_;