Changeset 17421

Timestamp:

03/12/14 10:45:59 (11 years ago)

Author:

Mathieu Morlighem

Message:

BUG: wrong gauss point for order 3 tetra

File:

• issm/trunk-jpl/src/c/shared/Numerics/GaussPoints.cpp (modified) (2 diffs)

issm/trunk-jpl/src/c/shared/Numerics/GaussPoints.cpp ¶

@@ -1253 +1253 @@
                 0.250000000000000, 0.500000000000000, 0.166666666666667, 
                 0.166666666666667, 0.166666666666667};
+        static IssmPDouble l23[]={
+                0.250000000000000, 0.166666666666667, 0.500000000000000, 
+                0.166666666666667, 0.166666666666667};
         static IssmPDouble l33[]={
                 0.250000000000000, 0.166666666666667, 0.166666666666667, 
@@ -1376 +1379 @@
         static IssmPDouble* wgtp[MAX_TETRA_SYM_ORD]={wgt1,wgt2,wgt3,wgt4,wgt5,wgt6};
         static IssmPDouble* l1p [MAX_TETRA_SYM_ORD]={l11 ,l12 ,l13 ,l14 ,l15 ,l16 };
-        static IssmPDouble* l2p [MAX_TETRA_SYM_ORD]={l21 ,l22 ,l32 ,l24 ,l25 ,l26 };
+        static IssmPDouble* l2p [MAX_TETRA_SYM_ORD]={l21 ,l22 ,l23 ,l24 ,l25 ,l26 };
         static IssmPDouble* l3p [MAX_TETRA_SYM_ORD]={l31 ,l32 ,l33 ,l34 ,l35 ,l36 };
         static IssmPDouble* l4p [MAX_TETRA_SYM_ORD]={l41 ,l42 ,l43 ,l44 ,l45 ,l46 };