Changeset 5739 for issm/trunk/src/c/objects/Gauss/GaussTria.cpp

Timestamp:

09/10/10 09:26:58 (15 years ago)

Author:

Mathieu Morlighem

Message:

NumericalFlux does not use beams anymore

File:

• issm/trunk/src/c/objects/Gauss/GaussTria.cpp (modified) (2 diffs)

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

@@ -56 +56 @@
         weights=(double*)xmalloc(numgauss*sizeof(double));
 
-        /*Reverse grid1 and 2 if necessary*/
+        /*Reverse index1 and 2 if necessary*/
         if (index1>index2){
                 index3=index1; index1=index2; index2=index3;
@@ -150 +150 @@
         coord2=ONETHIRD;
         coord3=ONETHIRD;
+
+}
+/*}}}*/
+/*FUNCTION GaussTria::GaussEdgeCenter{{{1*/
+void GaussTria::GaussEdgeCenter(int index1,int index2){
+
+        int     index3;
+
+        /*Reverse index1 and 2 if necessary*/
+        if (index1>index2){
+                index3=index1; index1=index2; index2=index3;
+        }
+
+        /*update static arrays*/
+        if (index1==0 && index2==1){
+                coord1=0.5;
+                coord2=0.5;
+                coord3=0.0;
+        }
+        else if (index1==0 && index2==2){
+                coord1=0.5;
+                coord2=0.0;
+                coord3=0.5;
+        }
+        else if (index1==1 && index2==2){
+                coord1=0.0;
+                coord2=0.5;
+                coord3=0.5;
+        }
+        else
+         ISSMERROR("The 2 indices provided are not supported yet (user provided %i and %i)",index1,index2);
 
 }