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 &)
 Construct given an LDU addressed mesh. More...
 
 LduMatrix (const LduMatrix< Type, DType, LUType > &)
 Copy constructor. More...
 
 LduMatrix (LduMatrix< Type, DType, LUType > &, bool reuse)
 Copy constructor or re-use as specified. More...
 
 LduMatrix (const lduMesh &, Istream &)
 Construct given an LDU addressed mesh and an Istream. More...
 
 ~LduMatrix ()
 
const lduMeshmesh () const
 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
 Return interfaces. More...
 
LduInterfaceFieldPtrsList< Type > & interfaces ()
 Return interfaces. More...
 
Field< DType > & diag ()
 
Field< LUType > & upper ()
 
Field< LUType > & lower ()
 
Field< Type > & source ()
 
FieldField< Field, LUType > & interfacesUpper ()
 
FieldField< Field, LUType > & interfacesLower ()
 
const Field< DType > & diag () const
 
const Field< LUType > & upper () const
 
const Field< LUType > & lower () const
 
const Field< Type > & source () const
 
const FieldField< Field, LUType > & interfacesUpper () const
 
const FieldField< Field, LUType > & interfacesLower () const
 
bool hasDiag () const
 
bool hasUpper () const
 
bool hasLower () const
 
bool hasSource () const
 
bool diagonal () const
 
bool symmetric () const
 
bool asymmetric () const
 
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 FieldField< Field, LUType > &interfaceCoeffs, const Field< Type > &psiif, Field< Type > &result) const
 Initialise the update of interfaced interfaces. More...
 
void updateMatrixInterfaces (const FieldField< Field, LUType > &interfaceCoeffs, const Field< Type > &psiif, Field< Type > &result) 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 > &)
 
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 69 of file LduMatrix.H.

Constructor & Destructor Documentation

◆ LduMatrix() [1/4]

LduMatrix ( const lduMesh mesh)

Construct given an LDU addressed mesh.

The coefficients are initially empty for subsequent setting.

Definition at line 32 of file LduMatrix.C.

Referenced by LduMatrix< Type, DType, LUType >::preconditioner::preconditionT().

Here is the caller graph for this function:

◆ LduMatrix() [2/4]

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

Copy constructor.

◆ LduMatrix() [3/4]

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

Copy constructor or re-use as specified.

◆ LduMatrix() [4/4]

LduMatrix ( const lduMesh mesh,
Istream is 
)

Construct given an LDU addressed mesh and an Istream.

from which the coefficients are read

Definition at line 144 of file LduMatrix.C.

◆ ~LduMatrix()

~LduMatrix ( )

Definition at line 163 of file LduMatrix.C.

Referenced by LduMatrix< Type, DType, LUType >::preconditioner::preconditionT().

Here is the caller graph for this function:

Member Function Documentation

◆ ClassName()

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

Referenced by LduMatrix< Type, DType, LUType >::preconditioner::preconditionT().

Here is the caller graph for this function:

◆ mesh()

const lduMesh& mesh ( ) const
inline

Return the LDU mesh from which the addressing is obtained.

Definition at line 492 of file LduMatrix.H.

◆ lduAddr()

◆ patchSchedule()

const lduSchedule& patchSchedule ( ) const
inline

Return the patch evaluation schedule.

Definition at line 504 of file LduMatrix.H.

References LduMatrix< Type, DType, LUType >::lduAddr(), and lduAddressing::patchSchedule().

Here is the call graph for this function:

◆ interfaces() [1/2]

const LduInterfaceFieldPtrsList<Type>& interfaces ( ) const
inline

Return interfaces.

Definition at line 510 of file LduMatrix.H.

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

Here is the caller graph for this function:

◆ interfaces() [2/2]

LduInterfaceFieldPtrsList<Type>& interfaces ( )
inline

◆ diag() [1/2]

◆ upper() [1/2]

◆ lower() [1/2]

◆ source() [1/2]

◆ interfacesUpper() [1/2]

FieldField<Field, LUType>& interfacesUpper ( )
inline

Definition at line 529 of file LduMatrix.H.

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

Here is the caller graph for this function:

◆ interfacesLower() [1/2]

FieldField<Field, LUType>& interfacesLower ( )
inline

◆ diag() [2/2]

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

Definition at line 260 of file LduMatrix.C.

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

Here is the call graph for this function:

◆ upper() [2/2]

const Foam::Field< LUType > & upper ( ) const

Definition at line 274 of file LduMatrix.C.

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

Here is the call graph for this function:

◆ lower() [2/2]

const Foam::Field< LUType > & lower ( ) const

Definition at line 295 of file LduMatrix.C.

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

Here is the call graph for this function:

◆ source() [2/2]

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

Definition at line 316 of file LduMatrix.C.

References Foam::abort(), IOstream::check(), Foam::endl(), Foam::FatalError, and FatalErrorInFunction.

Here is the call graph for this function:

◆ interfacesUpper() [2/2]

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

Definition at line 545 of file LduMatrix.H.

◆ interfacesLower() [2/2]

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

Definition at line 550 of file LduMatrix.H.

◆ hasDiag()

bool hasDiag ( ) const
inline

Definition at line 556 of file LduMatrix.H.

