invZGE.F90

Go to the documentation of this file.
00001 subroutine invZGE(nm,mat)
00002   implicit none 
00003   integer,intent(in)::nm
00004   complex(8),intent(inout)::mat(nm,nm)
00005   integer::ipiv(nm)
00006   integer::Lwork 
00007   complex(8),ALLOCATABLE::work(:)
00008   integer::info 
00009   Lwork=10*nm
00010   allocate(work(Lwork))
00011   info=0
00012   call zgetrf(nm,nm,mat,nm,ipiv,info)
00013 !  write(6,*)'finish zgetrf'
00014   call zgetri(nm,mat,nm,ipiv,work,Lwork,info)
00015 !  write(6,*)'finish zgetri'
00016   if(info/=0) then
00017     write(6,*)'info (subrouitine inv):',info
00018     stop
00019   end if 
00020   deallocate(work)
00021   return 
00022 end subroutine
00023