Changeset 3877 for issm/trunk/src/c/objects/Elements/Penta.cpp

Timestamp:

05/21/10 12:12:38 (15 years ago)

Author:

Mathieu Morlighem

Message:

some fixing

File:

• issm/trunk/src/c/objects/Elements/Penta.cpp (modified) (1 diff)

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

@@ -459 +459 @@
 
         return tria;
+}
+/*}}}*/
+/*FUNCTION Penta::GaussFromNode {{{1*/
+double* Penta::GaussFromNode(Node* node){
+
+        /*variable declaration*/
+        int i,pos;
+        double*  gauss=NULL;
+
+        /*Allocate gauss*/
+        gauss=(double*)xmalloc(4*sizeof(double));
+
+        for(i=0;i<6;i++){
+                if (node==nodes[i]){
+                        switch(i){
+                                case 0:
+                                        gauss[0]=1.0; gauss[1]=0.0; gauss[2]=0.0; gauss[3]= -1.0;
+                                        return gauss;
+                                case 1:
+                                        gauss[0]=0.0; gauss[1]=1.0; gauss[2]=0.0; gauss[3]= -1.0;
+                                        return gauss;
+                                case 2:
+                                        gauss[0]=0.0; gauss[1]=0.0; gauss[2]=1.0; gauss[3]= -1.0;
+                                        return gauss;
+                                case 3:
+                                        gauss[0]=1.0; gauss[1]=0.0; gauss[2]=0.0; gauss[3]=1.0;
+                                        return gauss;
+                                case 4:
+                                        gauss[0]=0.0; gauss[1]=1.0; gauss[2]=0.0; gauss[3]=1.0;
+                                        return gauss;
+                                case 5:
+                                        gauss[0]=0.0; gauss[1]=0.0; gauss[2]=1.0; gauss[3]=1.0;
+                                        return gauss;
+                        }
+                }
+        }
+
+        /*output error if not found*/
+        ISSMERROR("Node not found in Penta");
 }
 /*}}}*/