introduce minTau for adaptive time stepping

src/multi-body-problem/multi-body-problem.cc

@@ -362,9 +362,10 @@ int main(int argc, char *argv[]) {
         timeStepper(stepBase, contactNetwork, current,
                             programState.relativeTime, programState.relativeTau);*/
 
+    const auto minTau = parset.get<double>("initialTime.minRelativeTau");
     AdaptiveTimeStepper<NonlinearFactory, std::decay_t<decltype(contactNetwork)>, Updaters, std::decay_t<decltype(stateEnergyNorms)>>
         timeStepper(stepBase, contactNetwork, current,
-                            programState.relativeTime, programState.relativeTau,
+                            programState.relativeTime, programState.relativeTau, minTau,
                             mustRefine);
 
     size_t timeSteps = std::round(parset.get<double>("timeSteps.timeSteps"));

src/multi-body-problem/multi-body-problem.cfg

 # -*- mode:conf -*-
 [general]
-outPath = newmark # output written to ./output/outPath
+outPath = 5-body # output written to ./output/outPath
 gravity         = 9.81  # [m/s^2]
 
 
@@ -48,7 +48,7 @@ finalVelocity   = 2e-4      # [m/s]
 timeStep = 0
 relativeTime = 0.0
 relativeTau = 1e-4 # 1e-6
-
+minRelativeTau = 0.0
 
 [timeSteps]
 scheme = newmark
@@ -57,7 +57,7 @@ timeSteps = 1e6
 
 [problem]
 finalTime       = 100  # [s] #1000
-bodyCount       = 4
+bodyCount       = 5
 
 
 [io]
@@ -66,7 +66,7 @@ printProgress   = true
 restarts.first  = 0
 restarts.spacing= 50 #20
 restarts.write  = true #true
-vtk.write       = true
+vtk.write       = false
 
 
 [u0.solver]