diff --git a/src/postProcessing/functionObjects/field/fieldValues/cellSource/cellSourceTemplates.C b/src/postProcessing/functionObjects/field/fieldValues/cellSource/cellSourceTemplates.C index 372bd224b9..1acc0876b1 100644 --- a/src/postProcessing/functionObjects/field/fieldValues/cellSource/cellSourceTemplates.C +++ b/src/postProcessing/functionObjects/field/fieldValues/cellSource/cellSourceTemplates.C @@ -132,7 +132,7 @@ Type Foam::fieldValues::cellSource::processValues scalar mean = component(meanValue, d); scalar& res = setComponent(result, d); - res = sqrt(sum(V*sqr(vals - mean))/(V.size()*sum(V)))/mean; + res = sqrt(sum(V*sqr(vals - mean))/sum(V))/mean; } break; diff --git a/src/postProcessing/functionObjects/field/fieldValues/faceSource/faceSourceTemplates.C b/src/postProcessing/functionObjects/field/fieldValues/faceSource/faceSourceTemplates.C index c9e69771ff..584a0e03d6 100644 --- a/src/postProcessing/functionObjects/field/fieldValues/faceSource/faceSourceTemplates.C +++ b/src/postProcessing/functionObjects/field/fieldValues/faceSource/faceSourceTemplates.C @@ -236,9 +236,7 @@ Type Foam::fieldValues::faceSource::processSameTypeValues scalar mean = component(meanValue, d); scalar& res = setComponent(result, d); - res = - sqrt(sum(magSf*sqr(vals - mean))/(magSf.size()*sum(magSf))) - /mean; + res = sqrt(sum(magSf*sqr(vals - mean))/sum(magSf))/mean; } break;