如何在 Fortran 中正确调用 SGEMV？答案

【问题标题】：How to properly call the SGEMV in Fortran?如何在 Fortran 中正确调用 SGEMV？
【发布时间】：2019-01-14 01:20:45
【问题描述】：

我想使用 BLAS 的 SGEMV 子例程在 fortran 中执行矩阵向量乘积。我有一个类似的代码：

program test
integer, parameter :: DP = selected_real_kind(15)
real(kind=DP), dimension (3,3) :: A
real(kind=DP), dimension (3) :: XP,YP
call sgemv(A,XP,YP)

A 是 3x3 矩阵，XP 和 YP 是向量。在包含的模块中可以看到以下代码：

PURE SUBROUTINE SGEMV_F95(A,X,Y,ALPHA,BETA,TRANS)
    ! Fortran77 call:
    ! SGEMV(TRANS,M,N,ALPHA,A,LDA,X,INCX,BETA,Y,INCY)
    USE F95_PRECISION, ONLY: WP => SP
    REAL(WP), INTENT(IN), OPTIONAL :: ALPHA
    REAL(WP), INTENT(IN), OPTIONAL :: BETA
    CHARACTER(LEN=1), INTENT(IN), OPTIONAL :: TRANS
    REAL(WP), INTENT(IN) :: A(:,:)
    REAL(WP), INTENT(IN) :: X(:)
    REAL(WP), INTENT(INOUT) :: Y(:)
END SUBROUTINE SGEMV_F95

我知道有些参数是可选的，那么我在方法调用中哪里错了？

【问题讨论】：

请包含您在尝试编译时收到的错误消息。
您显示了子程序sgemv_f95 的接口，但call sgemv。您也不使用任何模块。因此，正如之前的评论所暗示的，我们需要更多关于您的问题的详细信息。

标签： visual-studio fortran blas intel-fortran intel-mkl

【解决方案1】：

当您查看 BLAS 或 LAPACK 例程时，您应该始终查看第一个字母：

S：单精度
D：双精度
C: 单精度复数
Z: 双精度复数

您使用以下语句将矩阵 A 以及向量 XP 和 YP 定义为双精度数：

integer, parameter :: DP = selected_real_kind(15)

为此，您需要使用dgemv 或将您的精度定义为单精度。

调用dgemv 和dgemv_f95 也有区别。 dgemv_f95 是英特尔 MKL 的一部分，并不是真正的通用命名。出于可移植性的原因，我不会使用该表示法，而是坚持使用经典的 dgemv，它也是英特尔 MKL 的一部分。

DGEMV 执行矩阵向量运算之一
y := alpha*A*x + beta*y,   or   y := alpha*A**T*x + beta*y,
其中alpha 和beta 是标量，x 和y 是向量，A 是 m by n 矩阵。

如果你想知道如何调用该函数，我建议看看here，但它最终应该是这样的：

call DGEMV('N',3,3,ALPHA,A,3,XP,1,BETA,YP,1)

【讨论】：

不仅解决了我的问题，还帮助我更好地理解了 BLAS 例程。非常感谢！

【解决方案2】：

精度不兼容。您正在调用 sgemv，它采用单精度参数，但您传递的是双精度数组和向量。

【讨论】：

使用模块有助于确保调用具有可比性。模块可用吗？

【解决方案3】：

也许trans 参数是必需的？

trans: Must be 'N', 'C', or 'T'.

（根据Developer Reference for Intel® Math Kernel Library - Fortran底部的注释。）

【讨论】：