running one

dG3D/benchmarks/multiScaleBodyForce_Num/BeamMultiscale.py

@@ -62,7 +62,7 @@ microSolver.addMicroBC(microBC)
 #microSolver.setPathFollowingControlType(1)
 #microSolver.setPathFollowingArcLengthStep(1e-3)
 
-microSolver.snlData(1,1.,2e-5,1e-10)
+microSolver.snlData(1,1.,1.e-4,1e-5)
 microSolver.setSystemType(1) 
 microSolver.Solver(2)
 microSolver.Scheme(1)
@@ -75,7 +75,7 @@ microSolver.setTangentAveragingMethod(2,1e-6)
 microSolver.stiffnessModification(bool(1));
 microSolver.iterativeProcedure(bool(1));
 
-MaxIter=15
+MaxIter=25
 StepIncrease= 3
 StepReducFactor=2.
 NumberReduction=10
@@ -90,7 +90,7 @@ sol = 2  # Gmm=0 (default) Taucs=1 PETsc=2
 soltype =1 # StaticLinear=0 (default) StaticNonLinear=1
 nstep = 50  # number of step (used only if soltype=1)
 ftime =1.   # Final time (used only if soltype=1)
-tol=2.e-5  # relative tolerance for NR scheme (used only if soltype=1)
+tol=1.e-4  # relative tolerance for NR scheme (used only if soltype=1)
 nstepArch=3 # Number of step between 2 archiving (used only if soltype=1)
 
 
@@ -115,14 +115,14 @@ mysolver.addMaterialLaw(macromat)
 #mysolver.setMultiscaleFlag(bool(1))
 mysolver.Scheme(soltype)
 mysolver.Solver(sol)
-mysolver.snlData(nstep,ftime,tol)
+mysolver.snlData(nstep,ftime,tol,tol/10.)
 
 #mysolver.pathFollowing(True,0)
 #mysolver.setPathFollowingControlType(1)
 #mysolver.setPathFollowingArcLengthStep(1e-3)
 #mysolver.setPathFollowingCorrectionMethod(2)
 
-MaxIter=15
+MaxIter=25
 StepIncrease= 2
 StepReducFactor=2.
 NumberReduction=10
@@ -188,6 +188,6 @@ except:
   os._exit(1)
 else:
   check = TestCheck()
-  check.equal(3.331748e-03,float(linedisp.split(';')[1]))
+  check.equal(-1.257474e+01,float(linedisp.split(';')[1]))
 
 

dG3D/benchmarks/multiScaleBodyForce_Num/Sample6_24.geo

@@ -29,7 +29,7 @@ Physical Line(3) = {3};
 Physical Line(4) = {4};
 Physical Point(5)={3};
 
-Transfinite Line {2,4} = 5 Using Progression 1;
+Transfinite Line {2,4} = 3 Using Progression 1;
 Transfinite Line {1,3} = 2 Using Progression 1;
 Transfinite Surface {100};
 

dG3D/benchmarks/multiScaleBodyForce_Num/Sample6_24_2.msh

@@ -21,20 +21,20 @@ $PartitionedEntities
 6 0 2 1 2 59.82000000000001 0 0 0 
 7 0 3 1 1 59.82000000000001 338.98 0 1 5 
 8 0 4 1 1 0 338.98 0 0 
-9 1 2 2 1 2 59.82000000000001 169.4900000000001 0 0 
-10 1 4 2 1 2 0 169.4899999999999 0 0 
+9 1 4 2 1 2 0 169.4899999999999 0 0 
+10 1 2 2 1 2 59.82000000000001 169.4900000000001 0 0 
 5 1 1 1 2 0 0 0 59.82000000000001 0 0 1 1 2 5 -6 
-6 1 2 1 2 59.82000000000001 0 0 59.82000000000001 169.4900000000001 0 1 2 2 6 -9 
-7 1 2 1 1 59.82000000000001 169.4900000000001 0 59.82000000000001 338.98 0 1 2 2 9 -7 
+6 1 2 1 2 59.82000000000001 0 0 59.82000000000001 169.4900000000001 0 1 2 2 6 -10 
+7 1 2 1 1 59.82000000000001 169.4900000000001 0 59.82000000000001 338.98 0 1 2 2 10 -7 
 8 1 3 1 1 0 338.98 0 59.82000000000001 338.98 0 1 3 2 7 -8 
