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...
Classes | |
class | fvSolver |
Public Types | |
typedef GeometricField< Type, fvPatchField, volMesh > | psiFieldType |
The geometric field type for psi. More... | |
typedef GeometricField< Type, fvsPatchField, surfaceMesh > | faceFluxFieldType |
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 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... | |
lduPrimitiveMeshAssembly * | lduMeshPtr () |
Access to lduPrimitiveMeshAssembly. More... | |
const lduPrimitiveMeshAssembly * | lduMeshPtr () 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 lduPrimitiveMeshAssembly & | lduMeshAssembly () |
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 dimensionSet & | dimensions () const noexcept |
Field< Type > & | source () noexcept |
const Field< Type > & | source () const noexcept |
const FieldField< Field, Type > & | internalCoeffs () const noexcept |
FieldField< Field, Type > & | internalCoeffs () noexcept |
const FieldField< Field, Type > & | boundaryCoeffs () const noexcept |
FieldField< Field, Type > & | boundaryCoeffs () noexcept |
faceFluxFieldPtrType & | faceFluxCorrectionPtr () |
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) |
void | setValues (const labelUList &cellLabels, const UList< Type > &values) |
void | setValues (const labelUList &cellLabels, const UIndirectList< Type > &values) |
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) |
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< fvSolver > | solver (const dictionary &) |
Construct and return the solver. More... | |
autoPtr< fvSolver > | solver () |
Construct and return the solver. More... | |
SolverPerformance< Type > | solveSegregatedOrCoupled (const dictionary &) |
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 () |
Solve returning the solution statistics. More... | |
tmp< Field< Type > > | residual () const |
Return the matrix residual. More... | |
tmp< scalarField > | D () const |
Return the matrix scalar diagonal. More... | |
tmp< Field< Type > > | DD () const |
Return the matrix Type diagonal. More... | |
tmp< volScalarField > | A () const |
Return the central coefficient. More... | |
tmp< GeometricField< Type, fvPatchField, volMesh > > | H () const |
Return the H operation source. More... | |
tmp< volScalarField > | H1 () const |
Return H(1) More... | |
tmp< GeometricField< Type, fvsPatchField, surfaceMesh > > | flux () const |
Return the face-flux field from the matrix. More... | |
const dictionary & | solverDict () const |
Return the solver dictionary 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 > &) |
void | setComponentReference (const label patchi, const label facei, const direction, const scalar value) |
Foam::autoPtr< Foam::fvMatrix< Foam::scalar >::fvSolver > | solver (const dictionary &solverControls) |
Foam::solverPerformance | solveSegregated (const dictionary &solverControls) |
Foam::tmp< Foam::scalarField > | residual () const |
Foam::tmp< Foam::volScalarField > | H () const |
Foam::tmp< Foam::volScalarField > | H1 () const |
void | setComponentReference (const label patchi, const label facei, const direction, const scalar value) |
autoPtr< fvMatrix< scalar >::fvSolver > | solver (const dictionary &) |
solverPerformance | solveSegregated (const dictionary &) |
tmp< scalarField > | residual () const |
tmp< volScalarField > | H () const |
tmp< volScalarField > | H1 () const |
Public Member Functions inherited from refCount | |
constexpr | refCount () noexcept |
Default construct, initializing count to 0. More... | |
int | 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 &) | |
~lduMatrix () | |
Destructor. More... | |
const lduMesh & | mesh () const |
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 lduAddressing & | lduAddr () const |
Return the LDU addressing. More... | |
const lduSchedule & | patchSchedule () const |
Return the patch evaluation schedule. More... | |
scalarField & | lower () |
scalarField & | diag () |
scalarField & | upper () |
scalarField & | lower (const label size) |
scalarField & | diag (const label nCoeffs) |
scalarField & | upper (const label nCoeffs) |
const scalarField & | lower () const |
const scalarField & | diag () const |
const scalarField & | upper () const |
bool | hasDiag () const |
bool | hasUpper () const |
bool | hasLower () const |
bool | diagonal () const |
bool | symmetric () const |
bool | asymmetric () const |
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< solveScalarField > | residual (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 |
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 |
template<class Type > | |
tmp< Field< Type > > | H (const Field< Type > &) const |
template<class Type > | |
tmp< Field< Type > > | H (const tmp< Field< Type > > &) const |
tmp< scalarField > | H1 () 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< lduMatrix > | info () const |
Return info proxy. 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 > > &) |
Ostream & | operator (Ostream &, const 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.
Definition at line 117 of file fvMatrix.H.
typedef GeometricField<Type, fvPatchField, volMesh> psiFieldType |
The geometric field type for psi.
Definition at line 129 of file fvMatrix.H.
typedef GeometricField<Type, fvsPatchField, surfaceMesh> faceFluxFieldType |
Field type for face flux (for non-orthogonal correction)
Definition at line 134 of file fvMatrix.H.
typedef GeometricField<Type, fvsPatchField, surfaceMesh>* faceFluxFieldPtrType |
Declare return type of the faceFluxCorrectionPtr() function.
Definition at line 498 of file fvMatrix.H.
fvMatrix | ( | const GeometricField< Type, fvPatchField, volMesh > & | psi, |
const dimensionSet & | ds | ||
) |
Construct given a field to solve for.
Definition at line 343 of file fvMatrix.C.
References faMesh::boundary(), fvMatrix< Type >::checkImplicit(), DebugInFunction, Foam::endl(), forAll, DimensionedField< Type, GeoMesh >::mesh(), IOobject::name(), psi, UPtrList< T >::size(), and Foam::Zero.
Copy construct.
Definition at line 388 of file fvMatrix.C.
References DebugInFunction, Foam::endl(), and IOobject::name().
Copy/move construct from tmp<fvMatrix<Type>>
Definition at line 416 of file fvMatrix.C.
References DebugInFunction, Foam::endl(), and IOobject::name().
|
inline |
Deprecated(2022-05) - construct with dimensionSet instead.
Definition at line 313 of file fvMatrix.H.
|
virtual |
Destructor.
Definition at line 456 of file fvMatrix.C.
References DebugInFunction, Foam::deleteDemandDrivenData(), and Foam::endl().
|
protected |
Add patch contribution to internal field.
Definition at line 48 of file fvMatrix.C.
References Foam::abort(), Foam::endl(), Foam::FatalError, FatalErrorInFunction, forAll, and UList< T >::size().
|
protected |
Definition at line 72 of file fvMatrix.C.
|
protected |
Subtract patch contribution from internal field.
Definition at line 86 of file fvMatrix.C.
References Foam::abort(), Foam::endl(), Foam::FatalError, FatalErrorInFunction, forAll, and UList< T >::size().
|
protected |
Definition at line 110 of file fvMatrix.C.
|
protected |
Name the implicit assembly addressing.
Definition at line 307 of file fvMatrix.C.
References Foam::endl(), forAll, Foam::name(), Foam::Pout, and psi.
Referenced by fvMatrix< Type >::fvMatrix().
|
protected |
Definition at line 123 of file fvMatrix.C.
References Foam::component(), Foam::diag(), forAll, and psi.
Referenced by fvMatrix< Type >::residual(), and fvMatrix< Type >::solveSegregated().
|
protected |
Definition at line 152 of file fvMatrix.C.
References Foam::cmptAv(), Foam::diag(), forAll, and psi.
Definition at line 176 of file fvMatrix.C.
References Foam::cmptMultiply(), fvPatchField< Type >::coupled(), forAll, fvPatchField< Type >::patchNeighbourField(), and psi.
Referenced by fvMatrix< Type >::H(), fvMatrix< Type >::residual(), and fvMatrix< Type >::solveSegregated().
|
protected |
Set solution in given cells to the specified values.
Definition at line 226 of file fvMatrix.C.
References cells, Foam::diag(), forAll, mesh, primitiveFieldRef(), and psi.
ClassName | ( | "fvMatrix< Type >" | ) |
Construct and return a clone.
Definition at line 323 of file fvMatrix.H.
References Foam::New().
Referenced by fvMatrix< Type >::addFvMatrix().
|
inline |
Definition at line 339 of file fvMatrix.H.
Referenced by fvMatrix< Type >::solveSegregated().
|
inline |
Definition at line 344 of file fvMatrix.H.
Referenced by fvMatrix< Type >::psi().
|
inline |
Definition at line 349 of file fvMatrix.H.
|
inline |
Definition at line 354 of file fvMatrix.H.
References fvMatrix< Type >::lduMeshPtr(), and lduPrimitiveMeshAssembly::patchMap().
void transferFvMatrixCoeffs |
Transfer lower, upper, diag and source to this fvMatrix.
Definition at line 813 of file fvMatrix.C.
References Foam::diag(), Foam::faceMap(), forAll, List< T >::setSize(), UList< T >::size(), and Foam::Zero.
Referenced by fvMatrix< Type >::solveSegregated().
void createOrUpdateLduPrimitiveAssembly |
Create or update ldu assembly.
Definition at line 906 of file fvMatrix.C.
References lduAddressing::clearOut(), Foam::endl(), Foam::Info, io(), lduPrimitiveMesh::lduAddr(), DimensionedField< Type, GeoMesh >::mesh(), IOobject::NO_READ, IOobject::NO_WRITE, psi, UPtrList< T >::set(), regIOobject::store(), and lduPrimitiveMeshAssembly::update().
Referenced by fvMatrix< Type >::solveSegregated().
Foam::lduPrimitiveMeshAssembly * lduMeshPtr |
Access to lduPrimitiveMeshAssembly.
Definition at line 880 of file fvMatrix.C.
References lduPrimitiveMesh::mesh(), and lduMesh::thisDb().
Referenced by fvMatrix< Type >::globalPatchID(), fvMatrix< Type >::lduMeshAssembly(), and fvMatrix< Type >::solveSegregated().
const Foam::lduPrimitiveMeshAssembly * lduMeshPtr |
Const Access to lduPrimitiveMeshAssembly.
Definition at line 893 of file fvMatrix.C.
void manipulateMatrix | ( | direction | cmp | ) |
Manipulate matrix.
Definition at line 790 of file fvMatrix.C.
References GeometricField< Type, PatchField, GeoMesh >::boundaryFieldRef(), forAll, and psi.
Referenced by fvMatrix< Type >::solveSegregated().
void setBounAndInterCoeffs |
Manipulate boundary/internal coeffs for coupling.
Definition at line 660 of file fvMatrix.C.
References boundary, GeometricField< Type, PatchField, GeoMesh >::boundaryField(), forAll, DimensionedField< Type, GeoMesh >::mesh(), psi, List< T >::set(), and Foam::Zero.
Referenced by fvMatrix< Type >::solveSegregated().
void setInterfaces | ( | lduInterfaceFieldPtrsList & | interfaces, |
PtrDynList< lduInterfaceField > & | newInterfaces | ||
) |
Set interfaces.
Definition at line 469 of file fvMatrix.C.
References PtrDynList< T, SizeMin >::append(), GeometricField< Type, PatchField, GeoMesh >::boundaryField(), forAll, UPtrList< T >::last(), psi, UPtrList< T >::set(), and UPtrList< T >::setSize().
Referenced by fvMatrix< Type >::solveSegregated().
void mapContributions | ( | label | fieldi, |
const FieldField< Field, Type > & | fluxContrib, | ||
FieldField< Field, Type > & | contrib, | ||
bool | internal | ||
) | const |
Add internal and boundary contribution to local patches.
Definition at line 547 of file fvMatrix.C.
References lduPrimitiveMeshAssembly::cellBoundMap(), cellId, Foam::cmptMultiply(), coupled(), faceId(), lduPrimitiveMeshAssembly::facePatchFaceMap(), forAll, polyPatch::masterImplicit(), DimensionedField< Type, GeoMesh >::mesh(), polyPatch::neighbPolyPatchID(), lduPrimitiveMeshAssembly::patchLocalToGlobalMap(), lduPrimitiveMeshAssembly::patchMap(), psi, List< T >::setSize(), and PtrList< T >::setSize().
|
inline |
Return optional lduAdressing.
Definition at line 405 of file fvMatrix.H.
References fvMatrix< Type >::lduMeshPtr().
Referenced by mixedEnergyFvPatchScalarField::manipulateMatrix(), cyclicFvPatchField< Type >::manipulateMatrix(), cyclicACMIFvPatchField< Type >::manipulateMatrix(), and cyclicAMIFvPatchField< Type >::manipulateMatrix().
|
inline |
Return psi.
Definition at line 411 of file fvMatrix.H.
References fvMatrix< Type >::matrix(), and psi.
Referenced by MRFZone::addCoriolis(), velocityDampingConstraint::addDamping(), dynamicOversetFvMesh::addInterpolation(), solver::addOptimisationTypeSource(), PhaseLimitStabilization< Type >::addSup(), jouleHeatingSource::addSup(), patchCellsSource::addSup(), radiation::addSup(), SemiImplicitSource< Type >::addSup(), multiphaseStabilizedTurbulence::addSup(), interRegionHeatTransferModel::addSup(), multiphaseMangrovesTurbulenceModel::addSup(), atmCoriolisUSource::addSup(), atmPlantCanopyUSource::addSup(), acousticDampingSource::addSup(), explicitPorositySource::addSup(), radialActuationDiskSource::addSup(), rotorDiskSource::addSup(), solidificationMeltingSource::addSup(), interRegionExplicitPorositySource::addSup(), multiphaseMangrovesSource::addSup(), Foam::checkMethod(), optionList::constrain(), cyclicFvPatchField< Type >::manipulateMatrix(), cyclicACMIFvPatchField< Type >::manipulateMatrix(), cyclicAMIFvPatchField< Type >::manipulateMatrix(), MomentumTransferPhaseSystem< BasePhaseSystem >::momentumTransfer(), dynamicOversetFvMesh::normalisation(), dynamicOversetFvMesh::solve(), and dynamicOversetFvMesh::write().
|
inline |
Definition at line 423 of file fvMatrix.H.
References fvMatrix< Type >::matrix(), and psi.
|
inline |
Clear multiple fvMatrices.
Definition at line 446 of file fvMatrix.H.
References PtrList< T >::clear().
|
inlinenoexcept |
Definition at line 453 of file fvMatrix.H.
Referenced by fvMatrix< Type >::addFvMatrix(), patchCellsSource::addSup(), SemiImplicitSource< Type >::addSup(), directionalPressureGradientExplicitSource::addSup(), explicitPorositySource::addSup(), meanVelocityForce::addSup(), rotorDiskSource::addSup(), interRegionExplicitPorositySource::addSup(), Foam::checkMethod(), and fvMatrix< Type >::operator*=().
|
inlinenoexcept |
Definition at line 458 of file fvMatrix.H.
Referenced by MRFZone::addCoriolis(), dynamicOversetFvMesh::addInterpolation(), effectivenessHeatExchangerSource::addSup(), radialActuationDiskSource::addSup(), solidificationMeltingSource::addSup(), interRegionExplicitPorositySource::addSup(), adjointkOmegaSST::addWallFunctionTerms(), EulerD2dt2Scheme< Type >::fvmD2dt2(), backwardDdtScheme< Type >::fvmDdt(), CoEulerDdtScheme< Type >::fvmDdt(), CrankNicolsonDdtScheme< Type >::fvmDdt(), EulerDdtScheme< Type >::fvmDdt(), localEulerDdtScheme< Type >::fvmDdt(), SLTSDdtScheme< Type >::fvmDdt(), gaussLaplacianScheme< Type, GType >::fvmLaplacian(), relaxedNonOrthoGaussLaplacianScheme< Type, GType >::fvmLaplacian(), ThermoCloud< Foam::DSMCCloud >::heatTransfer(), mixedEnergyFvPatchScalarField::manipulateMatrix(), kaqRWallFunctionFvPatchScalarField::manipulateMatrix(), adjointOutletVelocityFluxFvPatchVectorField::manipulateMatrix(), adjointWallVelocityFvPatchVectorField::manipulateMatrix(), ReactingCloud< CloudType >::SYi(), and dynamicOversetFvMesh::write().
|
inlinenoexcept |
Definition at line 463 of file fvMatrix.H.
|
inlinenoexcept |
fvBoundary scalar field containing pseudo-matrix coeffs for internal cells
Definition at line 470 of file fvMatrix.H.
Referenced by dynamicOversetFvMesh::addInterpolation(), gaussConvectionScheme< Type >::fvmDiv(), gaussLaplacianScheme< Type, GType >::fvmLaplacianUncorrected(), relaxedNonOrthoGaussLaplacianScheme< Type, GType >::fvmLaplacianUncorrected(), mixedEnergyFvPatchScalarField::manipulateMatrix(), waWallFunctionFvPatchScalarField::manipulateMatrix(), cyclicFvPatchField< Type >::manipulateMatrix(), cyclicACMIFvPatchField< Type >::manipulateMatrix(), cyclicAMIFvPatchField< Type >::manipulateMatrix(), dynamicOversetFvMesh::normalisation(), dynamicOversetFvMesh::solve(), and dynamicOversetFvMesh::write().
|
inlinenoexcept |
fvBoundary scalar field containing pseudo-matrix coeffs for internal cells
Definition at line 477 of file fvMatrix.H.
|
inlinenoexcept |
fvBoundary scalar field containing pseudo-matrix coeffs for boundary cells
Definition at line 484 of file fvMatrix.H.
Referenced by dynamicOversetFvMesh::addInterpolation(), gaussConvectionScheme< Type >::fvmDiv(), gaussLaplacianScheme< Type, GType >::fvmLaplacianUncorrected(), relaxedNonOrthoGaussLaplacianScheme< Type, GType >::fvmLaplacianUncorrected(), waWallFunctionFvPatchScalarField::manipulateMatrix(), cyclicFvPatchField< Type >::manipulateMatrix(), cyclicACMIFvPatchField< Type >::manipulateMatrix(), cyclicAMIFvPatchField< Type >::manipulateMatrix(), dynamicOversetFvMesh::solve(), and dynamicOversetFvMesh::write().
|
inlinenoexcept |
fvBoundary scalar field containing pseudo-matrix coeffs for boundary cells
Definition at line 491 of file fvMatrix.H.
|
inline |
Return pointer to face-flux non-orthogonal correction field.
Definition at line 501 of file fvMatrix.H.
Referenced by gaussLaplacianScheme< Type, GType >::fvmLaplacian(), and relaxedNonOrthoGaussLaplacianScheme< Type, GType >::fvmLaplacian().
|
inlinenoexcept |
True if face-flux non-orthogonal correction field exists.
Definition at line 507 of file fvMatrix.H.
References bool.
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.
Definition at line 968 of file fvMatrix.C.
Referenced by fixedTemperatureConstraint::constrain(), FixedValueConstraint< Type >::constrain(), epsilonWallFunctionFvPatchScalarField::manipulateMatrix(), omegaWallFunctionFvPatchScalarField::manipulateMatrix(), and fixedInternalValueFvPatchField< Type >::manipulateMatrix().
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 979 of file fvMatrix.C.
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 990 of file fvMatrix.C.
void setReference | ( | const label | celli, |
const Type & | value, | ||
const bool | forceReference = false |
||
) |
Set reference level for solution.
Definition at line 1001 of file fvMatrix.C.
References Foam::diag().
Referenced by Foam::CorrectPhi(), adjointSimple::mainIter(), and simple::mainIter().
void setReferences | ( | const labelUList & | cellLabels, |
const Type & | value, | ||
const bool | forceReference = false |
||
) |
Set reference level for solution.
Definition at line 1017 of file fvMatrix.C.
References cellId, Foam::diag(), and forAll.
void setReferences | ( | const labelUList & | cellLabels, |
const UList< Type > & | values, | ||
const bool | forceReference = false |
||
) |
Set reference level for solution.
Definition at line 1040 of file fvMatrix.C.
References cellId, Foam::diag(), and forAll.
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 37 of file fvMatrixSolve.C.
References Foam::diag().
void addFvMatrix | ( | fvMatrix< Type > & | matrix | ) |
Add fvMatrix.
Definition at line 1063 of file fvMatrix.C.
References Foam::abort(), fvMatrix< Type >::clone(), fvMatrix< Type >::dimensions(), Foam::endl(), Foam::FatalError, and FatalErrorInFunction.
void relax | ( | const scalar | alpha | ) |
Relax matrix (for steady-state solution).
alpha = 1 : diagonally equal alpha < 1 : diagonally dominant alpha = 0 : do nothing Note: Requires positive diagonal.
Definition at line 1092 of file fvMatrix.C.
References alpha, Foam::cmptMag(), Foam::cmptMax(), Foam::cmptMin(), Foam::component(), fvPatchField< Type >::coupled(), D, DebugInFunction, Foam::diag(), Foam::endl(), forAll, InfoInFunction, Foam::mag(), Foam::max(), UPstream::msgType(), Foam::nl, Foam::reduce(), Foam::returnReduce(), UList< T >::size(), and Foam::Zero.
Referenced by jouleHeatingSource::addSup(), hydrostaticPressure::calculateAndWrite(), kineticTheoryModel::correct(), IATE::correct(), radiativeIntensityRay::correct(), dynamicLagrangian< BasicTurbulenceModel >::correct(), thixotropicViscosity::correct(), advectionDiffusion::correct(), incompressiblePrimalSolver::correctBoundaryConditions(), waxSolventEvaporation::correctModel(), thermo::evolveRegion(), age::execute(), electricPotential::execute(), energyTransport::execute(), scalarTransport::execute(), adjointSimple::mainIter(), simple::mainIter(), adjointEikonalSolver::solve(), populationBalanceModel::solve(), twoPhaseSystem::solve(), thermalShell::solveEnergy(), reactingOneDim::solveEnergy(), thermalBaffle::solveEnergy(), and kinematicSingleLayer::solveMomentum().
void relax |
Relax matrix (for steady-state solution).
alpha is read from controlDict
Definition at line 1240 of file fvMatrix.C.
References Foam::name(), and relax().
void boundaryManipulate | ( | typename GeometricField< Type, fvPatchField, volMesh >::Boundary & | values | ) |
Manipulate based on a boundary field.
Definition at line 1256 of file fvMatrix.C.
References forAll.
Referenced by adjointSimple::mainIter(), adjointMeshMovementSolver::solve(), and multiphaseSystem::solveAlphas().
autoPtr< fvSolver > solver | ( | const dictionary & | ) |
Construct and return the solver.
Use the given solver controls
Referenced by radiativeIntensityRay::correct().
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.
Foam::SolverPerformance< Type > solveSegregatedOrCoupled | ( | const dictionary & | solverControls | ) |
Solve segregated or coupled returning the solution statistics.
Use the given solver controls
Definition at line 61 of file fvMatrixSolve.C.
References addProfiling, Foam::endl(), Foam::exit(), Foam::FatalIOError, FatalIOErrorInFunction, dictionary::getOrDefault(), Foam::Info, messageStream::masterStream(), mesh, regionName, and solve().
Referenced by displacementComponentLaplacianFvMotionSolver::solve(), velocityComponentLaplacianFvMotionSolver::solve(), displacementLaplacianFvMotionSolver::solve(), displacementSBRStressFvMotionSolver::solve(), solidBodyDisplacementLaplacianFvMotionSolver::solve(), surfaceAlignedSBRStressFvMotionSolver::solve(), velocityLaplacianFvMotionSolver::solve(), and fvMesh::solve().
Foam::SolverPerformance< Type > solveSegregated | ( | const dictionary & | solverControls | ) |
Solve segregated returning the solution statistics.
Use the given solver controls
Definition at line 111 of file fvMatrixSolve.C.
References Field< Type >::component(), Foam::diag(), Foam::endl(), Foam::exit(), Foam::FatalError, FatalErrorInFunction, dictionary::getOrDefault(), Foam::Info, messageStream::masterStream(), mesh, pTraits< bool >::nComponents, SolverPerformance< Type >::print(), psi, refPtr< T >::ref(), GeometricField< Type, PatchField, GeoMesh >::replace(), SolverPerformance< Type >::solverName(), and ensightPTraits< Type >::typeName.
Foam::SolverPerformance< Type > solveCoupled | ( | const dictionary & | solverControls | ) |
Solve coupled returning the solution statistics.
Use the given solver controls
Definition at line 251 of file fvMatrixSolve.C.
References GeometricField< Type, PatchField, GeoMesh >::component(), Foam::diag(), Foam::endl(), dictionary::getOrDefault(), Foam::Info, messageStream::masterStream(), mesh, SolverPerformance< Type >::print(), and psi.
Foam::SolverPerformance< Type > solve | ( | const dictionary & | solverControls | ) |
Solve returning the solution statistics.
Use the given solver controls
Definition at line 318 of file fvMatrixSolve.C.
Referenced by jouleHeatingSource::addSup(), Implicit< CloudType >::cacheFields(), hydrostaticPressure::calculateAndWrite(), kineticTheoryModel::correct(), IATE::correct(), dynamicLagrangian< BasicTurbulenceModel >::correct(), thixotropicViscosity::correct(), advectionDiffusion::correct(), waxSolventEvaporation::correctModel(), Foam::CorrectPhi(), thermo::evolveRegion(), age::execute(), electricPotential::execute(), energyTransport::execute(), scalarTransport::execute(), adjointSimple::mainIter(), simple::mainIter(), elasticityMotionSolver::solve(), laplacianMotionSolver::solve(), adjointEikonalSolver::solve(), populationBalanceModel::solve(), twoPhaseSystem::solve(), multiphaseSystem::solveAlphas(), reactingOneDim::solveContinuity(), thermalShell::solveEnergy(), reactingOneDim::solveEnergy(), thermalBaffle::solveEnergy(), sensitivityBezierFI::solveMeshMovementEqn(), reactingOneDim::solveSpeciesMass(), kinematicSingleLayer::solveThickness(), MultiComponentPhaseModel< BasePhaseModel, phaseThermo >::solveYi(), and Foam::fvc::spreadSource().
Foam::SolverPerformance< Type > solve |
Solve returning the solution statistics.
Solver controls read from fvSolution
Definition at line 343 of file fvMatrixSolve.C.
References fvMatrix< Type >::solve(), and fvMatrix< Type >::solverDict().
Referenced by fvMatrix< Type >::solve().
Foam::tmp< Foam::Field< Type > > residual |
Return the matrix residual.
Definition at line 350 of file fvMatrixSolve.C.
References fvMatrix< Type >::addBoundaryDiag(), fvMatrix< Type >::addBoundarySource(), GeometricField< Type, PatchField, GeoMesh >::boundaryField(), Field< Type >::component(), GeometricField< Type, PatchField, GeoMesh >::primitiveField(), tmp< T >::ref(), Field< Type >::replace(), faMatrix< Type >::residual(), GeometricBoundaryField< Type, PatchField, GeoMesh >::scalarInterfaces(), UList< T >::size(), and Foam::Zero.
Return the matrix scalar diagonal.
Definition at line 1270 of file fvMatrix.C.
References Foam::diag(), and tmp< T >::ref().
Foam::tmp< Foam::Field< Type > > DD |
Return the matrix Type diagonal.
Definition at line 1279 of file fvMatrix.C.
References fvPatchField< Type >::coupled(), Foam::diag(), forAll, tmp< T >::ref(), and UList< T >::size().
Return the central coefficient.
Definition at line 1303 of file fvMatrix.C.
References D, Foam::dimVol, IOobject::NO_READ, IOobject::NO_WRITE, and tmp< T >::ref().
Referenced by directionalPressureGradientExplicitSource::constrain(), meanVelocityForce::constrain(), incompressiblePrimalSolver::correctBoundaryConditions(), adjointSimple::mainIter(), simple::mainIter(), and kinematicSingleLayer::solveThickness().
Foam::tmp< Foam::GeometricField< Type, Foam::fvPatchField, Foam::volMesh > > H |
Return the H operation source.
Definition at line 1332 of file fvMatrix.C.
References GeometricField< Type, PatchField, GeoMesh >::correctBoundaryConditions(), DimensionedField< Type, GeoMesh >::dimensions(), Foam::dimVol, faMatrix< Type >::H(), Field< Type >::negate(), IOobject::NO_READ, IOobject::NO_WRITE, GeometricField< Type, PatchField, GeoMesh >::primitiveFieldRef(), tmp< T >::ref(), GeometricField< Type, PatchField, GeoMesh >::replace(), Field< Type >::replace(), and Foam::Zero.
Referenced by incompressiblePrimalSolver::correctBoundaryConditions(), adjointSimple::mainIter(), simple::mainIter(), and kinematicSingleLayer::solveThickness().
Return H(1)
Definition at line 1394 of file fvMatrix.C.
References Foam::component(), GeometricField< Type, PatchField, GeoMesh >::correctBoundaryConditions(), fvPatchField< Type >::coupled(), Foam::dimVol, forAll, lduMatrix::H1(), IOobject::NO_READ, IOobject::NO_WRITE, GeometricField< Type, PatchField, GeoMesh >::primitiveFieldRef(), tmp< T >::ref(), and UList< T >::size().
Referenced by adjointSimple::mainIter(), and simple::mainIter().
Return the face-flux field from the matrix.
Definition at line 1442 of file fvMatrix.C.
References Foam::abort(), GeometricField< Type, PatchField, GeoMesh >::boundaryFieldRef(), Foam::cmptMultiply(), lduMatrix::faceH(), Foam::FatalError, FatalErrorInFunction, forAll, IOobject::NO_READ, IOobject::NO_WRITE, GeometricField< Type, PatchField, GeoMesh >::primitiveFieldRef(), tmp< T >::ref(), Field< Type >::replace(), and DimensionedField< Type, GeoMesh >::setOriented().
Referenced by Implicit< CloudType >::cacheFields(), incompressiblePrimalSolver::correctBoundaryConditions(), Foam::CorrectPhi(), scalarTransport::execute(), adjointSimple::mainIter(), simple::mainIter(), twoPhaseSystem::solve(), multiphaseSystem::solveAlphas(), kinematicSingleLayer::solveThickness(), and MultiComponentPhaseModel< BasePhaseModel, phaseThermo >::solveYi().
const Foam::dictionary & solverDict |
Return the solver dictionary taking into account finalIteration.
Definition at line 1558 of file fvMatrix.C.
Referenced by fvMatrix< Type >::solve(), displacementComponentLaplacianFvMotionSolver::solve(), velocityComponentLaplacianFvMotionSolver::solve(), displacementLaplacianFvMotionSolver::solve(), displacementSBRStressFvMotionSolver::solve(), solidBodyDisplacementLaplacianFvMotionSolver::solve(), surfaceAlignedSBRStressFvMotionSolver::solve(), and velocityLaplacianFvMotionSolver::solve().
Definition at line 1574 of file fvMatrix.C.
References Foam::abort(), Foam::FatalError, FatalErrorInFunction, and lduMatrix::operator=().
Definition at line 1611 of file fvMatrix.C.
void negate |
Inplace negate.
Definition at line 1619 of file fvMatrix.C.
References lduMatrix::negate().
Definition at line 1634 of file fvMatrix.C.
References Foam::checkMethod(), and lduMatrix::operator+=().
Definition at line 1664 of file fvMatrix.C.
Definition at line 1672 of file fvMatrix.C.
References Foam::checkMethod(), and lduMatrix::operator-=().
Definition at line 1700 of file fvMatrix.C.
void operator+= | ( | const DimensionedField< Type, volMesh > & | su | ) |
Definition at line 1708 of file fvMatrix.C.
References Foam::checkMethod().
void operator+= | ( | const tmp< DimensionedField< Type, volMesh > > & | tsu | ) |
Definition at line 1719 of file fvMatrix.C.
void operator+= | ( | const tmp< GeometricField< Type, fvPatchField, volMesh > > & | tsu | ) |
Definition at line 1730 of file fvMatrix.C.
void operator-= | ( | const DimensionedField< Type, volMesh > & | su | ) |
Definition at line 1741 of file fvMatrix.C.
References Foam::checkMethod().
void operator-= | ( | const tmp< DimensionedField< Type, volMesh > > & | tsu | ) |
Definition at line 1752 of file fvMatrix.C.
void operator-= | ( | const tmp< GeometricField< Type, fvPatchField, volMesh > > & | tsu | ) |
Definition at line 1763 of file fvMatrix.C.
void operator+= | ( | const dimensioned< Type > & | su | ) |
Definition at line 1774 of file fvMatrix.C.
References DimensionedField< Type, GeoMesh >::mesh(), and psi.
void operator-= | ( | const dimensioned< Type > & | su | ) |
Definition at line 1784 of file fvMatrix.C.
References DimensionedField< Type, GeoMesh >::mesh(), and psi.
void operator+= | ( | const Foam::zero | ) |
Definition at line 1794 of file fvMatrix.C.
void operator-= | ( | const Foam::zero | ) |
Definition at line 1799 of file fvMatrix.C.
void operator*= | ( | const volScalarField::Internal & | dsf | ) |
Definition at line 1804 of file fvMatrix.C.
References Foam::abort(), fvMatrix< Type >::dimensions(), Foam::FatalError, FatalErrorInFunction, forAll, and lduMatrix::operator*=().
void operator*= | ( | const tmp< volScalarField::Internal > & | tfld | ) |
Definition at line 1834 of file fvMatrix.C.
void operator*= | ( | const tmp< volScalarField > & | tfld | ) |
Definition at line 1845 of file fvMatrix.C.
void operator*= | ( | const dimensioned< scalar > & | ds | ) |
Definition at line 1856 of file fvMatrix.C.
References fvMatrix< Type >::dimensions(), and lduMatrix::operator*=().
void setComponentReference | ( | const label | patchi, |
const label | facei, | ||
const | direction, | ||
const scalar | value | ||
) |
Definition at line 40 of file fvScalarMatrix.C.
References Foam::diag().
Foam::autoPtr< Foam::fvMatrix< Foam::scalar >::fvSolver > solver | ( | const dictionary & | solverControls | ) |
Definition at line 65 of file fvScalarMatrix.C.
References addProfiling, Foam::diag(), Foam::endl(), Foam::Info, messageStream::masterStream(), mesh, regionName, and solve().
Foam::solverPerformance solveSegregated | ( | const dictionary & | solverControls | ) |
Definition at line 165 of file fvScalarMatrix.C.
References fvMatrix< Type >::addBoundaryDiag(), fvMatrix< Type >::addBoundarySource(), lduMatrix::asymmetric(), GeometricField< Type, PatchField, GeoMesh >::boundaryField(), lduPrimitiveMeshAssembly::cellOffsets(), GeometricField< Type, PatchField, GeoMesh >::correctBoundaryConditions(), fvMatrix< Type >::createOrUpdateLduPrimitiveAssembly(), lduMatrix::diag(), Foam::endl(), schemesLookup::fluxRequired(), forAll, dictionary::getOrDefault(), Foam::Info, lduMatrix::lduAddr(), fvMatrix< Type >::lduMeshPtr(), lduMatrix::lower(), fvMatrix< Type >::manipulateMatrix(), messageStream::masterStream(), DimensionedField< Type, GeoMesh >::mesh(), lduMatrix::mesh(), IOobject::name(), Foam::name(), Time::New(), fvMatrix< Type >::nMatrices(), GeometricField< Type, PatchField, GeoMesh >::primitiveField(), primitiveFieldRef(), SolverPerformance< Type >::print(), psi, tmp< T >::ref(), GeometricBoundaryField< Type, PatchField, GeoMesh >::scalarInterfaces(), fvMatrix< Type >::setBounAndInterCoeffs(), fvMatrix< Type >::setInterfaces(), lduMatrix::setLduMesh(), List< T >::setSize(), UList< T >::size(), fvMatrix< Type >::transferFvMatrixCoeffs(), lduMatrix::upper(), and Foam::Zero.
Foam::tmp< Foam::scalarField > residual | ( | ) | const |
Definition at line 328 of file fvScalarMatrix.C.
References fvMatrix< Type >::addBoundaryDiag(), fvMatrix< Type >::addBoundarySource(), GeometricField< Type, PatchField, GeoMesh >::boundaryField(), GeometricField< Type, PatchField, GeoMesh >::primitiveField(), psi, refPtr< T >::ref(), faMatrix< Type >::residual(), GeometricBoundaryField< Type, PatchField, GeoMesh >::scalarInterfaces(), UList< T >::size(), and Foam::Zero.
Foam::tmp< Foam::volScalarField > H | ( | ) | const |
Definition at line 357 of file fvScalarMatrix.C.
References fvMatrix< Type >::addBoundarySource(), GeometricField< Type, PatchField, GeoMesh >::correctBoundaryConditions(), Foam::dimVol, faMatrix< Type >::H(), IOobject::instance(), DimensionedField< Type, GeoMesh >::mesh(), IOobject::name(), IOobject::NO_READ, IOobject::NO_WRITE, GeometricField< Type, PatchField, GeoMesh >::primitiveField(), GeometricField< Type, PatchField, GeoMesh >::primitiveFieldRef(), and tmp< T >::ref().
Foam::tmp< Foam::volScalarField > H1 | ( | ) | const |
Definition at line 389 of file fvScalarMatrix.C.
References GeometricField< Type, PatchField, GeoMesh >::correctBoundaryConditions(), DimensionedField< Type, GeoMesh >::dimensions(), Foam::dimVol, lduMatrix::H1(), IOobject::instance(), DimensionedField< Type, GeoMesh >::mesh(), IOobject::NO_READ, IOobject::NO_WRITE, GeometricField< Type, PatchField, GeoMesh >::primitiveFieldRef(), and tmp< T >::ref().
void setComponentReference | ( | const label | patchi, |
const label | facei, | ||
const | direction, | ||
const scalar | value | ||
) |
autoPtr< fvMatrix< scalar >::fvSolver > solver | ( | const dictionary & | ) |
solverPerformance solveSegregated | ( | const dictionary & | ) |
tmp< scalarField > residual | ( | ) | const |
tmp< volScalarField > H | ( | ) | const |
tmp< volScalarField > H1 | ( | ) | const |
|
friend |
Declare friendship with the fvSolver class.
Definition at line 181 of file fvMatrix.H.
|
friend |
|
friend |
|
friend |
|
friend |