|
From: Michael A. K. <ma...@ll...> - 2002-09-26 13:25:32
|
Nicolas,
1. (M/ A) eventually calls (GESV A I), where I is the identity
matrix, which intern calls the FORTRAN ?GESV which performs the LU
decomposition on the system of equations [A I] to form the inverse
of A.
2. Your statement concerning inverses is correct...and it is correctly
implimented in MATLISP.
> The point is that computing the inverse is often not recommended
> because of numerical stability. Usually one uses the LU decomposition
> directly which is done by the getrf/getrs-pair. Of course, if you
> only need to apply the inverse once you could be satisfied with m/.
>
mike
**************************************************
Dr Michael A. Koerber Micro$oft Free Zone
MIT/Lincoln Laboratory
|