add Hessian indices

dune/tectonic/globalnonlinearity.hh

@@ -27,6 +27,11 @@ class GlobalNonlinearity {
     }
   }
 
+  void addHessianIndices(Dune::MatrixIndexSet &indices) const {
+    for (size_t i = 0; i < indices.rows(); ++i)
+      indices.add(i, i);
+  }
+
   void addGradient(VectorType const &v, VectorType &gradient) const {
     for (size_t i = 0; i < v.size(); ++i) {
       auto res = restriction(i);

dune/tectonic/myblockproblem.hh

@@ -91,8 +91,7 @@ template <class MyConvexProblemTypeTEMPLATE> class MyBlockProblem {
     // construct sparsity pattern for linearization
     Dune::MatrixIndexSet indices(problem.A.N(), problem.A.M());
     indices.import(problem.A);
-    // TODO: handle smooth domain
-    // problem.phi.addHessianIndices(indices);
+    problem.phi.addHessianIndices(indices);
 
     // construct matrix from pattern and initialize it
     indices.exportIdx(linearization.A);