Changeset 16378
- Timestamp:
- 10/11/13 07:37:38 (11 years ago)
- Location:
- issm/trunk-jpl/src/c
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
TabularUnified issm/trunk-jpl/src/c/Makefile.am ¶
r16371 r16378 25 25 ./datastructures/datastructures.h\ 26 26 ./classes/classes.h\ 27 ./classes/gauss/GaussSeg.h\ 28 ./classes/gauss/GaussSeg.cpp\ 27 29 ./classes/gauss/GaussTria.h\ 28 30 ./classes/gauss/GaussTria.cpp\ -
TabularUnified issm/trunk-jpl/src/c/classes/Elements/SegRef.cpp ¶
r16377 r16378 76 76 IssmDouble z2=xyz_list[3*1+2]; 77 77 78 *Jdet=.5*sqrt(pow(x2-x1,2) + pow(y2-y1,2) + pow(z2-z1,2)); 79 if(*Jdet<0) _error_("negative jacobian determinant!"); 78 *J=.5*sqrt(pow(x2-x1,2) + pow(y2-y1,2) + pow(z2-z1,2)); 80 79 } 81 80 /*}}}*/ … … 86 85 87 86 /*Call Jacobian routine to get the jacobian:*/ 88 GetJacobian(&Jdet, xyz_list, gauss); 87 GetJacobian(Jdet, xyz_list, gauss); 88 if(*Jdet<0) _error_("negative jacobian determinant!"); 89 89 90 90 }
Note:
See TracChangeset
for help on using the changeset viewer.