A new problem

src/one-body-sample.cc

@@ -114,7 +114,7 @@ void assemble_neumann(GridView const &gridView, FEBasis const &feBasis,
   LocalVectorType SampleVector(0);
   // FIXME: random values (time-dependent)
   SampleVector[0] =
-      (time <= 0.5) ? sin(time * 2 * M_PI) * 1e3 : (time - 0.5) * 2 * 1e4;
+      (time <= 0.5) ? sin(time * 2 * M_PI) * 1e5 : (time - 0.5) * 2 * 1e5;
   SampleVector[1] = 0;
   ConstantFunction<LocalVectorType, LocalVectorType> const fNeumann(
       SampleVector);

src/one-body-sample.parset

 # -*- mode:conf -*-
-timesteps = 30
+timesteps = 300
 
 verbose = false
 printEvolution = true
@@ -7,10 +7,10 @@ printFrictionalBoundary = false
 printProgress = false
 printDifference = false
 
-writeVTK = false
+writeVTK = true
 
 [grid]
-refinements = 3
+refinements = 5
 
 [body]
 E = 1e6