OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-11 13:45:34 -06:00
Code Issues Packages Projects Releases Wiki Activity

Narrow down the possible Matrix types for ParallelOverlappingILU0

Browse Source 
It now has to be a BCRSMatrix, but the block type is flexible and
needs to this way. flow_legacy uses MatrixBlock, and flow_ebos uses
FieldMatrix.
This commit is contained in:
Markus Blatt 2017-07-03 13:50:21 +02:00
parent 1bb0968283
commit 542f181f7e
1 changed files with 10 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

16
opm/autodiff/ParallelOverlappingILU0.hpp
View File

@ -165,6 +165,7 @@ class ParallelOverlappingILU0
{
typedef ParallelInfoT ParallelInfo;
public:
//! \brief The matrix type the preconditioner is for.
typedef typename Dune::remove_const<Matrix>::type matrix_type;
@ -238,8 +239,9 @@ public:
\param n ILU fill in level (for testing). This does not work in parallel.
\param w The relaxation factor.
*/
template<class Matrix1>
ParallelOverlappingILU0 (const Matrix1& A, const int n, const field_type w )
template<class BlockType, class Alloc>
ParallelOverlappingILU0 (const Dune::BCRSMatrix<BlockType,Alloc>& A,
const int n, const field_type w )
: lower_(),
upper_(),
inv_(),
@ -257,8 +259,9 @@ public:
\param A The matrix to operate on.
\param w The relaxation factor.
*/
template<class Matrix1>
ParallelOverlappingILU0 (const Matrix1& A, const field_type w)
template<class BlockType, class Alloc>
ParallelOverlappingILU0 (const Dune::BCRSMatrix<BlockType,Alloc>& A,
const field_type w)
: ParallelOverlappingILU0( A, 0, w )
{
}
@ -270,8 +273,9 @@ public:
\param comm communication object, e.g. Dune::OwnerOverlapCopyCommunication
\param w The relaxation factor.
*/
template<class Matrix1>
ParallelOverlappingILU0 (const Matrix1& A, const ParallelInfo& comm, const field_type w)
template<class BlockType, class Alloc>
ParallelOverlappingILU0 (const Dune::BCRSMatrix<BlockType,Alloc>& A,
const ParallelInfo& comm, const field_type w)
: lower_(),
upper_(),
inv_(),