suppress condition warning on inv?

[Paul Kienzle]

Hi!  I'm working with some potentially ill-conditioned matrices
Which when I invert them may give a warning message:

warning: inverse: matrix singular to machine precision, rcond = 1.34804e-30

Is there a way to return rcond, or suppress the warning, or know
quickly ahead of time if the warning will be issued?  

I already have the schur decomposition and the eigenvalues/eigenvectors
of the matrix if that helps.  

Calling cond/svd takes 2x longer than the inverse itself, and if it
passes, I will still want to do the inverse, so that's not an ideal solution.

Shall I patch inv.cc so that it accepts the form: [X, rcond] = inv(A); ?

        ...

        retval(0) = m.inverse (info, rcond, 1);

        if (nargout > 1)
          retval(1) = rcond;
        else if (info == -1)
          warning ("inverse: matrix singular ...", rcond);

        ...

Paul Kienzle
address@hidden



-------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.

Octave's home on the web:  http://www.octave.org
How to fund new projects:  http://www.octave.org/funding.html
Subscription information:  http://www.octave.org/archive.html
-------------------------------------------------------------