[Octave-bug-tracker] [bug #29721] tests for function svds frequently fails

[Andreas Weber]

Follow-up Comment #5, bug #29721 (project octave):

Thank you for your reply.
I have removed:

libarpack++2-dev libarpack++2c2a libarpack2 libarpack2-dev liblapack-dev
liblapack3gf


re-ran configure && make but the test still fails. I have cropped the test to

n = 100;
k = 7;
A =
sparse([3:n,1:n,1:(n-2)],[1:(n-2),1:n,3:n],[ones(1,n-2),0.4*n*ones(1,n),ones(1,n-2)]);
[u,s,v] = svd(full(A));
s = diag(s);
[~, idx] = sort(abs(s));
s = s(idx);
randn('state',42);      % Initialize to make normest function reproducible
rand('state',42)
opts.v0 = rand (2*n,1);
[u2,s2,v2,flag] = svds(A,k,0,opts);
s2 = diag(s2);
assert(s2, s(k:-1:1), 1e-10);
error: assert (s2,s (k:-1:1),1e-10) expected
   38.060
   38.034
   38.034
   38.015
   38.015
   38.004
   38.004
but got
   38.034
   38.034
   38.015
   38.015
   38.004
   38.004
Dimensions don't match


If I change the rand state different to 42 (tried 41,43,44) the test passes:

octave:60> test svds
PASSES 5 out of 5 tests


please let me know if I can provide more informations. I also tried bisect but
found no version which passes the test beginning from


Änderung:        8417:654bcfb937bf
Nutzer:          David Bateman <address@hidden>
Datum:           Tue Dec 23 08:28:23 2008 +0100
Zusammenfassung: Add the eigs and svds functions
-verbatim 

Andy

    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?29721>

_______________________________________________
  Nachricht geschickt von/durch Savannah
  http://savannah.gnu.org/