Question

我正在尝试创建一个简单的示例，以便在MKL Intel库中使用Pardiso解算器。我一直在关注提供的示例，但是如果我将Pardiso的调用放在子例程中则不起作用。我担心这与INCLUDE语句或链接方面有关。

所有子程序使用的数组都包含在名为variables

的模块中

MODULE variables
INTEGER :: M        ! Lines
INTEGER :: N        ! Columns
REAL*8, DIMENSION(:,:), ALLOCATABLE :: MATRA ! original matrix
INTEGER, DIMENSION(:), ALLOCATABLE :: ROWSA,COLSA ! ia and ja in pardiso
REAL*8,    dimension(:), ALLOCATABLE :: VALSA, RHSVC, SOLVC ! a, b, x
END MODULE variables

包含该程序的文件如下：

INCLUDE 'mkl_pardiso.f90'
program PardisoFortran
use variables
use mkl_pardiso
implicit none

! do some stuff to create the matrices
call create_matrices

call INITPARDISO

end program

最后是initpardiso文件

subroutine INITPARDISO
USE VARIABLES
USE mkl_pardiso

! pardiso variable declaration    
TYPE(MKL_PARDISO_HANDLE), ALLOCATABLE  :: pt(:)
INTEGER maxfct, mnum, mtype, phase, nrhs, error, msglvl
INTEGER, ALLOCATABLE, DIMENSION(:) :: iparm
INTEGER i, idum
REAL*8 waltime1, waltime2, ddum

! --- then I allocate and fill the variables

! Finally I can call pardiso
phase = 11 ! only reordering and symbolic factorization phase

CALL pardiso_64 (PT, maxfct, mnum, mtype, phase, M, VALSA, COLSA, &
              ROWSA, idum, nrhs, iparm, msglvl, ddum, ddum, error)

end subroutine

现在，我还在Visual Studio项目配置中添加了附加包含目录（即

）

C:\Program Files (x86)\Intel\ComposerXE-2011\mkl\include;
C:\Program Files (x86)\Intel\ComposerXE-2011\mkl\lib\intel64;
C:\Program Files (x86)\Intel\ComposerXE-2011\mkl\lib\ia32)

如果我对pardiso的调用发表评论，一切都运行正常，否则会因此错误停止编译：

错误1错误＃6285：此通用子例程调用没有匹配的特定子例程。 [PARDISO_64] .... PardisoFortran \ initpardiso.f 144

有什么想法吗？这是INCLUDE声明的问题吗？我应该把它包括在哪里？

Answer 1

我相信pardiso_64是一个使用64位整数的子程序版本。您的代码段显示没有证据表明您已采取措施确保您的整数为64位。我想你可能已经使用编译器选项或其他方法告诉编译器默认为64位整数。

您发布的错误消息是在调用库通用过程时，在伪参数和实际参数之间存在不匹配时获得的错误消息的典型信息。假设您的整数只有32位，那么编译器正在寻找并且找不到一个名为pardiso_64的例程，该例程采用32位整数参数。

我认为该错误与include语句无关。

Answer 2

我在您的代码中看到一个问题。参数'ddum'是标量。它应该是一个维度数组（M，nrhs）。可能会有更多问题。我需要查看整个代码才能提供帮助。

与往常一样，询问英特尔MKL相关问题的最佳位置是MKL论坛：http://software.intel.com/en-us/forums/intel-math-kernel-library。在那里发布您的问题，附上您的整个程序。你会得到更快的答案！

英特尔MKL Pardiso解算器使用Visual Studio进行编译

2 个答案: