better convergence

dG3D/benchmarks/pathFollowing_TrueSnapback/twoHole.py

@@ -35,9 +35,9 @@ myfield1.stabilityParameters(beta1)
 # solver
 sol = 2  # Gmm=0 (default) Taucs=1 PETsc=2
 soltype =1 # StaticLinear=0 (default) StaticNonLinear=1
-nstep = 38  # number of step (used only if soltype=1)
+nstep = 30  # number of step (used only if soltype=1)
 ftime =1.   # Final time (used only if soltype=1)
-tol=1.e-4 # relative tolerance for NR scheme (used only if soltype=1)
+tol=1.e-5 # relative tolerance for NR scheme (used only if soltype=1)
 nstepArch=1 # Number of step between 2 archiving (used only if soltype=1)
 
 # creation of Solver
@@ -89,4 +89,4 @@ mysolver.archivingForceOnPhysicalGroup('Face',12,1)
 mysolver.solve()
 
 check = TestCheck()
-check.equal(-4.239069e+03,mysolver.getArchivedForceOnPhysicalGroup("Face", 12, 1),1.e-4)
+check.equal(-9.223652e+02,mysolver.getArchivedForceOnPhysicalGroup("Face", 12, 1),1.e-4)

dG3D/benchmarks/pathFollowing_snapback/twoHole.py

@@ -35,9 +35,9 @@ myfield1.stabilityParameters(beta1)
 # solver
 sol = 2  # Gmm=0 (default) Taucs=1 PETsc=2
 soltype =1 # StaticLinear=0 (default) StaticNonLinear=1
-nstep = 40  # number of step (used only if soltype=1)
+nstep = 35  # number of step (used only if soltype=1)
 ftime =1.   # Final time (used only if soltype=1)
-tol=1.e-3 # relative tolerance for NR scheme (used only if soltype=1)
+tol=1.e-5 # relative tolerance for NR scheme (used only if soltype=1)
 nstepArch=1 # Number of step between 2 archiving (used only if soltype=1)
 
 # creation of Solver
@@ -51,7 +51,7 @@ mysolver.snlData(nstep,ftime,tol)
 
 mysolver.pathFollowing(True,1)
 # time-step adaptation by number of NR iterations
-mysolver.setPathFollowingIncrementAdaptation(True,3)
+mysolver.setPathFollowingIncrementAdaptation(True,4)
 mysolver.setPathFollowingLocalSteps(1e-2,1e-6)
 mysolver.setPathFollowingLocalIncrementType(1); 
 
@@ -89,4 +89,4 @@ mysolver.archivingForceOnPhysicalGroup('Face',12,1)
 mysolver.solve()
 
 check = TestCheck()
-check.equal(-5.452552e+03,mysolver.getArchivedForceOnPhysicalGroup("Face", 12, 1),1.e-4)
+check.equal(-2.256506e+03,mysolver.getArchivedForceOnPhysicalGroup("Face", 12, 1),1.e-4)