diff --git a/src/averagedistanceassembler.hh b/src/averagedistanceassembler.hh
index 0f62f9846c5913e56f35b8948a31ee86e2a1c642..0dfd06d24e2f8ed1660cfdc00711d2f323ab0c71 100644
--- a/src/averagedistanceassembler.hh
+++ b/src/averagedistanceassembler.hh
@@ -45,8 +45,8 @@ public:
 #endif
     }
 
-    void assembleMatrix(const TargetSpace& x,
-                        Dune::FieldMatrix<double,size,size>& matrix) const
+    void assembleHessianApproximation(const TargetSpace& x,
+                                      Dune::FieldMatrix<double,size,size>& matrix) const
     {
         for (int i=0; i<size; i++)
             for (int j=0; j<size; j++)
diff --git a/src/targetspacertrsolver.cc b/src/targetspacertrsolver.cc
index f343b1ad9da1911f55913dfdc6dff32293cb1b5c..9cb11266fcb76b902f68d1139eafefb1b95b0a07 100644
--- a/src/targetspacertrsolver.cc
+++ b/src/targetspacertrsolver.cc
@@ -77,7 +77,7 @@ void TargetSpaceRiemannianTRSolver<TargetSpace>::solve()
         MatrixType hesseMatrix(1,1);
 
         assembler_->assembleGradient(x_, rhs[0]);
-        assembler_->assembleMatrix(x_, hesseMatrix[0][0]);
+        assembler_->assembleHessianApproximation(x_, hesseMatrix[0][0]);
 
         //gradientFDCheck(x_, rhs, *rodAssembler_);
         //hessianFDCheck(x_, *hessianMatrix_, *rodAssembler_);