My Project
|
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...
Public Types | |
using | AssembledLinearOperatorType = Dune::AssembledLinearOperator< Matrix, Vector, Vector > |
Public Member Functions | |
ISTLSolverEbos (const Simulator &simulator) | |
Construct a system solver. | |
void | eraseMatrix () |
void | prepare (const SparseMatrixAdapter &M, Vector &b) |
void | setResidual (Vector &) |
void | getResidual (Vector &b) const |
void | setMatrix (const SparseMatrixAdapter &) |
bool | solve (Vector &x) |
int | iterations () const |
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 |
Static Public Member Functions | |
static void | registerParameters () |
Protected Types | |
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::EclWellModel > |
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 | CommunicationType = Dune::CollectiveCommunication< int > |
Protected Member Functions | |
void | 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. | |
Vector | getTrueImpesWeights (int pressureVarIndex) const |
Matrix & | getMatrix () |
const Matrix & | getMatrix () const |
Protected Attributes | |
const Simulator & | simulator_ |
int | iterations_ |
int | calls_ |
bool | converged_ |
std::any | parallelInformation_ |
Matrix * | matrix_ |
Vector * | rhs_ |
detail::FlexibleSolverInfo< Matrix, Vector, CommunicationType > | flexibleSolver_ |
std::vector< int > | overlapRows_ |
std::vector< int > | interiorRows_ |
bool | useWellConn_ |
FlowLinearSolverParameters | parameters_ |
PropertyTree | prm_ |
std::shared_ptr< CommunicationType > | comm_ |
Static Protected Attributes | |
static constexpr std::size_t | pressureIndex = GetPropType<TypeTag, Properties::Indices>::pressureSwitchIdx |
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 .
|
inlineexplicit |
Construct a system solver.
[in] | parallelInformation | In the case of a parallel run with dune-istl the information about the parallelization. |
|
inline |
Solve the system of linear equations Ax = b, with A being the combined derivative matrix of the residual and b being the residual itself.
[in] | residual | residual object containing A and b. |
|
inline |