Re: Improved normest

[Jaroslav Hajek]

On Tue, Feb 2, 2010 at 4:28 PM, Marco Caliari <address@hidden> wrote:
> Dear Jaroslav,
>
> what about the following solution for the initial vector?
>
> 1) save the state of the random number generator with
>
> v = rand("state");
>
> 2) set the state to something fixed (eventually depending on the matrix,
> such as
>
> rand("state",trace(A))
>
> 3) "randomly" initialize the vector
>
> x = rand(size(A,1),1);
>
> 4) perform the powers method and finally, before endfunction, restore the
> state of the random number generator with
>
> rand("state",v).
>
> Cheers,
>
> Marco
>

I have no objections to this idea - if you wish, please submit a patch.

cheers

-- 
RNDr. Jaroslav Hajek, PhD
computing expert & GNU Octave developer
Aeronautical Research and Test Institute (VZLU)
Prague, Czech Republic
url: www.highegg.matfyz.cz