@ -210,12 +210,14 @@ void Foam::cyclicAMIFvPatchField<Type>::updateInterfaceMatrix
|
||||
|
||||
solveScalarField pnf(psiInternal, nbrFaceCells);
|
||||
|
||||
const labelUList& faceCells = lduAddr.patchAddr(patchId);
|
||||
|
||||
// Transform according to the transformation tensors
|
||||
transformCoupleField(pnf, cmpt);
|
||||
|
||||
if (cyclicAMIPatch_.applyLowWeightCorrection())
|
||||
{
|
||||
solveScalarField pif(psiInternal, cyclicAMIPatch_.faceCells());
|
||||
solveScalarField pif(psiInternal, faceCells);
|
||||
pnf = cyclicAMIPatch_.interpolate(pnf, pif);
|
||||
}
|
||||
else
|
||||
@ -223,8 +225,6 @@ void Foam::cyclicAMIFvPatchField<Type>::updateInterfaceMatrix
|
||||
pnf = cyclicAMIPatch_.interpolate(pnf);
|
||||
}
|
||||
|
||||
const labelUList& faceCells = lduAddr.patchAddr(patchId);
|
||||
|
||||
// Multiply the field by coefficients and add into the result
|
||||
this->addToInternalField(result, !add, faceCells, coeffs, pnf);
|
||||
}
|
||||
|
||||
@ -313,10 +313,21 @@ Foam::label Foam::fvMatrix<Type>::checkImplicit(const label fieldI)
|
||||
{
|
||||
if (bpsi[patchI].useImplicit())
|
||||
{
|
||||
if (debug)
|
||||
{
|
||||
Pout<< "fvMatrix<Type>::checkImplicit "
|
||||
<< " fieldi:" << fieldI
|
||||
<< " field:" << this->psi(fieldI).name()
|
||||
<< " on mesh:"
|
||||
<< this->psi(fieldI).mesh().name()
|
||||
<< " patch:" << bpsi[patchI].patch().name()
|
||||
<< endl;
|
||||
}
|
||||
|
||||
id += (label(2) << patchI);
|
||||
}
|
||||
}
|
||||
if (id > 0)
|
||||
if (id >= 0)
|
||||
{
|
||||
lduAssemblyName_ = word("lduAssembly") + name(id);
|
||||
useImplicit_ = true;
|
||||
|
||||
Reference in New Issue
Block a user