This class solves the fully implicit black-oil system by solving the reduced system (after eliminating well variables) as a block-structured matrix (one block for all cell variables) for a fixed number of cell variables np . More...

#include <ISTLSolver.hpp>

Inheritance diagram for Opm::ISTLSolver< TypeTag >:

Public Types
using	AssembledLinearOperatorType = Dune::AssembledLinearOperator< Matrix, Vector, Vector >
Public Types inherited from Opm::AbstractISTLSolver< GetPropType< TypeTag, Properties::SparseMatrixAdapter >, GetPropType< TypeTag, Properties::GlobalEqVector > >
using	CommunicationType
using	Matrix

Public Member Functions
	ISTLSolver (const Simulator &simulator, const FlowLinearSolverParameters &parameters, bool forceSerial=false)
	Construct a system solver.
	ISTLSolver (const Simulator &simulator)
	Construct a system solver.
void	initialize ()
void	eraseMatrix () override
	Signals that the memory for the matrix internally in the solver could be erased.
void	setActiveSolver (const int num) override
	Set the active solver by its index.
int	numAvailableSolvers () const override
	Get the number of available solvers.
void	initPrepare (const Matrix &M, Vector &b)
void	prepare (const SparseMatrixAdapter &M, Vector &b) override
void	prepare (const Matrix &M, Vector &b) override
void	setResidual (Vector &) override
void	getResidual (Vector &b) const override
void	setMatrix (const SparseMatrixAdapter &) override
int	getSolveCount () const override
	Get the count of how many times the solver has been called.
void	resetSolveCount ()
bool	solve (Vector &x) override
int	iterations () const override
	Solve the system of linear equations Ax = b, with A being the combined derivative matrix of the residual and b being the residual itself.
const std::any &	parallelInformation () const
const CommunicationType *	comm () const override
	Get the communication object used by the solver.
void	setDomainIndex (const int index)
bool	isNlddLocalSolver () const
Public Member Functions inherited from Opm::AbstractISTLSolver< GetPropType< TypeTag, Properties::SparseMatrixAdapter >, GetPropType< TypeTag, Properties::GlobalEqVector > >
virtual void	prepare (const Matrix &M, GetPropType< TypeTag, Properties::GlobalEqVector > &b)=0
	Prepare the solver with the given matrix and right-hand side vector.
virtual void	setResidual (GetPropType< TypeTag, Properties::GlobalEqVector > &b)=0
	Set the residual vector.
virtual void	getResidual (GetPropType< TypeTag, Properties::GlobalEqVector > &b) const=0
	Get the residual vector.
virtual void	setMatrix (const SparseMatrixAdapter &M)=0
	Set the matrix for the solver.
virtual bool	solve (GetPropType< TypeTag, Properties::GlobalEqVector > &x)=0
	Solve the system of equations Ax = b.

Static Public Member Functions
static void	registerParameters ()

Protected Types
enum	{ enablePolymerMolarWeight = getPropValue<TypeTag, Properties::EnablePolymerMW>() }
using	GridView = GetPropType<TypeTag, Properties::GridView>
using	Scalar = GetPropType<TypeTag, Properties::Scalar>
using	SparseMatrixAdapter = GetPropType<TypeTag, Properties::SparseMatrixAdapter>
using	Vector = GetPropType<TypeTag, Properties::GlobalEqVector>
using	Indices = GetPropType<TypeTag, Properties::Indices>
using	WellModel = GetPropType<TypeTag, Properties::WellModel>
using	Simulator = GetPropType<TypeTag, Properties::Simulator>
using	Matrix = typename SparseMatrixAdapter::IstlMatrix
using	ThreadManager = GetPropType<TypeTag, Properties::ThreadManager>
using	ElementContext = GetPropType<TypeTag, Properties::ElementContext>
using	AbstractSolverType = Dune::InverseOperator<Vector, Vector>
using	AbstractOperatorType = Dune::AssembledLinearOperator<Matrix, Vector, Vector>
using	AbstractPreconditionerType = Dune::PreconditionerWithUpdate<Vector, Vector>
using	WellModelOperator = WellModelAsLinearOperator<WellModel, Vector, Vector>
using	ElementMapper = GetPropType<TypeTag, Properties::ElementMapper>
using	ElementChunksType = ElementChunks<GridView, Dune::Partitions::All>
using	CommunicationType = Dune::Communication<int>

Protected Member Functions
bool	checkConvergence (const Dune::InverseOperatorResult &result) const
bool	isParallel () const
void	prepareFlexibleSolver ()
bool	shouldCreateSolver () const
	Return true if we should (re)create the whole solver, instead of just calling update() on the preconditioner.
std::function< Vector()>	getWeightsCalculator (const PropertyTree &prm, const Matrix &matrix, std::size_t pressIndex) const
Matrix &	getMatrix ()
const Matrix &	getMatrix () const

Protected Attributes
const Simulator &	simulator_
int	iterations_
int	solveCount_
std::any	parallelInformation_
Matrix *	matrix_
Vector *	rhs_
int	activeSolverNum_ = 0
std::vector< detail::FlexibleSolverInfo< Matrix, Vector, CommunicationType > >	flexibleSolver_
std::vector< int >	overlapRows_
std::vector< int >	interiorRows_
int	domainIndex_ = -1
bool	useWellConn_
std::vector< FlowLinearSolverParameters >	parameters_
bool	forceSerial_ = false
std::vector< PropertyTree >	prm_
std::shared_ptr< CommunicationType >	comm_
std::unique_ptr< ElementChunksType >	element_chunks_

