Changeset 398

Timestamp:

05/13/09 15:45:43 (16 years ago)

Author:

Mathieu Morlighem

Message:

bad gradient sign

File:

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

issm/trunk/src/c/objects/Tria.cpp

@@ -1674 +1674 @@
                 gauss_l1l2l3[2]=*(third_gauss_area_coord+ig);
 
-                
-                /*Compute velocities at gaussian point: */
-                GetParameterValue(&velocity_x, &vx_list[0],gauss_l1l2l3);
-                GetParameterValue(&velocity_y, &vy_list[0],gauss_l1l2l3);
-                #ifdef _DEBUG_ 
-                        printf("Velocity: %lf %lf\n", velocity_x,velocity_y);
-                #endif
-        
-                /*Compute obs_velocities at gaussian point: */
-                GetParameterValue(&obs_velocity_x, &obs_vx_list[0],gauss_l1l2l3);
-                GetParameterValue(&obs_velocity_y, &obs_vy_list[0],gauss_l1l2l3);
-                #ifdef _DEBUG_ 
-                        printf("Observed velocity: %g %g\n", obs_velocity_x,obs_velocity_y);
-                #endif
-
                 /* Get Jacobian determinant: */
                 GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss_l1l2l3);
@@ -1948 +1933 @@
                 /*Build gradje_g_gaussian vector (actually -dJ/ddrag): */
                 for (i=0;i<numgrids;i++){
-                        grade_g_gaussian[i]=2*drag*alpha_complement*( (lambda*vx+mu*vy))*Jdet*gauss_weight*l1l2l3[i]; 
+                        grade_g_gaussian[i]=-2*drag*alpha_complement*((lambda*vx+mu*vy))*Jdet*gauss_weight*l1l2l3[i]; 
                 }