Changeset 12493 for issm/trunk-jpl/src/c/objects/Elements/PentaRef.cpp
- Timestamp:
- 06/21/12 10:02:29 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/objects/Elements/PentaRef.cpp
r12471 r12493 949 949 /*Get Determinant*/ 950 950 Matrix3x3Determinant(Jdet,&J[0][0]); 951 if(*Jdet<0) _error _("negative jacobian determinant!");951 if(*Jdet<0) _error2_("negative jacobian determinant!"); 952 952 953 953 } … … 972 972 /*Jdet = norm( AB ^ AC ) / (2 * area of the reference triangle), with areaRef=sqrt(3) */ 973 973 *Jdet=SQRT3/6.0*pow(pow(((y2-y1)*(z3-z1)-(z2-z1)*(y3-y1)),2.0)+pow(((z2-z1)*(x3-x1)-(x2-x1)*(z3-z1)),2.0)+pow(((x2-x1)*(y3-y1)-(y2-y1)*(x3-x1)),2.0),0.5); 974 if(*Jdet<0) _error _("negative jacobian determinant!");974 if(*Jdet<0) _error2_("negative jacobian determinant!"); 975 975 } 976 976 /*}}}*/ … … 990 990 991 991 *Jdet=1.0/2.0*sqrt(pow(x2-x1,2.) + pow(y2-y1,2.) + pow(z2-z1,2.)); 992 if(*Jdet<0) _error _("negative jacobian determinant!");992 if(*Jdet<0) _error2_("negative jacobian determinant!"); 993 993 994 994 } … … 1226 1226 /*Area of a trabezoid = altitude * (base1 + base2)/2 */ 1227 1227 *Jdet= pow(pow(x2-x1,2.) + pow(y2-y1,2.),0.5) * (z4-z1 + z3-z2)/8; 1228 if(*Jdet<0) _error _("negative jacobian determinant!");1228 if(*Jdet<0) _error2_("negative jacobian determinant!"); 1229 1229 1230 1230 }
Note:
See TracChangeset
for help on using the changeset viewer.