Changeset 18151

Timestamp:

06/12/14 15:58:18 (11 years ago)

Author:

Mathieu Morlighem

Message:

BUG: fixed allocation required in 3d for four-field FS

File:

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

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

@@ -4410 +4410 @@
                 else{
                         _assert_(tnumnodes==4);
+                        d_zz = xNew<IssmDouble>(tnumnodes);
+                        d_xz = xNew<IssmDouble>(tnumnodes);
+                        d_yz = xNew<IssmDouble>(tnumnodes);
                         Matrix4x4Solve(&d_xx[0],Ke,pe_xx);
                         Matrix4x4Solve(&d_yy[0],Ke,pe_yy);
@@ -4495 +4498 @@
                 IssmDouble* tau_xz = NULL;
                 IssmDouble* tau_yz = NULL;
+                if(dim==3){
+                        tau_zz = xNew<IssmDouble>(tnumnodes);
+                        tau_xz = xNew<IssmDouble>(tnumnodes);
+                        tau_yz = xNew<IssmDouble>(tnumnodes);
+                }
                 Gauss* gauss = element->NewGauss();
                 for(int ig=0;ig<tnumnodes;ig++){