fvMatrix< Type > Class Template Reference

A special matrix type and solver, designed for finite volume solutions of scalar equations. Face addressing is used to make all matrix assembly and solution loops vectorise. More...

Inheritance diagram for fvMatrix< Type >:
Collaboration diagram for fvMatrix< Type >:

Classes

class  fvSolver
 Solver class returned by the solver function used for systems in which it is useful to cache the solver for reuse. More...
 

Public Types

typedef GeometricField< Type, fvPatchField, volMeshpsiFieldType
 The geometric field type for psi. More...
 
typedef GeometricField< Type, fvsPatchField, surfaceMeshfaceFluxFieldType
 Field type for face flux (for non-orthogonal correction) More...
 
typedef GeometricField< Type, fvsPatchField, surfaceMesh > * faceFluxFieldPtrType
 Declare return type of the faceFluxCorrectionPtr() function. More...
 
- Public Types inherited from lduMatrix
enum  normTypes : char { NO_NORM, DEFAULT_NORM, L1_SCALED_NORM }
 Enumerated matrix normalisation types. More...
 

Public Member Functions

 ClassName ("fvMatrix")
 
 fvMatrix (const GeometricField< Type, fvPatchField, volMesh > &psi, const dimensionSet &ds)
 Construct given a field to solve for. More...
 
 fvMatrix (const fvMatrix< Type > &)
 Copy construct. More...
 
 fvMatrix (const tmp< fvMatrix< Type >> &)
 Copy/move construct from tmp<fvMatrix<Type>> More...
 
 fvMatrix (const GeometricField< Type, fvPatchField, volMesh > &psi, Istream &is)
 Deprecated(2022-05) - construct with dimensionSet instead. More...
 
tmp< fvMatrix< Type > > clone () const
 Construct and return a clone. More...
 
virtual ~fvMatrix ()
 Destructor. More...
 
label nMatrices () const
 
const fvMatrix< Type > & matrix (const label i) const
 
fvMatrix< Type > & matrix (const label i)
 
label globalPatchID (const label fieldi, const label patchi) const
 
void transferFvMatrixCoeffs ()
 Transfer lower, upper, diag and source to this fvMatrix. More...
 
void createOrUpdateLduPrimitiveAssembly ()
 Create or update ldu assembly. More...
 
lduPrimitiveMeshAssemblylduMeshPtr ()
 Access to lduPrimitiveMeshAssembly. More...
 
const lduPrimitiveMeshAssemblylduMeshPtr () const
 Const Access to lduPrimitiveMeshAssembly. More...
 
void manipulateMatrix (direction cmp)
 Manipulate matrix. More...
 
void setBounAndInterCoeffs ()
 Manipulate boundary/internal coeffs for coupling. More...
 
void setInterfaces (lduInterfaceFieldPtrsList &, PtrDynList< lduInterfaceField > &newInterfaces)
 Set interfaces. More...
 
void mapContributions (label fieldi, const FieldField< Field, Type > &fluxContrib, FieldField< Field, Type > &contrib, bool internal) const
 Add internal and boundary contribution to local patches. More...
 
const lduPrimitiveMeshAssemblylduMeshAssembly ()
 Return optional lduAdressing. More...
 
const GeometricField< Type, fvPatchField, volMesh > & psi (const label i=0) const
 Return psi. More...
 
GeometricField< Type, fvPatchField, volMesh > & psi (const label i=0)
 
void clear ()
 Clear multiple fvMatrices. More...
 
const dimensionSetdimensions () const noexcept
 
Field< Type > & source () noexcept
 
const Field< Type > & source () const noexcept
 
const FieldField< Field, Type > & internalCoeffs () const noexcept
 fvBoundary scalar field containing pseudo-matrix coeffs for internal cells More...
 
FieldField< Field, Type > & internalCoeffs () noexcept
 fvBoundary scalar field containing pseudo-matrix coeffs for internal cells More...
 
const FieldField< Field, Type > & boundaryCoeffs () const noexcept
 fvBoundary scalar field containing pseudo-matrix coeffs for boundary cells More...
 
FieldField< Field, Type > & boundaryCoeffs () noexcept
 fvBoundary scalar field containing pseudo-matrix coeffs for boundary cells More...
 
faceFluxFieldPtrTypefaceFluxCorrectionPtr ()
 Return pointer to face-flux non-orthogonal correction field. More...
 
bool hasFaceFluxCorrection () const noexcept
 True if face-flux non-orthogonal correction field exists. More...
 
void setValues (const labelUList &cellLabels, const Type &value)
 Set solution in given cells to the specified value and eliminate the corresponding equations from the matrix. More...
 
void setValues (const labelUList &cellLabels, const UList< Type > &values)
 Set solution in given cells to the specified values and eliminate the corresponding equations from the matrix. More...
 
void setValues (const labelUList &cellLabels, const UIndirectList< Type > &values)
 Set solution in given cells to the specified values and eliminate the corresponding equations from the matrix. More...
 
void setReference (const label celli, const Type &value, const bool forceReference=false)
 Set reference level for solution. More...
 
void setReferences (const labelUList &cellLabels, const Type &value, const bool forceReference=false)
 Set reference level for solution. More...
 
void setReferences (const labelUList &cellLabels, const UList< Type > &values, const bool forceReference=false)
 Set reference level for solution. More...
 
void setComponentReference (const label patchi, const label facei, const direction cmpt, const scalar value)
 Set reference level for a component of the solution on a given patch face. More...
 
void addFvMatrix (fvMatrix< Type > &matrix)
 Add fvMatrix. More...
 
void relax (const scalar alpha)
 Relax matrix (for steady-state solution). More...
 
