Changeset 81

Timestamp:

04/28/09 13:56:46 (16 years ago)

Author:

seroussi

Message:

fixed overshooting in parallel

File:

• issm/trunk/src/c/parallel/diagnostic_core_nonlinear.cpp (modified) (2 diffs)

issm/trunk/src/c/parallel/diagnostic_core_nonlinear.cpp

@@ -63 +63 @@
                 if (debug) _printf_("   Updating inputs\n");
 
+                /*Set input parameters: */
+                if(old_ug)ParameterInputsAddFromVec(inputs,old_ug,"old_velocity");
+                if(ug)ParameterInputsAddFromVec(inputs,ug,"velocity");
+
                 //save pointer to old velocity
                 VecFree(&old_ug);old_ug=ug;
                 VecFree(&old_uf);old_uf=uf;
 
-                /*Set input parameters: */
-                if(ug)ParameterInputsAddFromVec(inputs,ug,"velocity");
 
                 /*Update parameters: */
@@ -111 +113 @@
                 if (debug) _printf_("   penalty constraints\n");
                 ParameterInputsAddFromVec(inputs,ug,"velocity");
+                ParameterInputsAddFromVec(inputs,old_ug,"old_velocity");
                 
                 PenaltyConstraintsx(&constraints_converged, &num_unstable_constraints, fem->elements,fem->nodes,loads,fem->materials,inputs,analysis_type);