-9 1 4 1 1 0 169.4899999999999 0 0 338.98 0 1 4 2 8 -10 
-10 1 4 1 2 0 0 0 0 169.4899999999999 0 1 4 2 10 -5 
-11 2 100 2 1 2 0 169.4899999999999 0 59.82000000000001 169.4900000000001 0 0 2 10 -9 
+9 1 4 1 1 0 169.4899999999999 0 0 338.98 0 1 4 2 8 -9 
+10 1 4 1 2 0 0 0 0 169.4899999999999 0 1 4 2 9 -5 
+11 2 100 2 1 2 0 169.4899999999999 0 59.82000000000001 169.4900000000001 0 0 2 9 -10 
 101 2 100 1 2 0 0 0 59.82000000000001 169.4900000000001 0 1 51 4 5 6 10 -11 
 102 2 100 1 1 0 169.4899999999999 0 59.82000000000001 338.98 0 1 51 4 7 8 9 11 
 $EndPartitionedEntities
 $Nodes
-15 27 1 27
+15 15 1 15
 0 5 0 1
 1
 0 0 0
@@ -48,87 +48,57 @@ $Nodes
 4
 0 338.98 0
 0 9 0 1
-7
-59.82000000000001 169.4900000000001 0
-0 10 0 1
-15
+10
 0 169.4899999999999 0
+0 10 0 1
+6
+59.82000000000001 169.4900000000001 0
 1 5 0 1
 5
 29.91 0 0
-1 6 0 3
-6
-9
-10
-59.82000000000001 84.7450000000001 0
-59.82000000000001 42.37250000000005 0
-59.82000000000001 127.1175000000001 0
-1 7 0 3
+1 6 0 1
+7
+59.82000000000001 84.74500000000006 0
+1 7 0 1
 8
-11
-12
 59.82000000000001 254.2350000000001 0
-59.82000000000001 211.8625000000001 0
-59.82000000000001 296.6075000000001 0
 1 8 0 1
-13
+9
 29.91 338.98 0
-1 9 0 3
-14
-17
-18
-0 254.2349999999999 0
-0 296.6075 0
-0 211.8624999999999 0
-1 10 0 3
-16
-19
-20
-0 84.74499999999992 0
-0 127.1174999999999 0
-0 42.37249999999995 0
+1 9 0 1
+11
+0 254.235 0
+1 10 0 1
+12
+0 84.74499999999995 0
 1 11 0 1
-23
+13
 29.91 169.49 0
-2 101 0 3
-21
-22
-24
+2 101 0 1
+14
 29.91 84.745 0
-29.91 42.3725 0
-29.91 127.1175 0
-2 102 0 3
-25
-26
-27
+2 102 0 1
+15
 29.91 254.235 0
-29.91 211.8625 0
-29.91 296.6075 0
 $EndNodes
 $Elements
-9 15 1 15
+9 9 1 9
 0 7 15 1
 1 3 
 1 5 8 1
 2 1 2 5 
-1 6 8 2
-3 2 6 9 
-4 6 7 10 
-1 7 8 2
-5 7 8 11 
-6 8 3 12 
+1 6 8 1
+3 2 6 7 
+1 7 8 1
+4 6 3 8 
 1 8 8 1
-7 3 4 13 
-1 9 8 2
-8 4 14 17 
-9 14 15 18 
-1 10 8 2
-10 15 16 19 
-11 16 1 20 
-2 101 10 2
-12 1 2 6 16 5 9 21 20 22 
-13 16 6 7 15 21 10 23 19 24 
-2 102 10 2
-14 15 7 8 14 23 11 25 18 26 
-15 14 8 3 4 25 12 13 17 27 
+5 3 4 9 
+1 9 8 1
+6 4 10 11 
+1 10 8 1
+7 10 1 12 
+2 101 10 1
+8 1 2 6 10 5 7 13 12 14 
+2 102 10 1
+9 10 6 3 4 13 8 9 11 15 
 $EndElements