diff options
| author |  Benoit Jacob <jacob.benoit.1@gmail.com> | 2008-08-09 15:50:07 +0000 |
|---|---|---|
| committer | Benoit Jacob <jacob.benoit.1@gmail.com> | 2008-08-09 15:50:07 +0000 |
| commit | becbeda50ac17288dba0a93c6adc67b663d32a7a (patch) | |
| tree | 0bc52ff80dff8efd259dc245edc3d7a5f9cd5d4e /test/submatrices.cpp | |
| parent | 681e9446705449a17d7c5ba669b91c068f9b98a0 (diff) | |
* reimplement the general case of inverse() on top of LU. Advantages:
- removes much code
- 2.5x faster (even though LU uses complete pivoting contrary to what
inverse used to do!)
- there _were_ numeric stability problems with the partial pivoting
approach of inverse(), with 200x200 matrices inversion failed almost
half of the times (overflow). Now these problems are solved thanks to
complete pivoting.
* remove some useless stuff in LU
Diffstat (limited to 'test/submatrices.cpp')
0 files changed, 0 insertions, 0 deletions