R: new test error for svds.m

[Marco Atzeri]

--- Gio 27/5/10, Ben Abbott ha scritto:

> Rik, 
> 
> With my last build, I encountered a test error for svds.m.
> My tip is
> 
> changeset:   10665:0f310fce905d
> tag:         tip
> user:        Jaroslav Hajek <address@hidden>
> date:        Thu May 27 08:54:44 2010
> +0200
> summary:     implement
> octave_base_matrix::matrix_ref const
> 
> The test error is …
> 
> >>>>> processing
> /Users/bpabbott/Development/mercurial/local_clone/scripts/sparse/svds.m
>   ***** testif HAVE_ARPACK
>  [u2,s2,v2,flag] = svds(a,k,0,opts);
>  s2 = diag(s2);
>  assert(flag,!1);
>  assert(s(k:-1:1), s2, 1e-10); 
> !!!!! test failed
> assert (s (k:-1:1),s2,1e-10) expected
>    38.060
>    38.060
>    38.034
>    38.034
>    38.015
>    38.015
>    38.004
> but got
>    38.060
>    38.034
>    38.034
>    38.015
>    38.015
>    38.004
>    38.004
> maximum absolute error 0.0263523 exceeds tolerance 1e-10
> 
> If I make the change below, the test passes.
> 
> %! [u2,s2,v2,flag] = svds(a,k,0);#,opts);
> 
> Are there more changes to come for svds.m which will
> resolve this?
> 
> Ben
> 

I guess we found that is an issue coming from eigs 
convergence and not really a fault of svds.
Repeating the test several times it fails around 30% of 
the time


Marco