void relax ()
 Relax matrix (for steady-state solution). More...
 
void boundaryManipulate (typename GeometricField< Type, fvPatchField, volMesh >::Boundary &values)
 Manipulate based on a boundary field. More...
 
autoPtr< fvSolversolver (const dictionary &)
 Construct and return the solver. More...
 
autoPtr< fvSolversolver ()
 Construct and return the solver. More...
 
SolverPerformance< Type > solveSegregatedOrCoupled (const dictionary &)
 Solve segregated or coupled returning the solution statistics. More...
 
SolverPerformance< Type > solveSegregatedOrCoupled ()
 Solve segregated or coupled returning the solution statistics. More...
 
SolverPerformance< Type > solveSegregated (const dictionary &)
 Solve segregated returning the solution statistics. More...
 
SolverPerformance< Type > solveCoupled (const dictionary &)
 Solve coupled returning the solution statistics. More...
 
SolverPerformance< Type > solve (const dictionary &)
 Solve returning the solution statistics. More...
 
SolverPerformance< Type > solve (const word &name)
 Solve returning the solution statistics. More...
 
SolverPerformance< Type > solve ()
 Solve returning the solution statistics. More...
 
tmp< Field< Type > > residual () const
 Return the matrix residual. More...
 
tmp< scalarFieldD () const
 Return the matrix scalar diagonal. More...
 
tmp< Field< Type > > DD () const
 Return the matrix Type diagonal. More...
 
tmp< volScalarFieldA () const
 Return the central coefficient. More...
 
tmp< GeometricField< Type, fvPatchField, volMesh > > H () const
 Return the H operation source. More...
 
tmp< volScalarFieldH1 () const
 Return H(1) More...
 
tmp< GeometricField< Type, fvsPatchField, surfaceMesh > > flux () const
 Return the face-flux field from the matrix. More...
 
const dictionarysolverDict (const word &name) const
 Return the solver dictionary (from fvSolution) for name. More...
 
const dictionarysolverDict () const
 Return the solver dictionary for psi, taking into account finalIteration. More...
 
void operator= (const fvMatrix< Type > &)
 
void operator= (const tmp< fvMatrix< Type >> &)
 
void negate ()
 Inplace negate. More...
 
void operator+= (const fvMatrix< Type > &)
 
void operator+= (const tmp< fvMatrix< Type >> &)
 
void operator-= (const fvMatrix< Type > &)
 
void operator-= (const tmp< fvMatrix< Type >> &)
 
void operator+= (const DimensionedField< Type, volMesh > &)
 
void operator+= (const tmp< DimensionedField< Type, volMesh >> &)
 
void operator+= (const tmp< GeometricField< Type, fvPatchField, volMesh >> &)
 
void operator-= (const DimensionedField< Type, volMesh > &)
 
void operator-= (const tmp< DimensionedField< Type, volMesh >> &)
 
void operator-= (const tmp< GeometricField< Type, fvPatchField, volMesh >> &)
 
void operator+= (const dimensioned< Type > &)
 
void operator-= (const dimensioned< Type > &)
 
void operator+= (const Foam::zero)
 
void operator-= (const Foam::zero)
 
void operator*= (const volScalarField::Internal &)
 
void operator*= (const tmp< volScalarField::Internal > &)
 
void operator*= (const tmp< volScalarField > &)
 
void operator*= (const dimensioned< scalar > &)
 
template<>
void setComponentReference (const label patchi, const label facei, const direction, const scalar value)
 
template<>
Foam::autoPtr< Foam::fvMatrix< Foam::scalar >::fvSolversolver (const dictionary &solverControls)
 
template<>
Foam::solverPerformance solveSegregated (const dictionary &solverControls)
 
template<>
Foam::tmp< Foam::scalarFieldresidual () const
 
template<>
Foam::tmp< Foam::volScalarFieldH () const
 
template<>
Foam::tmp< Foam::volScalarFieldH1 () const
 
template<>
void setComponentReference (const label patchi, const label facei, const direction, const scalar value)
 
template<>
autoPtr< fvMatrix< scalar >::fvSolversolver (const dictionary &)
 
template<>
solverPerformance solveSegregated (const dictionary &)
 
template<>
tmp< scalarFieldresidual () const
 
template<>
tmp< volScalarFieldH () const
 
template<>
tmp< volScalarFieldH1 () const
 
- Public Member Functions inherited from refCount
constexpr refCount () noexcept
 Default construct, initializing count to 0. More...
 
int use_count () const noexcept
 Return the current reference count. More...
 
bool unique () const noexcept
 Return true if the reference count is zero. More...
 
void operator++ () noexcept
 Increment the reference count. More...
 
void operator++ (int) noexcept
 Increment the reference count. More...
 
void operator-- () noexcept
 Decrement the reference count. More...
 
void operator-- (int) noexcept
 Decrement the reference count. More...
 
- Public Member Functions inherited from lduMatrix
 ClassName ("lduMatrix")
 
 lduMatrix (const lduMesh &)
 Construct given an LDU addressed mesh. More...
 
 lduMatrix (const lduMatrix &)
 Construct as copy. More...
 
 lduMatrix (lduMatrix &, bool reuse)
 Construct as copy or re-use as specified. More...
 
 lduMatrix (const lduMesh &, Istream &)
 Construct given an LDU addressed mesh and an Istream from which the coefficients are read. More...
 
 ~lduMatrix ()
 Destructor. More...
 
const lduMeshmesh () const noexcept
 Return the LDU mesh from which the addressing is obtained. More...
 
void setLduMesh (const lduMesh &m)
 Set the LDU mesh containing the addressing is obtained. More...
 
const lduAddressinglduAddr () const
 Return the LDU addressing. More...
 
