From 0286093a3d9384bdd3877bb3b2d62b4faf27d8ed Mon Sep 17 00:00:00 2001
From: Van Dung Nguyen <vdg.nguyen@gmail.com>
Date: Sat, 26 Aug 2017 22:06:30 +0200
Subject: [PATCH] better convergence
---
dG3D/benchmarks/pathFollowing_TrueSnapback/twoHole.py | 6 +++---
dG3D/benchmarks/pathFollowing_snapback/twoHole.py | 8 ++++----
2 files changed, 7 insertions(+), 7 deletions(-)
diff --git a/dG3D/benchmarks/pathFollowing_TrueSnapback/twoHole.py b/dG3D/benchmarks/pathFollowing_TrueSnapback/twoHole.py
index b2adab7c4..655b6439c 100644
--- a/dG3D/benchmarks/pathFollowing_TrueSnapback/twoHole.py
+++ b/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)
diff --git a/dG3D/benchmarks/pathFollowing_snapback/twoHole.py b/dG3D/benchmarks/pathFollowing_snapback/twoHole.py
index 52e5fcb0c..9b4dfe0f4 100644
--- a/dG3D/benchmarks/pathFollowing_snapback/twoHole.py
+++ b/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)
--
GitLab