Changeset 5625 for issm/trunk/src/c/objects/Elements/TriaRef.cpp

Timestamp:

08/30/10 17:01:51 (15 years ago)

Author:

Mathieu Morlighem

Message:

Added GaussTria for slope only for now

File:

• issm/trunk/src/c/objects/Elements/TriaRef.cpp (modified) (5 diffs)

issm/trunk/src/c/objects/Elements/TriaRef.cpp

@@ -180 +180 @@
 }
 /*}}}*/
-/*FUNCTION TriaRef::GetL {{{1*/
+/*FUNCTION TriaRef::GetL(double* L, double* xyz_list, double* gauss,int numdof) {{{1*/
 void TriaRef::GetL(double* L, double* xyz_list, double* gauss,int numdof){
         /*Compute L  matrix. L=[L1 L2 L3] where Li is square and of size numdof. 
@@ -221 +221 @@
 }
 /*}}}*/
-/*FUNCTION TriaRef::GetJacobian {{{1*/
+/*FUNCTION TriaRef::GetL(double* L, double* xyz_list,GaussTria* gauss,int numdof) {{{1*/
+void TriaRef::GetL(double* L, double* xyz_list,GaussTria* gauss,int numdof){
+        /*Compute L  matrix. L=[L1 L2 L3] where Li is square and of size numdof. 
+         * For grid i, Li can be expressed in the actual coordinate system
+         * by: 
+         *       numdof=1: 
+         *                 Li=h;
+         *       numdof=2:
+         *                 Li=[ h   0 ]
+         *                    [ 0   h ]
+         * where h is the interpolation function for grid i.
+         *
+         * We assume L has been allocated already, of size: numgrids (numdof=1), or numdofx(numdof*numgrids) (numdof=2)
+         */
+
+        int i;
+        const int NDOF2=2;
+        const int numgrids=3;
+        double l1l2l3[3];
+
+        /*Get l1l2l3 in actual coordinate system: */
+        GetNodalFunctions(l1l2l3,gauss);
+
+        /*Build L: */
+        if(numdof==1){
+                for (i=0;i<numgrids;i++){
+                        L[i]=l1l2l3[i]; 
+                }
+        }
+        else{
+                for (i=0;i<numgrids;i++){
+                        *(L+numdof*numgrids*0+numdof*i)=l1l2l3[i]; //L[0][NDOF2*i]=dh1dh3[0][i];
+                        *(L+numdof*numgrids*0+numdof*i+1)=0;
+                        *(L+numdof*numgrids*1+numdof*i)=0;
+                        *(L+numdof*numgrids*1+numdof*i+1)=l1l2l3[i];
+                }
+        }
+}
+/*}}}*/
+/*FUNCTION TriaRef::GetJacobian(double* J, double* xyz_list,double* gauss) {{{1*/
 void TriaRef::GetJacobian(double* J, double* xyz_list,double* gauss){
         /*The Jacobian is constant over the element, discard the gaussian points. 
@@ -244 +283 @@
 }
 /*}}}*/
-/*FUNCTION TriaRef::GetJacobianDeterminant2d {{{1*/
+/*FUNCTION TriaRef::GetJacobian(double* J, double* xyz_list,GaussTria* gauss) {{{1*/
+void TriaRef::GetJacobian(double* J, double* xyz_list,GaussTria* gauss){
+        /*The Jacobian is constant over the element, discard the gaussian points. 
+         * J is assumed to have been allocated of size NDOF2xNDOF2.*/
+
+        const int NDOF2=2;
+        const int numgrids=3;
+        double x1,y1,x2,y2,x3,y3;
+
+        x1=*(xyz_list+numgrids*0+0);
+        y1=*(xyz_list+numgrids*0+1);
+        x2=*(xyz_list+numgrids*1+0);
+        y2=*(xyz_list+numgrids*1+1);
+        x3=*(xyz_list+numgrids*2+0);
+        y3=*(xyz_list+numgrids*2+1);
+
+
+        *(J+NDOF2*0+0)=0.5*(x2-x1);
+        *(J+NDOF2*1+0)=SQRT3/6.0*(2*x3-x1-x2);
+        *(J+NDOF2*0+1)=0.5*(y2-y1);
+        *(J+NDOF2*1+1)=SQRT3/6.0*(2*y3-y1-y2);
+}
+/*}}}*/
+/*FUNCTION TriaRef::GetJacobianDeterminant2d(double* Jdet, double* xyz_list,double* gauss) {{{1*/
 void TriaRef::GetJacobianDeterminant2d(double* Jdet, double* xyz_list,double* gauss){
 
+        /*The Jacobian determinant is constant over the element, discard the gaussian points. 
+         * J is assumed to have been allocated of size NDOF2xNDOF2.*/
+        double J[2][2];
+
+        /*Get Jacobian*/
+        GetJacobian(&J[0][0],xyz_list,gauss);
+
+        /*Get Determinant*/
+        Matrix2x2Determinant(Jdet,&J[0][0]);
+        if(*Jdet<0) ISSMERROR("negative jacobian determinant!");
+
+}
+/*}}}*/
+/*FUNCTION TriaRef::GetJacobianDeterminant2d(double* Jdet, double* xyz_list,GaussTria* gauss) {{{1*/
+void TriaRef::GetJacobianDeterminant2d(double* Jdet, double* xyz_list,GaussTria* gauss){
         /*The Jacobian determinant is constant over the element, discard the gaussian points. 
          * J is assumed to have been allocated of size NDOF2xNDOF2.*/
@@ -296 +373 @@
 }
 /*}}}*/
-/*FUNCTION TriaRef::GetNodalFunctions {{{1*/
+/*FUNCTION TriaRef::GetNodalFunctions(double* l1l2l3, double* gauss) {{{1*/
 void TriaRef::GetNodalFunctions(double* l1l2l3, double* gauss){
         /*This routine returns the values of the nodal functions  at the gaussian point.*/
@@ -303 +380 @@
         l1l2l3[1]=gauss[1];
         l1l2l3[2]=gauss[2];
+
+}
+/*}}}*/
+/*FUNCTION TriaRef::GetNodalFunctions(double* l1l2l3,GaussTria* gauss) {{{1*/
+void TriaRef::GetNodalFunctions(double* l1l2l3,GaussTria* gauss){
+        /*This routine returns the values of the nodal functions  at the gaussian point.*/
+
+        l1l2l3[0]=gauss->coord1;
+        l1l2l3[1]=gauss->coord2;
+        l1l2l3[2]=gauss->coord3;
 
 }