const lduSchedulepatchSchedule () const
 Return the patch evaluation schedule. More...
 
scalarFieldlower ()
 
scalarFielddiag ()
 
scalarFieldupper ()
 
scalarFieldlower (const label size)
 
scalarFielddiag (const label nCoeffs)
 
scalarFieldupper (const label nCoeffs)
 
const scalarFieldlower () const
 
const scalarFielddiag () const
 
const scalarFieldupper () const
 
bool hasDiag () const noexcept
 
bool hasUpper () const noexcept
 
bool hasLower () const noexcept
 
bool diagonal () const noexcept
 
bool symmetric () const noexcept
 
bool asymmetric () const noexcept
 
void sumDiag ()
 
void negSumDiag ()
 
void sumMagOffDiag (scalarField &sumOff) const
 
void Amul (solveScalarField &, const tmp< solveScalarField > &, const FieldField< Field, scalar > &, const lduInterfaceFieldPtrsList &, const direction cmpt) const
 Matrix multiplication with updated interfaces. More...
 
void Tmul (solveScalarField &, const tmp< solveScalarField > &, const FieldField< Field, scalar > &, const lduInterfaceFieldPtrsList &, const direction cmpt) const
 Matrix transpose multiplication with updated interfaces. More...
 
void sumA (solveScalarField &, const FieldField< Field, scalar > &, const lduInterfaceFieldPtrsList &) const
 Sum the coefficients on each row of the matrix. More...
 
void residual (solveScalarField &rA, const solveScalarField &psi, const scalarField &source, const FieldField< Field, scalar > &interfaceBouCoeffs, const lduInterfaceFieldPtrsList &interfaces, const direction cmpt) const
 
tmp< solveScalarFieldresidual (const solveScalarField &psi, const scalarField &source, const FieldField< Field, scalar > &interfaceBouCoeffs, const lduInterfaceFieldPtrsList &interfaces, const direction cmpt) const
 
void initMatrixInterfaces (const bool add, const FieldField< Field, scalar > &interfaceCoeffs, const lduInterfaceFieldPtrsList &interfaces, const solveScalarField &psiif, solveScalarField &result, const direction cmpt) const
 Initialise the update of interfaced interfaces for matrix operations. More...
 
void updateMatrixInterfaces (const bool add, const FieldField< Field, scalar > &interfaceCoeffs, const lduInterfaceFieldPtrsList &interfaces, const solveScalarField &psiif, solveScalarField &result, const direction cmpt, const label startRequest) const
 Update interfaced interfaces for matrix operations. More...
 
void setResidualField (const scalarField &residual, const word &fieldName, const bool initial) const
 Set the residual field using an IOField on the object registry if it exists. More...
 
template<class Type >
tmp< Field< Type > > H (const Field< Type > &) const
 
template<class Type >
tmp< Field< Type > > H (const tmp< Field< Type >> &) const
 
tmp< scalarFieldH1 () const
 
template<class Type >
tmp< Field< Type > > faceH (const Field< Type > &) const
 
template<class Type >
tmp< Field< Type > > faceH (const tmp< Field< Type >> &) const
 
InfoProxy< lduMatrixinfo () const noexcept
 Return info proxy, used to print matrix information to a stream. More...
 
void operator= (const lduMatrix &)
 
void negate ()
 
void operator+= (const lduMatrix &)
 
void operator-= (const lduMatrix &)
 
void operator*= (const scalarField &)
 
void operator*= (scalar)
 
template<class Type >
Foam::tmp< Foam::Field< Type > > H (const Field< Type > &psi) const
 
template<class Type >
Foam::tmp< Foam::Field< Type > > H (const tmp< Field< Type >> &tpsi) const
 
template<class Type >
Foam::tmp< Foam::Field< Type > > faceH (const Field< Type > &psi) const
 
template<class Type >
Foam::tmp< Foam::Field< Type > > faceH (const tmp< Field< Type >> &tpsi) const
 

Protected Member Functions

template<class Type2 >
void addToInternalField (const labelUList &addr, const Field< Type2 > &pf, Field< Type2 > &intf) const
 Add patch contribution to internal field. More...
 
template<class Type2 >
void addToInternalField (const labelUList &addr, const tmp< Field< Type2 >> &tpf, Field< Type2 > &intf) const
 
template<class Type2 >
void subtractFromInternalField (const labelUList &addr, const Field< Type2 > &pf, Field< Type2 > &intf) const
 Subtract patch contribution from internal field. More...
 
template<class Type2 >
void subtractFromInternalField (const labelUList &addr, const tmp< Field< Type2 >> &tpf, Field< Type2 > &intf) const
 
bool checkImplicit (const label fieldi=0)
 Name the implicit assembly addressing. More...
 
void addBoundaryDiag (scalarField &diag, const direction cmpt) const
 
void addCmptAvBoundaryDiag (scalarField &diag) const
 
void addBoundarySource (Field< Type > &source, const bool couples=true) const
 
template<template< class > class ListType>
void setValuesFromList (const labelUList &cellLabels, const ListType< Type > &values)
 Set solution in given cells to the specified values. More...
 

Friends

class fvSolver
 Declare friendship with the fvSolver class. More...
 
tmp< GeometricField< Type, fvPatchField, volMesh > > operator & (const fvMatrix< Type > &, const DimensionedField< Type, volMesh > &)
 
tmp< GeometricField< Type, fvPatchField, volMesh > > operator & (const fvMatrix< Type > &, const tmp< GeometricField< Type, fvPatchField, volMesh >> &)
 
tmp< GeometricField< Type, fvPatchField, volMesh > > operator & (const tmp< fvMatrix< Type >> &, const DimensionedField< Type, volMesh > &)
 
