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 |