从 MKL 11.3 库调用 DPOSV 例程期间出现分段错误

Question

我正在尝试运行以下代码，使用 IFORT 2016 编译，链接 MKL 11.3 库：

program bug
implicit none
INCLUDE 'mkl.fi'
integer*4, parameter :: Npart=25
real(kind=8) :: u(1:Npart), lapackab(1:Npart,1:Npart)
integer*4 :: i
u=0.2d0
lapackab=0d0
do i=1,Npart
lapackab(i,i)=2d0
enddo
call DPOSV('L',Npart,1,lapackab,Npart,u,Npart,i)
write(*,*) "i=",i
end program

使用以下命令：

ifort -O0 -g -openmp -o file.o -c file.f90 -I/opt/share/INTEL/mkl/include
ifort -O0 -g -openmp file.o -o run -L/opt/share/INTEL/mkl/lib/  -I/opt/share/INTEL/mkl/include  -lmkl_intel_ilp64 -lmkl_intel_thread -lmkl_core -liomp5 -lpthread -lm -ldl

由于某种晦涩的原因，代码在 DPOSV 期间出现段错误。 我不知道为什么。说真的，即使是瓦尔格林德也没有报告任何奇怪的事情。 有没有人对这个特定的例程有同样的问题？

编辑：该代码与LAPACK 3.6.1完美配合，但不适用于MKL 11.3。我在英特尔支持论坛上发了一篇文章，但仍然需要得到版主的批准。

Answer 1

您正在链接 MKL 的ILP64版本，但该版本用于 64 位整数。您使用的是integer*4因此您的整数是 32 位的(以及默认整数(。使用 MKL 链接顾问正确链接 MKL。您可能需要LP64版本。