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

Less calls to abs()

This commit is contained in:
Andreas Fabri 2025-07-01 17:28:34 +01:00
parent f2e16bd0c9
commit 68d2be2798
1 changed files with 54 additions and 88 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

142
NewKernel_d/include/CGAL/NewKernel_d/orientationC4.h
View File

@ -115,94 +115,60 @@ inline int orientationC4( double p0, double p1, double p2, double p3, double q0,
det = determinant( m01, m02, m03, m04, m11, m12, m13, m14, m21, m22, m23, m24, m31, m32, m33, m34 ); det = determinant( m01, m02, m03, m04, m11, m12, m13, m14, m21, m22, m23, m24, m31, m32, m33, m34 );
int int_tmp_result; int int_tmp_result;
double eps; double eps;
double max1 = fabs(m01); double max1 = CGAL::abs(m01);
if( (max1 < fabs(m02)) ) double am = CGAL::abs(m02);
{ if( (max1 < am) ) { max1 = am; }
max1 = fabs(m02); am = CGAL::abs(m03);
} if( (max1 < am) ) { max1 = am; }
if( (max1 < fabs(m03)) ) am = CGAL::abs(m11);
{ if( (max1 < am) ) { max1 = am; }
max1 = fabs(m03); am = CGAL::abs(m12);
} if( (max1 < am) ) { max1 = am; }
if( (max1 < fabs(m11)) ) am = CGAL::abs(m13);
{ if( (max1 < am) ) { max1 = am; }
max1 = fabs(m11); am = CGAL::abs(m23);
} if( (max1 < am) ) { max1 = am; }
if( (max1 < fabs(m12)) )
{
max1 = fabs(m12); double max2 = CGAL::abs(m01);
} am = CGAL::abs(m02);
if( (max1 < fabs(m13)) ) if( (max2 < am) ) { max2 = am; }
{ am = CGAL::abs(m11);
max1 = fabs(m13); if( (max2 < am) ) { max2 = am; }
} am = CGAL::abs(m12);
if( (max1 < fabs(m23)) ) if( (max2 < am) ) { max2 = am; }
{ am = CGAL::abs(m21);
max1 = fabs(m23); if( (max2 < am) ) { max2 = am; }
} am = CGAL::abs(m22);
double max2 = fabs(m01); if( (max2 < am) ) { max2 = am; }
if( (max2 < fabs(m02)) ) am = CGAL::abs(m23);
{ if( (max2 < am) ) { max2 = am; }
max2 = fabs(m02); am = CGAL::abs(m33);
} if( (max2 < am) ) { max2 = am; }
if( (max2 < fabs(m11)) )
{
max2 = fabs(m11); double max3 = CGAL::abs(m04);
} am = CGAL::abs(m14);
if( (max2 < fabs(m12)) ) if( (max3 < am) ) { max3 = am; }
{ am = CGAL::abs(m24);
max2 = fabs(m12); if( (max3 < am) ) { max3 = am; }
} am = CGAL::abs(m34);
if( (max2 < fabs(m21)) ) if( (max3 < am) ) { max3 = am; }
{
max2 = fabs(m21);
} double max4 = CGAL::abs(m11);
if( (max2 < fabs(m22)) ) am = CGAL::abs(m12);
{ if( (max4 < am) ) { max4 = am; }
max2 = fabs(m22); am = CGAL::abs(m21);
} if( (max4 < am) ) { max4 = am; }
if( (max2 < fabs(m23)) ) am = CGAL::abs(m22);
{ if( (max4 < am) ) { max4 = am; }
max2 = fabs(m23); am = CGAL::abs(m31);
} if( (max4 < am) ) { max4 = am; }
if( (max2 < fabs(m33)) ) am = CGAL::abs(m32);
{ if( (max4 < am) ) { max4 = am; }
max2 = fabs(m33);
}
double max3 = fabs(m04);
if( (max3 < fabs(m14)) )
{
max3 = fabs(m14);
}
if( (max3 < fabs(m24)) )
{
max3 = fabs(m24);
}
if( (max3 < fabs(m34)) )
{
max3 = fabs(m34);
}
double max4 = fabs(m11);
if( (max4 < fabs(m12)) )
{
max4 = fabs(m12);
}
if( (max4 < fabs(m21)) )
{
max4 = fabs(m21);
}
if( (max4 < fabs(m22)) )
{
max4 = fabs(m22);
}
if( (max4 < fabs(m31)) )
{
max4 = fabs(m31);
}
if( (max4 < fabs(m32)) )
{
max4 = fabs(m32);
}
double lower_bound_1; double lower_bound_1;
double upper_bound_1; double upper_bound_1;
lower_bound_1 = max1; lower_bound_1 = max1;