Use Dune::ScaledIdentityMatrix<field_type, blocksize> instead of Dune::FieldMatrix<field_type,1,1>

dune/elasticity/common/trustregionsolver.hh

@@ -62,7 +62,7 @@ class TrustRegionSolver
 
     // Some types that I need
     typedef Dune::BCRSMatrix<Dune::FieldMatrix<field_type, blocksize, blocksize> > MatrixType;
-    typedef Dune::BCRSMatrix<Dune::FieldMatrix<field_type,1,1> >                   ScalarMatrixType;
+    typedef Dune::BCRSMatrix<Dune::ScaledIdentityMatrix<field_type, blocksize>>    ScalarMatrixType;
     typedef Dune::BlockVector<Dune::FieldVector<field_type, blocksize> >           CorrectionType;
     typedef VectorType                                                             SolutionType;