Upozornění kompilátoru (úroveň 4) C4938
'var' : Proměnná redukce s plovoucí desetinou čárkou může způsobit nekonzistentní výsledky v /fp:strict nebo #pragma fenv_access
Neměli byste používat /fp:strict nebo fenv_access se snížením desetinné čárky OpenMP, protože součet se vypočítá v jiném pořadí. Výsledky se tedy mohou lišit od výsledků bez /openmp.
Následující ukázka vygeneruje C4938:
// C4938.cpp
// compile with: /openmp /W4 /fp:strict /c
// #pragma fenv_access(on)
extern double *a;
double test(int first, int last) {
double sum = 0.0;
#pragma omp parallel for reduction(+: sum) // C4938
for (int i = first ; i <= last ; ++i)
sum += a[i];
return sum;
}
Bez explicitní paralelizace se součet vypočítá takto:
sum = a[first] + a[first + 1] + ... + a[last];
Při explicitní paralelizaci (a dvou vláknech) se součet vypočítá následujícím způsobem:
sum1 = a[first] + ... a[first + last / 2];
sum2 = a[(first + last / 2) + 1] + ... a[last];
sum = sum1 + sum2;