Changeset 16961

Timestamp:

11/26/13 17:15:43 (11 years ago)

Author:

seroussi

Message:

CHG: fixed leak

File:

• issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp (modified) (3 diffs)

issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp

@@ -3507 +3507 @@
         IssmDouble LprimeFSSSA[4][numdof2dm];
         IssmDouble DLFSSSA[4][4]={0.0};
-        IssmDouble Ke_drag_gaussian[numdof2dm][numdofs];
-        IssmDouble Ke_drag_gaussian2[numdof2d][numdof2dm];
+        IssmDouble Ke_drag[numdof2dm][numdofs];
+        IssmDouble Ke_drag2[numdof2d][numdof2dm];
         IssmDouble *xyz_list      = NULL;
         IssmDouble *xyz_list_tria = NULL;
@@ -3586 +3586 @@
                                         &DLSSAFS[0][0],8,8,0,
                                         &LprimeSSAFS[0][0],8,numdofs,0,
-                                        &Ke_drag_gaussian[0][0],0);
+                                        &Ke_drag[0][0],1);
 
                 TripleMultiply( &LFSSSA[0][0],4,numdof2d,1,
                                         &DLFSSSA[0][0],4,4,0,
                                         &LprimeFSSSA[0][0],4,numdof2dm,0,
-                                        &Ke_drag_gaussian2[0][0],0);
-                for(i=0;i<numdof2dm;i++) for(j=0;j<numdofs;j++) Ke->values[i*numdoftot+j+numdofm]+=Ke_drag_gaussian[i][j];
-                for(i=0;i<numdof2d;i++) for(j=0;j<numdof2dm;j++) Ke->values[(i+numdofm)*numdoftot+j]+=Ke_drag_gaussian2[i][j];
-        }
-
+                                        &Ke_drag2[0][0],1);
+        }
+
+        for(i=0;i<numdof2dm;i++) for(j=0;j<numdofs;j++) Ke->values[i*numdoftot+j+numdofm]+=Ke_drag[i][j];
+        for(i=0;i<numdof2d;i++) for(j=0;j<numdof2dm;j++) Ke->values[(i+numdofm)*numdoftot+j]+=Ke_drag2[i][j];
+                
         /*Transform Coordinate System*/
         element->TransformStiffnessMatrixCoord(Ke,node_list,numnodes,cs_list);
@@ -3917 +3918 @@
                 vzHO_input->GetInputDerivativeValue(&dw[0],xyz_list,gauss);
 
-                for(i=0;i<vnumnodes;i++){
+                for(i=0;i<6;i++){
                         pe->values[i*3+0]+=-Jdet*gauss->weight*viscosity*dw[0]*dbasis[2][i];
                         pe->values[i*3+1]+=-Jdet*gauss->weight*viscosity*dw[1]*dbasis[2][i];