The common code for the linearizers of non-linear systems of equations. More...

#include <tpfalinearizer.hh>

Public Member Functions

 TpfaLinearizer ()
 
 ~TpfaLinearizer ()
 
void init (Simulator &simulator)
 Initialize the linearizer. More...
 
void eraseMatrix ()
 Causes the Jacobian matrix to be recreated from scratch before the next iteration. More...
 
void linearize ()
 Linearize the full system of non-linear equations. More...
 
void linearizeDomain ()
 Linearize the part of the non-linear system of equations that is associated with the spatial domain. More...
 
template<class SubDomainType >
void linearizeDomain (const SubDomainType &domain)
 Linearize the part of the non-linear system of equations that is associated with a part of the spatial domain. More...
 
void finalize ()
 
void linearizeAuxiliaryEquations ()
 Linearize the part of the non-linear system of equations that is associated with the spatial domain. More...
 
const SparseMatrixAdapter & jacobian () const
 Return constant reference to global Jacobian matrix backend. More...
 
SparseMatrixAdapter & jacobian ()
 
const GlobalEqVector & residual () const
 Return constant reference to global residual vector. More...
 
GlobalEqVector & residual ()
 
void setLinearizationType (LinearizationType linearizationType)
 
const LinearizationTypegetLinearizationType () const
 
const auto & getFlowsInfo () const
 Return constant reference to the flowsInfo. More...
 
const auto & getFloresInfo () const
 Return constant reference to the floresInfo. More...
 
const auto & getVelocityInfo () const
 Return constant reference to the velocityInfo. More...
 
void updateDiscretizationParameters ()
 
void updateBoundaryConditionData ()
 
const std::map< unsigned, Constraints > constraintsMap () const
 Returns the map of constraint degrees of freedom. More...
 
template<class SubDomainType >
void resetSystem_ (const SubDomainType &domain)
 
void setResAndJacobi (VectorBlock &res, MatrixBlock &bMat, const ADVectorBlock &resid) const
 
void updateFlowsInfo ()
 

Static Public Member Functions

static void registerParameters ()
 Register all run-time parameters for the Jacobian linearizer. More...
 

Detailed Description

template<class TypeTag>
class Opm::TpfaLinearizer< TypeTag >

The common code for the linearizers of non-linear systems of equations.

This class assumes that these system of equations to be linearized are stemming from models that use an finite volume scheme for spatial discretization and an Euler scheme for time discretization.

Constructor & Destructor Documentation

◆ TpfaLinearizer()

template<class TypeTag >
Opm::TpfaLinearizer< TypeTag >::TpfaLinearizer ( )
inline

◆ ~TpfaLinearizer()

template<class TypeTag >
Opm::TpfaLinearizer< TypeTag >::~TpfaLinearizer ( )
inline

Member Function Documentation

◆ constraintsMap()

template<class TypeTag >
const std::map< unsigned, Constraints > Opm::TpfaLinearizer< TypeTag >::constraintsMap ( ) const
inline

Returns the map of constraint degrees of freedom.

(This object is only non-empty if the EnableConstraints property is true.)

◆ eraseMatrix()

template<class TypeTag >
void Opm::TpfaLinearizer< TypeTag >::eraseMatrix ( )
inline

Causes the Jacobian matrix to be recreated from scratch before the next iteration.

This method is usally called if the sparsity pattern has changed for some reason. (e.g. by modifications of the grid or changes of the auxiliary equations.)

Referenced by Opm::TpfaLinearizer< TypeTag >::init().

◆ finalize()

template<class TypeTag >
void Opm::TpfaLinearizer< TypeTag >::finalize ( )
inline

◆ getFloresInfo()

template<class TypeTag >
const auto & Opm::TpfaLinearizer< TypeTag >::getFloresInfo ( ) const
inline

Return constant reference to the floresInfo.

(This object is only non-empty if the FLORES keyword is true.)

◆ getFlowsInfo()

template<class TypeTag >
const auto & Opm::TpfaLinearizer< TypeTag >::getFlowsInfo ( ) const
inline

Return constant reference to the flowsInfo.

(This object is only non-empty if the FLOWS keyword is true.)

◆ getLinearizationType()

template<class TypeTag >
const LinearizationType & Opm::TpfaLinearizer< TypeTag >::getLinearizationType ( ) const
inline

◆ getVelocityInfo()

template<class TypeTag >
const auto & Opm::TpfaLinearizer< TypeTag >::getVelocityInfo ( ) const
inline

Return constant reference to the velocityInfo.

(This object is only non-empty if the DISPERC keyword is true.)

◆ init()

