Question

openmp+fortran,为什么串行和并行结果不同，怎么修改？其中c1是已经存在的矩阵，c2初始化为0

!$omp parallel do reduction(+:c2)
do i=1,5
do j=1,6
c=c1(i,j)
c2(i,j)=c2(i,j)+c
enddo
enddo
!$omp end parallel do

Answer 1

没看懂你要实现什么功能？只是想把c1加到c2上的话，去掉reduction，私有化j就行，如下：
!$omp parallel do private(j)
do i=1,5
do j=1,6
c2(i,j)=c2(i,j)+c1(i,j)
enddo
enddo
!$omp end parallel do