switch to multi threading

dune/tectonic/time-stepping/adaptivetimestepper.cc

@@ -78,7 +78,7 @@ int AdaptiveTimeStepper<Factory, ContactNetwork, Updaters, ErrorNorms>::coarsen(
 
       setDeformation(current_);
       auto C_Step = Step(stepBase_, current_, currentNBodyAssembler, relativeTime_, 2 * relativeTau_, iterationRegister_);
-      C_Step.run(Step::Mode::sameThread); //newThread
+      C_Step.run(Step::Mode::newThread); //newThread
 
       //updateReductionFactors(reductionFactors);
       //std::cout << "AdaptiveTimeStepper C computed!" << std::endl << std::endl;
@@ -92,7 +92,7 @@ int AdaptiveTimeStepper<Factory, ContactNetwork, Updaters, ErrorNorms>::coarsen(
       //auto R2_linearSolver = makeLinearSolver();
       auto&& nBodyAssembler = step(currentNBodyAssembler);
       auto R2_Step = Step(stepBase_, R1_.updaters, nBodyAssembler, relativeTime_ + relativeTau_, relativeTau_, iterationRegister_);
-      R2_Step.run(Step::Mode::sameThread); //newThread
+      R2_Step.run(Step::Mode::newThread); //newThread
 
       //updateReductionFactors(reductionFactors);
       //std::cout << "AdaptiveTimeStepper R2 computed!" << std::endl << std::endl;
@@ -195,7 +195,7 @@ int AdaptiveTimeStepper<Factory, ContactNetwork, Updaters, ErrorNorms>::determin
     setDeformation(current_);
     //auto C_linearSolver = makeLinearSolver();
     auto C_Step = Step(stepBase_, current_, currentNBodyAssembler, relativeTime_, 2 * relativeTau_, iterationRegister_);
-    C_Step.run(Step::Mode::sameThread); // newThread
+    C_Step.run(Step::Mode::newThread); // newThread
     //updateReductionFactors(reductionFactors);
     //std::cout << "AdaptiveTimeStepper C computed!" << std::endl << std::endl;
 
@@ -203,7 +203,7 @@ int AdaptiveTimeStepper<Factory, ContactNetwork, Updaters, ErrorNorms>::determin
     //auto R2_linearSolver = makeLinearSolver();
     auto&& nBodyAssembler = step(currentNBodyAssembler);
     auto R2_Step = Step(stepBase_, R1_.updaters, nBodyAssembler, relativeTime_ + relativeTau_, relativeTau_, iterationRegister_);
-    R2_Step.run(Step::Mode::sameThread); //newThread
+    R2_Step.run(Step::Mode::newThread); //newThread
 
     //updateReductionFactors(reductionFactors);
     //std::cout << "AdaptiveTimeStepper R2 computed!" << std::endl << std::endl;
@@ -221,7 +221,7 @@ int AdaptiveTimeStepper<Factory, ContactNetwork, Updaters, ErrorNorms>::determin
         setDeformation(current_);
         //auto F1_linearSolver = makeLinearSolver();
         auto F1_Step = Step(stepBase_, current_, currentNBodyAssembler, relativeTime_, relativeTau_ / 2.0, iterationRegister_);
-        F1_Step.run(Step::Mode::sameThread); //newThread
+        F1_Step.run(Step::Mode::newThread); //newThread
         //updateReductionFactors(reductionFactors);
         //std::cout << "AdaptiveTimeStepper F1 computed!" << std::endl << std::endl;
 

src/strikeslip/strikeslip-2D.cfg

@@ -6,7 +6,7 @@ smallestDiameter = 0.02 # 2e-3 [m]
 smallestDiameter = 0.02  # 2e-3 [m]
 
 [timeSteps]
-refinementTolerance = 1e-5 # 1e-5
+refinementTolerance = 1e-3# 1e-5
 
 [u0.solver]
 tolerance         = 1e-8
@@ -18,7 +18,7 @@ tolerance         = 1e-8
 tolerance         = 1e-8
 
 [v.fpi]
-tolerance         = 1e-5    # 1e-5
+tolerance         = 1e-3    # 1e-5
 
 [solver.tnnmg.preconditioner.basesolver]
 tolerance          = 1e-10

src/strikeslip/strikeslip.cfg

@@ -69,11 +69,11 @@ bodyCount       = 2
 [initialTime]
 timeStep = 0
 relativeTime = 0.0
-relativeTau = 2e-4 # 1e-6
+relativeTau = 2e-2 # 1e-6
 
 [timeSteps]
 scheme = newmark
-timeSteps = 5000
+timeSteps = 5
 
 [u0.solver]
 maximumIterations = 100