[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Improved normest
From: |
Marco Caliari |
Subject: |
Improved normest |
Date: |
Thu, 26 Nov 2009 12:19:05 +0100 (CET) |
User-agent: |
Alpine 1.00 (DEB 882 2007-12-20) |
Dear maintainers,
I have seen changeset
http://hg.savannah.gnu.org/hgweb/octave/rev/87595d714005 about normest, by
Jaroslav. As already discussed, the choice
x = norm (A, "columns").'
for the initial vector leads to wrong values (second largest eigenvalue)
for some "important" matrices, like
toeplitz([-2,1,zeros(1,2)]), toeplitz([-2,1,zeros(1,4)]), etc
This was the reason for the old random initial vector. I agree it is not
nice to have slightly different results at each run, but it is not nice to
have wrong values, too.
Best regards,
Marco
- Improved normest,
Marco Caliari <=