#include <stdio.h>
#define N 1000
int main(int argc, char **argv) {
  int b[N]={[0 ... N-1] = 1}, c[N]={[0 ... N-1] = 1};
  int i, sum=0, ps=0;

#pragma omp parallel for reduction(+:sum)
  for (i=0; i<N; i++)  
    sum = sum + b[i]*c[i];

  printf("Sum from reduction is %d\n",sum);

  sum = 0;
#pragma omp parallel firstprivate(ps)
  {
#pragma omp for
    for (i=0; i<N; i++)
      ps = ps + b[i]*c[i];
    printf("Partial sum is %d\n",ps);
#pragma omp critical
    sum = sum + ps;
  }

  printf("Sum from critical section is %d\n",sum);
}