OpenMP并行采用reduction，出现stack overflow问题
利用OpenMP进行并行运算2个矩阵P与Q的值，并利用reduction进行归约处理；先尝试了一个简单算例，计算没有问题；但是移至到另外一个大程序中时，出现了stack overflow问题。
以下是简单算例的程序部分：
! 主程序段
matrix = 0.d0
vector = 0.d0
!$omp parallel private(i, tmp_vector) &
!$omp reduction(+:matrix, vector)
!$omp do
do i =1, 48
call gemm(a, b, matrix, 'n', 'n', 1.d0, 1.d0)
call gemm(a, c, tmp_vector)
call gemm(b, tmp_vector, vector, 'n', 'n', 1.d0, 1.d0)
end do
!$omp end do
!$omp end parallel
在上述简单算例中运行正确，但移至到大型程序中出现以下问题：
OpenMP 可以一定程度解决效率问题。但并不能解决数组太大内存（堆栈）不足的问题。 好像并行编程有 OpenMP 和MPI两种，本人很想学习，但苦于资料比较少:-dizzy: 找到了
呵呵，http://fcode.cn/resource_ebook-list-1.html
慢慢啃，入手并行编程！
:-lol
页:
[1]