1  int
       2  foo (int *a, int *b)
       3  {
       4    int r = 0;
       5    #pragma omp parallel for reduction (inscan, +:r) default(none) firstprivate (a, b)
       6    for (int i = 0; i < 64; i++)
       7      {
       8        r += a[i];
       9        #pragma omp scan inclusive (r)
      10        b[i] = r;
      11      }
      12    return r;
      13  }