[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: svds.m error in maka check in recent 3.2.x (??????)
From: |
Tatsuro MATSUOKA |
Subject: |
Re: svds.m error in maka check in recent 3.2.x (??????) |
Date: |
Fri, 27 Nov 2009 16:56:48 +0900 (JST) |
Hello
No problem occurred in PC in the university.
Perhaps it relies on the ATLAS.
Regards
Tatsuro
--- Tatsuro MATSUOKA
> 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/
>
--------------------------------------
What is your No.1 Entertainment of 2009? -Yahoo! JAPAN Net BANZUKE 2009
http://pr.mail.yahoo.co.jp/banzuke/