OpenMP为啥这段程序是死锁？

Posted 2023-02-17

在学OpenMP，书中举了一个死锁的例子，但是不明白为什么会死锁：
...
#pragma omp critical
y=f(x);
...
double f(double x)
#pragma omp critical
z=g(x);
...
这段代码中，不是其中一个线程顺序执行下来就可以吗？为什么还会死锁呢？

参考技术A 才发现你这源程序混合了C、F77、F95三种语法…… 而且没有恰当使用private保护变量或者reduction。 正确的写法（F95+OpenMP）： program sct use omp_lib implicit none integer k,dct dct=0 !$OMP parallel do reduction(+:dct) do k=1,1000 dct=dct+1 enddo !$OMP end parallel do write(*,*)"dct=", dct pause end program sct