Changeset 15446

Timestamp:

07/05/13 14:51:53 (12 years ago)

Author:

Mathieu Morlighem

Message:

BUG: fixed P2 basis functions

Location:

issm/trunk-jpl/src/c/classes/Elements

Files:

• Tria.cpp (modified) (2 diffs)
• TriaRef.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/classes/Elements/Tria.cpp

@@ -992 +992 @@
 
         /* Start looping on the number of vertices: */
-        if (input){
+        if(input){
                 GaussTria* gauss=new GaussTria();
-                for (int iv=0;iv<this->NumberofNodes();iv++){
+                for(int iv=0;iv<this->NumberofNodes();iv++){
                         gauss->GaussNode(numnodes,iv);
                         input->GetInputValue(&pvalue[iv],gauss);
@@ -1001 +1001 @@
         }
         else{
-                for (int iv=0;iv<numnodes;iv++) pvalue[iv]=defaultvalue;
+                for(int iv=0;iv<numnodes;iv++) pvalue[iv]=defaultvalue;
         }
 }

issm/trunk-jpl/src/c/classes/Elements/TriaRef.cpp

@@ -541 +541 @@
                         dbasis[NUMNODESP2*1+0]=-2.*SQRT3/3.*gauss->coord1 + SQRT3/6.;
                         /*Nodal function 2*/
-                        dbasis[NUMNODESP2*0+1]=+2.*gauss->coord2 + 0.5;
+                        dbasis[NUMNODESP2*0+1]=+2.*gauss->coord2 - 0.5;
                         dbasis[NUMNODESP2*1+1]=-2.*SQRT3/3.*gauss->coord2 + SQRT3/6.;
                         /*Nodal function 3*/