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

Use SSE2 based versions of abs, min, max 

Their usage is controlled with CGAL_USE_SSE2_FABS CGAL_USE_SSE2_MAX
This commit is contained in:
Andreas Fabri 2011-04-01 14:02:58 +00:00
parent bce0e8c4b4
commit e7748c9dce
2 changed files with 49 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

31
Filtered_kernel/include/CGAL/internal/Static_filters/Orientation_3.h
View File

@ -28,6 +28,10 @@
namespace CGAL { namespace internal { namespace Static_filters_predicates {
#include <iostream>
template < typename K_base >
class Orientation_3
: public K_base::Orientation_3
@ -97,6 +101,7 @@ public:
// CGAL::abs uses fabs on platforms where it is faster than (a<0)?-a:a
// Then semi-static filter.
double maxx = CGAL::abs(pqx);
double maxy = CGAL::abs(pqy);
double maxz = CGAL::abs(pqz);
@ -109,18 +114,37 @@ public:
double aprz = CGAL::abs(prz);
double apsz = CGAL::abs(psz);
#ifdef CGAL_USE_SSE2_MAX
CGAL::Max<double> mmax;
maxx = mmax(maxx, aprx, apsx);
maxy = mmax(maxy, apry, apsy);
maxz = mmax(maxz, aprz, apsz);
#else
if (maxx < aprx) maxx = aprx;
if (maxx < apsx) maxx = apsx;
if (maxy < apry) maxy = apry;
if (maxy < apsy) maxy = apsy;
if (maxz < aprz) maxz = aprz;
if (maxz < apsz) maxz = apsz;
#endif
double det = CGAL::determinant(pqx, pqy, pqz,
prx, pry, prz,
psx, psy, psz);
double eps = 5.1107127829973299e-15 * maxx * maxy * maxz;
#ifdef CGAL_USE_SSE2_MAX
#if 0
CGAL::Min<double> mmin;
double tmp = mmin(maxx, maxy, maxz);
maxz = mmax(maxx, maxy, maxz);
maxx = tmp;
#else
sse2minmax(maxx,maxy,maxz);
// maxy can contain ANY element
#endif
#else
// Sort maxx < maxy < maxz.
if (maxx > maxz)
std::swap(maxx, maxz);
@ -128,7 +152,7 @@ public:
std::swap(maxy, maxz);
else if (maxy < maxx)
std::swap(maxx, maxy);
#endif
// Protect against underflow in the computation of eps.
if (maxx < 1e-97) /* cbrt(min_double/eps) */ {
if (maxx == 0)
@ -136,11 +160,10 @@ public:
}
// Protect against overflow in the computation of det.
else if (maxz < 1e102) /* cbrt(max_double [hadamard]/4) */ {
double eps = 5.1107127829973299e-15 * maxx * maxy * maxz;
if (det > eps) return POSITIVE;
if (det < -eps) return NEGATIVE;
}
CGAL_BRANCH_PROFILER_BRANCH_2(tmp);
}

25
Filtered_kernel/include/CGAL/internal/Static_filters/Side_of_oriented_sphere_3.h
View File

@ -98,6 +98,12 @@ public:
double artz = CGAL::abs(rtz);
double astz = CGAL::abs(stz);
#ifdef CGAL_USE_SSE2_MAX
CGAL::Max<double> mmax;
maxx = mmax(maxx, aqtx, artx, astx);
maxy = mmax(maxy, aqty, arty, asty);
maxz = mmax(maxz, aqtz, artz, astz);
#else
if (maxx < aqtx) maxx = aqtx;
if (maxx < artx) maxx = artx;
if (maxx < astx) maxx = astx;
@ -109,7 +115,21 @@ public:
if (maxz < aqtz) maxz = aqtz;
if (maxz < artz) maxz = artz;
if (maxz < astz) maxz = astz;
#endif
double eps = 1.2466136531027298e-13 * maxx * maxy * maxz;
#ifdef CGAL_USE_SSE2_MAX
/*
CGAL::Min<double> mmin;
double tmp = mmin(maxx, maxy, maxz);
maxz = mmax(maxx, maxy, maxz);
maxx = tmp;
*/
sse2minmax(maxx,maxy,maxz);
// maxy can contain ANY element
#else
// Sort maxx < maxy < maxz.
if (maxx > maxz)
std::swap(maxx, maxz);
@ -117,7 +137,7 @@ public:
std::swap(maxy, maxz);
else if (maxy < maxx)
std::swap(maxx, maxy);
#endif
double det = CGAL::determinant(ptx,pty,ptz,pt2,
rtx,rty,rtz,rt2,
qtx,qty,qtz,qt2,
@ -130,8 +150,7 @@ public:
}
// Protect against overflow in the computation of det.
else if (maxz < 1e61) /* sqrt^5(max_double/4 [hadamard]) */ {
double eps = 1.2466136531027298e-13 * maxx * maxy * maxz
* (maxz * maxz);
eps *= (maxz * maxz);
if (det > eps) return ON_POSITIVE_SIDE;
if (det < -eps) return ON_NEGATIVE_SIDE;
}