template<class TypeTag >
void Opm::TpfaLinearizer< TypeTag >::init ( Simulator simulator)
inline

Initialize the linearizer.

At this point we can assume that all objects in the simulator have been allocated. We cannot assume that they are fully initialized, though.

References Opm::TpfaLinearizer< TypeTag >::eraseMatrix().

◆ jacobian() [1/2]

template<class TypeTag >
SparseMatrixAdapter & Opm::TpfaLinearizer< TypeTag >::jacobian ( )
inline

◆ jacobian() [2/2]

template<class TypeTag >
const SparseMatrixAdapter & Opm::TpfaLinearizer< TypeTag >::jacobian ( ) const
inline

Return constant reference to global Jacobian matrix backend.

◆ linearize()

template<class TypeTag >
void Opm::TpfaLinearizer< TypeTag >::linearize ( )
inline

Linearize the full system of non-linear equations.

The linearizationType() controls the scheme used and the focus time index. The default is fully implicit scheme, and focus index equal to 0, i.e. current time (end of step).

This linearizes the spatial domain and all auxiliary equations.

References Opm::TpfaLinearizer< TypeTag >::linearizeAuxiliaryEquations(), and Opm::TpfaLinearizer< TypeTag >::linearizeDomain().

◆ linearizeAuxiliaryEquations()

template<class TypeTag >
void Opm::TpfaLinearizer< TypeTag >::linearizeAuxiliaryEquations ( )
inline

Linearize the part of the non-linear system of equations that is associated with the spatial domain.

References Opm::Simulator< TypeTag >::gridView().

Referenced by Opm::TpfaLinearizer< TypeTag >::linearize().

◆ linearizeDomain() [1/2]

template<class TypeTag >
void Opm::TpfaLinearizer< TypeTag >::linearizeDomain ( )
inline

Linearize the part of the non-linear system of equations that is associated with the spatial domain.

That means that the global Jacobian of the residual is assembled and the residual is evaluated for the current solution.

The current state of affairs (esp. the previous and the current solutions) is represented by the model object.

References Opm::Simulator< TypeTag >::gridView(), and Opm::TpfaLinearizer< TypeTag >::linearizeDomain().

Referenced by Opm::TpfaLinearizer< TypeTag >::linearize(), and Opm::TpfaLinearizer< TypeTag >::linearizeDomain().

◆ linearizeDomain() [2/2]

template<class TypeTag >
template<class SubDomainType >
void Opm::TpfaLinearizer< TypeTag >::linearizeDomain ( const SubDomainType &  domain)
inline

Linearize the part of the non-linear system of equations that is associated with a part of the spatial domain.

That means that the Jacobian of the residual is assembled and the residual is evaluated for the current solution, on the domain passed in as argument.

The current state of affairs (esp. the previous and the current solutions) is represented by the model object.

References Opm::TpfaLinearizer< TypeTag >::linearizeDomain(), and Opm::TpfaLinearizer< TypeTag >::resetSystem_().

◆ registerParameters()

template<class TypeTag >
static void Opm::TpfaLinearizer< TypeTag >::registerParameters ( )
inlinestatic

Register all run-time parameters for the Jacobian linearizer.

◆ resetSystem_()

template<class TypeTag >
template<class SubDomainType >
void Opm::TpfaLinearizer< TypeTag >::resetSystem_ ( const SubDomainType &  domain)
inline

◆ residual() [1/2]

template<class TypeTag >
GlobalEqVector & Opm::TpfaLinearizer< TypeTag >::residual ( )
inline

◆ residual() [2/2]

template<class TypeTag >
const GlobalEqVector & Opm::TpfaLinearizer< TypeTag >::residual ( ) const
inline

Return constant reference to global residual vector.

◆ setLinearizationType()

template<class TypeTag >
void Opm::TpfaLinearizer< TypeTag >::setLinearizationType ( LinearizationType  linearizationType)
inline

◆ setResAndJacobi()

template<class TypeTag >
void Opm::TpfaLinearizer< TypeTag >::setResAndJacobi ( VectorBlock &  res,
MatrixBlock bMat,
const ADVectorBlock &  resid 
) const
inline

◆ updateBoundaryConditionData()

template<class TypeTag >
void Opm::TpfaLinearizer< TypeTag >::updateBoundaryConditionData ( )
inline

◆ updateDiscretizationParameters()

template<class TypeTag >
void Opm::TpfaLinearizer< TypeTag >::updateDiscretizationParameters ( )
inline

◆ updateFlowsInfo()

template<class TypeTag >
void Opm::TpfaLinearizer< TypeTag >::updateFlowsInfo ( )
inline

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