[liuhuafei]

Maybe some bug in leastSquaresVectors.c

Description:
some error in calculation of dd in leastsquare vector

Source file:
leastSquaresVectors.c

function
void Foam::leastSquaresVectors::makeLeastSquaresVectors () const

if (p.coupled())
{
forAll(pd, patchFacei)
{
const vector& d = pd[patchFacei];

dd[faceCells[patchFacei]] +=
(pw[patchFacei]*pMagSf[patchFacei]/magSqr(d))*sqr(d);
}
}

dd maybe done by the following:
dd[faceCells[patchFacei]] +=
(1.-pw[patchFacei])*pMagSf[patchFacei]/magSqr(d))*sqr(d);

also a question:
for the couple boundary,dd is calculated only for P cell(owner),why not for N cell(neighbour)?