矩阵求逆问题

ffp · 发表于 2016-4-13 08:34:26

大家好，问下，我调用lapack里面的求逆函数，我的原矩阵中有A(i,i)=0的项，说LU分解和求逆有问题，对于这种中间有A(i,i)=0的这种应该怎么调用哪个求逆的函数命令，我用下面的这个是不对的

[Fortran] syntaxhighlighter_viewsource syntaxhighlighter_copycode

     using LAPACK
                call dgetrf(n,n,bitmp,ldb,ipiv,info)
                if(info.ne.0) then
                WRITE(6,*) 'problem in LU decomposition'
                WRITE(6,*) 'info=',info
                endif
                call dgetri(n,bitmp,ldb,ipiv,work,n,info)
                if(info.ne.0) then
                WRITE(6,*) 'problem in matrix inversion'
                WRITE(6,*) 'info=',info
                endif