Changeset 24391

Timestamp:

11/23/19 06:33:45 (5 years ago)

Author:

Mathieu Morlighem

Message:

BUG: cannot request a Gauss point of order 0

Location:

issm/trunk-jpl/src/c

Files:

• analyses/HydrologyShaktiAnalysis.cpp (modified) (1 diff)
• shared/Numerics/GaussPoints.cpp (modified) (4 diffs)

issm/trunk-jpl/src/c/analyses/HydrologyShaktiAnalysis.cpp ¶

@@ -412 +412 @@
 
         /*Get gap height derivatives at the center of the element*/
-        Gauss* gauss=element->NewGauss(0);
+        Gauss* gauss=element->NewGauss(1);
         head_input->GetInputDerivativeValue(&dh[0],xyz_list,gauss);
         delete gauss;

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

@@ -20 +20 @@
                 For degree p, the required number of Gauss-Legendre points is
                 n>=(p+1)/2.*/
+        _assert_(ngaus>0);
 
         /*Intermediaries*/
@@ -98 +99 @@
           Symmetrical Gaussian Quadrature Rules for the Triangle", IJNME,
           Vol. 21, pp. 1129-1148 (1985), as transcribed for Probe rules3.*/
+        _assert_(iord>0);
 
         /*Intermediaries*/
@@ -1218 +1220 @@
                 Quadrature Formulas", Computer Methods in Applied Mechanics and
                 Engineering, Vol. 55, pp. 339-348 (1986).*/
+        _assert_(iord>0);
 
         /*Intermediaries*/
@@ -1484 +1487 @@
           For degree p, the required number of Gauss-Lobatto points is
           n>=(p+1)/2+1 (one more than Gauss-Legendre).*/
+        _assert_(ngaus>0);
 
         int i;