From 281e3979d5a68ba0420dce49aa9572b52d892abd Mon Sep 17 00:00:00 2001 From: heckflosse Date: Tue, 23 Aug 2016 21:32:30 +0200 Subject: [PATCH] ImProcFunctions::moyeqt, simplified calculation --- rtengine/improcfun.cc | 11 ++++++----- 1 file changed, 6 insertions(+), 5 deletions(-) diff --git a/rtengine/improcfun.cc b/rtengine/improcfun.cc index f70341334..34e2eeb34 100644 --- a/rtengine/improcfun.cc +++ b/rtengine/improcfun.cc @@ -2867,7 +2867,8 @@ void ImProcFunctions::moyeqt (Imagefloat* working, float &moyS, float &eqty) int tWw = working->width; double moy = 0.0; double sqrs = 0.0; -#ifndef _DEBUG + +#ifdef _OPENMP #pragma omp parallel for reduction(+:moy,sqrs) schedule(dynamic,16) #endif @@ -2879,10 +2880,10 @@ void ImProcFunctions::moyeqt (Imagefloat* working, float &moyS, float &eqty) } } - double mo = moy / (tHh * tWw); - moyS = mo; - double eqt = (sqrs - 2.0 * mo * moy + tHh * tWw * SQR(mo)) / (tHh * tWw); - eqty = sqrt(eqt); + moy /= (tHh * tWw); + sqrs /= (tHh * tWw); + eqty = sqrt(sqrs - SQR(moy)); + moyS = moy; } static inline void