Fortran Coder

标题: 关于用mkl求特征值特征向量 [打印本页]

作者: 380623711 时间: 2024-11-26 19:16
标题: 关于用mkl求特征值特征向量
使用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

但是在另一个程序中加入这部分就出错

[attach]4077[/attach]
不知道是什么原因

作者: 380623711 时间: 2024-11-26 19:20
本帖最后由 380623711 于 2024-11-26 19:41 编辑

出错为floating divide by zero

我是在写一个FEM程序需要求解3*3对称矩阵的特征值特征向量，有没有大神能告诉一个最优的方法。

欢迎光临 Fortran Coder (http://bbs.fcode.cn/)