Why doesn't Fortran have a built-in matrix inverse function?

Actually it’s mostly solve that is used a lot, not matrix inversion, since it is not as accurate and has other issues as @ivanpribec pointed out.

That being said, inv(x) is useful for just experimenting or testing or debugging, I’ve used it a few times. It should go into stdlib into its linear algebra module.

1 Like