diff --git a/src/OpenFOAM/matrices/lduMatrix/lduAddressing/lduInterfaceFields/cyclicLduInterfaceField/cyclicLduInterfaceField.C b/src/OpenFOAM/matrices/lduMatrix/lduAddressing/lduInterfaceFields/cyclicLduInterfaceField/cyclicLduInterfaceField.C index e5da5b9e4dfc2db8e34dcc7440b6c78dfa1784cf..7c9d704616c3bd7b27179f97b3e8b5bb8f2e73d2 100644 --- a/src/OpenFOAM/matrices/lduMatrix/lduAddressing/lduInterfaceFields/cyclicLduInterfaceField/cyclicLduInterfaceField.C +++ b/src/OpenFOAM/matrices/lduMatrix/lduAddressing/lduInterfaceFields/cyclicLduInterfaceField/cyclicLduInterfaceField.C @@ -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()); + } } }