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

Rationalised loops 

git-svn-id: https://www.dynare.org/svn/dynare/trunk@2644 ac1d8469-bf42-47a9-8791-bf33cf982152
time-shift
george 2009-04-28 17:03:54 +00:00
parent 920e21561b
commit afb3cba759
1 changed files with 6 additions and 14 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

20
matlab/kordpert/dr1.m
View File

@ -151,27 +151,19 @@ function [dr,info,M_,options_,oo_] = dr1(dr,task,M_,options_,oo_)
ghs2=zeros(g_2rows, 1);
for i=1:g_2cols
if s0 < nspred & s1 < nspred
for j=1:g_2rows
ghxx(j,s0*nspred+s1+1) = 2*g_2(j,i);
ghxx(:,s0*nspred+s1+1) = 2*g_2(:,i);
if s1 > s0
ghxx(j,s1*nspred+s0+1) = 2*g_2(j,i);
ghxx(:,s1*nspred+s0+1) = 2*g_2(:,i);
end
end
elseif s0 < nspred & s1 < nspred+nExog
for j=1:g_2rows
ghxu(j,(s0*nExog+s1-nspred+1)) = 2*g_2(j,i);
end
ghxu(:,(s0*nExog+s1-nspred+1)) = 2*g_2(:,i);
elseif s0 < nspred+nExog & s1 < nspred+nExog
for j=1:g_2rows
ghuu(j,(s0-nspred)*nExog+s1-nspred +1) = 2*g_2(j,i);
ghuu(:,(s0-nspred)*nExog+s1-nspred +1) = 2*g_2(:,i);
if s1 > s0
ghuu(j,(s1-nspred)*nExog+s0-nspred+1) = 2*g_2(j,i);
ghuu(:,(s1-nspred)*nExog+s0-nspred+1) = 2*g_2(:,i);
end
end
else
for j=1:g_2rows
ghs2(j,1) = 2*g_2(j,end);
end
ghs2(:,1) = 2*g_2(:,end);
end
s1 = s1+1;
if s1 == nspred+nExog