tmp< GeometricField< Type, fvPatchField, volMesh > > operator & (const tmp< fvMatrix< Type >> &, const tmp< GeometricField< Type, fvPatchField, volMesh >> &)
 
Ostreamoperator (Ostream &, const fvMatrix< Type > &)
 

Additional Inherited Members

- Static Public Attributes inherited from lduMatrix
static const Enum< normTypesnormTypesNames_
 Names for the normTypes. More...
 
static constexpr const label defaultMaxIter = 1000
 Default maximum number of iterations for solvers (1000) More...
 
static const scalar defaultTolerance = 1e-6
 Default (absolute) tolerance (1e-6) More...
 

Detailed Description

template<class Type>
class Foam::fvMatrix< Type >

A special matrix type and solver, designed for finite volume solutions of scalar equations. Face addressing is used to make all matrix assembly and solution loops vectorise.

Source files

Definition at line 64 of file fvPatchField.H.

Member Typedef Documentation

◆ psiFieldType

The geometric field type for psi.

Definition at line 126 of file fvMatrix.H.

◆ faceFluxFieldType

Field type for face flux (for non-orthogonal correction)

Definition at line 133 of file fvMatrix.H.

◆ faceFluxFieldPtrType

Declare return type of the faceFluxCorrectionPtr() function.

Definition at line 583 of file fvMatrix.H.

Constructor & Destructor Documentation

◆ fvMatrix() [1/4]

fvMatrix ( const GeometricField< Type, fvPatchField, volMesh > &  psi,
const dimensionSet ds 
)

Construct given a field to solve for.

Definition at line 352 of file fvMatrix.C.

References DebugInFunction, Foam::endl(), forAll, psi, and Foam::Zero.

Here is the call graph for this function:

◆ fvMatrix() [2/4]

fvMatrix ( const fvMatrix< Type > &  fvm)

Copy construct.

Definition at line 396 of file fvMatrix.C.

References DebugInFunction, Foam::endl(), and IOobject::name().

Here is the call graph for this function:

◆ fvMatrix() [3/4]

fvMatrix ( const tmp< fvMatrix< Type >> &  tmat)

Copy/move construct from tmp<fvMatrix<Type>>

Definition at line 424 of file fvMatrix.C.

References DebugInFunction, Foam::endl(), and IOobject::name().

Here is the call graph for this function:

◆ fvMatrix() [4/4]

fvMatrix ( const GeometricField< Type, fvPatchField, volMesh > &  psi,
Istream is 
)
inline

Deprecated(2022-05) - construct with dimensionSet instead.

Deprecated:
(2022-05) - construct with dimensionSet instead

Definition at line 363 of file fvMatrix.H.

◆ ~fvMatrix()

~fvMatrix ( )
virtual

Destructor.

Definition at line 464 of file fvMatrix.C.

References DebugInFunction, Foam::deleteDemandDrivenData(), and Foam::endl().

Here is the call graph for this function:

Member Function Documentation

◆ addToInternalField() [1/2]

void addToInternalField ( const labelUList addr,
const Field< Type2 > &  pf,
Field< Type2 > &  intf 
) const
protected

Add patch contribution to internal field.

Definition at line 42 of file fvMatrix.C.

References Foam::abort(), Foam::endl(), Foam::FatalError, FatalErrorInFunction, forAll, UList< T >::size(), and UList< Type >::size().

Here is the call graph for this function:

◆ addToInternalField() [2/2]

void addToInternalField ( const labelUList addr,
const tmp< Field< Type2 >> &  tpf,
Field< Type2 > &  intf 
) const
protected

Definition at line 66 of file fvMatrix.C.

◆ subtractFromInternalField() [1/2]

void subtractFromInternalField ( const labelUList addr,
const Field< Type2 > &  pf,
Field< Type2 > &  intf 
) const
protected

Subtract patch contribution from internal field.

Definition at line 80 of file fvMatrix.C.

References Foam::abort(), Foam::endl(), Foam::FatalError, FatalErrorInFunction, forAll, UList< T >::size(), and UList< Type >::size().

Here is the call graph for this function:

◆ subtractFromInternalField() [2/2]

void subtractFromInternalField ( const labelUList addr,
const tmp< Field< Type2 >> &  tpf,
Field< Type2 > &  intf 
) const
protected

Definition at line 104 of file fvMatrix.C.

◆ checkImplicit()

bool checkImplicit ( const label  fieldi = 0)
protected

Name the implicit assembly addressing.

Returns
true if assembly is used

Definition at line 315 of file fvMatrix.C.

References Foam::ensightOutput::debug, Foam::endl(), forAll, Foam::name(), Foam::Pout, and psi.

Here is the call graph for this function:

◆ addBoundaryDiag()

void addBoundaryDiag ( scalarField diag,
const direction  cmpt 
) const
protected

Definition at line 117 of file fvMatrix.C.

References Foam::component(), Foam::diag(), forAll, and psi.

Referenced by fvMatrix< Type >::residual().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ addCmptAvBoundaryDiag()

void addCmptAvBoundaryDiag ( scalarField diag) const
protected

Definition at line 145 of file fvMatrix.C.

References Foam::cmptAv(), Foam::diag(), forAll, and psi.

Here is the call graph for this function:

◆ addBoundarySource()

void addBoundarySource ( Field< Type > &  source,
const bool  couples = true 
) const
protected

Definition at line 170 of file fvMatrix.C.

References Foam::cmptMultiply(), fvPatchFieldBase::coupled(), forAll, fvPatchField< Type >::patchNeighbourField(), and psi.

Referenced by fvMatrix< Type >::residual().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ setValuesFromList()

