BUG: cyclicLduInterfaceField : nonuniform transformation

This commit is contained in:
mattijs
2010-12-03 12:00:29 +00:00
parent 4e6f27634f
commit dac9c7c46d

View File

@ -41,16 +41,20 @@ Foam::cyclicLduInterfaceField::~cyclicLduInterfaceField()
void Foam::cyclicLduInterfaceField::transformCoupleField
(
scalarField& pnf,
scalarField& f,
const direction cmpt
) const
{
if (doTransform())
{
scalar forwardScale =
pow(diag(forwardT()[0]).component(cmpt), rank());
pnf *= forwardScale;
if (forwardT().size() == 1)
{
f *= pow(diag(forwardT()[0]).component(cmpt), rank());
}
else
{
f *= pow(diag(forwardT())().component(cmpt), rank());
}
}
}