stepan
/
dynare
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

The block using bespoke eigs() is commented out because eigs() intermitently returns different rts 

git-svn-id: https://www.dynare.org/svn/dynare/trunk@2527 ac1d8469-bf42-47a9-8791-bf33cf982152
time-shift
george 2009-03-26 14:51:47 +00:00
parent db899b7f6b
commit bab59b7972
1 changed files with 9 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

17
matlab/AIM/SPEigensystem.m
View File

@ -30,13 +30,14 @@
function [w,rts,lgroots,flag_trouble] = SPEigensystem(a,uprbnd,rowsLeft) function [w,rts,lgroots,flag_trouble] = SPEigensystem(a,uprbnd,rowsLeft)
opts.disp=0; opts.disp=0;
try % next block is commented out because eigs() intermitently returns different rts
[w,d] = eigs(a',rowsLeft,'LM',opts); %try
rts = diag(d); % [w,d] = eigs(a',rowsLeft,'LM',opts);
mag = abs(rts); % rts = diag(d);
[mag,k] = sort(-mag); % mag = abs(rts);
rts = rts(k); % [mag,k] = sort(-mag);
catch % rts = rts(k);
%catch
%disp('Catch in SPE'); %disp('Catch in SPE');
%pause(0.5); %pause(0.5);
%aStr=datestr(clock); %aStr=datestr(clock);
@ -53,7 +54,7 @@ catch
mag = abs(rts); mag = abs(rts);
[mag,k] = sort(-mag); [mag,k] = sort(-mag);
rts = rts(k); rts = rts(k);
end %end
flag_trouble=0; flag_trouble=0;
%ws=SPSparse(w); %ws=SPSparse(w);