void setValuesFromList ( const labelUList cellLabels,
const ListType< Type > &  values 
)
protected

Set solution in given cells to the specified values.

Definition at line 220 of file fvMatrix.C.

References cells, Foam::diag(), forAll, Foam::stringOps::lower(), mesh, primitiveFieldRef(), psi, Foam::stringOps::upper(), Foam::HashTableOps::values(), and Foam::Zero.

Here is the call graph for this function:

◆ ClassName()

ClassName ( "fvMatrix< Type >"  )

◆ clone()

tmp<fvMatrix<Type> > clone ( ) const
inline

Construct and return a clone.

Definition at line 374 of file fvMatrix.H.

References Foam::New().

Referenced by fvMatrix< Type >::addFvMatrix().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ nMatrices()

label nMatrices ( ) const
inline

Definition at line 392 of file fvMatrix.H.

◆ matrix() [1/2]

const fvMatrix<Type>& matrix ( const label  i) const
inline

Definition at line 397 of file fvMatrix.H.

Referenced by fvMatrix< Type >::psi().

Here is the caller graph for this function:

◆ matrix() [2/2]

fvMatrix<Type>& matrix ( const label  i)
inline

Definition at line 402 of file fvMatrix.H.

◆ globalPatchID()

label globalPatchID ( const label  fieldi,
const label  patchi 
) const
inline

Definition at line 408 of file fvMatrix.H.

References fvMatrix< Type >::lduMeshPtr(), and lduPrimitiveMeshAssembly::patchMap().

Here is the call graph for this function:

◆ transferFvMatrixCoeffs()

void transferFvMatrixCoeffs ( )

Transfer lower, upper, diag and source to this fvMatrix.

Definition at line 821 of file fvMatrix.C.

References Foam::diag(), Foam::faceMap(), forAll, Foam::stringOps::lower(), List< Type >::setSize(), UList< T >::size(), Foam::stringOps::upper(), and Foam::Zero.

Here is the call graph for this function:

◆ createOrUpdateLduPrimitiveAssembly()

void createOrUpdateLduPrimitiveAssembly ( )

◆ lduMeshPtr() [1/2]

Access to lduPrimitiveMeshAssembly.

Definition at line 888 of file fvMatrix.C.

References lduPrimitiveMesh::mesh(), and lduMesh::thisDb().

Referenced by fvMatrix< Type >::globalPatchID(), and fvMatrix< Type >::lduMeshAssembly().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ lduMeshPtr() [2/2]

const Foam::lduPrimitiveMeshAssembly * lduMeshPtr ( ) const

Const Access to lduPrimitiveMeshAssembly.

Definition at line 901 of file fvMatrix.C.

◆ manipulateMatrix()

void manipulateMatrix ( direction  cmp)

Manipulate matrix.

Definition at line 798 of file fvMatrix.C.

References GeometricField< Type, PatchField, GeoMesh >::boundaryFieldRef(), forAll, and psi.

Here is the call graph for this function:

◆ setBounAndInterCoeffs()

void setBounAndInterCoeffs ( )

Manipulate boundary/internal coeffs for coupling.

Definition at line 668 of file fvMatrix.C.

References boundary, GeometricField< Type, PatchField, GeoMesh >::boundaryField(), forAll, DimensionedField< Type, GeoMesh >::mesh(), psi, List< T >::set(), and Foam::Zero.

Here is the call graph for this function:

◆ setInterfaces()

void setInterfaces ( lduInterfaceFieldPtrsList interfaces,
PtrDynList< lduInterfaceField > &  newInterfaces 
)

◆ mapContributions()

void mapContributions ( label  fieldi,
const FieldField< Field, Type > &  fluxContrib,
FieldField< Field, Type > &  contrib,
bool  internal 
) const

◆ lduMeshAssembly()

const lduPrimitiveMeshAssembly& lduMeshAssembly ( )
inline

Return optional lduAdressing.

Definition at line 476 of file fvMatrix.H.

References fvMatrix< Type >::lduMeshPtr().

Referenced by mixedEnergyFvPatchScalarField::manipulateMatrix(), cyclicFvPatchField< vector >::manipulateMatrix(), cyclicACMIFvPatchField< Type >::manipulateMatrix(), and cyclicAMIFvPatchField< scalar >::manipulateMatrix().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ psi() [1/2]

const GeometricField<Type, fvPatchField, volMesh>& psi ( const label  i = 0) const
inline

Return psi.

Definition at line 485 of file fvMatrix.H.

References fvMatrix< Type >::matrix().

Referenced by MRFZone::addCoriolis(), velocityDampingConstraint::addDamping(), oversetFvMeshBase::addInterpolation(), isotropic::addRegularisationTerm(), objectivePowerDissipation::addSource(), objectiveNutSqr::addSource(), objectiveIncompressible::addSource(), radiation::addSup(), topOSource::addSup(), multiphaseMangrovesSource::addSup(), patchCellsSource::addSup(), multiphaseMangrovesTurbulenceModel::addSup(), atmPlantCanopyUSource::addSup(), PhaseLimitStabilization< Type >::addSup(), radialActuationDiskSource::addSup(), explicitPorositySource::addSup(), interRegionExplicitPorositySource::addSup(), atmCoriolisUSource::addSup(), acousticDampingSource::addSup(), multiphaseStabilizedTurbulence::addSup(), interRegionHeatTransferModel::addSup(), jouleHeatingSource::addSup(), SemiImplicitSource< Type >::addSup(), solidificationMeltingSource::addSup(), rotorDiskSource::addSup(), Foam::checkMethod(), optionList::constrain(), oversetFvPatchField< Type >::fringeFlux(), cyclicFvPatchField< vector >::manipulateMatrix(), cyclicACMIFvPatchField< Type >::manipulateMatrix(), cyclicAMIFvPatchField< scalar >::manipulateMatrix(), MomentumTransferPhaseSystem< BasePhaseSystem >::momentumTransfer(), oversetFvMeshBase::normalisation(), Foam::oversetPatchPhiErr(), fvMatrix< Type >::psi(), regularisationPDE::setValues(), oversetFvMeshBase::solveOverset(), and oversetFvMeshBase::write().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ psi() [2/2]

