Matrix inverse

Hi,
I need to know if there exist a method for the Matrix and IntervalMatrix to compute the inverse.
Thanks a lot !
Etienne

Dear Etienne

Yes. Given a Matrix A of size n*n:

Matrix invA(n,n);

real_inverse(A,invA)

will store in invA the inverse of A.

For IntervalMatrix, you can obtain an enclosure of the inverse using backward multiplication. Given an IntervalMatrix A of size n*n, try:

IntervalMatrix invA(n,n);

bwd_mul(Matrix::eye(n), invA, A);

Gilles