Skip to content
Snippets Groups Projects
Commit 4fff09cd authored by Klaus Mosthaf's avatar Klaus Mosthaf
Browse files

Added the BCRSMatrix type as template parameter to the call of SuperLU...

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
parent 9ec29767
No related merge requests found
......@@ -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_);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment