From 474ccfb8c3331f6f2337d8c587558b95c3c90d33 Mon Sep 17 00:00:00 2001
From: Lisa Julia Nebel <lisa_julia.nebel@tu-dresden.de>
Date: Wed, 18 Jan 2023 11:55:36 +0100
Subject: [PATCH] Use Dune::ScaledIdentityMatrix<field_type, blocksize> instead
of Dune::FieldMatrix<field_type,1,1>
---
dune/elasticity/common/trustregionsolver.hh | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/dune/elasticity/common/trustregionsolver.hh b/dune/elasticity/common/trustregionsolver.hh
index c80b181..ef4ef03 100644
--- a/dune/elasticity/common/trustregionsolver.hh
+++ b/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;
--
GitLab