[Changeset] Add the rcond function

View: New views
2 Messages — Rating Filter:   Alert me  

[Changeset] Add the rcond function

by David Bateman :: Rate this Message:

Reply to Author | View Threaded | Show Only this Message

Although this changeset is large, most of its just replacing "rcond"
with "rcon" to avoid shadow warnings for the rcond methods that the
changeset introduces in Matrix and ComplexMatrix classes. The oct-file
itself is then trivial.

Note that although I started doing the same methods for sparse matrices,
matlab doesn't support this and so I didn't do it in the end. Use
condest(a) if the matrix "a" is sparse.. Also it appears that matlab
only uses the LU factorization for their rcond function and not the
hermitian or triangular matrix versions of the LAPACK factorizations.
Therefore the Octave code will be faster for these cases, but won't
necessarily return the same result as matlab. This is a 1-norm estimate
of the reciprocal condition number in any case and so I believe the
difference can be ignore as its not the exact value we are interested by
rather its magnitude.

D.

--
David Bateman                                David.Bateman@...
Motorola Labs - Paris                        +33 1 69 35 48 04 (Ph)
Parc Les Algorithmes, Commune de St Aubin    +33 6 72 01 06 33 (Mob)
91193 Gif-Sur-Yvette FRANCE                  +33 1 69 35 77 01 (Fax)

The information contained in this communication has been classified as:

[x] General Business Information
[ ] Motorola Internal Use Only
[ ] Motorola Confidential Proprietary



patch8090.bz2 (11K) Download Attachment

[Changeset] Add the rcond function

by John W. Eaton :: Rate this Message:

Reply to Author | View Threaded | Show Only this Message

On 14-May-2008, David Bateman wrote:

| Although this changeset is large, most of its just replacing "rcond"
| with "rcon" to avoid shadow warnings for the rcond methods that the
| changeset introduces in Matrix and ComplexMatrix classes. The oct-file
| itself is then trivial.

I applied this changeset.

| Note that although I started doing the same methods for sparse matrices,
| matlab doesn't support this and so I didn't do it in the end. Use
| condest(a) if the matrix "a" is sparse.. Also it appears that matlab
| only uses the LU factorization for their rcond function and not the
| hermitian or triangular matrix versions of the LAPACK factorizations.
| Therefore the Octave code will be faster for these cases, but won't
| necessarily return the same result as matlab. This is a 1-norm estimate
| of the reciprocal condition number in any case and so I believe the
| difference can be ignore as its not the exact value we are interested by
| rather its magnitude.

OK.

Thanks,

jwe