FIXMEs/TODOs for the Laursen nonlinearity

dune/tectonic/globallaursennonlinearity.hh

@@ -43,11 +43,7 @@ class GlobalLaursenNonlinearity
 
     sigma_n [(1/eta \bar Gamma)* + mu id]
 
-    TODO: We chose Gamma = id, so that (\bar Gamma)* = \Gamma^{-1}
-    = id^{-1} = id. The factor 1/eta cancels in this special case, leaving us
-    with
-
-    sigma_n [id + mu id] = sigma_n (1 + mu) id
+    TODO: We currently assume eta = 1
   */
   virtual shared_ptr<LocalNonlinearity<dim> const> restriction(int i) const {
     if ((*nodalIntegrals)[i] == 0)

src/assemblers.cc

@@ -88,6 +88,8 @@ assemble_nonlinearity(
   } else if (friction_model == std::string("Laursen")) {
     return
         // TODO: take state and h into account
+        // FIXME: We should be using a quadratic rather than a linear function
+        // here!
         Dune::make_shared<
             Dune::GlobalLaursenNonlinearity<Dune::LinearFunction, VectorType,
                                             MatrixType> const>(mu, normalStress,