Commit 16cb03a9 authored by Valentin Nigolian's avatar Valentin Nigolian

might have done the trick

parent f112bd83
Pipeline #14860 passed with stage
in 15 minutes and 3 seconds
......@@ -160,10 +160,10 @@ private:
static Scalar tet_weight(Scalar alpha, Scalar beta, Scalar theta) {
if(abs(sin(alpha)) > std::numeric_limits<Scalar>::epsilon() &&
abs(sin(beta)) > std::numeric_limits<Scalar>::epsilon() &&
abs(sin(theta)) > std::numeric_limits<Scalar>::epsilon() &&
abs(cos(theta)) > std::numeric_limits<Scalar>::epsilon()){
if(fabs(sin(alpha)) > std::numeric_limits<Scalar>::epsilon() &&
fabs(sin(beta)) > std::numeric_limits<Scalar>::epsilon() &&
fabs(sin(theta)) > std::numeric_limits<Scalar>::epsilon() &&
fabs(cos(theta)) > std::numeric_limits<Scalar>::epsilon()){
double cotan_alpha = cot(alpha);
double cotan_beta = cot(beta);
......@@ -182,10 +182,10 @@ private:
std::cout<<" abs sin theta : "<<std::setprecision(20)<<abs(sin(theta))<<std::endl;
std::cout<<" abs cos theta : "<<std::setprecision(20)<<abs(cos(theta))<<std::endl;
std::cout<<" comparisons : "<<
(abs(sin(alpha)) > std::numeric_limits<Scalar>::epsilon()) << " "<<
(abs(sin(beta)) > std::numeric_limits<Scalar>::epsilon()) << " "<<
(abs(sin(theta)) > std::numeric_limits<Scalar>::epsilon()) << " "<<
(abs(cos(theta)) > std::numeric_limits<Scalar>::epsilon()) << std::endl;
(fabs(sin(alpha)) > std::numeric_limits<Scalar>::epsilon()) << " " <<
(fabs(sin(beta)) > std::numeric_limits<Scalar>::epsilon()) << " "<<
(fabs(sin(theta)) > std::numeric_limits<Scalar>::epsilon()) << " "<<
(fabs(cos(theta)) > std::numeric_limits<Scalar>::epsilon()) << std::endl;
return 0.f;
}
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment