octave-maintainers
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

svds.m error in maka check in recent 3.2.x (??????)


From: Tatsuro MATSUOKA
Subject: svds.m error in maka check in recent 3.2.x (??????)
Date: Fri, 27 Nov 2009 05:41:12 +0900 (JST)

Hello

I have build octave from 3.2.x source updated yesterday night.

New error occurred in 'make check' results

>>>>> processing 
>>>>> c:\usr\Tatsu\mingwhome\octaves\hg\release-3-2-x-work\scripts/sparse\svds.m
  ***** testif HAVE_ARPACK
 [u2,s2,v2,flag] = svds(a,k,0);
 s2 = diag(s2);
 assert(flag,!1);
 assert(s(k:-1:1), s2, 1e-10); 
!!!!! test failed
assert (flag,!1) expected
0
but got
 1
values do not match>>>>> processing
c:\usr\Tatsu\mingwhome\octaves\hg\release-3-2-x-work\scripts/specfun\beta.m

make check on 3.2.3 does not show errors

However on octave 3.2.x, 
./run-octave

octave:1> n = 100;
octave:2> k = 7;
octave:3> 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)]);
octave:4> %%a = 
sparse([3:n,1:n,1:(n-2)],[1:(n-2),1:n,3:n],[ones(1,n-2),1:n,-ones(1,n-2)]);
octave:4> [u,s,v] = svd(full(a));
octave:5> s = diag(s);
octave:6> [dum, idx] = sort(abs(s));
octave:7> s = s(idx);
octave:8> u = u(:,idx);
octave:9> v = v(:,idx);
octave:10> [u2,s2,v2,flag] = svds(a,k,0);
octave:11> s2 = diag(s2);
octave:12> assert(flag,!1);
octave:13> assert(s(k:-1:1), s2, 1e-10);
octave:14>

No error happened!

??????

It is too mysterious to understand.

I will check at PC in university.

Regards

Tatsuro

--------------------------------------
What is your No.1 Entertainment of 2009? -Yahoo! JAPAN Net BANZUKE 2009
http://pr.mail.yahoo.co.jp/banzuke/


reply via email to

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