Changeset 26533

Timestamp:

11/04/21 08:28:38 (3 years ago)

Author:

Mathieu Morlighem

Message:

CHG: prevent CoDiPack from complaining about power of a negative number

File:

• issm/trunk-jpl/src/c/classes/Elements/TriaRef.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/classes/Elements/TriaRef.cpp

@@ -236 +236 @@
         IssmDouble y2 = xyz_list[3*1+1];
 
-        *Jdet = .5*sqrt(pow(x2-x1,2) + pow(y2-y1,2));
+        *Jdet = .5*sqrt((x2-x1)*(x2-x1) + (y2-y1)*(y2-y1));
         if(*Jdet<0) _error_("negative jacobian determinant!");