urbach/tmLQCD — pull request #7 by kostrzewa (2012-10-29 11:37:08+00:00)
Title: do manual reduction in assign_mul_add_r_and_square

OpenMP reductions add the reduction variable in an undefined order,
resulting in differing rounding errors from call to call. To alleviate this
we use the omp accumulator arrays and collect the result manually.