Changeset 11558

Timestamp:

02/28/12 09:11:09 (13 years ago)

Author:

Mathieu Morlighem

Message:

Fixed streamline upwinding in case vel=0

File:

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

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

@@ -524 +524 @@
                 if(stabilization==2){
                         /*Streamline upwinding*/
-                        vel=sqrt(pow(vx,2.)+pow(vy,2.));
+                        vel=sqrt(pow(vx,2.)+pow(vy,2.))+1.e-8;
                         K[0][0]=h/(2*vel)*vx*vx;
                         K[1][0]=h/(2*vel)*vy*vx;
@@ -3667 +3667 @@
         //      vy_input->GetInputValue(&vy,gauss);
         //      grade_g[iv]=-2*drag*alpha_complement*((lambda*vx+mu*vy));
-        //      VecSetValues(gradient,NUMVERTICES,doflist1,(const double*)grade_g,INSERT_VALUES);
         //}
+        //VecSetValues(gradient,NUMVERTICES,doflist1,(const double*)grade_g,INSERT_VALUES);
         /*End Analytical gradient*/