|
|
|
|
|
by beagle3
5025 days ago
|
|
The reason you want to avoid matrix inversion is not computational performance - it is precision. You can generally "multiply by the inverse" without actually computing the inverse, in a way that needs less intermediate floating point precision. If your matrices have a large spread of eigenvalues, it makes a lot of difference - double precision often doesn't have enough precision for direct inversion in the real world. (And even if you do want to invert, it is more numerically stable to do that as a solution of Ax=e, for each 'e' a basis element, as long as you compute A^-1 \* e indirectly using a numerically stable method) |
|