LduMatrix< Type, DType, LUType > Class Template Reference

LduMatrix is a general matrix class in which the coefficients are stored as three arrays, one for the upper triangle, one for the lower triangle and a third for the diagonal. More...

Classes

class  preconditioner
 Abstract base-class for LduMatrix preconditioners. More...
 
class  smoother
 Abstract base-class for LduMatrix smoothers. More...
 
class  solver
 Abstract base-class for LduMatrix solvers. More...
 

Public Member Functions

 ClassName ("LduMatrix")
 
 LduMatrix (const lduMesh &mesh)
 Construct given an LDU addressed mesh. More...
 
 LduMatrix (const LduMatrix< Type, DType, LUType > &)
 Copy construct. More...
 
 LduMatrix (LduMatrix< Type, DType, LUType > &&)
 Move construct. More...
 
 LduMatrix (LduMatrix< Type, DType, LUType > &, bool reuse)
 Construct as copy or re-use as specified. More...
 
 LduMatrix (const lduMesh &mesh, Istream &is)
 Construct given an LDU addressed mesh and an Istream from which the coefficients are read. More...
 
 ~LduMatrix ()=default
 Destructor. More...
 
const lduMeshmesh () const noexcept
 Return the LDU mesh from which the addressing is obtained. More...
 
const lduAddressinglduAddr () const
 Return the LDU addressing. More...
 
const lduSchedulepatchSchedule () const
 Return the patch evaluation schedule. More...
 
const LduInterfaceFieldPtrsList< Type > & interfaces () const noexcept
 Const access to the interfaces. More...
 
LduInterfaceFieldPtrsList< Type > & interfaces () noexcept
 Non-const access to the interfaces. More...
 
const Field< DType > & diag () const
 
const Field< LUType > & upper () const
 
const Field< LUType > & lower () const
 
const Field< Type > & source () const
 
Field< DType > & diag ()
 
Field< LUType > & upper ()
 
Field< LUType > & lower ()
 
Field< Type > & source ()
 
const FieldField< Field, LUType > & interfacesUpper () const noexcept
 
const FieldField< Field, LUType > & interfacesLower () const noexcept
 
FieldField< Field, LUType > & interfacesUpper () noexcept
 
FieldField< Field, LUType > & interfacesLower () noexcept
 
word matrixTypeName () const
 The matrix type (empty, diagonal, symmetric, ...) More...
 
bool hasDiag () const noexcept
 
bool hasUpper () const noexcept
 
bool hasLower () const noexcept
 
bool hasSource () const noexcept
 
bool diagonal () const noexcept
 Matrix has diagonal only. More...
 
bool symmetric () const noexcept
 Matrix is symmetric. More...
 
bool asymmetric () const noexcept
 Matrix is asymmetric (ie, full) More...
 
void sumDiag ()
 
void negSumDiag ()
 
void sumMagOffDiag (Field< LUType > &sumOff) const
 
void Amul (Field< Type > &, const tmp< Field< Type >> &) const
 Matrix multiplication. More...
 
void Tmul (Field< Type > &, const tmp< Field< Type >> &) const
 Matrix transpose multiplication. More...
 
void sumA (Field< Type > &) const
 Sum the coefficients on each row of the matrix. More...
 
void residual (Field< Type > &rA, const Field< Type > &psi) const
 
tmp< Field< Type > > residual (const Field< Type > &psi) const
 
void initMatrixInterfaces (const bool add, const FieldField< Field, LUType > &interfaceCoeffs, const Field< Type > &psiif, Field< Type > &result) const
 Initialise the update of interfaced interfaces. More...
 
void updateMatrixInterfaces (const bool add, const FieldField< Field, LUType > &interfaceCoeffs, const Field< Type > &psiif, Field< Type > &result, const label startRequest) const
 Update interfaced interfaces for matrix operations. More...
 
tmp< Field< Type > > H (const Field< Type > &) const
 
tmp< Field< Type > > H (const tmp< Field< Type >> &) const
 
tmp< Field< Type > > faceH (const Field< Type > &) const
 
tmp< Field< Type > > faceH (const tmp< Field< Type >> &) const
 
void operator= (const LduMatrix< Type, DType, LUType > &)
 Copy assignment. More...
 
void operator= (LduMatrix< Type, DType, LUType > &&)
 Move assignment. More...
 
void negate ()
 
void operator+= (const LduMatrix< Type, DType, LUType > &)
 
