Changeset 22983

Timestamp:

07/19/18 16:37:45 (7 years ago)

Author:

Eric.Larour

Message:

CHG: added new JacobianDeterminant3D, which does not error out if jacobian is negative. Only used for domain 3d surface.

Location:

issm/trunk-jpl/src/c/classes/Elements

Files:

• TriaRef.cpp (modified) (1 diff)
• TriaRef.h (modified) (1 diff)

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

@@ -116 +116 @@
         Matrix2x2Determinant(Jdet,&J[0][0]);
         if(*Jdet<0) _error_("negative jacobian determinant!");
+
+}
+/*}}}*/
+void TriaRef::GetJacobianDeterminant3D(IssmDouble* Jdet, IssmDouble* xyz_list,Gauss* gauss){/*{{{*/
+        /*The Jacobian determinant is constant over the element, discard the gaussian points. 
+         * J is assumed to have been allocated of size NDOF2xNDOF2.*/
+        IssmDouble J[2][2];
+
+        /*Get Jacobian*/
+        GetJacobian(&J[0][0],xyz_list,gauss);
+
+        /*Get Determinant*/
+        Matrix2x2Determinant(Jdet,&J[0][0]);
 
 }

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

@@ -20 +20 @@
                 void GetJacobian(IssmDouble* J, IssmDouble* xyz_list,Gauss* gauss);
                 void GetJacobianDeterminant(IssmDouble* Jdet, IssmDouble* xyz_list,Gauss* gauss);
+                void GetJacobianDeterminant3D(IssmDouble* Jdet, IssmDouble* xyz_list,Gauss* gauss);
                 void GetJacobianInvert(IssmDouble*  Jinv, IssmDouble* xyz_list,Gauss* gauss);
                 void GetNodalFunctions(IssmDouble* basis,Gauss* gauss,int finiteelement);