inv.F90

Go to the documentation of this file.

subroutine invmat(nm,mat)
  implicit none 
  integer,intent(in)::nm
  real(8),intent(inout)::mat(nm,nm)
  integer::ipiv(nm)
  integer::Lwork 
  real(8),allocatable::work(:)
  integer::info 

  Lwork = 10*nm
  allocate(work(Lwork))
  info = 0
  call dgetrf(nm,nm,mat,nm,ipiv,info)
  call dgetri(nm,mat,nm,ipiv,work,Lwork,info)
  if(info /= 0) then
   write(6,*)'info (subrouitine inv):',info
   stop
  end if 
  deallocate(work)
  return 
end subroutine