0

这是本书中实现的矩阵到向量乘法的并行(逐字)版本

使用 OpenMP

subroutine matvecmul(mat,vec,res,m,n)
use omp_lib
implicit none
integer::m,n,i,j
real*8,dimension(m,n),intent(in)::mat
real*8,dimension(n),intent(in)::vec
real*8,dimension(m)::res 
res(:)=0.0
 !$omp parallel do default(none)
 !$omp shared(mat,vec,res,m,n) private(i,j)
  do i=1,m
   do j=1,n
    res(i)=res(i)+mat(i,j)*vec(j)
   end do
  end do
 !$omp end parallel do
 return
end subroutine matvecmul

我得到这个语法错误“未分类的 openMP 子句共享”,“未分类的 openMP 子句私有”用 gfortran 4.4.5 和 gfortran 4.6.0 编译

共享子句是否在“gomp”中被删除/弃用或未实现,或者我犯了一些愚蠢的错误或书的勘误表......使用默认值(一)我没有得到任何语法错误我该怎么说哪些变量是私有的和要分享哪个?

这是我如何实现没有错误(并且正确)

subroutine matvecmul(mat,vec,res,m,n)
use omp_lib
implicit none
integer::m,n,i,j
real*8,dimension(m,n),intent(in)::mat
real*8,dimension(n),intent(in)::vec
real*8,dimension(m)::res
res(:)=0.0
 !$omp parallel do private(i,j)
  do i=1,m
   do j=1,n
    res(i)=res(i)+mat(i,j)*vec(j)
   end do
  end do
 !$omp end parallel do
return
end subroutine matvecmul

我仍然是openmp的新手(但我需要快速学习它......我的工作的一部分!)这里的一些解释将不胜感激......是否有任何只使用gomp的openMP教程(比如gcc/gfortran 4.4.5) ?

4

1 回答 1

1

这与 gcc 的 OpenMP 实现无关;intel fortran 也拒绝编译它。这只是关于如何在 OpenMP + Fortran 中继续行。

你想要做的,这是完全正确的,是:

!$omp parallel do default(none) shared(mat,vec,res,m,n) private(i,j)

我衷心赞同使用默认无并明确指定共享和私有变量作为最佳实践的方法。如果您使用上面的行,您的例程将编译

为清楚起见,代码将其分成多行;这并没有什么问题,除了语法是,与其他现代 Fortran 续行一样,在续行的末尾需要一个 & 符号:

 !$omp parallel do default(none)  &
 !$omp shared(mat,vec,res,m,n) private(i,j)

这大概是书中的一个错字。因为 OpenMP 是一种标准,所以 OpenMP 没有特殊的 gomp(或 iomp)版本。兼容的 OpenMP 程序应该在任何兼容的编译器上编译。

于 2011-08-05T16:01:09.803 回答