Changeset 17371

Timestamp:

03/03/14 10:30:35 (11 years ago)

Author:

jbondzio

Message:

BUG: replaced static arrays by dynamic ones

File:

• issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp (modified) (4 diffs)

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

@@ -98 +98 @@
         IssmDouble Jdet, dt, D_scalar;
         IssmDouble h,hx,hy,hz;
-        IssmDouble vel;
-        IssmDouble calvingrate;
-        IssmDouble dlsf[dim], norm_dlsf, normal[dim];
+        IssmDouble vel, calvingrate;
+        IssmDouble norm_dlsf;
         IssmDouble* xyz_list = NULL;
 
@@ -115 +114 @@
         IssmDouble*    w        = xNew<IssmDouble>(dim);
         IssmDouble*    c        = xNew<IssmDouble>(dim);
+        IssmDouble*    dlsf     = xNew<IssmDouble>(dim);
 
         /*Retrieve all inputs and parameters*/
@@ -156 +156 @@
                 norm_dlsf=sqrt(norm_dlsf);
 
-                if(norm_dlsf>1.e-10){
-                        for(i=0;i<dim;i++) normal[i]=dlsf[i]/norm_dlsf;
-                        for(i=0;i<dim;i++) c[i]=calvingrate*normal[i];
-                }
+                if(norm_dlsf>1.e-10)
+                        for(i=0;i<dim;i++) c[i]=calvingrate*dlsf[i]/norm_dlsf;
                 else
                         for(i=0;i<dim;i++) c[i]=0.;
@@ -230 +228 @@
         xDelete<IssmDouble>(w);
         xDelete<IssmDouble>(c);
+        xDelete<IssmDouble>(dlsf);
         delete gauss;
         return Ke;