diff --git a/dune/tectonic/globallaursennonlinearity.hh b/dune/tectonic/globallaursennonlinearity.hh
index 1dfc677590ca1e1f41aed25f14a30fda09bfcf91..932c6152b6eed46c6c3b81c32c05d7b752d765fb 100644
--- a/dune/tectonic/globallaursennonlinearity.hh
+++ b/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)
diff --git a/src/assemblers.cc b/src/assemblers.cc
index 383279bb56a41994e337463a0ad6131569739408..df8cffe18bd24686d71518024a9924ea963499c5 100644
--- a/src/assemblers.cc
+++ b/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,