openmp parallel for loop with two or more reductions

pyCthon picture pyCthon · Feb 22, 2012 · Viewed 34.5k times · Source

Hi just wondering if this is the right way to go going about having a regular for loop but with two reductions , is this the right approach below? Would this work with more then two reductions as well. Is there a better way to do this? also is there any chance to integrate this with an MPI_ALLREDUCE command?

heres the psuedo code

      #pragma omp parallel for \
      default(shared) private(i) \
      //todo first  reduction(+:sum)
      //todo second reduction(+:result)

      for loop i < n; i ++; {
        y = fun(x,z,i)
        sum += fun2(y,x)
        result += fun3(y,z)
      }

Answer

devil picture devil · Feb 22, 2012

You can do reduction by specifying more than one variable separated by a comma, i.e. a list:

#pragma omp parallel for default(shared) reduction(+:sum,result) ...

Private thread variables will be created for sum and result that will be combined using + and assigned to the original global variables at the end of the thread block.

Also, variable y should be marked private.

See https://computing.llnl.gov/tutorials/openMP/#REDUCTION