49 boundaryCoeffs_[patchI][edgeI] = value;
60 <<
"solving faMatrix<scalar>" 74 addBoundaryDiag(
diag(), 0);
77 addBoundarySource(totalSource, 0);
86 psi_.boundaryField().scalarInterfaces(),
97 psi.correctBoundaryConditions();
99 psi.mesh().setSolverPerformance(
psi.name(), solverPerf);
109 addBoundaryDiag(boundaryDiag, 0);
120 source_ - boundaryDiag*psif,
122 psi_.boundaryField().scalarInterfaces(),
128 addBoundarySource(tres_s.
ref());
139 "H(" + psi_.name() +
')',
144 auto& Hphi = tHphi.ref();
146 Hphi.primitiveFieldRef() = (
lduMatrix::H(psi_.primitiveField()) + source_);
147 addBoundarySource(Hphi.primitiveFieldRef());
149 Hphi.ref() /= psi_.mesh().S();
150 Hphi.correctBoundaryConditions();
A list of keyword definitions, which are a keyword followed by a number of values (eg...
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.
Ostream & endl(Ostream &os)
Add newline and flush stream.
tmp< DimensionedField< TypeR, GeoMesh > > New(const tmp< DimensionedField< TypeR, GeoMesh >> &tf1, const word &name, const dimensionSet &dimensions, const bool initCopy=false)
Global function forwards to reuseTmpDimensionedField::New.
labelList faceLabels(nFaceLabels)
volScalarField H(IOobject("H", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), mesh, dimensionedScalar(dimLength, Zero))
void diag(pointPatchField< vector > &, const pointPatchField< tensor > &)
SolverPerformance< Type > solve()
Solve returning the solution statistics.
#define DebugInFunction
Report an information message using Foam::Info.
A const Field/List wrapper with possible data conversion.
tmp< GeometricField< Type, faPatchField, areaMesh > > H() const
Return the H operation source.
int debug
Static debugging option.
Container< Type > & ref() const
Return non-const reference to the contents of a non-null managed pointer.
tmp< Field< Type > > residual() const
Return the matrix residual.
messageStream Info
Information stream (stdout output on master, null elsewhere)
const volScalarField & psi
T getOrDefault(const word &keyword, const T &deflt, enum keyType::option matchOpt=keyType::REGEX) const
Find and return a T, or return the given default value. FatalIOError if it is found and the number of...
A class for managing temporary objects.
const word extrapolatedCalculatedType
A combined zero-gradient and calculated patch field type.
const dimensionSet dimArea(sqr(dimLength))
static constexpr const zero Zero
Global zero (0)