mirror of https://github.com/CGAL/cgal
- Changed .value() -> .to_double(), and .to_double() -> .dbl().
This commit is contained in:
Sylvain Pion
parent
8e15655d50
commit
a37e58669e
|
|
@ -33,16 +33,6 @@
|
|||
#include <CGAL/Arithmetic_filter/predicates_on_ftC3.h>
|
||||
#endif
|
||||
|
||||
// #if defined( CGAL_PREDICATES_ON_RTH2_H ) && \
|
||||
// !defined( CGAL_ARITHMETIC_FILTER_PREDICATES_ON_RTH2_H )
|
||||
// #include <CGAL/Arithmetic_filter/predicates_on_rtH2.h>
|
||||
// #endif
|
||||
//
|
||||
// #if defined( CGAL_PREDICATES_ON_RTH3_H ) && \
|
||||
// !defined( CGAL_ARITHMETIC_FILTER_PREDICATES_ON_RTH3_H )
|
||||
// #include <CGAL/Arithmetic_filter/predicates_on_rtH3.h>
|
||||
// #endif
|
||||
|
||||
#if defined( CGAL_PREDICATES_SIGN_OF_DETERMINANT_H ) && \
|
||||
!defined( CGAL_ARITHMETIC_FILTER_PREDICATES_SIGN_OF_DETERMINANT_H )
|
||||
#include <CGAL/Arithmetic_filter/predicates/sign_of_determinant.h>
|
||||
|
|
|
|||
|
|
@ -197,34 +197,34 @@ power_testC2(
|
|||
|
||||
// Check the bounds. All arguments must be <= SAF_bound.
|
||||
if (
|
||||
fabs(px.value()) > SAF_bound ||
|
||||
fabs(py.value()) > SAF_bound ||
|
||||
fabs(pwt.value()) > SAF_bound ||
|
||||
fabs(qx.value()) > SAF_bound ||
|
||||
fabs(qy.value()) > SAF_bound ||
|
||||
fabs(qwt.value()) > SAF_bound ||
|
||||
fabs(rx.value()) > SAF_bound ||
|
||||
fabs(ry.value()) > SAF_bound ||
|
||||
fabs(rwt.value()) > SAF_bound ||
|
||||
fabs(tx.value()) > SAF_bound ||
|
||||
fabs(ty.value()) > SAF_bound ||
|
||||
fabs(twt.value()) > SAF_bound)
|
||||
fabs(px.to_double()) > SAF_bound ||
|
||||
fabs(py.to_double()) > SAF_bound ||
|
||||
fabs(pwt.to_double()) > SAF_bound ||
|
||||
fabs(qx.to_double()) > SAF_bound ||
|
||||
fabs(qy.to_double()) > SAF_bound ||
|
||||
fabs(qwt.to_double()) > SAF_bound ||
|
||||
fabs(rx.to_double()) > SAF_bound ||
|
||||
fabs(ry.to_double()) > SAF_bound ||
|
||||
fabs(rwt.to_double()) > SAF_bound ||
|
||||
fabs(tx.to_double()) > SAF_bound ||
|
||||
fabs(ty.to_double()) > SAF_bound ||
|
||||
fabs(twt.to_double()) > SAF_bound)
|
||||
{
|
||||
re_adjust:
|
||||
// Compute the new bound.
|
||||
double NEW_bound = 0.0;
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pwt.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qwt.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ry.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rwt.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(tx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ty.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(twt.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pwt.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qwt.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ry.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rwt.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(tx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ty.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(twt.to_double()));
|
||||
// Re-adjust the context.
|
||||
Static_Filtered_power_testC2_12::new_bound(NEW_bound);
|
||||
}
|
||||
|
|
@ -232,18 +232,18 @@ re_adjust:
|
|||
try
|
||||
{
|
||||
return Static_Filtered_power_testC2_12::epsilon_variant(
|
||||
px.to_double(),
|
||||
py.to_double(),
|
||||
pwt.to_double(),
|
||||
qx.to_double(),
|
||||
qy.to_double(),
|
||||
qwt.to_double(),
|
||||
rx.to_double(),
|
||||
ry.to_double(),
|
||||
rwt.to_double(),
|
||||
tx.to_double(),
|
||||
ty.to_double(),
|
||||
twt.to_double(),
|
||||
px.dbl(),
|
||||
py.dbl(),
|
||||
pwt.dbl(),
|
||||
qx.dbl(),
|
||||
qy.dbl(),
|
||||
qwt.dbl(),
|
||||
rx.dbl(),
|
||||
ry.dbl(),
|
||||
rwt.dbl(),
|
||||
tx.dbl(),
|
||||
ty.dbl(),
|
||||
twt.dbl(),
|
||||
SAF_epsilon_0);
|
||||
}
|
||||
catch (Restricted_double::unsafe_comparison)
|
||||
|
|
@ -435,28 +435,28 @@ power_testC2(
|
|||
|
||||
// Check the bounds. All arguments must be <= SAF_bound.
|
||||
if (
|
||||
fabs(px.value()) > SAF_bound ||
|
||||
fabs(py.value()) > SAF_bound ||
|
||||
fabs(pwt.value()) > SAF_bound ||
|
||||
fabs(qx.value()) > SAF_bound ||
|
||||
fabs(qy.value()) > SAF_bound ||
|
||||
fabs(qwt.value()) > SAF_bound ||
|
||||
fabs(tx.value()) > SAF_bound ||
|
||||
fabs(ty.value()) > SAF_bound ||
|
||||
fabs(twt.value()) > SAF_bound)
|
||||
fabs(px.to_double()) > SAF_bound ||
|
||||
fabs(py.to_double()) > SAF_bound ||
|
||||
fabs(pwt.to_double()) > SAF_bound ||
|
||||
fabs(qx.to_double()) > SAF_bound ||
|
||||
fabs(qy.to_double()) > SAF_bound ||
|
||||
fabs(qwt.to_double()) > SAF_bound ||
|
||||
fabs(tx.to_double()) > SAF_bound ||
|
||||
fabs(ty.to_double()) > SAF_bound ||
|
||||
fabs(twt.to_double()) > SAF_bound)
|
||||
{
|
||||
re_adjust:
|
||||
// Compute the new bound.
|
||||
double NEW_bound = 0.0;
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pwt.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qwt.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(tx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ty.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(twt.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pwt.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qwt.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(tx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ty.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(twt.to_double()));
|
||||
// Re-adjust the context.
|
||||
Static_Filtered_power_testC2_9::new_bound(NEW_bound);
|
||||
}
|
||||
|
|
@ -464,15 +464,15 @@ re_adjust:
|
|||
try
|
||||
{
|
||||
return Static_Filtered_power_testC2_9::epsilon_variant(
|
||||
px.to_double(),
|
||||
py.to_double(),
|
||||
pwt.to_double(),
|
||||
qx.to_double(),
|
||||
qy.to_double(),
|
||||
qwt.to_double(),
|
||||
tx.to_double(),
|
||||
ty.to_double(),
|
||||
twt.to_double(),
|
||||
px.dbl(),
|
||||
py.dbl(),
|
||||
pwt.dbl(),
|
||||
qx.dbl(),
|
||||
qy.dbl(),
|
||||
qwt.dbl(),
|
||||
tx.dbl(),
|
||||
ty.dbl(),
|
||||
twt.dbl(),
|
||||
SAF_epsilon_0,
|
||||
SAF_epsilon_1,
|
||||
SAF_epsilon_2,
|
||||
|
|
|
|||
|
|
@ -257,50 +257,50 @@ power_testC3(
|
|||
|
||||
// Check the bounds. All arguments must be <= SAF_bound.
|
||||
if (
|
||||
fabs(px.value()) > SAF_bound ||
|
||||
fabs(py.value()) > SAF_bound ||
|
||||
fabs(pz.value()) > SAF_bound ||
|
||||
fabs(pwt.value()) > SAF_bound ||
|
||||
fabs(qx.value()) > SAF_bound ||
|
||||
fabs(qy.value()) > SAF_bound ||
|
||||
fabs(qz.value()) > SAF_bound ||
|
||||
fabs(qwt.value()) > SAF_bound ||
|
||||
fabs(rx.value()) > SAF_bound ||
|
||||
fabs(ry.value()) > SAF_bound ||
|
||||
fabs(rz.value()) > SAF_bound ||
|
||||
fabs(rwt.value()) > SAF_bound ||
|
||||
fabs(sx.value()) > SAF_bound ||
|
||||
fabs(sy.value()) > SAF_bound ||
|
||||
fabs(sz.value()) > SAF_bound ||
|
||||
fabs(swt.value()) > SAF_bound ||
|
||||
fabs(tx.value()) > SAF_bound ||
|
||||
fabs(ty.value()) > SAF_bound ||
|
||||
fabs(tz.value()) > SAF_bound ||
|
||||
fabs(twt.value()) > SAF_bound)
|
||||
fabs(px.to_double()) > SAF_bound ||
|
||||
fabs(py.to_double()) > SAF_bound ||
|
||||
fabs(pz.to_double()) > SAF_bound ||
|
||||
fabs(pwt.to_double()) > SAF_bound ||
|
||||
fabs(qx.to_double()) > SAF_bound ||
|
||||
fabs(qy.to_double()) > SAF_bound ||
|
||||
fabs(qz.to_double()) > SAF_bound ||
|
||||
fabs(qwt.to_double()) > SAF_bound ||
|
||||
fabs(rx.to_double()) > SAF_bound ||
|
||||
fabs(ry.to_double()) > SAF_bound ||
|
||||
fabs(rz.to_double()) > SAF_bound ||
|
||||
fabs(rwt.to_double()) > SAF_bound ||
|
||||
fabs(sx.to_double()) > SAF_bound ||
|
||||
fabs(sy.to_double()) > SAF_bound ||
|
||||
fabs(sz.to_double()) > SAF_bound ||
|
||||
fabs(swt.to_double()) > SAF_bound ||
|
||||
fabs(tx.to_double()) > SAF_bound ||
|
||||
fabs(ty.to_double()) > SAF_bound ||
|
||||
fabs(tz.to_double()) > SAF_bound ||
|
||||
fabs(twt.to_double()) > SAF_bound)
|
||||
{
|
||||
re_adjust:
|
||||
// Compute the new bound.
|
||||
double NEW_bound = 0.0;
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pwt.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qwt.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ry.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rwt.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(sx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(sy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(sz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(swt.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(tx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ty.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(tz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(twt.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pz.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pwt.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qz.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qwt.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ry.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rz.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rwt.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(sx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(sy.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(sz.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(swt.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(tx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ty.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(tz.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(twt.to_double()));
|
||||
// Re-adjust the context.
|
||||
Static_Filtered_power_testC3_20::new_bound(NEW_bound);
|
||||
}
|
||||
|
|
@ -308,26 +308,26 @@ re_adjust:
|
|||
try
|
||||
{
|
||||
return Static_Filtered_power_testC3_20::epsilon_variant(
|
||||
px.to_double(),
|
||||
py.to_double(),
|
||||
pz.to_double(),
|
||||
pwt.to_double(),
|
||||
qx.to_double(),
|
||||
qy.to_double(),
|
||||
qz.to_double(),
|
||||
qwt.to_double(),
|
||||
rx.to_double(),
|
||||
ry.to_double(),
|
||||
rz.to_double(),
|
||||
rwt.to_double(),
|
||||
sx.to_double(),
|
||||
sy.to_double(),
|
||||
sz.to_double(),
|
||||
swt.to_double(),
|
||||
tx.to_double(),
|
||||
ty.to_double(),
|
||||
tz.to_double(),
|
||||
twt.to_double(),
|
||||
px.dbl(),
|
||||
py.dbl(),
|
||||
pz.dbl(),
|
||||
pwt.dbl(),
|
||||
qx.dbl(),
|
||||
qy.dbl(),
|
||||
qz.dbl(),
|
||||
qwt.dbl(),
|
||||
rx.dbl(),
|
||||
ry.dbl(),
|
||||
rz.dbl(),
|
||||
rwt.dbl(),
|
||||
sx.dbl(),
|
||||
sy.dbl(),
|
||||
sz.dbl(),
|
||||
swt.dbl(),
|
||||
tx.dbl(),
|
||||
ty.dbl(),
|
||||
tz.dbl(),
|
||||
twt.dbl(),
|
||||
SAF_epsilon_0);
|
||||
}
|
||||
catch (Restricted_double::unsafe_comparison)
|
||||
|
|
@ -619,42 +619,42 @@ power_testC3(
|
|||
|
||||
// Check the bounds. All arguments must be <= SAF_bound.
|
||||
if (
|
||||
fabs(px.value()) > SAF_bound ||
|
||||
fabs(py.value()) > SAF_bound ||
|
||||
fabs(pz.value()) > SAF_bound ||
|
||||
fabs(pwt.value()) > SAF_bound ||
|
||||
fabs(qx.value()) > SAF_bound ||
|
||||
fabs(qy.value()) > SAF_bound ||
|
||||
fabs(qz.value()) > SAF_bound ||
|
||||
fabs(qwt.value()) > SAF_bound ||
|
||||
fabs(rx.value()) > SAF_bound ||
|
||||
fabs(ry.value()) > SAF_bound ||
|
||||
fabs(rz.value()) > SAF_bound ||
|
||||
fabs(rwt.value()) > SAF_bound ||
|
||||
fabs(tx.value()) > SAF_bound ||
|
||||
fabs(ty.value()) > SAF_bound ||
|
||||
fabs(tz.value()) > SAF_bound ||
|
||||
fabs(twt.value()) > SAF_bound)
|
||||
fabs(px.to_double()) > SAF_bound ||
|
||||
fabs(py.to_double()) > SAF_bound ||
|
||||
fabs(pz.to_double()) > SAF_bound ||
|
||||
fabs(pwt.to_double()) > SAF_bound ||
|
||||
fabs(qx.to_double()) > SAF_bound ||
|
||||
fabs(qy.to_double()) > SAF_bound ||
|
||||
fabs(qz.to_double()) > SAF_bound ||
|
||||
fabs(qwt.to_double()) > SAF_bound ||
|
||||
fabs(rx.to_double()) > SAF_bound ||
|
||||
fabs(ry.to_double()) > SAF_bound ||
|
||||
fabs(rz.to_double()) > SAF_bound ||
|
||||
fabs(rwt.to_double()) > SAF_bound ||
|
||||
fabs(tx.to_double()) > SAF_bound ||
|
||||
fabs(ty.to_double()) > SAF_bound ||
|
||||
fabs(tz.to_double()) > SAF_bound ||
|
||||
fabs(twt.to_double()) > SAF_bound)
|
||||
{
|
||||
re_adjust:
|
||||
// Compute the new bound.
|
||||
double NEW_bound = 0.0;
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pwt.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qwt.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ry.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rwt.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(tx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ty.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(tz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(twt.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pz.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pwt.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qz.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qwt.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ry.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rz.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rwt.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(tx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ty.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(tz.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(twt.to_double()));
|
||||
// Re-adjust the context.
|
||||
Static_Filtered_power_testC3_16::new_bound(NEW_bound);
|
||||
}
|
||||
|
|
@ -662,22 +662,22 @@ re_adjust:
|
|||
try
|
||||
{
|
||||
return Static_Filtered_power_testC3_16::epsilon_variant(
|
||||
px.to_double(),
|
||||
py.to_double(),
|
||||
pz.to_double(),
|
||||
pwt.to_double(),
|
||||
qx.to_double(),
|
||||
qy.to_double(),
|
||||
qz.to_double(),
|
||||
qwt.to_double(),
|
||||
rx.to_double(),
|
||||
ry.to_double(),
|
||||
rz.to_double(),
|
||||
rwt.to_double(),
|
||||
tx.to_double(),
|
||||
ty.to_double(),
|
||||
tz.to_double(),
|
||||
twt.to_double(),
|
||||
px.dbl(),
|
||||
py.dbl(),
|
||||
pz.dbl(),
|
||||
pwt.dbl(),
|
||||
qx.dbl(),
|
||||
qy.dbl(),
|
||||
qz.dbl(),
|
||||
qwt.dbl(),
|
||||
rx.dbl(),
|
||||
ry.dbl(),
|
||||
rz.dbl(),
|
||||
rwt.dbl(),
|
||||
tx.dbl(),
|
||||
ty.dbl(),
|
||||
tz.dbl(),
|
||||
twt.dbl(),
|
||||
SAF_epsilon_0,
|
||||
SAF_epsilon_1,
|
||||
SAF_epsilon_2,
|
||||
|
|
@ -920,34 +920,34 @@ power_testC3(
|
|||
|
||||
// Check the bounds. All arguments must be <= SAF_bound.
|
||||
if (
|
||||
fabs(px.value()) > SAF_bound ||
|
||||
fabs(py.value()) > SAF_bound ||
|
||||
fabs(pz.value()) > SAF_bound ||
|
||||
fabs(pwt.value()) > SAF_bound ||
|
||||
fabs(qx.value()) > SAF_bound ||
|
||||
fabs(qy.value()) > SAF_bound ||
|
||||
fabs(qz.value()) > SAF_bound ||
|
||||
fabs(qwt.value()) > SAF_bound ||
|
||||
fabs(tx.value()) > SAF_bound ||
|
||||
fabs(ty.value()) > SAF_bound ||
|
||||
fabs(tz.value()) > SAF_bound ||
|
||||
fabs(twt.value()) > SAF_bound)
|
||||
fabs(px.to_double()) > SAF_bound ||
|
||||
fabs(py.to_double()) > SAF_bound ||
|
||||
fabs(pz.to_double()) > SAF_bound ||
|
||||
fabs(pwt.to_double()) > SAF_bound ||
|
||||
fabs(qx.to_double()) > SAF_bound ||
|
||||
fabs(qy.to_double()) > SAF_bound ||
|
||||
fabs(qz.to_double()) > SAF_bound ||
|
||||
fabs(qwt.to_double()) > SAF_bound ||
|
||||
fabs(tx.to_double()) > SAF_bound ||
|
||||
fabs(ty.to_double()) > SAF_bound ||
|
||||
fabs(tz.to_double()) > SAF_bound ||
|
||||
fabs(twt.to_double()) > SAF_bound)
|
||||
{
|
||||
re_adjust:
|
||||
// Compute the new bound.
|
||||
double NEW_bound = 0.0;
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pwt.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qwt.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(tx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ty.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(tz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(twt.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pz.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pwt.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qz.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qwt.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(tx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ty.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(tz.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(twt.to_double()));
|
||||
// Re-adjust the context.
|
||||
Static_Filtered_power_testC3_12::new_bound(NEW_bound);
|
||||
}
|
||||
|
|
@ -955,18 +955,18 @@ re_adjust:
|
|||
try
|
||||
{
|
||||
return Static_Filtered_power_testC3_12::epsilon_variant(
|
||||
px.to_double(),
|
||||
py.to_double(),
|
||||
pz.to_double(),
|
||||
pwt.to_double(),
|
||||
qx.to_double(),
|
||||
qy.to_double(),
|
||||
qz.to_double(),
|
||||
qwt.to_double(),
|
||||
tx.to_double(),
|
||||
ty.to_double(),
|
||||
tz.to_double(),
|
||||
twt.to_double(),
|
||||
px.dbl(),
|
||||
py.dbl(),
|
||||
pz.dbl(),
|
||||
pwt.dbl(),
|
||||
qx.dbl(),
|
||||
qy.dbl(),
|
||||
qz.dbl(),
|
||||
qwt.dbl(),
|
||||
tx.dbl(),
|
||||
ty.dbl(),
|
||||
tz.dbl(),
|
||||
twt.dbl(),
|
||||
SAF_epsilon_0,
|
||||
SAF_epsilon_1,
|
||||
SAF_epsilon_2,
|
||||
|
|
|
|||
|
|
@ -237,42 +237,42 @@ power_testH2(
|
|||
|
||||
// Check the bounds. All arguments must be <= SAF_bound.
|
||||
if (
|
||||
fabs(phx.value()) > SAF_bound ||
|
||||
fabs(phy.value()) > SAF_bound ||
|
||||
fabs(phw.value()) > SAF_bound ||
|
||||
fabs(pwt.value()) > SAF_bound ||
|
||||
fabs(qhx.value()) > SAF_bound ||
|
||||
fabs(qhy.value()) > SAF_bound ||
|
||||
fabs(qhw.value()) > SAF_bound ||
|
||||
fabs(qwt.value()) > SAF_bound ||
|
||||
fabs(rhx.value()) > SAF_bound ||
|
||||
fabs(rhy.value()) > SAF_bound ||
|
||||
fabs(rhw.value()) > SAF_bound ||
|
||||
fabs(rwt.value()) > SAF_bound ||
|
||||
fabs(thx.value()) > SAF_bound ||
|
||||
fabs(thy.value()) > SAF_bound ||
|
||||
fabs(thw.value()) > SAF_bound ||
|
||||
fabs(twt.value()) > SAF_bound)
|
||||
fabs(phx.to_double()) > SAF_bound ||
|
||||
fabs(phy.to_double()) > SAF_bound ||
|
||||
fabs(phw.to_double()) > SAF_bound ||
|
||||
fabs(pwt.to_double()) > SAF_bound ||
|
||||
fabs(qhx.to_double()) > SAF_bound ||
|
||||
fabs(qhy.to_double()) > SAF_bound ||
|
||||
fabs(qhw.to_double()) > SAF_bound ||
|
||||
fabs(qwt.to_double()) > SAF_bound ||
|
||||
fabs(rhx.to_double()) > SAF_bound ||
|
||||
fabs(rhy.to_double()) > SAF_bound ||
|
||||
fabs(rhw.to_double()) > SAF_bound ||
|
||||
fabs(rwt.to_double()) > SAF_bound ||
|
||||
fabs(thx.to_double()) > SAF_bound ||
|
||||
fabs(thy.to_double()) > SAF_bound ||
|
||||
fabs(thw.to_double()) > SAF_bound ||
|
||||
fabs(twt.to_double()) > SAF_bound)
|
||||
{
|
||||
re_adjust:
|
||||
// Compute the new bound.
|
||||
double NEW_bound = 0.0;
|
||||
NEW_bound = std::max(NEW_bound, fabs(phx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(phy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(phw.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pwt.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qhx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qhy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qhw.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qwt.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rhx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rhy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rhw.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rwt.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(thx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(thy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(thw.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(twt.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(phx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(phy.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(phw.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pwt.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qhx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qhy.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qhw.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qwt.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rhx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rhy.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rhw.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rwt.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(thx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(thy.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(thw.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(twt.to_double()));
|
||||
// Re-adjust the context.
|
||||
Static_Filtered_power_testH2_16::new_bound(NEW_bound);
|
||||
}
|
||||
|
|
@ -280,22 +280,22 @@ re_adjust:
|
|||
try
|
||||
{
|
||||
return Static_Filtered_power_testH2_16::epsilon_variant(
|
||||
phx.to_double(),
|
||||
phy.to_double(),
|
||||
phw.to_double(),
|
||||
pwt.to_double(),
|
||||
qhx.to_double(),
|
||||
qhy.to_double(),
|
||||
qhw.to_double(),
|
||||
qwt.to_double(),
|
||||
rhx.to_double(),
|
||||
rhy.to_double(),
|
||||
rhw.to_double(),
|
||||
rwt.to_double(),
|
||||
thx.to_double(),
|
||||
thy.to_double(),
|
||||
thw.to_double(),
|
||||
twt.to_double(),
|
||||
phx.dbl(),
|
||||
phy.dbl(),
|
||||
phw.dbl(),
|
||||
pwt.dbl(),
|
||||
qhx.dbl(),
|
||||
qhy.dbl(),
|
||||
qhw.dbl(),
|
||||
qwt.dbl(),
|
||||
rhx.dbl(),
|
||||
rhy.dbl(),
|
||||
rhw.dbl(),
|
||||
rwt.dbl(),
|
||||
thx.dbl(),
|
||||
thy.dbl(),
|
||||
thw.dbl(),
|
||||
twt.dbl(),
|
||||
SAF_epsilon_0);
|
||||
}
|
||||
catch (Restricted_double::unsafe_comparison)
|
||||
|
|
@ -529,34 +529,34 @@ power_testH2(
|
|||
|
||||
// Check the bounds. All arguments must be <= SAF_bound.
|
||||
if (
|
||||
fabs(phx.value()) > SAF_bound ||
|
||||
fabs(phy.value()) > SAF_bound ||
|
||||
fabs(phw.value()) > SAF_bound ||
|
||||
fabs(pwt.value()) > SAF_bound ||
|
||||
fabs(qhx.value()) > SAF_bound ||
|
||||
fabs(qhy.value()) > SAF_bound ||
|
||||
fabs(qhw.value()) > SAF_bound ||
|
||||
fabs(qwt.value()) > SAF_bound ||
|
||||
fabs(thx.value()) > SAF_bound ||
|
||||
fabs(thy.value()) > SAF_bound ||
|
||||
fabs(thw.value()) > SAF_bound ||
|
||||
fabs(twt.value()) > SAF_bound)
|
||||
fabs(phx.to_double()) > SAF_bound ||
|
||||
fabs(phy.to_double()) > SAF_bound ||
|
||||
fabs(phw.to_double()) > SAF_bound ||
|
||||
fabs(pwt.to_double()) > SAF_bound ||
|
||||
fabs(qhx.to_double()) > SAF_bound ||
|
||||
fabs(qhy.to_double()) > SAF_bound ||
|
||||
fabs(qhw.to_double()) > SAF_bound ||
|
||||
fabs(qwt.to_double()) > SAF_bound ||
|
||||
fabs(thx.to_double()) > SAF_bound ||
|
||||
fabs(thy.to_double()) > SAF_bound ||
|
||||
fabs(thw.to_double()) > SAF_bound ||
|
||||
fabs(twt.to_double()) > SAF_bound)
|
||||
{
|
||||
re_adjust:
|
||||
// Compute the new bound.
|
||||
double NEW_bound = 0.0;
|
||||
NEW_bound = std::max(NEW_bound, fabs(phx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(phy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(phw.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pwt.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qhx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qhy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qhw.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qwt.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(thx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(thy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(thw.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(twt.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(phx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(phy.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(phw.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pwt.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qhx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qhy.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qhw.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qwt.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(thx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(thy.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(thw.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(twt.to_double()));
|
||||
// Re-adjust the context.
|
||||
Static_Filtered_power_testH2_12::new_bound(NEW_bound);
|
||||
}
|
||||
|
|
@ -564,18 +564,18 @@ re_adjust:
|
|||
try
|
||||
{
|
||||
return Static_Filtered_power_testH2_12::epsilon_variant(
|
||||
phx.to_double(),
|
||||
phy.to_double(),
|
||||
phw.to_double(),
|
||||
pwt.to_double(),
|
||||
qhx.to_double(),
|
||||
qhy.to_double(),
|
||||
qhw.to_double(),
|
||||
qwt.to_double(),
|
||||
thx.to_double(),
|
||||
thy.to_double(),
|
||||
thw.to_double(),
|
||||
twt.to_double(),
|
||||
phx.dbl(),
|
||||
phy.dbl(),
|
||||
phw.dbl(),
|
||||
pwt.dbl(),
|
||||
qhx.dbl(),
|
||||
qhy.dbl(),
|
||||
qhw.dbl(),
|
||||
qwt.dbl(),
|
||||
thx.dbl(),
|
||||
thy.dbl(),
|
||||
thw.dbl(),
|
||||
twt.dbl(),
|
||||
SAF_epsilon_0,
|
||||
SAF_epsilon_1);
|
||||
}
|
||||
|
|
|
|||
|
|
@ -313,60 +313,60 @@ power_testH3(
|
|||
|
||||
// Check the bounds. All arguments must be <= SAF_bound.
|
||||
if (
|
||||
fabs(phx.value()) > SAF_bound ||
|
||||
fabs(phy.value()) > SAF_bound ||
|
||||
fabs(phz.value()) > SAF_bound ||
|
||||
fabs(phw.value()) > SAF_bound ||
|
||||
fabs(pwt.value()) > SAF_bound ||
|
||||
fabs(qhx.value()) > SAF_bound ||
|
||||
fabs(qhy.value()) > SAF_bound ||
|
||||
fabs(qhz.value()) > SAF_bound ||
|
||||
fabs(qhw.value()) > SAF_bound ||
|
||||
fabs(qwt.value()) > SAF_bound ||
|
||||
fabs(rhx.value()) > SAF_bound ||
|
||||
fabs(rhy.value()) > SAF_bound ||
|
||||
fabs(rhz.value()) > SAF_bound ||
|
||||
fabs(rhw.value()) > SAF_bound ||
|
||||
fabs(rwt.value()) > SAF_bound ||
|
||||
fabs(shx.value()) > SAF_bound ||
|
||||
fabs(shy.value()) > SAF_bound ||
|
||||
fabs(shz.value()) > SAF_bound ||
|
||||
fabs(shw.value()) > SAF_bound ||
|
||||
fabs(swt.value()) > SAF_bound ||
|
||||
fabs(thx.value()) > SAF_bound ||
|
||||
fabs(thy.value()) > SAF_bound ||
|
||||
fabs(thz.value()) > SAF_bound ||
|
||||
fabs(thw.value()) > SAF_bound ||
|
||||
fabs(twt.value()) > SAF_bound)
|
||||
fabs(phx.to_double()) > SAF_bound ||
|
||||
fabs(phy.to_double()) > SAF_bound ||
|
||||
fabs(phz.to_double()) > SAF_bound ||
|
||||
fabs(phw.to_double()) > SAF_bound ||
|
||||
fabs(pwt.to_double()) > SAF_bound ||
|
||||
fabs(qhx.to_double()) > SAF_bound ||
|
||||
fabs(qhy.to_double()) > SAF_bound ||
|
||||
fabs(qhz.to_double()) > SAF_bound ||
|
||||
fabs(qhw.to_double()) > SAF_bound ||
|
||||
fabs(qwt.to_double()) > SAF_bound ||
|
||||
fabs(rhx.to_double()) > SAF_bound ||
|
||||
fabs(rhy.to_double()) > SAF_bound ||
|
||||
fabs(rhz.to_double()) > SAF_bound ||
|
||||
fabs(rhw.to_double()) > SAF_bound ||
|
||||
fabs(rwt.to_double()) > SAF_bound ||
|
||||
fabs(shx.to_double()) > SAF_bound ||
|
||||
fabs(shy.to_double()) > SAF_bound ||
|
||||
fabs(shz.to_double()) > SAF_bound ||
|
||||
fabs(shw.to_double()) > SAF_bound ||
|
||||
fabs(swt.to_double()) > SAF_bound ||
|
||||
fabs(thx.to_double()) > SAF_bound ||
|
||||
fabs(thy.to_double()) > SAF_bound ||
|
||||
fabs(thz.to_double()) > SAF_bound ||
|
||||
fabs(thw.to_double()) > SAF_bound ||
|
||||
fabs(twt.to_double()) > SAF_bound)
|
||||
{
|
||||
re_adjust:
|
||||
// Compute the new bound.
|
||||
double NEW_bound = 0.0;
|
||||
NEW_bound = std::max(NEW_bound, fabs(phx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(phy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(phz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(phw.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pwt.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qhx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qhy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qhz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qhw.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qwt.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rhx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rhy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rhz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rhw.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rwt.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(shx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(shy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(shz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(shw.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(swt.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(thx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(thy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(thz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(thw.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(twt.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(phx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(phy.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(phz.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(phw.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pwt.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qhx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qhy.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qhz.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qhw.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qwt.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rhx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rhy.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rhz.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rhw.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rwt.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(shx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(shy.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(shz.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(shw.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(swt.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(thx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(thy.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(thz.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(thw.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(twt.to_double()));
|
||||
// Re-adjust the context.
|
||||
Static_Filtered_power_testH3_25::new_bound(NEW_bound);
|
||||
}
|
||||
|
|
@ -374,31 +374,31 @@ re_adjust:
|
|||
try
|
||||
{
|
||||
return Static_Filtered_power_testH3_25::epsilon_variant(
|
||||
phx.to_double(),
|
||||
phy.to_double(),
|
||||
phz.to_double(),
|
||||
phw.to_double(),
|
||||
pwt.to_double(),
|
||||
qhx.to_double(),
|
||||
qhy.to_double(),
|
||||
qhz.to_double(),
|
||||
qhw.to_double(),
|
||||
qwt.to_double(),
|
||||
rhx.to_double(),
|
||||
rhy.to_double(),
|
||||
rhz.to_double(),
|
||||
rhw.to_double(),
|
||||
rwt.to_double(),
|
||||
shx.to_double(),
|
||||
shy.to_double(),
|
||||
shz.to_double(),
|
||||
shw.to_double(),
|
||||
swt.to_double(),
|
||||
thx.to_double(),
|
||||
thy.to_double(),
|
||||
thz.to_double(),
|
||||
thw.to_double(),
|
||||
twt.to_double(),
|
||||
phx.dbl(),
|
||||
phy.dbl(),
|
||||
phz.dbl(),
|
||||
phw.dbl(),
|
||||
pwt.dbl(),
|
||||
qhx.dbl(),
|
||||
qhy.dbl(),
|
||||
qhz.dbl(),
|
||||
qhw.dbl(),
|
||||
qwt.dbl(),
|
||||
rhx.dbl(),
|
||||
rhy.dbl(),
|
||||
rhz.dbl(),
|
||||
rhw.dbl(),
|
||||
rwt.dbl(),
|
||||
shx.dbl(),
|
||||
shy.dbl(),
|
||||
shz.dbl(),
|
||||
shw.dbl(),
|
||||
swt.dbl(),
|
||||
thx.dbl(),
|
||||
thy.dbl(),
|
||||
thz.dbl(),
|
||||
thw.dbl(),
|
||||
twt.dbl(),
|
||||
SAF_epsilon_0);
|
||||
}
|
||||
catch (Restricted_double::unsafe_comparison)
|
||||
|
|
|
|||
|
|
@ -163,28 +163,28 @@ in_smallest_orthogonalcircleC2(
|
|||
|
||||
// Check the bounds. All arguments must be <= SAF_bound.
|
||||
if (
|
||||
fabs(px.value()) > SAF_bound ||
|
||||
fabs(py.value()) > SAF_bound ||
|
||||
fabs(pw.value()) > SAF_bound ||
|
||||
fabs(qx.value()) > SAF_bound ||
|
||||
fabs(qy.value()) > SAF_bound ||
|
||||
fabs(qw.value()) > SAF_bound ||
|
||||
fabs(tx.value()) > SAF_bound ||
|
||||
fabs(ty.value()) > SAF_bound ||
|
||||
fabs(tw.value()) > SAF_bound)
|
||||
fabs(px.to_double()) > SAF_bound ||
|
||||
fabs(py.to_double()) > SAF_bound ||
|
||||
fabs(pw.to_double()) > SAF_bound ||
|
||||
fabs(qx.to_double()) > SAF_bound ||
|
||||
fabs(qy.to_double()) > SAF_bound ||
|
||||
fabs(qw.to_double()) > SAF_bound ||
|
||||
fabs(tx.to_double()) > SAF_bound ||
|
||||
fabs(ty.to_double()) > SAF_bound ||
|
||||
fabs(tw.to_double()) > SAF_bound)
|
||||
{
|
||||
re_adjust:
|
||||
// Compute the new bound.
|
||||
double NEW_bound = 0.0;
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pw.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qw.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(tx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ty.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(tw.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pw.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qw.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(tx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ty.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(tw.to_double()));
|
||||
// Re-adjust the context.
|
||||
Static_Filtered_in_smallest_orthogonalcircleC2_9::new_bound(NEW_bound);
|
||||
}
|
||||
|
|
@ -192,15 +192,15 @@ re_adjust:
|
|||
try
|
||||
{
|
||||
return Static_Filtered_in_smallest_orthogonalcircleC2_9::epsilon_variant(
|
||||
px.to_double(),
|
||||
py.to_double(),
|
||||
pw.to_double(),
|
||||
qx.to_double(),
|
||||
qy.to_double(),
|
||||
qw.to_double(),
|
||||
tx.to_double(),
|
||||
ty.to_double(),
|
||||
tw.to_double(),
|
||||
px.dbl(),
|
||||
py.dbl(),
|
||||
pw.dbl(),
|
||||
qx.dbl(),
|
||||
qy.dbl(),
|
||||
qw.dbl(),
|
||||
tx.dbl(),
|
||||
ty.dbl(),
|
||||
tw.dbl(),
|
||||
SAF_epsilon_0);
|
||||
}
|
||||
catch (Restricted_double::unsafe_comparison)
|
||||
|
|
|
|||
|
|
@ -115,18 +115,18 @@ sign_of_determinant2x2(
|
|||
|
||||
// Check the bounds. All arguments must be <= SAF_bound.
|
||||
if (
|
||||
fabs(a00.value()) > SAF_bound ||
|
||||
fabs(a01.value()) > SAF_bound ||
|
||||
fabs(a10.value()) > SAF_bound ||
|
||||
fabs(a11.value()) > SAF_bound)
|
||||
fabs(a00.to_double()) > SAF_bound ||
|
||||
fabs(a01.to_double()) > SAF_bound ||
|
||||
fabs(a10.to_double()) > SAF_bound ||
|
||||
fabs(a11.to_double()) > SAF_bound)
|
||||
{
|
||||
re_adjust:
|
||||
// Compute the new bound.
|
||||
double NEW_bound = 0.0;
|
||||
NEW_bound = std::max(NEW_bound, fabs(a00.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a01.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a10.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a11.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a00.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a01.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a10.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a11.to_double()));
|
||||
// Re-adjust the context.
|
||||
Static_Filtered_sign_of_determinant2x2_4::new_bound(NEW_bound);
|
||||
}
|
||||
|
|
@ -134,10 +134,10 @@ re_adjust:
|
|||
try
|
||||
{
|
||||
return Static_Filtered_sign_of_determinant2x2_4::epsilon_variant(
|
||||
a00.to_double(),
|
||||
a01.to_double(),
|
||||
a10.to_double(),
|
||||
a11.to_double(),
|
||||
a00.dbl(),
|
||||
a01.dbl(),
|
||||
a10.dbl(),
|
||||
a11.dbl(),
|
||||
SAF_epsilon_0);
|
||||
}
|
||||
catch (Restricted_double::unsafe_comparison)
|
||||
|
|
@ -283,28 +283,28 @@ sign_of_determinant3x3(
|
|||
|
||||
// Check the bounds. All arguments must be <= SAF_bound.
|
||||
if (
|
||||
fabs(a00.value()) > SAF_bound ||
|
||||
fabs(a01.value()) > SAF_bound ||
|
||||
fabs(a02.value()) > SAF_bound ||
|
||||
fabs(a10.value()) > SAF_bound ||
|
||||
fabs(a11.value()) > SAF_bound ||
|
||||
fabs(a12.value()) > SAF_bound ||
|
||||
fabs(a20.value()) > SAF_bound ||
|
||||
fabs(a21.value()) > SAF_bound ||
|
||||
fabs(a22.value()) > SAF_bound)
|
||||
fabs(a00.to_double()) > SAF_bound ||
|
||||
fabs(a01.to_double()) > SAF_bound ||
|
||||
fabs(a02.to_double()) > SAF_bound ||
|
||||
fabs(a10.to_double()) > SAF_bound ||
|
||||
fabs(a11.to_double()) > SAF_bound ||
|
||||
fabs(a12.to_double()) > SAF_bound ||
|
||||
fabs(a20.to_double()) > SAF_bound ||
|
||||
fabs(a21.to_double()) > SAF_bound ||
|
||||
fabs(a22.to_double()) > SAF_bound)
|
||||
{
|
||||
re_adjust:
|
||||
// Compute the new bound.
|
||||
double NEW_bound = 0.0;
|
||||
NEW_bound = std::max(NEW_bound, fabs(a00.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a01.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a02.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a10.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a11.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a12.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a20.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a21.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a22.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a00.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a01.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a02.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a10.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a11.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a12.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a20.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a21.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a22.to_double()));
|
||||
// Re-adjust the context.
|
||||
Static_Filtered_sign_of_determinant3x3_9::new_bound(NEW_bound);
|
||||
}
|
||||
|
|
@ -312,15 +312,15 @@ re_adjust:
|
|||
try
|
||||
{
|
||||
return Static_Filtered_sign_of_determinant3x3_9::epsilon_variant(
|
||||
a00.to_double(),
|
||||
a01.to_double(),
|
||||
a02.to_double(),
|
||||
a10.to_double(),
|
||||
a11.to_double(),
|
||||
a12.to_double(),
|
||||
a20.to_double(),
|
||||
a21.to_double(),
|
||||
a22.to_double(),
|
||||
a00.dbl(),
|
||||
a01.dbl(),
|
||||
a02.dbl(),
|
||||
a10.dbl(),
|
||||
a11.dbl(),
|
||||
a12.dbl(),
|
||||
a20.dbl(),
|
||||
a21.dbl(),
|
||||
a22.dbl(),
|
||||
SAF_epsilon_0);
|
||||
}
|
||||
catch (Restricted_double::unsafe_comparison)
|
||||
|
|
@ -515,42 +515,42 @@ sign_of_determinant4x4(
|
|||
|
||||
// Check the bounds. All arguments must be <= SAF_bound.
|
||||
if (
|
||||
fabs(a00.value()) > SAF_bound ||
|
||||
fabs(a01.value()) > SAF_bound ||
|
||||
fabs(a02.value()) > SAF_bound ||
|
||||
fabs(a03.value()) > SAF_bound ||
|
||||
fabs(a10.value()) > SAF_bound ||
|
||||
fabs(a11.value()) > SAF_bound ||
|
||||
fabs(a12.value()) > SAF_bound ||
|
||||
fabs(a13.value()) > SAF_bound ||
|
||||
fabs(a20.value()) > SAF_bound ||
|
||||
fabs(a21.value()) > SAF_bound ||
|
||||
fabs(a22.value()) > SAF_bound ||
|
||||
fabs(a23.value()) > SAF_bound ||
|
||||
fabs(a30.value()) > SAF_bound ||
|
||||
fabs(a31.value()) > SAF_bound ||
|
||||
fabs(a32.value()) > SAF_bound ||
|
||||
fabs(a33.value()) > SAF_bound)
|
||||
fabs(a00.to_double()) > SAF_bound ||
|
||||
fabs(a01.to_double()) > SAF_bound ||
|
||||
fabs(a02.to_double()) > SAF_bound ||
|
||||
fabs(a03.to_double()) > SAF_bound ||
|
||||
fabs(a10.to_double()) > SAF_bound ||
|
||||
fabs(a11.to_double()) > SAF_bound ||
|
||||
fabs(a12.to_double()) > SAF_bound ||
|
||||
fabs(a13.to_double()) > SAF_bound ||
|
||||
fabs(a20.to_double()) > SAF_bound ||
|
||||
fabs(a21.to_double()) > SAF_bound ||
|
||||
fabs(a22.to_double()) > SAF_bound ||
|
||||
fabs(a23.to_double()) > SAF_bound ||
|
||||
fabs(a30.to_double()) > SAF_bound ||
|
||||
fabs(a31.to_double()) > SAF_bound ||
|
||||
fabs(a32.to_double()) > SAF_bound ||
|
||||
fabs(a33.to_double()) > SAF_bound)
|
||||
{
|
||||
re_adjust:
|
||||
// Compute the new bound.
|
||||
double NEW_bound = 0.0;
|
||||
NEW_bound = std::max(NEW_bound, fabs(a00.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a01.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a02.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a03.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a10.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a11.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a12.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a13.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a20.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a21.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a22.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a23.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a30.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a31.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a32.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a33.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a00.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a01.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a02.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a03.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a10.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a11.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a12.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a13.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a20.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a21.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a22.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a23.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a30.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a31.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a32.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a33.to_double()));
|
||||
// Re-adjust the context.
|
||||
Static_Filtered_sign_of_determinant4x4_16::new_bound(NEW_bound);
|
||||
}
|
||||
|
|
@ -558,22 +558,22 @@ re_adjust:
|
|||
try
|
||||
{
|
||||
return Static_Filtered_sign_of_determinant4x4_16::epsilon_variant(
|
||||
a00.to_double(),
|
||||
a01.to_double(),
|
||||
a02.to_double(),
|
||||
a03.to_double(),
|
||||
a10.to_double(),
|
||||
a11.to_double(),
|
||||
a12.to_double(),
|
||||
a13.to_double(),
|
||||
a20.to_double(),
|
||||
a21.to_double(),
|
||||
a22.to_double(),
|
||||
a23.to_double(),
|
||||
a30.to_double(),
|
||||
a31.to_double(),
|
||||
a32.to_double(),
|
||||
a33.to_double(),
|
||||
a00.dbl(),
|
||||
a01.dbl(),
|
||||
a02.dbl(),
|
||||
a03.dbl(),
|
||||
a10.dbl(),
|
||||
a11.dbl(),
|
||||
a12.dbl(),
|
||||
a13.dbl(),
|
||||
a20.dbl(),
|
||||
a21.dbl(),
|
||||
a22.dbl(),
|
||||
a23.dbl(),
|
||||
a30.dbl(),
|
||||
a31.dbl(),
|
||||
a32.dbl(),
|
||||
a33.dbl(),
|
||||
SAF_epsilon_0);
|
||||
}
|
||||
catch (Restricted_double::unsafe_comparison)
|
||||
|
|
@ -831,60 +831,60 @@ sign_of_determinant5x5(
|
|||
|
||||
// Check the bounds. All arguments must be <= SAF_bound.
|
||||
if (
|
||||
fabs(a00.value()) > SAF_bound ||
|
||||
fabs(a01.value()) > SAF_bound ||
|
||||
fabs(a02.value()) > SAF_bound ||
|
||||
fabs(a03.value()) > SAF_bound ||
|
||||
fabs(a04.value()) > SAF_bound ||
|
||||
fabs(a10.value()) > SAF_bound ||
|
||||
fabs(a11.value()) > SAF_bound ||
|
||||
fabs(a12.value()) > SAF_bound ||
|
||||
fabs(a13.value()) > SAF_bound ||
|
||||
fabs(a14.value()) > SAF_bound ||
|
||||
fabs(a20.value()) > SAF_bound ||
|
||||
fabs(a21.value()) > SAF_bound ||
|
||||
fabs(a22.value()) > SAF_bound ||
|
||||
fabs(a23.value()) > SAF_bound ||
|
||||
fabs(a24.value()) > SAF_bound ||
|
||||
fabs(a30.value()) > SAF_bound ||
|
||||
fabs(a31.value()) > SAF_bound ||
|
||||
fabs(a32.value()) > SAF_bound ||
|
||||
fabs(a33.value()) > SAF_bound ||
|
||||
fabs(a34.value()) > SAF_bound ||
|
||||
fabs(a40.value()) > SAF_bound ||
|
||||
fabs(a41.value()) > SAF_bound ||
|
||||
fabs(a42.value()) > SAF_bound ||
|
||||
fabs(a43.value()) > SAF_bound ||
|
||||
fabs(a44.value()) > SAF_bound)
|
||||
fabs(a00.to_double()) > SAF_bound ||
|
||||
fabs(a01.to_double()) > SAF_bound ||
|
||||
fabs(a02.to_double()) > SAF_bound ||
|
||||
fabs(a03.to_double()) > SAF_bound ||
|
||||
fabs(a04.to_double()) > SAF_bound ||
|
||||
fabs(a10.to_double()) > SAF_bound ||
|
||||
fabs(a11.to_double()) > SAF_bound ||
|
||||
fabs(a12.to_double()) > SAF_bound ||
|
||||
fabs(a13.to_double()) > SAF_bound ||
|
||||
fabs(a14.to_double()) > SAF_bound ||
|
||||
fabs(a20.to_double()) > SAF_bound ||
|
||||
fabs(a21.to_double()) > SAF_bound ||
|
||||
fabs(a22.to_double()) > SAF_bound ||
|
||||
fabs(a23.to_double()) > SAF_bound ||
|
||||
fabs(a24.to_double()) > SAF_bound ||
|
||||
fabs(a30.to_double()) > SAF_bound ||
|
||||
fabs(a31.to_double()) > SAF_bound ||
|
||||
fabs(a32.to_double()) > SAF_bound ||
|
||||
fabs(a33.to_double()) > SAF_bound ||
|
||||
fabs(a34.to_double()) > SAF_bound ||
|
||||
fabs(a40.to_double()) > SAF_bound ||
|
||||
fabs(a41.to_double()) > SAF_bound ||
|
||||
fabs(a42.to_double()) > SAF_bound ||
|
||||
fabs(a43.to_double()) > SAF_bound ||
|
||||
fabs(a44.to_double()) > SAF_bound)
|
||||
{
|
||||
re_adjust:
|
||||
// Compute the new bound.
|
||||
double NEW_bound = 0.0;
|
||||
NEW_bound = std::max(NEW_bound, fabs(a00.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a01.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a02.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a03.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a04.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a10.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a11.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a12.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a13.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a14.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a20.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a21.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a22.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a23.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a24.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a30.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a31.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a32.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a33.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a34.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a40.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a41.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a42.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a43.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a44.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a00.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a01.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a02.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a03.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a04.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a10.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a11.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a12.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a13.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a14.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a20.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a21.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a22.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a23.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a24.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a30.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a31.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a32.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a33.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a34.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a40.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a41.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a42.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a43.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a44.to_double()));
|
||||
// Re-adjust the context.
|
||||
Static_Filtered_sign_of_determinant5x5_25::new_bound(NEW_bound);
|
||||
}
|
||||
|
|
@ -892,31 +892,31 @@ re_adjust:
|
|||
try
|
||||
{
|
||||
return Static_Filtered_sign_of_determinant5x5_25::epsilon_variant(
|
||||
a00.to_double(),
|
||||
a01.to_double(),
|
||||
a02.to_double(),
|
||||
a03.to_double(),
|
||||
a04.to_double(),
|
||||
a10.to_double(),
|
||||
a11.to_double(),
|
||||
a12.to_double(),
|
||||
a13.to_double(),
|
||||
a14.to_double(),
|
||||
a20.to_double(),
|
||||
a21.to_double(),
|
||||
a22.to_double(),
|
||||
a23.to_double(),
|
||||
a24.to_double(),
|
||||
a30.to_double(),
|
||||
a31.to_double(),
|
||||
a32.to_double(),
|
||||
a33.to_double(),
|
||||
a34.to_double(),
|
||||
a40.to_double(),
|
||||
a41.to_double(),
|
||||
a42.to_double(),
|
||||
a43.to_double(),
|
||||
a44.to_double(),
|
||||
a00.dbl(),
|
||||
a01.dbl(),
|
||||
a02.dbl(),
|
||||
a03.dbl(),
|
||||
a04.dbl(),
|
||||
a10.dbl(),
|
||||
a11.dbl(),
|
||||
a12.dbl(),
|
||||
a13.dbl(),
|
||||
a14.dbl(),
|
||||
a20.dbl(),
|
||||
a21.dbl(),
|
||||
a22.dbl(),
|
||||
a23.dbl(),
|
||||
a24.dbl(),
|
||||
a30.dbl(),
|
||||
a31.dbl(),
|
||||
a32.dbl(),
|
||||
a33.dbl(),
|
||||
a34.dbl(),
|
||||
a40.dbl(),
|
||||
a41.dbl(),
|
||||
a42.dbl(),
|
||||
a43.dbl(),
|
||||
a44.dbl(),
|
||||
SAF_epsilon_0);
|
||||
}
|
||||
catch (Restricted_double::unsafe_comparison)
|
||||
|
|
@ -1251,82 +1251,82 @@ sign_of_determinant6x6(
|
|||
|
||||
// Check the bounds. All arguments must be <= SAF_bound.
|
||||
if (
|
||||
fabs(a00.value()) > SAF_bound ||
|
||||
fabs(a01.value()) > SAF_bound ||
|
||||
fabs(a02.value()) > SAF_bound ||
|
||||
fabs(a03.value()) > SAF_bound ||
|
||||
fabs(a04.value()) > SAF_bound ||
|
||||
fabs(a05.value()) > SAF_bound ||
|
||||
fabs(a10.value()) > SAF_bound ||
|
||||
fabs(a11.value()) > SAF_bound ||
|
||||
fabs(a12.value()) > SAF_bound ||
|
||||
fabs(a13.value()) > SAF_bound ||
|
||||
fabs(a14.value()) > SAF_bound ||
|
||||
fabs(a15.value()) > SAF_bound ||
|
||||
fabs(a20.value()) > SAF_bound ||
|
||||
fabs(a21.value()) > SAF_bound ||
|
||||
fabs(a22.value()) > SAF_bound ||
|
||||
fabs(a23.value()) > SAF_bound ||
|
||||
fabs(a24.value()) > SAF_bound ||
|
||||
fabs(a25.value()) > SAF_bound ||
|
||||
fabs(a30.value()) > SAF_bound ||
|
||||
fabs(a31.value()) > SAF_bound ||
|
||||
fabs(a32.value()) > SAF_bound ||
|
||||
fabs(a33.value()) > SAF_bound ||
|
||||
fabs(a34.value()) > SAF_bound ||
|
||||
fabs(a35.value()) > SAF_bound ||
|
||||
fabs(a40.value()) > SAF_bound ||
|
||||
fabs(a41.value()) > SAF_bound ||
|
||||
fabs(a42.value()) > SAF_bound ||
|
||||
fabs(a43.value()) > SAF_bound ||
|
||||
fabs(a44.value()) > SAF_bound ||
|
||||
fabs(a45.value()) > SAF_bound ||
|
||||
fabs(a50.value()) > SAF_bound ||
|
||||
fabs(a51.value()) > SAF_bound ||
|
||||
fabs(a52.value()) > SAF_bound ||
|
||||
fabs(a53.value()) > SAF_bound ||
|
||||
fabs(a54.value()) > SAF_bound ||
|
||||
fabs(a55.value()) > SAF_bound)
|
||||
fabs(a00.to_double()) > SAF_bound ||
|
||||
fabs(a01.to_double()) > SAF_bound ||
|
||||
fabs(a02.to_double()) > SAF_bound ||
|
||||
fabs(a03.to_double()) > SAF_bound ||
|
||||
fabs(a04.to_double()) > SAF_bound ||
|
||||
fabs(a05.to_double()) > SAF_bound ||
|
||||
fabs(a10.to_double()) > SAF_bound ||
|
||||
fabs(a11.to_double()) > SAF_bound ||
|
||||
fabs(a12.to_double()) > SAF_bound ||
|
||||
fabs(a13.to_double()) > SAF_bound ||
|
||||
fabs(a14.to_double()) > SAF_bound ||
|
||||
fabs(a15.to_double()) > SAF_bound ||
|
||||
fabs(a20.to_double()) > SAF_bound ||
|
||||
fabs(a21.to_double()) > SAF_bound ||
|
||||
fabs(a22.to_double()) > SAF_bound ||
|
||||
fabs(a23.to_double()) > SAF_bound ||
|
||||
fabs(a24.to_double()) > SAF_bound ||
|
||||
fabs(a25.to_double()) > SAF_bound ||
|
||||
fabs(a30.to_double()) > SAF_bound ||
|
||||
fabs(a31.to_double()) > SAF_bound ||
|
||||
fabs(a32.to_double()) > SAF_bound ||
|
||||
fabs(a33.to_double()) > SAF_bound ||
|
||||
fabs(a34.to_double()) > SAF_bound ||
|
||||
fabs(a35.to_double()) > SAF_bound ||
|
||||
fabs(a40.to_double()) > SAF_bound ||
|
||||
fabs(a41.to_double()) > SAF_bound ||
|
||||
fabs(a42.to_double()) > SAF_bound ||
|
||||
fabs(a43.to_double()) > SAF_bound ||
|
||||
fabs(a44.to_double()) > SAF_bound ||
|
||||
fabs(a45.to_double()) > SAF_bound ||
|
||||
fabs(a50.to_double()) > SAF_bound ||
|
||||
fabs(a51.to_double()) > SAF_bound ||
|
||||
fabs(a52.to_double()) > SAF_bound ||
|
||||
fabs(a53.to_double()) > SAF_bound ||
|
||||
fabs(a54.to_double()) > SAF_bound ||
|
||||
fabs(a55.to_double()) > SAF_bound)
|
||||
{
|
||||
re_adjust:
|
||||
// Compute the new bound.
|
||||
double NEW_bound = 0.0;
|
||||
NEW_bound = std::max(NEW_bound, fabs(a00.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a01.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a02.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a03.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a04.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a05.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a10.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a11.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a12.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a13.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a14.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a15.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a20.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a21.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a22.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a23.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a24.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a25.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a30.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a31.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a32.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a33.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a34.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a35.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a40.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a41.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a42.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a43.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a44.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a45.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a50.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a51.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a52.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a53.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a54.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a55.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a00.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a01.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a02.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a03.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a04.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a05.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a10.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a11.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a12.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a13.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a14.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a15.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a20.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a21.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a22.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a23.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a24.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a25.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a30.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a31.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a32.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a33.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a34.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a35.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a40.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a41.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a42.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a43.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a44.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a45.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a50.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a51.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a52.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a53.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a54.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(a55.to_double()));
|
||||
// Re-adjust the context.
|
||||
Static_Filtered_sign_of_determinant6x6_36::new_bound(NEW_bound);
|
||||
}
|
||||
|
|
@ -1334,42 +1334,42 @@ re_adjust:
|
|||
try
|
||||
{
|
||||
return Static_Filtered_sign_of_determinant6x6_36::epsilon_variant(
|
||||
a00.to_double(),
|
||||
a01.to_double(),
|
||||
a02.to_double(),
|
||||
a03.to_double(),
|
||||
a04.to_double(),
|
||||
a05.to_double(),
|
||||
a10.to_double(),
|
||||
a11.to_double(),
|
||||
a12.to_double(),
|
||||
a13.to_double(),
|
||||
a14.to_double(),
|
||||
a15.to_double(),
|
||||
a20.to_double(),
|
||||
a21.to_double(),
|
||||
a22.to_double(),
|
||||
a23.to_double(),
|
||||
a24.to_double(),
|
||||
a25.to_double(),
|
||||
a30.to_double(),
|
||||
a31.to_double(),
|
||||
a32.to_double(),
|
||||
a33.to_double(),
|
||||
a34.to_double(),
|
||||
a35.to_double(),
|
||||
a40.to_double(),
|
||||
a41.to_double(),
|
||||
a42.to_double(),
|
||||
a43.to_double(),
|
||||
a44.to_double(),
|
||||
a45.to_double(),
|
||||
a50.to_double(),
|
||||
a51.to_double(),
|
||||
a52.to_double(),
|
||||
a53.to_double(),
|
||||
a54.to_double(),
|
||||
a55.to_double(),
|
||||
a00.dbl(),
|
||||
a01.dbl(),
|
||||
a02.dbl(),
|
||||
a03.dbl(),
|
||||
a04.dbl(),
|
||||
a05.dbl(),
|
||||
a10.dbl(),
|
||||
a11.dbl(),
|
||||
a12.dbl(),
|
||||
a13.dbl(),
|
||||
a14.dbl(),
|
||||
a15.dbl(),
|
||||
a20.dbl(),
|
||||
a21.dbl(),
|
||||
a22.dbl(),
|
||||
a23.dbl(),
|
||||
a24.dbl(),
|
||||
a25.dbl(),
|
||||
a30.dbl(),
|
||||
a31.dbl(),
|
||||
a32.dbl(),
|
||||
a33.dbl(),
|
||||
a34.dbl(),
|
||||
a35.dbl(),
|
||||
a40.dbl(),
|
||||
a41.dbl(),
|
||||
a42.dbl(),
|
||||
a43.dbl(),
|
||||
a44.dbl(),
|
||||
a45.dbl(),
|
||||
a50.dbl(),
|
||||
a51.dbl(),
|
||||
a52.dbl(),
|
||||
a53.dbl(),
|
||||
a54.dbl(),
|
||||
a55.dbl(),
|
||||
SAF_epsilon_0);
|
||||
}
|
||||
catch (Restricted_double::unsafe_comparison)
|
||||
|
|
|
|||
|
|
@ -151,24 +151,24 @@ compare_xC2(
|
|||
|
||||
// Check the bounds. All arguments must be <= SAF_bound.
|
||||
if (
|
||||
fabs(px.value()) > SAF_bound ||
|
||||
fabs(l1a.value()) > SAF_bound ||
|
||||
fabs(l1b.value()) > SAF_bound ||
|
||||
fabs(l1c.value()) > SAF_bound ||
|
||||
fabs(l2a.value()) > SAF_bound ||
|
||||
fabs(l2b.value()) > SAF_bound ||
|
||||
fabs(l2c.value()) > SAF_bound)
|
||||
fabs(px.to_double()) > SAF_bound ||
|
||||
fabs(l1a.to_double()) > SAF_bound ||
|
||||
fabs(l1b.to_double()) > SAF_bound ||
|
||||
fabs(l1c.to_double()) > SAF_bound ||
|
||||
fabs(l2a.to_double()) > SAF_bound ||
|
||||
fabs(l2b.to_double()) > SAF_bound ||
|
||||
fabs(l2c.to_double()) > SAF_bound)
|
||||
{
|
||||
re_adjust:
|
||||
// Compute the new bound.
|
||||
double NEW_bound = 0.0;
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l1a.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l1b.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l1c.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l2a.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l2b.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l2c.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l1a.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l1b.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l1c.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l2a.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l2b.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l2c.to_double()));
|
||||
// Re-adjust the context.
|
||||
Static_Filtered_compare_xC2_7::new_bound(NEW_bound);
|
||||
}
|
||||
|
|
@ -176,13 +176,13 @@ re_adjust:
|
|||
try
|
||||
{
|
||||
return Static_Filtered_compare_xC2_7::epsilon_variant(
|
||||
px.to_double(),
|
||||
l1a.to_double(),
|
||||
l1b.to_double(),
|
||||
l1c.to_double(),
|
||||
l2a.to_double(),
|
||||
l2b.to_double(),
|
||||
l2c.to_double(),
|
||||
px.dbl(),
|
||||
l1a.dbl(),
|
||||
l1b.dbl(),
|
||||
l1c.dbl(),
|
||||
l2a.dbl(),
|
||||
l2b.dbl(),
|
||||
l2c.dbl(),
|
||||
SAF_epsilon_0,
|
||||
SAF_epsilon_1);
|
||||
}
|
||||
|
|
@ -374,34 +374,34 @@ compare_xC2(
|
|||
|
||||
// Check the bounds. All arguments must be <= SAF_bound.
|
||||
if (
|
||||
fabs(l1a.value()) > SAF_bound ||
|
||||
fabs(l1b.value()) > SAF_bound ||
|
||||
fabs(l1c.value()) > SAF_bound ||
|
||||
fabs(l2a.value()) > SAF_bound ||
|
||||
fabs(l2b.value()) > SAF_bound ||
|
||||
fabs(l2c.value()) > SAF_bound ||
|
||||
fabs(h1a.value()) > SAF_bound ||
|
||||
fabs(h1b.value()) > SAF_bound ||
|
||||
fabs(h1c.value()) > SAF_bound ||
|
||||
fabs(h2a.value()) > SAF_bound ||
|
||||
fabs(h2b.value()) > SAF_bound ||
|
||||
fabs(h2c.value()) > SAF_bound)
|
||||
fabs(l1a.to_double()) > SAF_bound ||
|
||||
fabs(l1b.to_double()) > SAF_bound ||
|
||||
fabs(l1c.to_double()) > SAF_bound ||
|
||||
fabs(l2a.to_double()) > SAF_bound ||
|
||||
fabs(l2b.to_double()) > SAF_bound ||
|
||||
fabs(l2c.to_double()) > SAF_bound ||
|
||||
fabs(h1a.to_double()) > SAF_bound ||
|
||||
fabs(h1b.to_double()) > SAF_bound ||
|
||||
fabs(h1c.to_double()) > SAF_bound ||
|
||||
fabs(h2a.to_double()) > SAF_bound ||
|
||||
fabs(h2b.to_double()) > SAF_bound ||
|
||||
fabs(h2c.to_double()) > SAF_bound)
|
||||
{
|
||||
re_adjust:
|
||||
// Compute the new bound.
|
||||
double NEW_bound = 0.0;
|
||||
NEW_bound = std::max(NEW_bound, fabs(l1a.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l1b.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l1c.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l2a.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l2b.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l2c.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(h1a.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(h1b.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(h1c.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(h2a.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(h2b.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(h2c.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l1a.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l1b.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l1c.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l2a.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l2b.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l2c.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(h1a.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(h1b.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(h1c.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(h2a.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(h2b.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(h2c.to_double()));
|
||||
// Re-adjust the context.
|
||||
Static_Filtered_compare_xC2_12::new_bound(NEW_bound);
|
||||
}
|
||||
|
|
@ -409,18 +409,18 @@ re_adjust:
|
|||
try
|
||||
{
|
||||
return Static_Filtered_compare_xC2_12::epsilon_variant(
|
||||
l1a.to_double(),
|
||||
l1b.to_double(),
|
||||
l1c.to_double(),
|
||||
l2a.to_double(),
|
||||
l2b.to_double(),
|
||||
l2c.to_double(),
|
||||
h1a.to_double(),
|
||||
h1b.to_double(),
|
||||
h1c.to_double(),
|
||||
h2a.to_double(),
|
||||
h2b.to_double(),
|
||||
h2c.to_double(),
|
||||
l1a.dbl(),
|
||||
l1b.dbl(),
|
||||
l1c.dbl(),
|
||||
l2a.dbl(),
|
||||
l2b.dbl(),
|
||||
l2c.dbl(),
|
||||
h1a.dbl(),
|
||||
h1b.dbl(),
|
||||
h1c.dbl(),
|
||||
h2a.dbl(),
|
||||
h2b.dbl(),
|
||||
h2c.dbl(),
|
||||
SAF_epsilon_0,
|
||||
SAF_epsilon_1,
|
||||
SAF_epsilon_2);
|
||||
|
|
@ -558,20 +558,20 @@ compare_y_at_xC2(
|
|||
|
||||
// Check the bounds. All arguments must be <= SAF_bound.
|
||||
if (
|
||||
fabs(px.value()) > SAF_bound ||
|
||||
fabs(py.value()) > SAF_bound ||
|
||||
fabs(la.value()) > SAF_bound ||
|
||||
fabs(lb.value()) > SAF_bound ||
|
||||
fabs(lc.value()) > SAF_bound)
|
||||
fabs(px.to_double()) > SAF_bound ||
|
||||
fabs(py.to_double()) > SAF_bound ||
|
||||
fabs(la.to_double()) > SAF_bound ||
|
||||
fabs(lb.to_double()) > SAF_bound ||
|
||||
fabs(lc.to_double()) > SAF_bound)
|
||||
{
|
||||
re_adjust:
|
||||
// Compute the new bound.
|
||||
double NEW_bound = 0.0;
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(la.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(lb.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(lc.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(la.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(lb.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(lc.to_double()));
|
||||
// Re-adjust the context.
|
||||
Static_Filtered_compare_y_at_xC2_5::new_bound(NEW_bound);
|
||||
}
|
||||
|
|
@ -579,11 +579,11 @@ re_adjust:
|
|||
try
|
||||
{
|
||||
return Static_Filtered_compare_y_at_xC2_5::epsilon_variant(
|
||||
px.to_double(),
|
||||
py.to_double(),
|
||||
la.to_double(),
|
||||
lb.to_double(),
|
||||
lc.to_double(),
|
||||
px.dbl(),
|
||||
py.dbl(),
|
||||
la.dbl(),
|
||||
lb.dbl(),
|
||||
lc.dbl(),
|
||||
SAF_epsilon_0,
|
||||
SAF_epsilon_1);
|
||||
}
|
||||
|
|
@ -731,24 +731,24 @@ compare_y_at_xC2(
|
|||
|
||||
// Check the bounds. All arguments must be <= SAF_bound.
|
||||
if (
|
||||
fabs(px.value()) > SAF_bound ||
|
||||
fabs(l1a.value()) > SAF_bound ||
|
||||
fabs(l1b.value()) > SAF_bound ||
|
||||
fabs(l1c.value()) > SAF_bound ||
|
||||
fabs(l2a.value()) > SAF_bound ||
|
||||
fabs(l2b.value()) > SAF_bound ||
|
||||
fabs(l2c.value()) > SAF_bound)
|
||||
fabs(px.to_double()) > SAF_bound ||
|
||||
fabs(l1a.to_double()) > SAF_bound ||
|
||||
fabs(l1b.to_double()) > SAF_bound ||
|
||||
fabs(l1c.to_double()) > SAF_bound ||
|
||||
fabs(l2a.to_double()) > SAF_bound ||
|
||||
fabs(l2b.to_double()) > SAF_bound ||
|
||||
fabs(l2c.to_double()) > SAF_bound)
|
||||
{
|
||||
re_adjust:
|
||||
// Compute the new bound.
|
||||
double NEW_bound = 0.0;
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l1a.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l1b.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l1c.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l2a.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l2b.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l2c.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l1a.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l1b.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l1c.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l2a.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l2b.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l2c.to_double()));
|
||||
// Re-adjust the context.
|
||||
Static_Filtered_compare_y_at_xC2_7::new_bound(NEW_bound);
|
||||
}
|
||||
|
|
@ -756,13 +756,13 @@ re_adjust:
|
|||
try
|
||||
{
|
||||
return Static_Filtered_compare_y_at_xC2_7::epsilon_variant(
|
||||
px.to_double(),
|
||||
l1a.to_double(),
|
||||
l1b.to_double(),
|
||||
l1c.to_double(),
|
||||
l2a.to_double(),
|
||||
l2b.to_double(),
|
||||
l2c.to_double(),
|
||||
px.dbl(),
|
||||
l1a.dbl(),
|
||||
l1b.dbl(),
|
||||
l1c.dbl(),
|
||||
l2a.dbl(),
|
||||
l2b.dbl(),
|
||||
l2c.dbl(),
|
||||
SAF_epsilon_0,
|
||||
SAF_epsilon_1,
|
||||
SAF_epsilon_2);
|
||||
|
|
@ -927,28 +927,28 @@ compare_y_at_xC2(
|
|||
|
||||
// Check the bounds. All arguments must be <= SAF_bound.
|
||||
if (
|
||||
fabs(l1a.value()) > SAF_bound ||
|
||||
fabs(l1b.value()) > SAF_bound ||
|
||||
fabs(l1c.value()) > SAF_bound ||
|
||||
fabs(l2a.value()) > SAF_bound ||
|
||||
fabs(l2b.value()) > SAF_bound ||
|
||||
fabs(l2c.value()) > SAF_bound ||
|
||||
fabs(ha.value()) > SAF_bound ||
|
||||
fabs(hb.value()) > SAF_bound ||
|
||||
fabs(hc.value()) > SAF_bound)
|
||||
fabs(l1a.to_double()) > SAF_bound ||
|
||||
fabs(l1b.to_double()) > SAF_bound ||
|
||||
fabs(l1c.to_double()) > SAF_bound ||
|
||||
fabs(l2a.to_double()) > SAF_bound ||
|
||||
fabs(l2b.to_double()) > SAF_bound ||
|
||||
fabs(l2c.to_double()) > SAF_bound ||
|
||||
fabs(ha.to_double()) > SAF_bound ||
|
||||
fabs(hb.to_double()) > SAF_bound ||
|
||||
fabs(hc.to_double()) > SAF_bound)
|
||||
{
|
||||
re_adjust:
|
||||
// Compute the new bound.
|
||||
double NEW_bound = 0.0;
|
||||
NEW_bound = std::max(NEW_bound, fabs(l1a.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l1b.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l1c.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l2a.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l2b.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l2c.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ha.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(hb.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(hc.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l1a.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l1b.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l1c.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l2a.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l2b.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l2c.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ha.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(hb.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(hc.to_double()));
|
||||
// Re-adjust the context.
|
||||
Static_Filtered_compare_y_at_xC2_9::new_bound(NEW_bound);
|
||||
}
|
||||
|
|
@ -956,15 +956,15 @@ re_adjust:
|
|||
try
|
||||
{
|
||||
return Static_Filtered_compare_y_at_xC2_9::epsilon_variant(
|
||||
l1a.to_double(),
|
||||
l1b.to_double(),
|
||||
l1c.to_double(),
|
||||
l2a.to_double(),
|
||||
l2b.to_double(),
|
||||
l2c.to_double(),
|
||||
ha.to_double(),
|
||||
hb.to_double(),
|
||||
hc.to_double(),
|
||||
l1a.dbl(),
|
||||
l1b.dbl(),
|
||||
l1c.dbl(),
|
||||
l2a.dbl(),
|
||||
l2b.dbl(),
|
||||
l2c.dbl(),
|
||||
ha.dbl(),
|
||||
hb.dbl(),
|
||||
hc.dbl(),
|
||||
SAF_epsilon_0,
|
||||
SAF_epsilon_1,
|
||||
SAF_epsilon_2,
|
||||
|
|
@ -1160,34 +1160,34 @@ compare_y_at_xC2(
|
|||
|
||||
// Check the bounds. All arguments must be <= SAF_bound.
|
||||
if (
|
||||
fabs(l1a.value()) > SAF_bound ||
|
||||
fabs(l1b.value()) > SAF_bound ||
|
||||
fabs(l1c.value()) > SAF_bound ||
|
||||
fabs(l2a.value()) > SAF_bound ||
|
||||
fabs(l2b.value()) > SAF_bound ||
|
||||
fabs(l2c.value()) > SAF_bound ||
|
||||
fabs(h1a.value()) > SAF_bound ||
|
||||
fabs(h1b.value()) > SAF_bound ||
|
||||
fabs(h1c.value()) > SAF_bound ||
|
||||
fabs(h2a.value()) > SAF_bound ||
|
||||
fabs(h2b.value()) > SAF_bound ||
|
||||
fabs(h2c.value()) > SAF_bound)
|
||||
fabs(l1a.to_double()) > SAF_bound ||
|
||||
fabs(l1b.to_double()) > SAF_bound ||
|
||||
fabs(l1c.to_double()) > SAF_bound ||
|
||||
fabs(l2a.to_double()) > SAF_bound ||
|
||||
fabs(l2b.to_double()) > SAF_bound ||
|
||||
fabs(l2c.to_double()) > SAF_bound ||
|
||||
fabs(h1a.to_double()) > SAF_bound ||
|
||||
fabs(h1b.to_double()) > SAF_bound ||
|
||||
fabs(h1c.to_double()) > SAF_bound ||
|
||||
fabs(h2a.to_double()) > SAF_bound ||
|
||||
fabs(h2b.to_double()) > SAF_bound ||
|
||||
fabs(h2c.to_double()) > SAF_bound)
|
||||
{
|
||||
re_adjust:
|
||||
// Compute the new bound.
|
||||
double NEW_bound = 0.0;
|
||||
NEW_bound = std::max(NEW_bound, fabs(l1a.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l1b.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l1c.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l2a.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l2b.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l2c.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(h1a.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(h1b.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(h1c.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(h2a.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(h2b.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(h2c.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l1a.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l1b.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l1c.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l2a.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l2b.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(l2c.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(h1a.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(h1b.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(h1c.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(h2a.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(h2b.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(h2c.to_double()));
|
||||
// Re-adjust the context.
|
||||
Static_Filtered_compare_y_at_xC2_12::new_bound(NEW_bound);
|
||||
}
|
||||
|
|
@ -1195,18 +1195,18 @@ re_adjust:
|
|||
try
|
||||
{
|
||||
return Static_Filtered_compare_y_at_xC2_12::epsilon_variant(
|
||||
l1a.to_double(),
|
||||
l1b.to_double(),
|
||||
l1c.to_double(),
|
||||
l2a.to_double(),
|
||||
l2b.to_double(),
|
||||
l2c.to_double(),
|
||||
h1a.to_double(),
|
||||
h1b.to_double(),
|
||||
h1c.to_double(),
|
||||
h2a.to_double(),
|
||||
h2b.to_double(),
|
||||
h2c.to_double(),
|
||||
l1a.dbl(),
|
||||
l1b.dbl(),
|
||||
l1c.dbl(),
|
||||
l2a.dbl(),
|
||||
l2b.dbl(),
|
||||
l2c.dbl(),
|
||||
h1a.dbl(),
|
||||
h1b.dbl(),
|
||||
h1c.dbl(),
|
||||
h2a.dbl(),
|
||||
h2b.dbl(),
|
||||
h2c.dbl(),
|
||||
SAF_epsilon_0,
|
||||
SAF_epsilon_1,
|
||||
SAF_epsilon_2,
|
||||
|
|
@ -1329,18 +1329,18 @@ compare_deltax_deltayC2(
|
|||
|
||||
// Check the bounds. All arguments must be <= SAF_bound.
|
||||
if (
|
||||
fabs(px.value()) > SAF_bound ||
|
||||
fabs(qx.value()) > SAF_bound ||
|
||||
fabs(ry.value()) > SAF_bound ||
|
||||
fabs(sy.value()) > SAF_bound)
|
||||
fabs(px.to_double()) > SAF_bound ||
|
||||
fabs(qx.to_double()) > SAF_bound ||
|
||||
fabs(ry.to_double()) > SAF_bound ||
|
||||
fabs(sy.to_double()) > SAF_bound)
|
||||
{
|
||||
re_adjust:
|
||||
// Compute the new bound.
|
||||
double NEW_bound = 0.0;
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ry.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(sy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ry.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(sy.to_double()));
|
||||
// Re-adjust the context.
|
||||
Static_Filtered_compare_deltax_deltayC2_4::new_bound(NEW_bound);
|
||||
}
|
||||
|
|
@ -1348,10 +1348,10 @@ re_adjust:
|
|||
try
|
||||
{
|
||||
return Static_Filtered_compare_deltax_deltayC2_4::epsilon_variant(
|
||||
px.to_double(),
|
||||
qx.to_double(),
|
||||
ry.to_double(),
|
||||
sy.to_double(),
|
||||
px.dbl(),
|
||||
qx.dbl(),
|
||||
ry.dbl(),
|
||||
sy.dbl(),
|
||||
SAF_epsilon_0);
|
||||
}
|
||||
catch (Restricted_double::unsafe_comparison)
|
||||
|
|
@ -1475,22 +1475,22 @@ orientationC2(
|
|||
|
||||
// Check the bounds. All arguments must be <= SAF_bound.
|
||||
if (
|
||||
fabs(px.value()) > SAF_bound ||
|
||||
fabs(py.value()) > SAF_bound ||
|
||||
fabs(qx.value()) > SAF_bound ||
|
||||
fabs(qy.value()) > SAF_bound ||
|
||||
fabs(rx.value()) > SAF_bound ||
|
||||
fabs(ry.value()) > SAF_bound)
|
||||
fabs(px.to_double()) > SAF_bound ||
|
||||
fabs(py.to_double()) > SAF_bound ||
|
||||
fabs(qx.to_double()) > SAF_bound ||
|
||||
fabs(qy.to_double()) > SAF_bound ||
|
||||
fabs(rx.to_double()) > SAF_bound ||
|
||||
fabs(ry.to_double()) > SAF_bound)
|
||||
{
|
||||
re_adjust:
|
||||
// Compute the new bound.
|
||||
double NEW_bound = 0.0;
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ry.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ry.to_double()));
|
||||
// Re-adjust the context.
|
||||
Static_Filtered_orientationC2_6::new_bound(NEW_bound);
|
||||
}
|
||||
|
|
@ -1498,12 +1498,12 @@ re_adjust:
|
|||
try
|
||||
{
|
||||
return Static_Filtered_orientationC2_6::epsilon_variant(
|
||||
px.to_double(),
|
||||
py.to_double(),
|
||||
qx.to_double(),
|
||||
qy.to_double(),
|
||||
rx.to_double(),
|
||||
ry.to_double(),
|
||||
px.dbl(),
|
||||
py.dbl(),
|
||||
qx.dbl(),
|
||||
qy.dbl(),
|
||||
rx.dbl(),
|
||||
ry.dbl(),
|
||||
SAF_epsilon_0);
|
||||
}
|
||||
catch (Restricted_double::unsafe_comparison)
|
||||
|
|
@ -1671,26 +1671,26 @@ side_of_oriented_circleC2(
|
|||
|
||||
// Check the bounds. All arguments must be <= SAF_bound.
|
||||
if (
|
||||
fabs(px.value()) > SAF_bound ||
|
||||
fabs(py.value()) > SAF_bound ||
|
||||
fabs(qx.value()) > SAF_bound ||
|
||||
fabs(qy.value()) > SAF_bound ||
|
||||
fabs(rx.value()) > SAF_bound ||
|
||||
fabs(ry.value()) > SAF_bound ||
|
||||
fabs(tx.value()) > SAF_bound ||
|
||||
fabs(ty.value()) > SAF_bound)
|
||||
fabs(px.to_double()) > SAF_bound ||
|
||||
fabs(py.to_double()) > SAF_bound ||
|
||||
fabs(qx.to_double()) > SAF_bound ||
|
||||
fabs(qy.to_double()) > SAF_bound ||
|
||||
fabs(rx.to_double()) > SAF_bound ||
|
||||
fabs(ry.to_double()) > SAF_bound ||
|
||||
fabs(tx.to_double()) > SAF_bound ||
|
||||
fabs(ty.to_double()) > SAF_bound)
|
||||
{
|
||||
re_adjust:
|
||||
// Compute the new bound.
|
||||
double NEW_bound = 0.0;
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ry.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(tx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ty.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ry.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(tx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ty.to_double()));
|
||||
// Re-adjust the context.
|
||||
Static_Filtered_side_of_oriented_circleC2_8::new_bound(NEW_bound);
|
||||
}
|
||||
|
|
@ -1698,14 +1698,14 @@ re_adjust:
|
|||
try
|
||||
{
|
||||
return Static_Filtered_side_of_oriented_circleC2_8::epsilon_variant(
|
||||
px.to_double(),
|
||||
py.to_double(),
|
||||
qx.to_double(),
|
||||
qy.to_double(),
|
||||
rx.to_double(),
|
||||
ry.to_double(),
|
||||
tx.to_double(),
|
||||
ty.to_double(),
|
||||
px.dbl(),
|
||||
py.dbl(),
|
||||
qx.dbl(),
|
||||
qy.dbl(),
|
||||
rx.dbl(),
|
||||
ry.dbl(),
|
||||
tx.dbl(),
|
||||
ty.dbl(),
|
||||
SAF_epsilon_0);
|
||||
}
|
||||
catch (Restricted_double::unsafe_comparison)
|
||||
|
|
@ -1855,26 +1855,26 @@ side_of_bounded_circleC2(
|
|||
|
||||
// Check the bounds. All arguments must be <= SAF_bound.
|
||||
if (
|
||||
fabs(px.value()) > SAF_bound ||
|
||||
fabs(py.value()) > SAF_bound ||
|
||||
fabs(qx.value()) > SAF_bound ||
|
||||
fabs(qy.value()) > SAF_bound ||
|
||||
fabs(rx.value()) > SAF_bound ||
|
||||
fabs(ry.value()) > SAF_bound ||
|
||||
fabs(tx.value()) > SAF_bound ||
|
||||
fabs(ty.value()) > SAF_bound)
|
||||
fabs(px.to_double()) > SAF_bound ||
|
||||
fabs(py.to_double()) > SAF_bound ||
|
||||
fabs(qx.to_double()) > SAF_bound ||
|
||||
fabs(qy.to_double()) > SAF_bound ||
|
||||
fabs(rx.to_double()) > SAF_bound ||
|
||||
fabs(ry.to_double()) > SAF_bound ||
|
||||
fabs(tx.to_double()) > SAF_bound ||
|
||||
fabs(ty.to_double()) > SAF_bound)
|
||||
{
|
||||
re_adjust:
|
||||
// Compute the new bound.
|
||||
double NEW_bound = 0.0;
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ry.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(tx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ty.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ry.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(tx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ty.to_double()));
|
||||
// Re-adjust the context.
|
||||
Static_Filtered_side_of_bounded_circleC2_8::new_bound(NEW_bound);
|
||||
}
|
||||
|
|
@ -1882,14 +1882,14 @@ re_adjust:
|
|||
try
|
||||
{
|
||||
return Static_Filtered_side_of_bounded_circleC2_8::epsilon_variant(
|
||||
px.to_double(),
|
||||
py.to_double(),
|
||||
qx.to_double(),
|
||||
qy.to_double(),
|
||||
rx.to_double(),
|
||||
ry.to_double(),
|
||||
tx.to_double(),
|
||||
ty.to_double(),
|
||||
px.dbl(),
|
||||
py.dbl(),
|
||||
qx.dbl(),
|
||||
qy.dbl(),
|
||||
rx.dbl(),
|
||||
ry.dbl(),
|
||||
tx.dbl(),
|
||||
ty.dbl(),
|
||||
SAF_epsilon_0,
|
||||
SAF_epsilon_1);
|
||||
}
|
||||
|
|
@ -2020,22 +2020,22 @@ cmp_dist_to_pointC2(
|
|||
|
||||
// Check the bounds. All arguments must be <= SAF_bound.
|
||||
if (
|
||||
fabs(px.value()) > SAF_bound ||
|
||||
fabs(py.value()) > SAF_bound ||
|
||||
fabs(qx.value()) > SAF_bound ||
|
||||
fabs(qy.value()) > SAF_bound ||
|
||||
fabs(rx.value()) > SAF_bound ||
|
||||
fabs(ry.value()) > SAF_bound)
|
||||
fabs(px.to_double()) > SAF_bound ||
|
||||
fabs(py.to_double()) > SAF_bound ||
|
||||
fabs(qx.to_double()) > SAF_bound ||
|
||||
fabs(qy.to_double()) > SAF_bound ||
|
||||
fabs(rx.to_double()) > SAF_bound ||
|
||||
fabs(ry.to_double()) > SAF_bound)
|
||||
{
|
||||
re_adjust:
|
||||
// Compute the new bound.
|
||||
double NEW_bound = 0.0;
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ry.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ry.to_double()));
|
||||
// Re-adjust the context.
|
||||
Static_Filtered_cmp_dist_to_pointC2_6::new_bound(NEW_bound);
|
||||
}
|
||||
|
|
@ -2043,12 +2043,12 @@ re_adjust:
|
|||
try
|
||||
{
|
||||
return Static_Filtered_cmp_dist_to_pointC2_6::epsilon_variant(
|
||||
px.to_double(),
|
||||
py.to_double(),
|
||||
qx.to_double(),
|
||||
qy.to_double(),
|
||||
rx.to_double(),
|
||||
ry.to_double(),
|
||||
px.dbl(),
|
||||
py.dbl(),
|
||||
qx.dbl(),
|
||||
qy.dbl(),
|
||||
rx.dbl(),
|
||||
ry.dbl(),
|
||||
SAF_epsilon_0);
|
||||
}
|
||||
catch (Restricted_double::unsafe_comparison)
|
||||
|
|
@ -2182,24 +2182,24 @@ cmp_signed_dist_to_lineC2(
|
|||
|
||||
// Check the bounds. All arguments must be <= SAF_bound.
|
||||
if (
|
||||
fabs(la.value()) > SAF_bound ||
|
||||
fabs(lb.value()) > SAF_bound ||
|
||||
fabs(lc.value()) > SAF_bound ||
|
||||
fabs(px.value()) > SAF_bound ||
|
||||
fabs(py.value()) > SAF_bound ||
|
||||
fabs(qx.value()) > SAF_bound ||
|
||||
fabs(qy.value()) > SAF_bound)
|
||||
fabs(la.to_double()) > SAF_bound ||
|
||||
fabs(lb.to_double()) > SAF_bound ||
|
||||
fabs(lc.to_double()) > SAF_bound ||
|
||||
fabs(px.to_double()) > SAF_bound ||
|
||||
fabs(py.to_double()) > SAF_bound ||
|
||||
fabs(qx.to_double()) > SAF_bound ||
|
||||
fabs(qy.to_double()) > SAF_bound)
|
||||
{
|
||||
re_adjust:
|
||||
// Compute the new bound.
|
||||
double NEW_bound = 0.0;
|
||||
NEW_bound = std::max(NEW_bound, fabs(la.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(lb.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(lc.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(la.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(lb.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(lc.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.to_double()));
|
||||
// Re-adjust the context.
|
||||
Static_Filtered_cmp_signed_dist_to_lineC2_7::new_bound(NEW_bound);
|
||||
}
|
||||
|
|
@ -2207,13 +2207,13 @@ re_adjust:
|
|||
try
|
||||
{
|
||||
return Static_Filtered_cmp_signed_dist_to_lineC2_7::epsilon_variant(
|
||||
la.to_double(),
|
||||
lb.to_double(),
|
||||
lc.to_double(),
|
||||
px.to_double(),
|
||||
py.to_double(),
|
||||
qx.to_double(),
|
||||
qy.to_double(),
|
||||
la.dbl(),
|
||||
lb.dbl(),
|
||||
lc.dbl(),
|
||||
px.dbl(),
|
||||
py.dbl(),
|
||||
qx.dbl(),
|
||||
qy.dbl(),
|
||||
SAF_epsilon_0);
|
||||
}
|
||||
catch (Restricted_double::unsafe_comparison)
|
||||
|
|
@ -2354,26 +2354,26 @@ cmp_signed_dist_to_lineC2(
|
|||
|
||||
// Check the bounds. All arguments must be <= SAF_bound.
|
||||
if (
|
||||
fabs(px.value()) > SAF_bound ||
|
||||
fabs(py.value()) > SAF_bound ||
|
||||
fabs(qx.value()) > SAF_bound ||
|
||||
fabs(qy.value()) > SAF_bound ||
|
||||
fabs(rx.value()) > SAF_bound ||
|
||||
fabs(ry.value()) > SAF_bound ||
|
||||
fabs(sx.value()) > SAF_bound ||
|
||||
fabs(sy.value()) > SAF_bound)
|
||||
fabs(px.to_double()) > SAF_bound ||
|
||||
fabs(py.to_double()) > SAF_bound ||
|
||||
fabs(qx.to_double()) > SAF_bound ||
|
||||
fabs(qy.to_double()) > SAF_bound ||
|
||||
fabs(rx.to_double()) > SAF_bound ||
|
||||
fabs(ry.to_double()) > SAF_bound ||
|
||||
fabs(sx.to_double()) > SAF_bound ||
|
||||
fabs(sy.to_double()) > SAF_bound)
|
||||
{
|
||||
re_adjust:
|
||||
// Compute the new bound.
|
||||
double NEW_bound = 0.0;
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ry.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(sx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(sy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ry.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(sx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(sy.to_double()));
|
||||
// Re-adjust the context.
|
||||
Static_Filtered_cmp_signed_dist_to_lineC2_8::new_bound(NEW_bound);
|
||||
}
|
||||
|
|
@ -2381,14 +2381,14 @@ re_adjust:
|
|||
try
|
||||
{
|
||||
return Static_Filtered_cmp_signed_dist_to_lineC2_8::epsilon_variant(
|
||||
px.to_double(),
|
||||
py.to_double(),
|
||||
qx.to_double(),
|
||||
qy.to_double(),
|
||||
rx.to_double(),
|
||||
ry.to_double(),
|
||||
sx.to_double(),
|
||||
sy.to_double(),
|
||||
px.dbl(),
|
||||
py.dbl(),
|
||||
qx.dbl(),
|
||||
qy.dbl(),
|
||||
rx.dbl(),
|
||||
ry.dbl(),
|
||||
sx.dbl(),
|
||||
sy.dbl(),
|
||||
SAF_epsilon_0);
|
||||
}
|
||||
catch (Restricted_double::unsafe_comparison)
|
||||
|
|
|
|||
|
|
@ -173,28 +173,28 @@ collinearC3(
|
|||
|
||||
// Check the bounds. All arguments must be <= SAF_bound.
|
||||
if (
|
||||
fabs(px.value()) > SAF_bound ||
|
||||
fabs(py.value()) > SAF_bound ||
|
||||
fabs(pz.value()) > SAF_bound ||
|
||||
fabs(qx.value()) > SAF_bound ||
|
||||
fabs(qy.value()) > SAF_bound ||
|
||||
fabs(qz.value()) > SAF_bound ||
|
||||
fabs(rx.value()) > SAF_bound ||
|
||||
fabs(ry.value()) > SAF_bound ||
|
||||
fabs(rz.value()) > SAF_bound)
|
||||
fabs(px.to_double()) > SAF_bound ||
|
||||
fabs(py.to_double()) > SAF_bound ||
|
||||
fabs(pz.to_double()) > SAF_bound ||
|
||||
fabs(qx.to_double()) > SAF_bound ||
|
||||
fabs(qy.to_double()) > SAF_bound ||
|
||||
fabs(qz.to_double()) > SAF_bound ||
|
||||
fabs(rx.to_double()) > SAF_bound ||
|
||||
fabs(ry.to_double()) > SAF_bound ||
|
||||
fabs(rz.to_double()) > SAF_bound)
|
||||
{
|
||||
re_adjust:
|
||||
// Compute the new bound.
|
||||
double NEW_bound = 0.0;
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ry.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pz.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qz.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ry.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rz.to_double()));
|
||||
// Re-adjust the context.
|
||||
Static_Filtered_collinearC3_9::new_bound(NEW_bound);
|
||||
}
|
||||
|
|
@ -202,15 +202,15 @@ re_adjust:
|
|||
try
|
||||
{
|
||||
return Static_Filtered_collinearC3_9::epsilon_variant(
|
||||
px.to_double(),
|
||||
py.to_double(),
|
||||
pz.to_double(),
|
||||
qx.to_double(),
|
||||
qy.to_double(),
|
||||
qz.to_double(),
|
||||
rx.to_double(),
|
||||
ry.to_double(),
|
||||
rz.to_double(),
|
||||
px.dbl(),
|
||||
py.dbl(),
|
||||
pz.dbl(),
|
||||
qx.dbl(),
|
||||
qy.dbl(),
|
||||
qz.dbl(),
|
||||
rx.dbl(),
|
||||
ry.dbl(),
|
||||
rz.dbl(),
|
||||
SAF_epsilon_0,
|
||||
SAF_epsilon_1,
|
||||
SAF_epsilon_2);
|
||||
|
|
@ -381,34 +381,34 @@ orientationC3(
|
|||
|
||||
// Check the bounds. All arguments must be <= SAF_bound.
|
||||
if (
|
||||
fabs(px.value()) > SAF_bound ||
|
||||
fabs(py.value()) > SAF_bound ||
|
||||
fabs(pz.value()) > SAF_bound ||
|
||||
fabs(qx.value()) > SAF_bound ||
|
||||
fabs(qy.value()) > SAF_bound ||
|
||||
fabs(qz.value()) > SAF_bound ||
|
||||
fabs(rx.value()) > SAF_bound ||
|
||||
fabs(ry.value()) > SAF_bound ||
|
||||
fabs(rz.value()) > SAF_bound ||
|
||||
fabs(sx.value()) > SAF_bound ||
|
||||
fabs(sy.value()) > SAF_bound ||
|
||||
fabs(sz.value()) > SAF_bound)
|
||||
fabs(px.to_double()) > SAF_bound ||
|
||||
fabs(py.to_double()) > SAF_bound ||
|
||||
fabs(pz.to_double()) > SAF_bound ||
|
||||
fabs(qx.to_double()) > SAF_bound ||
|
||||
fabs(qy.to_double()) > SAF_bound ||
|
||||
fabs(qz.to_double()) > SAF_bound ||
|
||||
fabs(rx.to_double()) > SAF_bound ||
|
||||
fabs(ry.to_double()) > SAF_bound ||
|
||||
fabs(rz.to_double()) > SAF_bound ||
|
||||
fabs(sx.to_double()) > SAF_bound ||
|
||||
fabs(sy.to_double()) > SAF_bound ||
|
||||
fabs(sz.to_double()) > SAF_bound)
|
||||
{
|
||||
re_adjust:
|
||||
// Compute the new bound.
|
||||
double NEW_bound = 0.0;
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ry.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(sx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(sy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(sz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pz.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qz.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ry.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rz.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(sx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(sy.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(sz.to_double()));
|
||||
// Re-adjust the context.
|
||||
Static_Filtered_orientationC3_12::new_bound(NEW_bound);
|
||||
}
|
||||
|
|
@ -416,18 +416,18 @@ re_adjust:
|
|||
try
|
||||
{
|
||||
return Static_Filtered_orientationC3_12::epsilon_variant(
|
||||
px.to_double(),
|
||||
py.to_double(),
|
||||
pz.to_double(),
|
||||
qx.to_double(),
|
||||
qy.to_double(),
|
||||
qz.to_double(),
|
||||
rx.to_double(),
|
||||
ry.to_double(),
|
||||
rz.to_double(),
|
||||
sx.to_double(),
|
||||
sy.to_double(),
|
||||
sz.to_double(),
|
||||
px.dbl(),
|
||||
py.dbl(),
|
||||
pz.dbl(),
|
||||
qx.dbl(),
|
||||
qy.dbl(),
|
||||
qz.dbl(),
|
||||
rx.dbl(),
|
||||
ry.dbl(),
|
||||
rz.dbl(),
|
||||
sx.dbl(),
|
||||
sy.dbl(),
|
||||
sz.dbl(),
|
||||
SAF_epsilon_0);
|
||||
}
|
||||
catch (Restricted_double::unsafe_comparison)
|
||||
|
|
@ -651,40 +651,40 @@ side_of_oriented_sphereC3(
|
|||
|
||||
// Check the bounds. All arguments must be <= SAF_bound.
|
||||
if (
|
||||
fabs(px.value()) > SAF_bound ||
|
||||
fabs(py.value()) > SAF_bound ||
|
||||
fabs(pz.value()) > SAF_bound ||
|
||||
fabs(qx.value()) > SAF_bound ||
|
||||
fabs(qy.value()) > SAF_bound ||
|
||||
fabs(qz.value()) > SAF_bound ||
|
||||
fabs(rx.value()) > SAF_bound ||
|
||||
fabs(ry.value()) > SAF_bound ||
|
||||
fabs(rz.value()) > SAF_bound ||
|
||||
fabs(sx.value()) > SAF_bound ||
|
||||
fabs(sy.value()) > SAF_bound ||
|
||||
fabs(sz.value()) > SAF_bound ||
|
||||
fabs(tx.value()) > SAF_bound ||
|
||||
fabs(ty.value()) > SAF_bound ||
|
||||
fabs(tz.value()) > SAF_bound)
|
||||
fabs(px.to_double()) > SAF_bound ||
|
||||
fabs(py.to_double()) > SAF_bound ||
|
||||
fabs(pz.to_double()) > SAF_bound ||
|
||||
fabs(qx.to_double()) > SAF_bound ||
|
||||
fabs(qy.to_double()) > SAF_bound ||
|
||||
fabs(qz.to_double()) > SAF_bound ||
|
||||
fabs(rx.to_double()) > SAF_bound ||
|
||||
fabs(ry.to_double()) > SAF_bound ||
|
||||
fabs(rz.to_double()) > SAF_bound ||
|
||||
fabs(sx.to_double()) > SAF_bound ||
|
||||
fabs(sy.to_double()) > SAF_bound ||
|
||||
fabs(sz.to_double()) > SAF_bound ||
|
||||
fabs(tx.to_double()) > SAF_bound ||
|
||||
fabs(ty.to_double()) > SAF_bound ||
|
||||
fabs(tz.to_double()) > SAF_bound)
|
||||
{
|
||||
re_adjust:
|
||||
// Compute the new bound.
|
||||
double NEW_bound = 0.0;
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ry.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(sx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(sy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(sz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(tx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ty.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(tz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pz.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qz.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ry.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rz.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(sx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(sy.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(sz.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(tx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ty.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(tz.to_double()));
|
||||
// Re-adjust the context.
|
||||
Static_Filtered_side_of_oriented_sphereC3_15::new_bound(NEW_bound);
|
||||
}
|
||||
|
|
@ -692,21 +692,21 @@ re_adjust:
|
|||
try
|
||||
{
|
||||
return Static_Filtered_side_of_oriented_sphereC3_15::epsilon_variant(
|
||||
px.to_double(),
|
||||
py.to_double(),
|
||||
pz.to_double(),
|
||||
qx.to_double(),
|
||||
qy.to_double(),
|
||||
qz.to_double(),
|
||||
rx.to_double(),
|
||||
ry.to_double(),
|
||||
rz.to_double(),
|
||||
sx.to_double(),
|
||||
sy.to_double(),
|
||||
sz.to_double(),
|
||||
tx.to_double(),
|
||||
ty.to_double(),
|
||||
tz.to_double(),
|
||||
px.dbl(),
|
||||
py.dbl(),
|
||||
pz.dbl(),
|
||||
qx.dbl(),
|
||||
qy.dbl(),
|
||||
qz.dbl(),
|
||||
rx.dbl(),
|
||||
ry.dbl(),
|
||||
rz.dbl(),
|
||||
sx.dbl(),
|
||||
sy.dbl(),
|
||||
sz.dbl(),
|
||||
tx.dbl(),
|
||||
ty.dbl(),
|
||||
tz.dbl(),
|
||||
SAF_epsilon_0);
|
||||
}
|
||||
catch (Restricted_double::unsafe_comparison)
|
||||
|
|
@ -917,40 +917,40 @@ side_of_bounded_sphereC3(
|
|||
|
||||
// Check the bounds. All arguments must be <= SAF_bound.
|
||||
if (
|
||||
fabs(px.value()) > SAF_bound ||
|
||||
fabs(py.value()) > SAF_bound ||
|
||||
fabs(pz.value()) > SAF_bound ||
|
||||
fabs(qx.value()) > SAF_bound ||
|
||||
fabs(qy.value()) > SAF_bound ||
|
||||
fabs(qz.value()) > SAF_bound ||
|
||||
fabs(rx.value()) > SAF_bound ||
|
||||
fabs(ry.value()) > SAF_bound ||
|
||||
fabs(rz.value()) > SAF_bound ||
|
||||
fabs(sx.value()) > SAF_bound ||
|
||||
fabs(sy.value()) > SAF_bound ||
|
||||
fabs(sz.value()) > SAF_bound ||
|
||||
fabs(tx.value()) > SAF_bound ||
|
||||
fabs(ty.value()) > SAF_bound ||
|
||||
fabs(tz.value()) > SAF_bound)
|
||||
fabs(px.to_double()) > SAF_bound ||
|
||||
fabs(py.to_double()) > SAF_bound ||
|
||||
fabs(pz.to_double()) > SAF_bound ||
|
||||
fabs(qx.to_double()) > SAF_bound ||
|
||||
fabs(qy.to_double()) > SAF_bound ||
|
||||
fabs(qz.to_double()) > SAF_bound ||
|
||||
fabs(rx.to_double()) > SAF_bound ||
|
||||
fabs(ry.to_double()) > SAF_bound ||
|
||||
fabs(rz.to_double()) > SAF_bound ||
|
||||
fabs(sx.to_double()) > SAF_bound ||
|
||||
fabs(sy.to_double()) > SAF_bound ||
|
||||
fabs(sz.to_double()) > SAF_bound ||
|
||||
fabs(tx.to_double()) > SAF_bound ||
|
||||
fabs(ty.to_double()) > SAF_bound ||
|
||||
fabs(tz.to_double()) > SAF_bound)
|
||||
{
|
||||
re_adjust:
|
||||
// Compute the new bound.
|
||||
double NEW_bound = 0.0;
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ry.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(sx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(sy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(sz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(tx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ty.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(tz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pz.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qz.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ry.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rz.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(sx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(sy.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(sz.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(tx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ty.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(tz.to_double()));
|
||||
// Re-adjust the context.
|
||||
Static_Filtered_side_of_bounded_sphereC3_15::new_bound(NEW_bound);
|
||||
}
|
||||
|
|
@ -958,21 +958,21 @@ re_adjust:
|
|||
try
|
||||
{
|
||||
return Static_Filtered_side_of_bounded_sphereC3_15::epsilon_variant(
|
||||
px.to_double(),
|
||||
py.to_double(),
|
||||
pz.to_double(),
|
||||
qx.to_double(),
|
||||
qy.to_double(),
|
||||
qz.to_double(),
|
||||
rx.to_double(),
|
||||
ry.to_double(),
|
||||
rz.to_double(),
|
||||
sx.to_double(),
|
||||
sy.to_double(),
|
||||
sz.to_double(),
|
||||
tx.to_double(),
|
||||
ty.to_double(),
|
||||
tz.to_double(),
|
||||
px.dbl(),
|
||||
py.dbl(),
|
||||
pz.dbl(),
|
||||
qx.dbl(),
|
||||
qy.dbl(),
|
||||
qz.dbl(),
|
||||
rx.dbl(),
|
||||
ry.dbl(),
|
||||
rz.dbl(),
|
||||
sx.dbl(),
|
||||
sy.dbl(),
|
||||
sz.dbl(),
|
||||
tx.dbl(),
|
||||
ty.dbl(),
|
||||
tz.dbl(),
|
||||
SAF_epsilon_0,
|
||||
SAF_epsilon_1);
|
||||
}
|
||||
|
|
@ -1128,28 +1128,28 @@ cmp_dist_to_pointC3(
|
|||
|
||||
// Check the bounds. All arguments must be <= SAF_bound.
|
||||
if (
|
||||
fabs(px.value()) > SAF_bound ||
|
||||
fabs(py.value()) > SAF_bound ||
|
||||
fabs(pz.value()) > SAF_bound ||
|
||||
fabs(qx.value()) > SAF_bound ||
|
||||
fabs(qy.value()) > SAF_bound ||
|
||||
fabs(qz.value()) > SAF_bound ||
|
||||
fabs(rx.value()) > SAF_bound ||
|
||||
fabs(ry.value()) > SAF_bound ||
|
||||
fabs(rz.value()) > SAF_bound)
|
||||
fabs(px.to_double()) > SAF_bound ||
|
||||
fabs(py.to_double()) > SAF_bound ||
|
||||
fabs(pz.to_double()) > SAF_bound ||
|
||||
fabs(qx.to_double()) > SAF_bound ||
|
||||
fabs(qy.to_double()) > SAF_bound ||
|
||||
fabs(qz.to_double()) > SAF_bound ||
|
||||
fabs(rx.to_double()) > SAF_bound ||
|
||||
fabs(ry.to_double()) > SAF_bound ||
|
||||
fabs(rz.to_double()) > SAF_bound)
|
||||
{
|
||||
re_adjust:
|
||||
// Compute the new bound.
|
||||
double NEW_bound = 0.0;
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ry.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pz.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qz.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ry.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(rz.to_double()));
|
||||
// Re-adjust the context.
|
||||
Static_Filtered_cmp_dist_to_pointC3_9::new_bound(NEW_bound);
|
||||
}
|
||||
|
|
@ -1157,15 +1157,15 @@ re_adjust:
|
|||
try
|
||||
{
|
||||
return Static_Filtered_cmp_dist_to_pointC3_9::epsilon_variant(
|
||||
px.to_double(),
|
||||
py.to_double(),
|
||||
pz.to_double(),
|
||||
qx.to_double(),
|
||||
qy.to_double(),
|
||||
qz.to_double(),
|
||||
rx.to_double(),
|
||||
ry.to_double(),
|
||||
rz.to_double(),
|
||||
px.dbl(),
|
||||
py.dbl(),
|
||||
pz.dbl(),
|
||||
qx.dbl(),
|
||||
qy.dbl(),
|
||||
qz.dbl(),
|
||||
rx.dbl(),
|
||||
ry.dbl(),
|
||||
rz.dbl(),
|
||||
SAF_epsilon_0);
|
||||
}
|
||||
catch (Restricted_double::unsafe_comparison)
|
||||
|
|
@ -1320,30 +1320,30 @@ cmp_signed_dist_to_planeC3(
|
|||
|
||||
// Check the bounds. All arguments must be <= SAF_bound.
|
||||
if (
|
||||
fabs(pa.value()) > SAF_bound ||
|
||||
fabs(pb.value()) > SAF_bound ||
|
||||
fabs(pc.value()) > SAF_bound ||
|
||||
fabs(pd.value()) > SAF_bound ||
|
||||
fabs(px.value()) > SAF_bound ||
|
||||
fabs(py.value()) > SAF_bound ||
|
||||
fabs(pz.value()) > SAF_bound ||
|
||||
fabs(qx.value()) > SAF_bound ||
|
||||
fabs(qy.value()) > SAF_bound ||
|
||||
fabs(qz.value()) > SAF_bound)
|
||||
fabs(pa.to_double()) > SAF_bound ||
|
||||
fabs(pb.to_double()) > SAF_bound ||
|
||||
fabs(pc.to_double()) > SAF_bound ||
|
||||
fabs(pd.to_double()) > SAF_bound ||
|
||||
fabs(px.to_double()) > SAF_bound ||
|
||||
fabs(py.to_double()) > SAF_bound ||
|
||||
fabs(pz.to_double()) > SAF_bound ||
|
||||
fabs(qx.to_double()) > SAF_bound ||
|
||||
fabs(qy.to_double()) > SAF_bound ||
|
||||
fabs(qz.to_double()) > SAF_bound)
|
||||
{
|
||||
re_adjust:
|
||||
// Compute the new bound.
|
||||
double NEW_bound = 0.0;
|
||||
NEW_bound = std::max(NEW_bound, fabs(pa.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pb.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pc.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pd.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pa.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pb.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pc.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pd.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pz.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qz.to_double()));
|
||||
// Re-adjust the context.
|
||||
Static_Filtered_cmp_signed_dist_to_planeC3_10::new_bound(NEW_bound);
|
||||
}
|
||||
|
|
@ -1351,16 +1351,16 @@ re_adjust:
|
|||
try
|
||||
{
|
||||
return Static_Filtered_cmp_signed_dist_to_planeC3_10::epsilon_variant(
|
||||
pa.to_double(),
|
||||
pb.to_double(),
|
||||
pc.to_double(),
|
||||
pd.to_double(),
|
||||
px.to_double(),
|
||||
py.to_double(),
|
||||
pz.to_double(),
|
||||
qx.to_double(),
|
||||
qy.to_double(),
|
||||
qz.to_double(),
|
||||
pa.dbl(),
|
||||
pb.dbl(),
|
||||
pc.dbl(),
|
||||
pd.dbl(),
|
||||
px.dbl(),
|
||||
py.dbl(),
|
||||
pz.dbl(),
|
||||
qx.dbl(),
|
||||
qy.dbl(),
|
||||
qz.dbl(),
|
||||
SAF_epsilon_0);
|
||||
}
|
||||
catch (Restricted_double::unsafe_comparison)
|
||||
|
|
@ -1556,40 +1556,40 @@ cmp_signed_dist_to_planeC3(
|
|||
|
||||
// Check the bounds. All arguments must be <= SAF_bound.
|
||||
if (
|
||||
fabs(ppx.value()) > SAF_bound ||
|
||||
fabs(ppy.value()) > SAF_bound ||
|
||||
fabs(ppz.value()) > SAF_bound ||
|
||||
fabs(pqx.value()) > SAF_bound ||
|
||||
fabs(pqy.value()) > SAF_bound ||
|
||||
fabs(pqz.value()) > SAF_bound ||
|
||||
fabs(prx.value()) > SAF_bound ||
|
||||
fabs(pry.value()) > SAF_bound ||
|
||||
fabs(prz.value()) > SAF_bound ||
|
||||
fabs(px.value()) > SAF_bound ||
|
||||
fabs(py.value()) > SAF_bound ||
|
||||
fabs(pz.value()) > SAF_bound ||
|
||||
fabs(qx.value()) > SAF_bound ||
|
||||
fabs(qy.value()) > SAF_bound ||
|
||||
fabs(qz.value()) > SAF_bound)
|
||||
fabs(ppx.to_double()) > SAF_bound ||
|
||||
fabs(ppy.to_double()) > SAF_bound ||
|
||||
fabs(ppz.to_double()) > SAF_bound ||
|
||||
fabs(pqx.to_double()) > SAF_bound ||
|
||||
fabs(pqy.to_double()) > SAF_bound ||
|
||||
fabs(pqz.to_double()) > SAF_bound ||
|
||||
fabs(prx.to_double()) > SAF_bound ||
|
||||
fabs(pry.to_double()) > SAF_bound ||
|
||||
fabs(prz.to_double()) > SAF_bound ||
|
||||
fabs(px.to_double()) > SAF_bound ||
|
||||
fabs(py.to_double()) > SAF_bound ||
|
||||
fabs(pz.to_double()) > SAF_bound ||
|
||||
fabs(qx.to_double()) > SAF_bound ||
|
||||
fabs(qy.to_double()) > SAF_bound ||
|
||||
fabs(qz.to_double()) > SAF_bound)
|
||||
{
|
||||
re_adjust:
|
||||
// Compute the new bound.
|
||||
double NEW_bound = 0.0;
|
||||
NEW_bound = std::max(NEW_bound, fabs(ppx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ppy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ppz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pqx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pqy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pqz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(prx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pry.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(prz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qz.value()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ppx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ppy.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(ppz.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pqx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pqy.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pqz.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(prx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pry.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(prz.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(px.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(py.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(pz.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qx.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qy.to_double()));
|
||||
NEW_bound = std::max(NEW_bound, fabs(qz.to_double()));
|
||||
// Re-adjust the context.
|
||||
Static_Filtered_cmp_signed_dist_to_planeC3_15::new_bound(NEW_bound);
|
||||
}
|
||||
|
|
@ -1597,21 +1597,21 @@ re_adjust:
|
|||
try
|
||||
{
|
||||
return Static_Filtered_cmp_signed_dist_to_planeC3_15::epsilon_variant(
|
||||
ppx.to_double(),
|
||||
ppy.to_double(),
|
||||
ppz.to_double(),
|
||||
pqx.to_double(),
|
||||
pqy.to_double(),
|
||||
pqz.to_double(),
|
||||
prx.to_double(),
|
||||
pry.to_double(),
|
||||
prz.to_double(),
|
||||
px.to_double(),
|
||||
py.to_double(),
|
||||
pz.to_double(),
|
||||
qx.to_double(),
|
||||
qy.to_double(),
|
||||
qz.to_double(),
|
||||
ppx.dbl(),
|
||||
ppy.dbl(),
|
||||
ppz.dbl(),
|
||||
pqx.dbl(),
|
||||
pqy.dbl(),
|
||||
pqz.dbl(),
|
||||
prx.dbl(),
|
||||
pry.dbl(),
|
||||
prz.dbl(),
|
||||
px.dbl(),
|
||||
py.dbl(),
|
||||
pz.dbl(),
|
||||
qx.dbl(),
|
||||
qy.dbl(),
|
||||
qz.dbl(),
|
||||
SAF_epsilon_0);
|
||||
}
|
||||
catch (Restricted_double::unsafe_comparison)
|
||||
|
|
|
|||
Loading…
Reference in New Issue