void operator-= (const LduMatrix< Type, DType, LUType > &)
 
void operator*= (const scalarField &)
 
void operator*= (scalar)
 

Friends

class SolverPerformance< Type >
 
Ostreamoperator (Ostream &, const LduMatrix< Type, DType, LUType > &)
 

Detailed Description

template<class Type, class DType, class LUType>
class Foam::LduMatrix< Type, DType, LUType >

LduMatrix is a general matrix class in which the coefficients are stored as three arrays, one for the upper triangle, one for the lower triangle and a third for the diagonal.

Addressing arrays must be supplied for the upper and lower triangles.

Note
It might be better if this class were organised as a hierarchy starting from an empty matrix, then deriving diagonal, symmetric and asymmetric matrices.
Source files

Definition at line 68 of file LduMatrix.H.

Constructor & Destructor Documentation

◆ LduMatrix() [1/5]

LduMatrix ( const lduMesh mesh)

Construct given an LDU addressed mesh.

The coefficients are initially empty for subsequent setting. Not yet 'explicit' (legacy code may rely on implicit construct)

Definition at line 28 of file LduMatrix.C.

◆ LduMatrix() [2/5]

LduMatrix ( const LduMatrix< Type, DType, LUType > &  )

Copy construct.

◆ LduMatrix() [3/5]

LduMatrix ( LduMatrix< Type, DType, LUType > &&  )

Move construct.

◆ LduMatrix() [4/5]

LduMatrix ( LduMatrix< Type, DType, LUType > &  ,
bool  reuse 
)

Construct as copy or re-use as specified.

◆ LduMatrix() [5/5]

LduMatrix ( const lduMesh mesh,
Istream is 
)

Construct given an LDU addressed mesh and an Istream from which the coefficients are read.

Definition at line 117 of file LduMatrix.C.

◆ ~LduMatrix()

~LduMatrix ( )
default

Destructor.

Member Function Documentation

◆ ClassName()

ClassName ( "LduMatrix< Type, DType, LUType >"  )

◆ mesh()

const lduMesh& mesh ( ) const
inlinenoexcept

Return the LDU mesh from which the addressing is obtained.

Definition at line 608 of file LduMatrix.H.

◆ lduAddr()

const lduAddressing& lduAddr ( ) const
inline

Return the LDU addressing.

Definition at line 616 of file LduMatrix.H.

Referenced by TDILUPreconditioner< Type, DType, LUType >::calcInvD(), and TGaussSeidelSmoother< Type, DType, LUType >::smooth().

Here is the caller graph for this function:

◆ patchSchedule()

const lduSchedule& patchSchedule ( ) const
inline

Return the patch evaluation schedule.

Definition at line 624 of file LduMatrix.H.

◆ interfaces() [1/2]

const LduInterfaceFieldPtrsList<Type>& interfaces ( ) const
inlinenoexcept

Const access to the interfaces.

Definition at line 632 of file LduMatrix.H.

◆ interfaces() [2/2]

LduInterfaceFieldPtrsList<Type>& interfaces ( )
inlinenoexcept

Non-const access to the interfaces.

Definition at line 640 of file LduMatrix.H.

◆ diag() [1/2]

const Foam::Field< DType > & diag ( ) const

Definition at line 151 of file LduMatrix.C.

References Foam::abort(), Foam::FatalError, and FatalErrorInFunction.

Referenced by fvMatrix< Type >::solveCoupled(), and TGaussSeidelSmoother< Type, DType, LUType >::TGaussSeidelSmoother().

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

◆ upper() [1/2]

◆ lower() [1/2]

◆ source() [1/2]

const Foam::Field< Type > & source ( ) const

Definition at line 267 of file LduMatrix.C.

References Foam::abort(), Foam::FatalError, and FatalErrorInFunction.

Referenced by TGaussSeidelSmoother< Type, DType, LUType >::smooth().

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

◆ diag() [2/2]

Foam::Field< DType > & diag ( )

Definition at line 165 of file LduMatrix.C.

◆ upper() [2/2]

Foam::Field< LUType > & upper ( )

Definition at line 199 of file LduMatrix.C.

◆ lower() [2/2]

Foam::Field< LUType > & lower ( )

Definition at line 244 of file LduMatrix.C.

◆ source() [2/2]

Foam::Field< Type > & source ( )

Definition at line 281 of file LduMatrix.C.

◆ interfacesUpper() [1/2]

const FieldField<Field, LUType>& interfacesUpper ( ) const
inlinenoexcept

