diff --git a/dumux/linear/pdesolver.hh b/dumux/linear/pdesolver.hh index dca3534564750aef506958e9e40d5404ccaa2ca3..8a6181153ab3b0fc658b54332dfafae7d74b0fe2 100644 --- a/dumux/linear/pdesolver.hh +++ b/dumux/linear/pdesolver.hh @@ -260,7 +260,6 @@ private: * Specialization for linear solvers that can handle MultiType matrices. * */ - template<class LS = LinearSolver, class V = SolutionVector> typename std::enable_if_t<linearSolverAcceptsMultiTypeMatrix<LS>() && isMultiTypeBlockVector<V>(), bool> @@ -270,9 +269,7 @@ private: SolutionVector& b) { assert(this->checkSizesOfSubMatrices(A) && "Sub-blocks of MultiTypeBlockMatrix have wrong sizes!"); - - // TODO: automatically derive the precondBlockLevel - return ls.template solve</*precondBlockLevel=*/2>(A, x, b); + return ls.solve(A, x, b); } /*!