Changeset 4339

Timestamp:

06/30/10 10:56:05 (15 years ago)

Author:

seroussi

Message:

fixed problem with transient when overshoot:
make sure that vx and vxold are similar once the solution has converged

Location:

issm/trunk/src

Files:

• c/solvers/solver_diagnostic_nonlinear.cpp (modified) (1 diff)
• m/solvers/solver_diagnostic_nonlinear.m (modified) (1 diff)

issm/trunk/src/c/solvers/solver_diagnostic_nonlinear.cpp

@@ -125 +125 @@
         }
         
+        //Now that the solution has converged, update once again the solution to make sure that vx and vxold are similar (for next step in transient or steadystate)
+        InputUpdateFromSolutionx( femmodel->elements,femmodel->nodes, femmodel->vertices, femmodel->loads, femmodel->materials, femmodel->parameters,ug);
+
         /*extrude if we are in 3D: */
         if (dim==3){

issm/trunk/src/m/solvers/solver_diagnostic_nonlinear.m

@@ -79 +79 @@
         end
         
+        %Now that the solution has converged, update once again the solution to make sure that vx and vxold are similar
+        femmodel.elements=InputUpdateFromSolution(femmodel.elements,femmodel.nodes,femmodel.vertices,loads,femmodel.materials,femmodel.parameters,ug);
+
         %extrude if we are in 3D: 
         if dim==3,