Changeset 17750 for issm/trunk-jpl/src/c/classes/gauss/GaussTetra.cpp

Timestamp:

04/16/14 10:31:31 (11 years ago)

Author:

Mathieu Morlighem

Message:

CHG: fixing indices of tetra base that need to be direct

File:

• issm/trunk-jpl/src/c/classes/gauss/GaussTetra.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/classes/gauss/GaussTetra.cpp

@@ -53 +53 @@
                 for(int i=0;i<numgauss;i++) coords4[i]=0.;
         }
-        else if(index1==0 && index2==1 && index3==3){
+        else if(index1==0 && index2==3 && index3==1){
                 GaussLegendreTria(&numgauss,&coords1,&coords2,&coords4,&weights,order);
                 coords3=xNew<IssmDouble>(numgauss);
                 for(int i=0;i<numgauss;i++) coords3[i]=0.;
         }
-        else if(index1==1 && index2==2 && index3==3){
+        else if(index1==1 && index2==3 && index3==2){
                 GaussLegendreTria(&numgauss,&coords2,&coords3,&coords4,&weights,order);
                 coords1=xNew<IssmDouble>(numgauss);
                 for(int i=0;i<numgauss;i++) coords1[i]=0.;
         }
-        else if(index1==2 && index2==0 && index3==3){
+        else if(index1==0 && index2==2 && index3==3){
                 GaussLegendreTria(&numgauss,&coords1,&coords3,&coords4,&weights,order);
                 coords2=xNew<IssmDouble>(numgauss);