octave-maintainers
[Top][All Lists]
Advanced

[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


reply via email to

[Prev in Thread] Current Thread [Next in Thread]