Fortran Coder

查看: 79|回复: 4

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

[复制链接]

6

帖子

2

主题

0

精华

入门

F 币
29 元
贡献
19 点
发表于 2020-2-14 10:59:57 | 显示全部楼层 |阅读模式
5F 币
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的函数,求指点。

回复

使用道具 举报

444

帖子

1

主题

0

精华

大宗师

F 币
2297 元
贡献
1479 点
发表于 2020-2-14 21:22:38 | 显示全部楼层
use lapack95,并需要连接 mkl_lapack95.lib 。
1.png
回复

使用道具 举报

6

帖子

2

主题

0

精华

入门

F 币
29 元
贡献
19 点
 楼主| 发表于 2020-2-15 12:26:55 | 显示全部楼层
li913 发表于 2020-2-14 21:22
use lapack95,并需要连接 mkl_lapack95.lib 。

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

使用道具 举报

444

帖子

1

主题

0

精华

大宗师

F 币
2297 元
贡献
1479 点
发表于 2020-2-16 14:36:45 | 显示全部楼层
linux的话,文件名可能不同。比如 ifort -o lapack.o test.o -O test.out
回复

使用道具 举报

6

帖子

2

主题

0

精华

入门

F 币
29 元
贡献
19 点
 楼主| 发表于 昨天 10:01 | 显示全部楼层
li913 发表于 2020-2-16 14:36
linux的话,文件名可能不同。比如 ifort -o lapack.o test.o -O test.out

感谢,问题我已经解决了
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

QQ|捐赠本站|Archiver|关于我们 About Us|群聊|Fcode

GMT+8, 2020-2-28 03:20

Powered by Discuz! X3.2

© 2001-2017 Comsenz Inc.

快速回复 返回顶部 返回列表