Changeset 5675

Timestamp:

09/03/10 11:22:19 (15 years ago)

Author:

Mathieu Morlighem

Message:

Added GaussFaceTria

Location:

issm/trunk/src/c/objects/Gauss

Files:

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

issm/trunk/src/c/objects/Gauss/GaussPenta.cpp

@@ -195 +195 @@
 }
 /*}}}*/
+/*FUNCTION GaussPenta::GaussFaceTria{{{1*/
+void GaussPenta::GaussFaceTria(int index1, int index2, int index3, int order){
+
+        /*in debugging mode: check that the default constructor has been called*/
+        ISSMASSERT(numgauss==-1);
+
+        /*Basal Tria*/
+        if(index1==0 && index2==1 && index3==2){
+                GaussLegendreTria(&numgauss,&coords1,&coords2,&coords3,&weights,order);
+                coords4=(double*)xmalloc(numgauss*sizeof(double));
+                for(int i=0;i<numgauss;i++) coords4[i]=-1.0;
+        }
+        else{
+                ISSMERROR("Tria not supported yet");
+        }
+
+}
+/*}}}*/
 /*FUNCTION GaussPenta::begin{{{1*/
 int GaussPenta::begin(void){

issm/trunk/src/c/objects/Gauss/GaussPenta.h

@@ -42 +42 @@
                 void GaussPoint(int ig);
                 void GaussVertex(int iv);
+                void GaussFaceTria(int index1, int index2, int index3, int order);
                 void GaussCenter(void);
                 void SynchronizeGaussTria(GaussTria* gauss_tria);