Definition at line 661 of file LduMatrix.H.

Referenced by TGaussSeidelSmoother< Type, DType, LUType >::smooth().

Here is the caller graph for this function:

◆ interfacesLower() [1/2]

const FieldField<Field, LUType>& interfacesLower ( ) const
inlinenoexcept

Definition at line 666 of file LduMatrix.H.

◆ interfacesUpper() [2/2]

FieldField<Field, LUType>& interfacesUpper ( )
inlinenoexcept

Definition at line 671 of file LduMatrix.H.

◆ interfacesLower() [2/2]

FieldField<Field, LUType>& interfacesLower ( )
inlinenoexcept

Definition at line 676 of file LduMatrix.H.

◆ matrixTypeName()

Foam::word matrixTypeName ( ) const

The matrix type (empty, diagonal, symmetric, ...)

Definition at line 133 of file LduMatrix.C.

◆ hasDiag()

bool hasDiag ( ) const
inlinenoexcept

Definition at line 689 of file LduMatrix.H.

◆ hasUpper()

bool hasUpper ( ) const
inlinenoexcept

Definition at line 690 of file LduMatrix.H.

◆ hasLower()

bool hasLower ( ) const
inlinenoexcept

Definition at line 691 of file LduMatrix.H.

◆ hasSource()

bool hasSource ( ) const
inlinenoexcept

Definition at line 692 of file LduMatrix.H.

◆ diagonal()

bool diagonal ( ) const
inlinenoexcept

Matrix has diagonal only.

Definition at line 697 of file LduMatrix.H.

Referenced by LduMatrix< Type, DType, LUType >::solver::New().

Here is the caller graph for this function:

◆ symmetric()

bool symmetric ( ) const
inlinenoexcept

Matrix is symmetric.

Definition at line 705 of file LduMatrix.H.

Referenced by LduMatrix< Type, DType, LUType >::solver::New(), and LduMatrix< Type, DType, LUType >::smoother::New().

Here is the caller graph for this function:

◆ asymmetric()

bool asymmetric ( ) const
inlinenoexcept

Matrix is asymmetric (ie, full)

Definition at line 713 of file LduMatrix.H.

Referenced by LduMatrix< Type, DType, LUType >::solver::New(), and LduMatrix< Type, DType, LUType >::smoother::New().

Here is the caller graph for this function:

◆ sumDiag()

void sumDiag ( )

Definition at line 27 of file LduMatrixOperations.C.

References Foam::diag(), LduMatrix< Type, DType, LUType >::lower(), UList< T >::size(), and LduMatrix< Type, DType, LUType >::upper().

Here is the call graph for this function:

◆ negSumDiag()

void negSumDiag ( )

Definition at line 45 of file LduMatrixOperations.C.

References Foam::diag(), LduMatrix< Type, DType, LUType >::lower(), UList< T >::size(), and LduMatrix< Type, DType, LUType >::upper().

Here is the call graph for this function:

◆ sumMagOffDiag()

void sumMagOffDiag ( Field< LUType > &  sumOff) const

Definition at line 64 of file LduMatrixOperations.C.

References Foam::cmptMag(), LduMatrix< Type, DType, LUType >::lower(), UList< T >::size(), and LduMatrix< Type, DType, LUType >::upper().

Here is the call graph for this function:

◆ Amul()

void Amul ( Field< Type > &  Apsi,
const tmp< Field< Type >> &  tpsi 
) const

Matrix multiplication.

Definition at line 28 of file LduMatrixATmul.C.

References UList< Type >::begin(), Foam::diag(), Foam::dot(), Foam::stringOps::lower(), psi, and Foam::stringOps::upper().

Here is the call graph for this function:

◆ Tmul()

void Tmul ( Field< Type > &  Tpsi,
const tmp< Field< Type >> &  tpsi 
) const

Matrix transpose multiplication.

Definition at line 87 of file LduMatrixATmul.C.

References UList< Type >::begin(), Foam::diag(), Foam::dot(), Foam::stringOps::lower(), psi, and Foam::stringOps::upper().

Here is the call graph for this function:

◆ sumA()

void sumA ( Field< Type > &  sumA) const

Sum the coefficients on each row of the matrix.

Definition at line 145 of file LduMatrixATmul.C.

References UList< Type >::begin(), Foam::diag(), Foam::dot(), forAll, Foam::stringOps::lower(), and Foam::stringOps::upper().

Here is the call graph for this function:

