有没有大神用过intel的mkl库求解大型矩阵的特征值
intel的mkl库求解大型矩阵(20000*20000)的特征值。系统是ubuntu,用到的函数是mkl库的syevd。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的函数,求指点。
use lapack95,并需要连接 mkl_lapack95.lib 。
li913 发表于 2020-2-14 21:22
use lapack95,并需要连接 mkl_lapack95.lib 。
请问linux编译时怎么链接mkl_lapack95.lib? linux的话,文件名可能不同。比如 ifort -o lapack.o test.o -O test.out li913 发表于 2020-2-16 14:36
linux的话,文件名可能不同。比如 ifort -o lapack.o test.o -O test.out
感谢,问题我已经解决了 li913 发表于 2020-2-14 21:22
use lapack95,并需要连接 mkl_lapack95.lib 。
大神,请问如果我只需要计算出最大特征值,不需要计算特征向量和其他特征值,应该是用什么mkl的函数? 本帖最后由 风平老涡 于 2020-8-29 22:57 编辑
ssyevd是F77版, syevd是F95版。可能是系统设置问题。 本帖最后由 风平老涡 于 2020-8-29 22:58 编辑
mmSCI 发表于 2020-8-29 15:58
大神,请问如果我只需要计算出最大特征值,不需要计算特征向量和其他特征值,应该是用什么mkl的函数? ...
call syevd(c,w,'N')
or
call syevd(c,w)
风平老涡 发表于 2020-8-29 22:57
call syevd(c,w,'N')
or
call syevd(c,w)
请问你有没有对比过,部分特征值计算的函数和这个syevd全部特征值计算函数,那个效率更高?我试了试,20000*20000的特征值计算,虽然syevd是计算全部特征值,但是还是它比较快一些,其他计算部分特征值的函数,时间会稍长。
页:
[1]