source: issm/oecreview/Archive/16133-16554/ISSM-16377-16378.diff

../trunk-jpl/src/c/Makefile.am

@@ -24 +24 @@
                                         ./datastructures/Object.h\
                                         ./datastructures/datastructures.h\
                                         ./classes/classes.h\
+                                        ./classes/gauss/GaussSeg.h\
+                                        ./classes/gauss/GaussSeg.cpp\
                                         ./classes/gauss/GaussTria.h\
                                         ./classes/gauss/GaussTria.cpp\
                                         ./classes/Update.h\

../trunk-jpl/src/c/classes/Elements/SegRef.cpp

@@ -75 +75 @@
         IssmDouble y2=xyz_list[3*1+1];
         IssmDouble z2=xyz_list[3*1+2];
 
-        *Jdet=.5*sqrt(pow(x2-x1,2) + pow(y2-y1,2) + pow(z2-z1,2));
-        if(*Jdet<0) _error_("negative jacobian determinant!");
+        *J=.5*sqrt(pow(x2-x1,2) + pow(y2-y1,2) + pow(z2-z1,2));
 }
 /*}}}*/
 /*FUNCTION SegRef::GetJacobianDeterminant{{{*/
@@ -85 +84 @@
          * J is assumed to have been allocated of size NDOF2xNDOF2.*/
 
         /*Call Jacobian routine to get the jacobian:*/
-        GetJacobian(&Jdet, xyz_list, gauss);
+        GetJacobian(Jdet, xyz_list, gauss);
+        if(*Jdet<0) _error_("negative jacobian determinant!");
 
 }
 /*}}}*/