Changeset 5630

Timestamp:

08/31/10 08:53:35 (15 years ago)

Author:

Mathieu Morlighem

Message:

Added GaussVertex method of GaussTria, very useful

Location:

issm/trunk/src/c/objects

Files:

• Elements/Tria.cpp (modified) (1 diff)
• Gauss/GaussTria.cpp (modified) (1 diff)
• Gauss/GaussTria.h (modified) (1 diff)

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

@@ -851 +851 @@
 /*FUNCTION Tria::GetBedList {{{1*/
 void  Tria::GetBedList(double* bedlist){
-        
-        const int numvertices=3;
-        double  gaussgrids[numvertices][numvertices]={{1,0,0},{0,1,0},{0,0,1}};
-        
-        inputs->GetParameterValues(bedlist,&gaussgrids[0][0],3,BedEnum);
+
+        Input     *bed_input = NULL;
+        GaussTria *gauss     = NULL;
+
+        /*Recover input*/
+        bed_input=inputs->GetInput(BedEnum);
+
+        /*Checks in debugging mode*/
+        ISSMASSERT(bedlist);
+        ISSMASSERT(bed_input);
+
+        /* Start looping on the number of vertices: */
+        gauss=new GaussTria();
+        for (int iv=0;iv<3;iv++){
+                gauss->GaussVertex(iv);
+                bed_input->GetParameterValue(&bedlist[iv],gauss);
+        }
+
+        /*clean-up*/
+        delete gauss;
 
 }

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

@@ -116 +116 @@
 }
 /*}}}*/
+/*FUNCTION GaussTria::GaussVertex{{{1*/
+void GaussTria::GaussVertex(int iv){
+
+        /*in debugging mode: check that the default constructor has been called*/
+        ISSMASSERT(numgauss==-1);
+
+        /*update static arrays*/
+        switch(iv){
+                case 0:
+                        coord1=1; coord2=0; coord3=0;
+                        break;
+                case 1:
+                        coord1=0; coord2=1; coord3=0;
+                        break;
+                case 2:
+                        coord1=0; coord2=0; coord3=1;
+                        break;
+                default:
+                        ISSMERROR("vertex index should be in [0 2]");
+
+        }
+
+}
+/*}}}*/
 /*FUNCTION GaussTria::begin{{{1*/
 int GaussTria::begin(void){

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

@@ -39 +39 @@
                 void Echo(void);
                 void GaussPoint(int ig);
+                void GaussVertex(int iv);
 };
 #endif  /* _GAUSSTRIA_H_ */