GeometricField<Type, fvPatchField, volMesh>& psi ( const label  i = 0)
inline

Definition at line 497 of file fvMatrix.H.

References fvMatrix< Type >::matrix(), and fvMatrix< Type >::psi().

Here is the call graph for this function:

◆ clear()

void clear ( )
inline

Clear multiple fvMatrices.

Definition at line 521 of file fvMatrix.H.

References PtrList< T >::clear().

Here is the call graph for this function:

◆ dimensions()

◆ source() [1/2]

◆ source() [2/2]

const Field<Type>& source ( ) const
inlinenoexcept

Definition at line 538 of file fvMatrix.H.

◆ internalCoeffs() [1/2]

◆ internalCoeffs() [2/2]

FieldField<Field, Type>& internalCoeffs ( )
inlinenoexcept

fvBoundary scalar field containing pseudo-matrix coeffs for internal cells

Definition at line 556 of file fvMatrix.H.

◆ boundaryCoeffs() [1/2]

◆ boundaryCoeffs() [2/2]

FieldField<Field, Type>& boundaryCoeffs ( )
inlinenoexcept

fvBoundary scalar field containing pseudo-matrix coeffs for boundary cells

Definition at line 574 of file fvMatrix.H.

◆ faceFluxCorrectionPtr()

faceFluxFieldPtrType& faceFluxCorrectionPtr ( )
inline

Return pointer to face-flux non-orthogonal correction field.

Definition at line 588 of file fvMatrix.H.

Referenced by gaussLaplacianScheme< Type, GType >::fvmLaplacian(), and relaxedNonOrthoGaussLaplacianScheme< Type, GType >::fvmLaplacian().

Here is the caller graph for this function:

◆ hasFaceFluxCorrection()

bool hasFaceFluxCorrection ( ) const
inlinenoexcept

True if face-flux non-orthogonal correction field exists.

Definition at line 596 of file fvMatrix.H.

◆ setValues() [1/3]

void setValues ( const labelUList cellLabels,
const Type &  value 
)

◆ setValues() [2/3]

void setValues ( const labelUList cellLabels,
const UList< Type > &  values 
)

Set solution in given cells to the specified values and eliminate the corresponding equations from the matrix.

Definition at line 989 of file fvMatrix.C.

References Foam::HashTableOps::values().

Here is the call graph for this function:

◆ setValues() [3/3]

void setValues ( const labelUList cellLabels,
const UIndirectList< Type > &  values 
)

Set solution in given cells to the specified values and eliminate the corresponding equations from the matrix.

Definition at line 1000 of file fvMatrix.C.

References Foam::HashTableOps::values().

Here is the call graph for this function:

◆ setReference()

void setReference ( const label  celli,
const Type &  value,
const bool  forceReference = false 
)

Set reference level for solution.

Definition at line 1011 of file fvMatrix.C.

References Foam::diag().

Referenced by Foam::CorrectPhi(), adjointSimple::mainIter(), and simple::mainIter().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ setReferences() [1/2]

void setReferences ( const labelUList cellLabels,
const Type &  value,
const bool  forceReference = false 
)

Set reference level for solution.

Definition at line 1027 of file fvMatrix.C.

References cellId, Foam::diag(), and forAll.

Here is the call graph for this function:

◆ setReferences() [2/2]

void setReferences ( const labelUList cellLabels,
const UList< Type > &  values,
const bool  forceReference = false 
)

Set reference level for solution.

Definition at line 1050 of file fvMatrix.C.

References cellId, Foam::diag(), forAll, and Foam::HashTableOps::values().

Here is the call graph for this function:

◆ setComponentReference() [1/3]

void setComponentReference ( const label  patchi,
const label  facei,
const direction  cmpt,
const scalar  value 
)

Set reference level for a component of the solution on a given patch face.

Definition at line 31 of file fvMatrixSolve.C.

References Foam::diag().

Here is the call graph for this function:

◆ addFvMatrix()

void addFvMatrix ( fvMatrix< Type > &  matrix)

Add fvMatrix.

Definition at line 1072 of file fvMatrix.C.

References Foam::abort(), fvMatrix< Type >::clone(), fvMatrix< Type >::dimensions(), Foam::endl(), Foam::FatalError, and FatalErrorInFunction.

Here is the call graph for this function:

◆ relax() [1/2]

◆ relax() [2/2]

void relax ( )

Relax matrix (for steady-state solution).

alpha is read from controlDict

Definition at line 1249 of file fvMatrix.C.

References Foam::name(), and relax().

Here is the call graph for this function:

◆ boundaryManipulate()

void boundaryManipulate ( typename GeometricField< Type, fvPatchField, volMesh >::Boundary &  values)

Manipulate based on a boundary field.

Definition at line 1267 of file fvMatrix.C.

References forAll.

Referenced by kOmegaSSTBase< eddyViscosity< RASModel< BasicTurbulenceModel > > >::correct(), adjointSimple::mainIter(), adjointMeshMovementSolver::solve(), multiphaseSystem::solveAlphas(), and oversetFvMeshBase::solveOverset().

Here is the caller graph for this function:

◆ solver() [1/4]

Construct and return the solver.

