Only assemble nonlinearity once

src/one-body-sample.cc

@@ -264,6 +264,10 @@ int main(int argc, char *argv[]) {
     assemble_frictional<GridType, GridType::LeafGridView, SmallVector, P1Basis>(
         leafView, p1Basis, frictionalNodes, nodalIntegrals);
 
+    Dune::shared_ptr<Dune::GlobalNonlinearity<dim> const> myGlobalNonlinearity;
+    assemble_nonlinearity(grid.size(grid.maxLevel(), dim), parset,
+                          myGlobalNonlinearity, nodalIntegrals);
+
     for (size_t run = 1; run <= parset.get<size_t>("timesteps"); ++run) {
       std::cout << ".";
       std::cout.flush();
@@ -279,12 +283,6 @@ int main(int argc, char *argv[]) {
       stiffnessMatrix.mmv(u2, b2);
       stiffnessMatrix.mmv(u3, b3);
 
-      Dune::shared_ptr<Dune::GlobalNonlinearity<dim> const>
-      myGlobalNonlinearity;
-      assemble_nonlinearity(grid.size(grid.maxLevel(), dim), parset,
-                            myGlobalNonlinearity, nodalIntegrals);
-      // }}}
-
       if (parset.get<bool>("useNonlinearGS")) {
         MyConvexProblemType myConvexProblem(
             stiffnessMatrix, *myGlobalNonlinearity, b1, u1_diff_new);