[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: Octave QR factorization
From: |
taltman |
Subject: |
Re: Octave QR factorization |
Date: |
Wed, 17 Dec 2003 21:00:47 +0000 (UTC) |
On Dec 17, 2003 at 9:53pm, Bart Vandewoestyne wrote:
Bart.V >address@hidden wrote:
Bart.V >So this is what you mean i should use instead of Q*R==A:
Bart.V >
Bart.V >if ( norm(my_Q*my_R-A)/norm(A) < eps )
Bart.V > fprintf('Our calculations are correct!\n');
Bart.V >else
Bart.V > fprintf('Our calculations are wrong!\n');
Bart.V >end
Bart.V >
Bart.V >?
Well, in a particular case, this may work. Not always.
Usually, the '< eps' is written '< tol', where 'tol' is some function
of the dimensions of the matrices involved, their condition number,
their rank, and the eps ( not necessarily all of those, but some ). As
I mentioned before, there's no catch-all way. Every algorithm &
application will have its own desired precision, and is own precision
limitations due to round-off error.
Is there a numerical analyst in the house who can be more precise?
Thanks,
~Tomer
-------------------------------------------------------------
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
-------------------------------------------------------------
Re: Octave QR factorization, A S Hodel, 2003/12/17