Fortran Coder

查看: 105|回复: 1

[数学库] 关于用mkl求特征值特征向量

[复制链接]

2

帖子

1

主题

0

精华

新人

F 币
12 元
贡献
5 点
发表于 2024-11-26 19:16:57 | 显示全部楼层 |阅读模式
使用dgeev函数求特征值特征向量,在独立程序中可以正常求解
program test
    use lapack95
    implicit none
    integer  info, lda, ldvl, ldvr, lwork, rwork,n
    REAL*8 vl(3,3), vr(3,3),wr(3),wi(3),arr(3,3),work(20)
    arr(1,:) = [1.0, 1.0 ,0.0]
    arr(2,:) = [ 0.0 ,1.0 ,0.0]
    arr(3,:) = [0.0,0.0,1.0]
    n=3
    lda=3
    ldvl=3
    ldvr=3
    lwork=12
    rwork=12
    call dgeev('N','V',n,arr,lda,wr,wi,vl,ldvl,vr,ldvr,work,lwork,info)
    print*, 'eigenvalue',wr
    print*, '----------------------------------------------------'
    print*, 'eigenvector',vr
    print*, '----------------------------------------------------'
    pause
    end program test
Screenshot 2024_11_26 19_06_56.png
但是在另一个程序中加入这部分就出错
str_model (正在调试) - Microsoft Visual Studio 2024_11_26 19_03_19 (2).png

不知道是什么原因


2

帖子

1

主题

0

精华

新人

F 币
12 元
贡献
5 点
 楼主| 发表于 2024-11-26 19:20:56 | 显示全部楼层
本帖最后由 380623711 于 2024-11-26 19:41 编辑

出错为floating divide by zero Screenshot 2024_11_26 19_00_54.png

我是在写一个FEM程序需要求解3*3对称矩阵的特征值特征向量,有没有大神能告诉一个最优的方法。
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2024-12-11 12:05

Powered by Tencent X3.4

© 2013-2024 Tencent

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