From 4a2d1514741801a0ab2cc371469e13b86c921de5 Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Fri, 25 Jan 2013 14:33:16 +0100
Subject: [PATCH] Use a ScaledIdentityMatrix
---
src/one-body-sample.cc | 4 +++-
1 file changed, 3 insertions(+), 1 deletion(-)
diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index 72f47a00..b959bfa3 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -44,6 +44,7 @@
#include <dune/grid/utility/structuredgridfactory.hh>
#include <dune/istl/bcrsmatrix.hh>
#include <dune/istl/bvector.hh>
+#include <dune/istl/scaledidmatrix.hh>
#include <dune/fufem/assemblers/functionalassembler.hh>
#include <dune/fufem/assemblers/localassemblers/l2functionalassembler.hh>
@@ -240,7 +241,8 @@ int main(int argc, char *argv[]) {
{
MassAssembler<GridType, P1Basis::LocalFiniteElement,
- P1Basis::LocalFiniteElement, SmallMatrix> const localMass;
+ P1Basis::LocalFiniteElement,
+ Dune::ScaledIdentityMatrix<double, dims>> const localMass;
OperatorAssembler<P1Basis, P1Basis>(p1Basis, p1Basis)
.assemble(localMass, massMatrix);
massMatrix *= density;
--
GitLab