Changeset 25441

Timestamp:

08/21/20 11:31:20 (5 years ago)

Author:

Mathieu Morlighem

Message:

BUG: missing initialization of ig

File:

• issm/trunk-jpl/src/c/classes/gauss/GaussTria.cpp (modified) (3 diffs)

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

@@ -309 +309 @@
          *
          */
-        int         ig;
         IssmDouble x,y;
         IssmDouble xy_list[3][2];
@@ -325 +324 @@
         this->weights=xNew<IssmDouble>(this->numgauss);
 
-        for(ig=0;ig<gauss1->numgauss;ig++){ // Add the first triangle gauss points
-                this->coords1[ig]=gauss1->coords1[ig];
-                this->coords2[ig]=gauss1->coords2[ig];
-                this->coords3[ig]=gauss1->coords3[ig];
-                this->weights[ig]=gauss1->weights[ig];
-        }
-        for(ig=0;ig<gauss2->numgauss;ig++){ // Add the second triangle gauss points
-                this->coords1[gauss1->numgauss+ig]=gauss2->coords1[ig];
-                this->coords2[gauss1->numgauss+ig]=gauss2->coords2[ig];
-                this->coords3[gauss1->numgauss+ig]=gauss2->coords3[ig];
-                this->weights[gauss1->numgauss+ig]=gauss2->weights[ig];
+        for(int ii=0;ii<gauss1->numgauss;ii++){ // Add the first triangle gauss points
+                this->coords1[ii]=gauss1->coords1[ii];
+                this->coords2[ii]=gauss1->coords2[ii];
+                this->coords3[ii]=gauss1->coords3[ii];
+                this->weights[ii]=gauss1->weights[ii];
+        }
+        for(int ii=0;ii<gauss2->numgauss;ii++){ // Add the second triangle gauss points
+                this->coords1[gauss1->numgauss+ii]=gauss2->coords1[ii];
+                this->coords2[gauss1->numgauss+ii]=gauss2->coords2[ii];
+                this->coords3[gauss1->numgauss+ii]=gauss2->coords3[ii];
+                this->weights[gauss1->numgauss+ii]=gauss2->weights[ii];
         }
 
@@ -343 +342 @@
 
         /*Initialize static fields as undefined*/
+        ig    = -1;
         weight=UNDEF;
         coord1=UNDEF;