OpenMP中运行后数值问题

978142355

练习OpenMP在Fortran中的应用，主要是代码中Single那一部分，传入的a值应该是0，为什么运行后a值变成1了？谢谢各位。

[Fortran] 纯文本查看 复制代码

查看代码

打印?

01	program test17

02	use omp_lib

03	implicit none

04	integer tid,cpus

05

integer a

06

07	call omp_set_num_threads(3)

08

09	!$omp parallel private(a,tid)

10	tid=omp_get_thread_num()

11

a=tid

12	write(*,*) 'a=',a,'tid=',tid,' before single'

13	!$omp single

14	write(*,*)

15	cpus=omp_get_num_threads()

16	tid=omp_get_thread_num()

17

18	a=a+2**tid

19

20	write(*,*) 'single threads:', cpus

21	write(*,*) 'a=',a,'tid=',tid,' during single'

22	write(*,*)

23	!$omp end single copyprivate(a)

24	tid=omp_get_thread_num()

25	write(*,*) 'a=',a,'tid=',tid,'after single'

26	!$omp end parallel

27

read(*,*)

28	end program test17

978142355

自评一下：脑子抽筋了，那里的符号是2**tid，当成了2*tid。当tid=0，2**tid结果是1，a传入的是0，所以最终值肯定是1。。。。。。。。。。。。。。。。。。