Use the given solver controls

◆ solver() [2/4]

Foam::autoPtr< typename Foam::fvMatrix< Type >::fvSolver > solver ( )

Construct and return the solver.

Solver controls read from fvSolution

Definition at line 329 of file fvMatrixSolve.C.

◆ solveSegregatedOrCoupled() [1/2]

◆ solveSegregatedOrCoupled() [2/2]

Foam::SolverPerformance< Type > solveSegregatedOrCoupled ( )

Solve segregated or coupled returning the solution statistics.

Solver controls read from fvSolution

Definition at line 311 of file fvMatrixSolve.C.

◆ solveSegregated() [1/3]

◆ solveCoupled()

Foam::SolverPerformance< Type > solveCoupled ( const dictionary solverControls)

Solve coupled returning the solution statistics.

Use the given solver controls

Definition at line 245 of file fvMatrixSolve.C.

References Foam::ensightOutput::debug, Foam::diag(), LduMatrix< Type, DType, LUType >::diag(), Foam::endl(), dictionary::getOrDefault(), Foam::Info, Foam::stringOps::lower(), messageStream::masterStream(), mesh, SolverPerformance< Type >::print(), psi, and Foam::stringOps::upper().

Here is the call graph for this function:

◆ solve() [1/3]

◆ solve() [2/3]

Foam::SolverPerformance< Type > solve ( const word name)

Solve returning the solution statistics.

Uses name solver controls from fvSolution

Definition at line 343 of file fvMatrixSolve.C.

References fvMatrix< Type >::solve(), and fvMatrix< Type >::solverDict().

Here is the call graph for this function:

◆ solve() [3/3]

Foam::SolverPerformance< Type > solve ( )

Solve returning the solution statistics.

Solver controls read from fvSolution

Definition at line 350 of file fvMatrixSolve.C.

References fvMatrix< Type >::solverDict().

Referenced by fvMatrix< Type >::solve().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ residual() [1/3]

◆ D()

Foam::tmp< Foam::scalarField > D ( ) const

Return the matrix scalar diagonal.

Definition at line 1280 of file fvMatrix.C.

References Foam::diag(), and tmp< T >::New().

Here is the call graph for this function:

◆ DD()

Foam::tmp< Foam::Field< Type > > DD ( ) const

Return the matrix Type diagonal.

Definition at line 1289 of file fvMatrix.C.

References fvPatchFieldBase::coupled(), Foam::diag(), forAll, tmp< T >::ref(), and UList< Type >::size().

Here is the call graph for this function:

◆ A()

Return the central coefficient.

Definition at line 1313 of file fvMatrix.C.

References D, Foam::dimVol, fvPatchFieldBase::extrapolatedCalculatedType(), and GeometricField< scalar, fvPatchField, volMesh >::New().

Referenced by meanVelocityForce::constrain(), directionalPressureGradientExplicitSource::constrain(), incompressiblePrimalSolver::correctBoundaryConditions(), and simple::mainIter().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ H() [1/3]

◆ H1() [1/3]

Return H(1)

Definition at line 1384 of file fvMatrix.C.

References Foam::component(), fvPatchFieldBase::coupled(), Foam::dimVol, fvPatchFieldBase::extrapolatedCalculatedType(), forAll, lduMatrix::H1(), GeometricField< scalar, fvPatchField, volMesh >::New(), and UList< Type >::size().

Referenced by simple::mainIter().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ flux()

◆ solverDict() [1/2]

const Foam::dictionary & solverDict ( const word name) const

Return the solver dictionary (from fvSolution) for name.

Definition at line 1526 of file fvMatrix.C.

References Foam::name().

Here is the call graph for this function:

◆ solverDict() [2/2]

const Foam::dictionary & solverDict ( ) const

Return the solver dictionary for psi, taking into account finalIteration.

Definition at line 1535 of file fvMatrix.C.

Referenced by fvMatrix< Type >::solve().

Here is the caller graph for this function:

◆ operator=() [1/2]

void operator= ( const fvMatrix< Type > &  fvmv)

Definition at line 1547 of file fvMatrix.C.

References Foam::abort(), Foam::FatalError, FatalErrorInFunction, and lduMatrix::operator=().

Here is the call graph for this function:

◆ operator=() [2/2]

void operator= ( const tmp< fvMatrix< Type >> &  tfvmv)

Definition at line 1584 of file fvMatrix.C.

◆ negate()

void negate ( )

Inplace negate.

Definition at line 1592 of file fvMatrix.C.

References lduMatrix::negate().

Here is the call graph for this function:

◆ operator+=() [1/7]

void operator+= ( const fvMatrix< Type > &  fvmv)

Definition at line 1607 of file fvMatrix.C.

References Foam::checkMethod(), and lduMatrix::operator+=().

Here is the call graph for this function:

◆ operator+=() [2/7]

void operator+= ( const tmp< fvMatrix< Type >> &  tfvmv)

Definition at line 1637 of file fvMatrix.C.

◆ operator-=() [1/7]

void operator-= ( const fvMatrix< Type > &  fvmv)

Definition at line 1645 of file fvMatrix.C.

References Foam::checkMethod(), and lduMatrix::operator-=().

Here is the call graph for this function:

◆ operator-=() [2/7]

void operator-= ( const tmp< fvMatrix< Type >> &  tfvmv)

Definition at line 1673 of file fvMatrix.C.

◆ operator+=() [3/7]

void operator+= ( const DimensionedField< Type, volMesh > &  su)

Definition at line 1682 of file fvMatrix.C.

References Foam::checkMethod().

Here is the call graph for this function:

◆ operator+=() [4/7]

void operator+= ( const tmp< DimensionedField< Type, volMesh >> &  tsu)

