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;