◆ residual() [1/2]

void residual ( Field< Type > &  rA,
const Field< Type > &  psi 
) const

Definition at line 193 of file LduMatrixATmul.C.

References UList< Type >::begin(), Foam::diag(), Foam::dot(), Foam::stringOps::lower(), psi, and Foam::stringOps::upper().

Here is the call graph for this function:

◆ residual() [2/2]

Foam::tmp< Foam::Field< Type > > residual ( const Field< Type > &  psi) const

Definition at line 253 of file LduMatrixATmul.C.

References Foam::New(), and psi.

Here is the call graph for this function:

◆ initMatrixInterfaces()

void initMatrixInterfaces ( const bool  add,
const FieldField< Field, LUType > &  interfaceCoeffs,
const Field< Type > &  psiif,
Field< Type > &  result 
) const

Initialise the update of interfaced interfaces.

for matrix operations

Definition at line 29 of file LduMatrixUpdateMatrixInterfaces.C.

References Foam::add(), Foam::exit(), Foam::FatalError, FatalErrorInFunction, forAll, and UList< T >::size().

Referenced by TGaussSeidelSmoother< Type, DType, LUType >::smooth().

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

◆ updateMatrixInterfaces()

void updateMatrixInterfaces ( const bool  add,
const FieldField< Field, LUType > &  interfaceCoeffs,
const Field< Type > &  psiif,
Field< Type > &  result,
const label  startRequest 
) const

Update interfaced interfaces for matrix operations.

Definition at line 101 of file LduMatrixUpdateMatrixInterfaces.C.

References Foam::add(), Foam::exit(), Foam::FatalError, FatalErrorInFunction, and forAll.

Referenced by TGaussSeidelSmoother< Type, DType, LUType >::smooth().

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

◆ H() [1/2]

Foam::tmp< Foam::Field< Type > > H ( const Field< Type > &  psi) const

Definition at line 84 of file LduMatrixOperations.C.

References Foam::stringOps::lower(), Foam::New(), psi, and Foam::stringOps::upper().

Here is the call graph for this function:

◆ H() [2/2]

Foam::tmp< Foam::Field< Type > > H ( const tmp< Field< Type >> &  tpsi) const

Definition at line 114 of file LduMatrixOperations.C.

References H().

Here is the call graph for this function:

◆ faceH() [1/2]

Foam::tmp< Foam::Field< Type > > faceH ( const Field< Type > &  psi) const

◆ faceH() [2/2]

Foam::tmp< Foam::Field< Type > > faceH ( const tmp< Field< Type >> &  tpsi) const

Definition at line 147 of file LduMatrixOperations.C.

◆ operator=() [1/2]

void operator= ( const LduMatrix< Type, DType, LUType > &  )

Copy assignment.

◆ operator=() [2/2]

void operator= ( LduMatrix< Type, DType, LUType > &&  )

Move assignment.

◆ negate()

void negate ( )

Definition at line 217 of file LduMatrixOperations.C.

References Foam::negate().

Here is the call graph for this function:

◆ operator+=()

void operator+= ( const LduMatrix< Type, DType, LUType > &  )

Definition at line 245 of file LduMatrixOperations.C.

References A, Foam::abort(), Foam::diag(), Foam::FatalError, FatalErrorInFunction, Foam::stringOps::lower(), and Foam::stringOps::upper().

Here is the call graph for this function:

◆ operator-=()

void operator-= ( const LduMatrix< Type, DType, LUType > &  )

Definition at line 322 of file LduMatrixOperations.C.

References A, Foam::abort(), Foam::diag(), Foam::FatalError, FatalErrorInFunction, Foam::stringOps::lower(), and Foam::stringOps::upper().

Here is the call graph for this function:

◆ operator*=() [1/2]

void operator*= ( const scalarField sf)

Definition at line 400 of file LduMatrixOperations.C.

References Foam::abort(), Foam::FatalError, FatalErrorInFunction, Foam::stringOps::lower(), and Foam::stringOps::upper().

Here is the call graph for this function:

◆ operator*=() [2/2]

void operator*= ( scalar  s)

Definition at line 447 of file LduMatrixOperations.C.

References s.

Friends And Related Function Documentation

◆ SolverPerformance< Type >

friend class SolverPerformance< Type >
friend

Definition at line 130 of file LduMatrix.H.

◆ operator

Ostream& operator ( Ostream ,
const LduMatrix< Type, DType, LUType > &   
)
friend

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