TODO

src/one-body-sample.cc

@@ -147,7 +147,6 @@ int main(int argc, char *argv[]) {
         verbose ? Solver::FULL : Solver::QUIET;
 
     // {{{ Set up grid
-
     typedef Dune::ALUGrid<dim, dim, Dune::simplex, Dune::nonconforming>
     GridType;
     Dune::FieldVector<typename GridType::ctype, dim> lowerLeft(0);
@@ -274,6 +273,8 @@ int main(int argc, char *argv[]) {
                 : 1;
         for (size_t state_fpi = 0; state_fpi < fixed_point_iterations;
              ++state_fpi) {
+          // TODO: We should not have to reassemble the nonlinearity anew
+          // everytime
           auto myGlobalNonlinearity =
               assemble_nonlinearity<VectorType, MatrixType>(
                   parset.sub("boundary.friction"), nodalIntegrals, s4_new, h);