Changeset 11563

Timestamp:

02/28/12 10:54:09 (14 years ago)

Author:

Mathieu Morlighem

Message:

Fixing analtical gradient

File:

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

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

@@ -3597 +3597 @@
         int        analysis_type;
         int        doflist1[NUMVERTICES];
+        int        connectivity[NUMVERTICES];
         double     vx,vy,lambda,mu,alpha_complement,Jdet;
         double     bed,thickness,Neff,drag;
@@ -3614 +3615 @@
         GetVerticesCoordinates(&xyz_list[0][0], nodes, NUMVERTICES);
         GradientIndexing(&doflist1[0],control_index);
+        this->GetConnectivityList(&connectivity[0]);
 
         /*Build frictoin element, needed later: */
@@ -3666 +3668 @@
         //      vx_input->GetInputValue(&vx,gauss);
         //      vy_input->GetInputValue(&vy,gauss);
-        //      grade_g[iv]=-2*drag*alpha_complement*((lambda*vx+mu*vy));
+        //      grade_g[iv] = -2*drag*alpha_complement*(lambda*vx+mu*vy)/((double)connectivity[iv]);
         //}
-        //VecSetValues(gradient,NUMVERTICES,doflist1,(const double*)grade_g,INSERT_VALUES);
         /*End Analytical gradient*/