Re: Willing to work for octave.

[msander]

Jordi Gutiérrez Hermoso-2 wrote
> Compile Octave and run the test suite with "make check". There are
> still some things that don't have tests. You can also look at the
> tests that are failing, maybe you can make some of them should pass!
> 
> - Jordi G. H.

 
i have built & run make check.   i looked into the very first failure.  
failing in svds.  output of "test svds verbose" below.  I have looked at the
code & tests.  it seems to be failing for technical reasons beyond my
expertise.  I'm not a mathematician.    Presumably this test passes, or used
to pass in certain environments.   Relaxing the pass criterion does not seem
like the correct approach.

Is there someone i could punt this to?    the file lists an author, but no
other contact info.   any ideas?




>> test svds verbose
>>>>> /home/mike/Documents/octave/.build/../scripts/sparse/svds.m
  ***** shared n, k, A, u, s, v, opts, rand_state, randn_state, tol
 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);
 u = u(:, idx);
 v = v(:, idx);
 randn_state = randn ("state");
 rand_state = rand ("state");
 randn ("state", 42);      % Initialize to make normest function
reproducible
 rand ("state", 42);
 opts.v0 = rand (2*n,1); % Initialize eigs ARPACK starting vector
                         % to guarantee reproducible results

  ***** testif HAVE_ARPACK
 [u2,s2,v2,flag] = svds (A,k);
 s2 = diag (s2);
 assert (flag, !1);
 tol = 10 * eps() * norm(s2, 1);
 assert (s2, s(end:-1:end-k+1), tol);

  ***** testif HAVE_ARPACK, HAVE_UMFPACK
 [u2,s2,v2,flag] = svds (A,k,0,opts);
 s2 = diag (s2);
 assert (flag, !1);
 tol = 10 * eps() * norm(s2, 1);
 assert (s2, s(k:-1:1), tol);

warning: returning fewer singular values than requested
warning: try increasing the value of sigma
!!!!! test failed
ASSERT errors for:  assert (s2,s (k:-1:1),tol)

  Location  |  Observed  |  Expected  |  Reason
     .          O(6x1)       E(7x1)      Dimensions don't match
>>



--
View this message in context: 
http://octave.1599824.n4.nabble.com/Willing-to-work-for-octave-tp4660415p4660946.html
Sent from the Octave - Maintainers mailing list archive at Nabble.com.