ENH: face source: updated sumDirection op for vector

This commit is contained in:
andy
2013-01-17 11:22:58 +00:00
parent 5e13cdf1ff
commit d4a937aafc

View File

@ -540,23 +540,7 @@ Foam::vector Foam::fieldValues::faceSource::processValues
case opSumDirection:
{
const vector direction(dict_.lookup("direction"));
vector v(vector::zero);
forAll(Sf, i)
{
scalar d = Sf[i] & direction;
if (d > 0)
{
v += pos(values[i] & direction)*values[i];
}
else
{
v += neg(values[i] & direction)*values[i];
}
}
return v;
return sum(pos(values & direction)*values);
}
case opAreaNormalAverage:
{