有没有大神用过intel的mkl库求解大型矩阵的特征值

mmSCI

intel的mkl库求解大型矩阵（20000*20000）的特征值。系统是ubuntu，用到的函数是mkl库的syevd。

[Fortran] 纯文本查看 复制代码

program main
      implicit none
      integer,parameter :: M=20000,N=500,WORKL=M,iWORKL=M
      integer i,j,info
      real dataM(M,N),C(M,M),T(N,M),W(M),WORK(WORKL)
      integer iWORK(iWORKL)
      open(5,file='output_1.47_0.000_1.dat')
      do i=1,M
          read(5,*) (dataM(i,j),j=1,N)
      end do
      write(*,*) "load the dat file!"
      do i=1,M
          dataM(i,:)=dataM(i,:)/sqrt(sum(dataM(i,:)**2))
      end do
      write(*,*) "dataM has been unitary"
      T=transpose(dataM)
      C=matmul(dataM,T)
      write(*,*) "C has been done!"
      call syevd(C,W,'V')

编译ifort -mkl test.f -o test,文件名test.f.但是系统找不到这个函数，可是我用同样的操作ssyevd函数就可以编译，但是ssyecd函数输入的参数太多，执行起来也很困难。有没有大神用过这个mkl的函数，求指点。

li913

use lapack95,并需要连接 mkl_lapack95.lib 。

mmSCI

li913 发表于 2020-2-14 21:22
use lapack95,并需要连接 mkl_lapack95.lib 。

请问linux编译时怎么链接mkl_lapack95.lib?

li913

linux的话，文件名可能不同。比如 ifort -o lapack.o test.o -O test.out

mmSCI

li913 发表于 2020-2-16 14:36
linux的话，文件名可能不同。比如 ifort -o lapack.o test.o -O test.out

感谢，问题我已经解决了

mmSCI

li913 发表于 2020-2-14 21:22
use lapack95,并需要连接 mkl_lapack95.lib 。

大神，请问如果我只需要计算出最大特征值，不需要计算特征向量和其他特征值，应该是用什么mkl的函数？

风平老涡

ssyevd是F77版， syevd是F95版。可能是系统设置问题。

风平老涡

mmSCI 发表于 2020-8-29 15:58
大神，请问如果我只需要计算出最大特征值，不需要计算特征向量和其他特征值，应该是用什么mkl的函数？ ...

call syevd(c,w,'N')
or
call syevd(c,w)

mmSCI

风平老涡 发表于 2020-8-29 22:57
call syevd(c,w,'N')
or
call syevd(c,w)

请问你有没有对比过，部分特征值计算的函数和这个syevd全部特征值计算函数，那个效率更高？我试了试，20000*20000的特征值计算，虽然syevd是计算全部特征值，但是还是它比较快一些，其他计算部分特征值的函数，时间会稍长。