Definition at line 1693 of file fvMatrix.C.

◆ operator+=() [5/7]

void operator+= ( const tmp< GeometricField< Type, fvPatchField, volMesh >> &  tsu)

Definition at line 1704 of file fvMatrix.C.

◆ operator-=() [3/7]

void operator-= ( const DimensionedField< Type, volMesh > &  su)

Definition at line 1715 of file fvMatrix.C.

References Foam::checkMethod().

Here is the call graph for this function:

◆ operator-=() [4/7]

void operator-= ( const tmp< DimensionedField< Type, volMesh >> &  tsu)

Definition at line 1726 of file fvMatrix.C.

◆ operator-=() [5/7]

void operator-= ( const tmp< GeometricField< Type, fvPatchField, volMesh >> &  tsu)

Definition at line 1737 of file fvMatrix.C.

◆ operator+=() [6/7]

void operator+= ( const dimensioned< Type > &  su)

Definition at line 1748 of file fvMatrix.C.

References DimensionedField< Type, GeoMesh >::mesh(), and psi.

Here is the call graph for this function:

◆ operator-=() [6/7]

void operator-= ( const dimensioned< Type > &  su)

Definition at line 1758 of file fvMatrix.C.

References DimensionedField< Type, GeoMesh >::mesh(), and psi.

Here is the call graph for this function:

◆ operator+=() [7/7]

void operator+= ( const Foam::zero  )
inline

Definition at line 851 of file fvMatrix.H.

◆ operator-=() [7/7]

void operator-= ( const Foam::zero  )
inline

Definition at line 852 of file fvMatrix.H.

◆ operator*=() [1/4]

void operator*= ( const volScalarField::Internal dsf)

Definition at line 1768 of file fvMatrix.C.

References Foam::abort(), fvMatrix< Type >::dimensions(), Foam::FatalError, FatalErrorInFunction, forAll, and lduMatrix::operator*=().

Here is the call graph for this function:

◆ operator*=() [2/4]

void operator*= ( const tmp< volScalarField::Internal > &  tfld)

Definition at line 1798 of file fvMatrix.C.

◆ operator*=() [3/4]

void operator*= ( const tmp< volScalarField > &  tfld)

Definition at line 1809 of file fvMatrix.C.

◆ operator*=() [4/4]

void operator*= ( const dimensioned< scalar > &  ds)

Definition at line 1820 of file fvMatrix.C.

References lduMatrix::operator*=().

Here is the call graph for this function:

◆ setComponentReference() [2/3]

void setComponentReference ( const label  patchi,
const label  facei,
const direction  ,
const scalar  value 
)

Definition at line 34 of file fvScalarMatrix.C.

References Foam::diag().

Here is the call graph for this function:

◆ solver() [3/4]

Foam::autoPtr< Foam::fvMatrix< Foam::scalar >::fvSolver > solver ( const dictionary solverControls)

Definition at line 59 of file fvScalarMatrix.C.

References addProfiling, Foam::ensightOutput::debug, Foam::diag(), Foam::endl(), Foam::Info, messageStream::masterStream(), mesh, Foam::New(), regionName, and solve().

Here is the call graph for this function:

◆ solveSegregated() [2/3]

◆ residual() [2/3]

Foam::tmp< Foam::scalarField > residual ( ) const

Definition at line 325 of file fvScalarMatrix.C.

References psi, refPtr< Container< Type > >::ref(), and Foam::Zero.

Here is the call graph for this function:

◆ H() [2/3]

Definition at line 354 of file fvScalarMatrix.C.

References Foam::dimVol, Foam::fieldTypes::extrapolatedCalculatedType, H(), and Foam::New().

Here is the call graph for this function:

◆ H1() [2/3]

Definition at line 376 of file fvScalarMatrix.C.

References Foam::dimVol, Foam::fieldTypes::extrapolatedCalculatedType, and Foam::New().

Here is the call graph for this function:

◆ setComponentReference() [3/3]

void setComponentReference ( const label  patchi,
const label  facei,
const direction  ,
const scalar  value 
)

◆ solver() [4/4]

autoPtr< fvMatrix< scalar >::fvSolver > solver ( const dictionary )

◆ solveSegregated() [3/3]

solverPerformance solveSegregated ( const dictionary )

◆ residual() [3/3]

tmp< scalarField > residual ( ) const

◆ H() [3/3]

tmp< volScalarField > H ( ) const

◆ H1() [3/3]

tmp< volScalarField > H1 ( ) const

Friends And Related Function Documentation

◆ fvSolver

friend class fvSolver
friend

Declare friendship with the fvSolver class.

Definition at line 203 of file fvMatrix.H.

◆ operator & [1/4]

tmp<GeometricField<Type, fvPatchField, volMesh> > operator& ( const fvMatrix< Type > &  ,
const DimensionedField< Type, volMesh > &   
)
friend

◆ operator & [2/4]

tmp<GeometricField<Type, fvPatchField, volMesh> > operator& ( const fvMatrix< Type > &  ,
const tmp< GeometricField< Type, fvPatchField, volMesh >> &   
)
friend

◆ operator & [3/4]

tmp<GeometricField<Type, fvPatchField, volMesh> > operator& ( const tmp< fvMatrix< Type >> &  ,
const DimensionedField< Type, volMesh > &   
)
friend

◆ operator & [4/4]

tmp<GeometricField<Type, fvPatchField, volMesh> > operator& ( const tmp< fvMatrix< Type >> &  ,
const tmp< GeometricField< Type, fvPatchField, volMesh >> &   
)
friend

◆ operator

Ostream& operator ( Ostream ,
const fvMatrix< Type > &   
)
friend

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