Question about thread save subroutines for OpenMP

[guangye_li]

Hi, I am a customer at Shell. I am writing a code using Fotran with OpenMP and ifort. To have high level parallelizm, each thread needs to call a subroutine to work on its own private column. After the the subroutine is done, I made sure that all threads return correct results. But after "!$omp end parallel", the values become all zeros. Does the Intel ifort compiler support thread-save? How should I do it?
Thanks.
Guangye Li

[roddur]

ifort do support threading and openmp.
what result do you get when you run with OMP_NUM_THREADS=1?

[TimP]

Quoting - guangye.li

after "!$omp end parallel", the values become all zeros. Does the Intel ifort compiler support thread-save? How should I do it?

Do you mean lastprivate?

[jimdempseyatthecove]

Quoting - guangye.li

Hi, I am a customer at Shell. I am writing a code using Fotran with OpenMP and ifort. To have high level parallelizm, each thread needs to call a subroutine to work on its own private column. After the the subroutine is done, I made sure that all threads return correct results. But after "!$omp end parallel", the values become all zeros. Does the Intel ifort compiler support thread-save? How should I do it?
Thanks.
Guangye Li

Please supply a sampe of your code (or mock-up of code).

This sounds like a case where each thread works independently on a column of a application-wide array .AND. you are declaring that column private in the !$omp paralllel statement when it (the entire array)needs to be declared shared (but where each thread works on seperate columns).

Jim Dempsey