Static Protected Attributes
static constexpr std::size_t	pressureIndex = GetPropType<TypeTag, Properties::Indices>::pressureSwitchIdx
static constexpr bool	isIncompatibleWithCprw = enablePolymerMolarWeight

Additional Inherited Members
Static Protected Member Functions inherited from Opm::AbstractISTLSolver< GetPropType< TypeTag, Properties::SparseMatrixAdapter >, GetPropType< TypeTag, Properties::GlobalEqVector > >
static bool	checkConvergence (const Dune::InverseOperatorResult &result, const FlowLinearSolverParameters &parameters)
	Check the convergence of the linear solver.

Detailed Description

template<class TypeTag>
class Opm::ISTLSolver< TypeTag >

This class solves the fully implicit black-oil system by solving the reduced system (after eliminating well variables) as a block-structured matrix (one block for all cell variables) for a fixed number of cell variables np .

Constructor & Destructor Documentation

◆ ISTLSolver() [1/2]

template<class TypeTag>

Opm::ISTLSolver< TypeTag >::ISTLSolver	(	const Simulator &	simulator,
		const FlowLinearSolverParameters &	parameters,
		bool	forceSerial = false )

inline

Construct a system solver.

Parameters

[in]	simulator	The opm-models simulator object
[in]	parameters	Explicit parameters for solver setup, do not read them from command line parameters.
[in]	forceSerial	If true, will set up a serial linear solver only, local to the current rank, instead of creating a parallel (MPI distributed) linear solver.

◆ ISTLSolver() [2/2]

template<class TypeTag>

Opm::ISTLSolver< TypeTag >::ISTLSolver ( const Simulator & simulator )

inlineexplicit

Construct a system solver.

Parameters

[in] simulator The opm-models simulator object

Member Function Documentation

◆ comm()

template<class TypeTag>

const CommunicationType * Opm::ISTLSolver< TypeTag >::comm ( ) const

inlineoverridevirtual

Get the communication object used by the solver.

This method returns a pointer to the communication object used by the solver.

Returns: A pointer to the communication object.

Implements Opm::AbstractISTLSolver< GetPropType< TypeTag, Properties::SparseMatrixAdapter >, GetPropType< TypeTag, Properties::GlobalEqVector > >.

◆ eraseMatrix()

template<class TypeTag>

void Opm::ISTLSolver< TypeTag >::eraseMatrix ( )

inlineoverridevirtual

Signals that the memory for the matrix internally in the solver could be erased.

Note: This call could be ignored by the solver, but it is a hint that the solver does not need the matrix anymore.

Implements Opm::AbstractISTLSolver< GetPropType< TypeTag, Properties::SparseMatrixAdapter >, GetPropType< TypeTag, Properties::GlobalEqVector > >.

◆ getSolveCount()

template<class TypeTag>

int Opm::ISTLSolver< TypeTag >::getSolveCount ( ) const

inlineoverridevirtual

Get the count of how many times the solver has been called.

This method returns the number of times the solve() method has been called.

Returns: The count of solve calls.

Implements Opm::AbstractISTLSolver< GetPropType< TypeTag, Properties::SparseMatrixAdapter >, GetPropType< TypeTag, Properties::GlobalEqVector > >.

◆ iterations()

template<class TypeTag>

int Opm::ISTLSolver< TypeTag >::iterations ( ) const

inlineoverridevirtual

Solve the system of linear equations Ax = b, with A being the combined derivative matrix of the residual and b being the residual itself.

Parameters

[in] residual residual object containing A and b.

Returns: the solution x

Implements Opm::AbstractISTLSolver< GetPropType< TypeTag, Properties::SparseMatrixAdapter >, GetPropType< TypeTag, Properties::GlobalEqVector > >.

◆ numAvailableSolvers()

template<class TypeTag>

int Opm::ISTLSolver< TypeTag >::numAvailableSolvers ( ) const

inlineoverridevirtual

Get the number of available solvers.

Returns: The number of solvers that can be used.

Implements Opm::AbstractISTLSolver< GetPropType< TypeTag, Properties::SparseMatrixAdapter >, GetPropType< TypeTag, Properties::GlobalEqVector > >.

◆ parallelInformation()

template<class TypeTag>

const std::any & Opm::ISTLSolver< TypeTag >::parallelInformation ( ) const

inline

◆ setActiveSolver()

template<class TypeTag>

void Opm::ISTLSolver< TypeTag >::setActiveSolver ( const int num )

inlineoverridevirtual

Set the active solver by its index.

Parameters

num	The index of the solver to set as active.

Note: The index corresponds to the order in which solvers are registered.

Implements Opm::AbstractISTLSolver< GetPropType< TypeTag, Properties::SparseMatrixAdapter >, GetPropType< TypeTag, Properties::GlobalEqVector > >.

The documentation for this class was generated from the following files:

opm/simulators/linalg/FlowLinearSolverParameters.hpp
opm/simulators/linalg/ISTLSolver.hpp

Public Types

Public Member Functions

Static Public Member Functions

Protected Types

Protected Member Functions

Protected Attributes

Static Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ ISTLSolver() [1/2]

◆ ISTLSolver() [2/2]

Member Function Documentation

◆ comm()

◆ eraseMatrix()

◆ getSolveCount()

◆ iterations()

◆ numAvailableSolvers()

◆ parallelInformation()

◆ setActiveSolver()