◆ hasUpper()

bool hasUpper ( ) const
inline

Definition at line 561 of file LduMatrix.H.

◆ hasLower()

bool hasLower ( ) const
inline

Definition at line 566 of file LduMatrix.H.

◆ hasSource()

bool hasSource ( ) const
inline

Definition at line 571 of file LduMatrix.H.

◆ diagonal()

bool diagonal ( ) const
inline

Definition at line 576 of file LduMatrix.H.

References UPstream::parRun(), and List< T >::size().

Referenced by LduMatrix< Type, DType, LUType >::solver::New(), LduMatrix< Type, DType, LUType >::operator+=(), and LduMatrix< Type, DType, LUType >::operator-=().

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

◆ symmetric()

bool symmetric ( ) const
inline

Definition at line 590 of file LduMatrix.H.

References UPstream::parRun(), and List< T >::size().

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

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

◆ asymmetric()

◆ sumDiag()

void sumDiag ( )

Definition at line 31 of file LduMatrixOperations.C.

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

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

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

◆ negSumDiag()

void negSumDiag ( )

Definition at line 54 of file LduMatrixOperations.C.

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

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

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

◆ sumMagOffDiag()

void sumMagOffDiag ( Field< LUType > &  sumOff) const

Definition at line 78 of file LduMatrixOperations.C.

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

Referenced by LduMatrix< Type, DType, LUType >::asymmetric(), and LduMatrix< Type, DType, LUType >::negSumDiag().

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

◆ Amul()

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

Matrix multiplication.

Definition at line 64 of file LduMatrixATmul.C.

References UList< Type >::begin(), Foam::diag(), Foam::dot(), psi, and LduMatrix< Type, DType, LUType >::Tmul().

Referenced by Amultiplier< Type, LUType >::addAmul(), and LduMatrix< Type, DType, LUType >::asymmetric().

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

◆ Tmul()

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

Matrix transpose multiplication.

Definition at line 118 of file LduMatrixATmul.C.

References UList< Type >::begin(), Foam::diag(), Foam::dot(), psi, and LduMatrix< Type, DType, LUType >::sumA().

Referenced by LduMatrix< Type, DType, LUType >::Amul(), and LduMatrix< Type, DType, LUType >::asymmetric().

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

◆ sumA()

void sumA ( Field< Type > &  sumA) const

Sum the coefficients on each row of the matrix.

Definition at line 171 of file LduMatrixATmul.C.

References UList< Type >::begin(), Foam::diag(), Foam::dot(), forAll, patchi, and LduMatrix< Type, DType, LUType >::residual().

Referenced by LduMatrix< Type, DType, LUType >::asymmetric(), and LduMatrix< Type, DType, LUType >::Tmul().

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

◆ residual() [1/2]

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

Definition at line 219 of file LduMatrixATmul.C.

References UList< Type >::begin(), Foam::diag(), Foam::dot(), forAll, patchi, and source().

Referenced by LduMatrix< Type, DType, LUType >::asymmetric(), and LduMatrix< Type, DType, LUType >::sumA().

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

◆ residual() [2/2]

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

Definition at line 284 of file LduMatrixATmul.C.

References psi, tmp< T >::ref(), and List< Type >::size().

Here is the call graph for this function:

◆ initMatrixInterfaces()

void initMatrixInterfaces ( 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 33 of file LduMatrixUpdateMatrixInterfaces.C.

References Foam::exit(), Foam::FatalError, FatalErrorInFunction, forAll, List< T >::size(), and LduMatrix< Type, DType, LUType >::updateMatrixInterfaces().

Referenced by LduMatrix< Type, DType, LUType >::asymmetric(), and 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 FieldField< Field, LUType > &  interfaceCoeffs,
const Field< Type > &  psiif,
Field< Type > &  result 
) const

Update interfaced interfaces for matrix operations.

Definition at line 98 of file LduMatrixUpdateMatrixInterfaces.C.

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

Referenced by LduMatrix< Type, DType, LUType >::asymmetric(), LduMatrix< Type, DType, LUType >::initMatrixInterfaces(), and 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 103 of file LduMatrixOperations.C.

References UList< Type >::begin(), and Foam::Zero.

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

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

◆ H() [2/2]

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

Definition at line 138 of file LduMatrixOperations.C.

◆ faceH() [1/2]

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

Definition at line 148 of file LduMatrixOperations.C.

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

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

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

◆ faceH() [2/2]

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

Definition at line 171 of file LduMatrixOperations.C.

◆ operator=()

◆ negate()

void negate ( )

Definition at line 227 of file LduMatrixOperations.C.

References Foam::negate().

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

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

◆ operator+=()

◆ operator-=()

◆ operator*=() [1/2]

void operator*= ( const scalarField sf)

Definition at line 410 of file LduMatrixOperations.C.

References Foam::abort(), Foam::FatalError, FatalErrorInFunction, sf(), and List< T >::size().

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

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

◆ operator*=() [2/2]

void operator*= ( scalar  s)

Definition at line 457 of file LduMatrixOperations.C.

References s().

Here is the call graph for this function:

Friends And Related Function Documentation

◆ SolverPerformance< Type >

friend class SolverPerformance< Type >
friend

Definition at line 109 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: