Множественные сокращения OpenMP SIMD (сумма, минимум и максимум) в одиночном цикле For ⇐ C++

[Гость]

У меня есть следующий цикл для расчета базовой сводной статистики (среднее значение, стандартное отклонение, минимум и максимум) в C++, пропуская пропущенные значения (x — двойной вектор):

int k = 0; длинная двойная сумма = 0,0, sq_sum = 0,0; двойной минв = 1,0/0,0, максв = -1,0/0,0; #pragma omp simd сокращение(+:sum,sq_sum,k) сокращение(max:maxv) сокращение(min:minv) for (int i = 0; i < n; ++i) { двойной xi = x; int tmp = xi == xi; сумма += ТМП ? хи: 0,0; sq_sum += tmp ? кси * кси: 0,0; к += тмп ? 1:0; минв = минв > кси? xi: минв; maxv = maxv