diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc index ca38283c038986ecdf086eed0127b9f6edace010..9ea4016b8468f7fe5c0d8d97fc9179433dc95aaa 100644 --- a/src/one-body-sample.cc +++ b/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);