octave-maintainers
[Top][All Lists]
Advanced

[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/


reply via email to

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