diff --git a/src/functionObjects/field/Curle/Curle.C b/src/functionObjects/field/Curle/Curle.C index 39f38def4e..1a8133a83d 100644 --- a/src/functionObjects/field/Curle/Curle.C +++ b/src/functionObjects/field/Curle/Curle.C @@ -213,11 +213,11 @@ bool Foam::functionObjects::Curle::execute() forAll(observerPositions_, pointi) { - const vectorField r(Cfp - observerPositions_[pointi]); + const vectorField r(observerPositions_[pointi] - Cfp); const scalarField invMagR(1/(mag(r) + ROOTVSMALL)); pDash[pointi] += - sum((pp*sqr(invMagR) + invMagR/c0_*dpdtp)*(Sfp & r)); + sum((pp*sqr(invMagR) + invMagR/c0_*dpdtp)*(Sfp & (r*invMagR))); } }