Back to index

inverse

function inverse(mat) result(imat)
  implicit none
  real(kind=8) :: mat(:,:)
  real(kind=8) :: imat(size(m,1),size(m,1))

This function returns the inverse of a squared matrix mat implementing a LU decomposition.

Internal dependencies: error