Commit c18c7752 authored by Valentin Nigolian's avatar Valentin Nigolian

fixed fabs issue

parent 9da61db0
Pipeline #14862 passed with stage
in 4 minutes and 6 seconds
......@@ -47,7 +47,7 @@ public:
Scalar weight(0);
double alpha, beta, theta;
double alpha(0), beta(0), theta(0);
auto hehf_iter = mesh.hehf_iter(edge);
......@@ -126,10 +126,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);
......
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