Parallel For Loop with OpenMP

Question

I have problem to parallel for-loop code in OpenMP, result of parallel for-loop is different with a sequential for-loop. How to make this code parallel with same result as sequential code.

counter = 0;
#pragma omp parallel for
for(i=0; i<L; i++) {
    int sum_found = 0;
    for(j=0; j<M; j++) {
        int found = 0;
        for(k=0; k<N_SUBSET; k++) {
            if (i==0 && unsorted_arr[j*N_SUBSET + k] < intervals[i]) {
                s_prime[counter] = unsorted_arr[j*N_SUBSET + k];
                counter++;
                found++;
            }
            else if (i!=0 && unsorted_arr[j*N_SUBSET + k] >= intervals[i-1] && unsorted_arr[j*N_SUBSET + k] < intervals[i]) {
                s_prime[counter] = unsorted_arr[j*N_SUBSET + k];
                counter++;
                found++;
            }
            else if (i==L-1 && unsorted_arr[j*N_SUBSET + k] >= intervals[i-1]) {
                s_prime[counter] = unsorted_arr[j*N_SUBSET + k];
                counter++;
                found++;
            }
        }
        C[i][j] = found;
        sum_found += found;
    }
    n_intervals_len[i] = sum_found;
}

Try to add -fsanitize=thread. It will tell you if and where you have a data race — sv90
– sv90, Commented Apr 4, 2019 at 18:21

Michael Klemm · Accepted Answer · 2019-04-03 14:50:32Z

1

As some of the context of the above code is missing (esp. the variable declarations), it's rather hard to tell what is going wrong. But here are some ideas:

the loop counters j and k should be declared private, because otherwise you might have a race condition on them
the counter++ and found++ statements are also races for updates from different threads on the same variable. You would either have to use an atomic construct for them or use a proper OpenMP reduction to get rid of the race condition.

answered Apr 3, 2019 at 14:50

Michael Klemm

2,9531 gold badge16 silver badges16 bronze badges

Sign up to request clarification or add additional context in comments.

Collectives™ on Stack Overflow

Parallel For Loop with OpenMP

1 Answer 1

Comments

Your Answer

Hot Network Questions

Collectives™ on Stack Overflow

1 Answer 1

Comments

Your Answer

Sign up or log in

Post as a guest

Related