[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/
- svds.m error in maka check in recent 3.2.x (??????),
Tatsuro MATSUOKA <=