From 4fff09cd11bdb40db0622b299809786ecc5ea7a5 Mon Sep 17 00:00:00 2001
From: Klaus Mosthaf <klmos@env.dtu.dk>
Date: Thu, 28 Jul 2011 14:18:20 +0000
Subject: [PATCH] Added the BCRSMatrix type as template parameter to the call
of SuperLU (required for the coupled models)
git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@6324 2fb0f335-1f38-0410-981e-8018bf24f1b0
---
dumux/linear/seqsolverbackend.hh | 8 +++++++-
1 file changed, 7 insertions(+), 1 deletion(-)
diff --git a/dumux/linear/seqsolverbackend.hh b/dumux/linear/seqsolverbackend.hh
index 7b2a65ae0c..06d0417e32 100644
--- a/dumux/linear/seqsolverbackend.hh
+++ b/dumux/linear/seqsolverbackend.hh
@@ -377,7 +377,13 @@ public:
static const double residReduction = GET_PROP_VALUE(TypeTag, PTAG(LSResidualReduction));
Vector bTmp(b);
- Dune::SuperLU<Matrix> solver(A);
+
+ typedef typename GET_PROP_TYPE(TypeTag, PTAG(Scalar)) Scalar;
+ enum {numEq = GET_PROP_VALUE(TypeTag, PTAG(NumEq))};
+ typedef typename Dune::FieldMatrix<Scalar, numEq, numEq> MatrixBlock;
+ typedef typename Dune::BCRSMatrix<MatrixBlock> ISTLMatrix;
+
+ Dune::SuperLU<ISTLMatrix> solver(A);
solver.apply(x, bTmp, result_);
--
GitLab