Changeset 16785

Timestamp:

11/15/13 11:43:44 (11 years ago)

Author:

Mathieu Morlighem

Message:

BUG: fixed bug in balance velocities

File:

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

issm/trunk-jpl/src/c/classes/Elements/Tria.cpp

@@ -8611 +8611 @@
         IssmDouble*    basis  = xNew<IssmDouble>(numnodes);
         IssmDouble*    dbasis = xNew<IssmDouble>(numnodes*2);
-        IssmDouble*    HNx    = xNew<IssmDouble>(numnodes);
-        IssmDouble*    HNy    = xNew<IssmDouble>(numnodes);
         IssmDouble*    H      = xNew<IssmDouble>(numnodes);
         IssmDouble*    Nx     = xNew<IssmDouble>(numnodes);
@@ -8651 +8649 @@
                 GetNodalFunctionsDerivatives(dbasis,&xyz_list[0][0],gauss);
 
-                TriaRef::GetInputDerivativeValue(&dhnx[0],HNx,&xyz_list[0][0],gauss);
-                TriaRef::GetInputDerivativeValue(&dhny[0],HNy,&xyz_list[0][0],gauss);
-                TriaRef::GetInputValue(&hnx,HNx,gauss);
-                TriaRef::GetInputValue(&hny,HNy,gauss);
+                TriaRef::GetInputDerivativeValue(&dhnx[0],Nx,&xyz_list[0][0],gauss);
+                TriaRef::GetInputDerivativeValue(&dhny[0],Ny,&xyz_list[0][0],gauss);
+                TriaRef::GetInputValue(&hnx,Nx,gauss);
+                TriaRef::GetInputValue(&hny,Ny,gauss);
 
                 gamma=h/(2.*thickness+1.e-10);
@@ -8671 +8669 @@
         xDelete<IssmDouble>(Nx);
         xDelete<IssmDouble>(Ny);
-        xDelete<IssmDouble>(HNx);
-        xDelete<IssmDouble>(HNy);
         delete gauss;
         return pe;