关于用mkl求特征值特征向量
使用dgeev函数求特征值特征向量，在独立程序中可以正常求解program test
use lapack95
implicit none
integerinfo, 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,:) =
arr(2,:) = [ 0.0 ,1.0 ,0.0]
arr(3,:) =
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
但是在另一个程序中加入这部分就出错
不知道是什么原因
本帖最后由 380623711 于 2024-11-26 19:41 编辑
出错为floating divide by zero
我是在写一个FEM程序需要求解3*3对称矩阵的特征值特征向量，有没有大神能告诉一个最优的方法。
页:
[1]