关于用mkl求特征值特征向量

380623711

使用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

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


不知道是什么原因

380623711

出错为floating divide by zero

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