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

Last change on this file since 16556 was 16556, checked in by Mathieu Morlighem, 11 years ago

NEW: added Archive/16133-16554
File size: 1.3 KB
RevLine 
[16556]1Index: ../trunk-jpl/src/c/Makefile.am
2===================================================================
3--- ../trunk-jpl/src/c/Makefile.am (revision 16377)
4+++ ../trunk-jpl/src/c/Makefile.am (revision 16378)
5@@ -24,6 +24,8 @@
6 ./datastructures/Object.h\
7 ./datastructures/datastructures.h\
8 ./classes/classes.h\
9+ ./classes/gauss/GaussSeg.h\
10+ ./classes/gauss/GaussSeg.cpp\
11 ./classes/gauss/GaussTria.h\
12 ./classes/gauss/GaussTria.cpp\
13 ./classes/Update.h\
14Index: ../trunk-jpl/src/c/classes/Elements/SegRef.cpp
15===================================================================
16--- ../trunk-jpl/src/c/classes/Elements/SegRef.cpp (revision 16377)
17+++ ../trunk-jpl/src/c/classes/Elements/SegRef.cpp (revision 16378)
18@@ -75,8 +75,7 @@
19 IssmDouble y2=xyz_list[3*1+1];
20 IssmDouble z2=xyz_list[3*1+2];
21
22- *Jdet=.5*sqrt(pow(x2-x1,2) + pow(y2-y1,2) + pow(z2-z1,2));
23- if(*Jdet<0) _error_("negative jacobian determinant!");
24+ *J=.5*sqrt(pow(x2-x1,2) + pow(y2-y1,2) + pow(z2-z1,2));
25 }
26 /*}}}*/
27 /*FUNCTION SegRef::GetJacobianDeterminant{{{*/
28@@ -85,7 +84,8 @@
29 * J is assumed to have been allocated of size NDOF2xNDOF2.*/
30
31 /*Call Jacobian routine to get the jacobian:*/
32- GetJacobian(&Jdet, xyz_list, gauss);
33+ GetJacobian(Jdet, xyz_list, gauss);
34+ if(*Jdet<0) _error_("negative jacobian determinant!");
35
36 }
37 /*}}}*/
Note: See TracBrowser for help on using the repository browser.