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

fixed bug: transpose after squeeze

time-shift
Marco Ratto 2010-02-12 19:30:27 +01:00
parent cfda6bf6bb
commit c8568a7fb8
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
matlab/DiffuseKalmanSmoother1.m
View File

@ -131,7 +131,7 @@ while notsteady & t<smpl
a(:,t+1) = T*atilde(:,t);
aK(1,:,t+1) = a(:,t+1);
for jnk=2:nk,
aK(jnk,:,t+jnk) = T*squeeze(aK(jnk-1,:,t+jnk-1));
aK(jnk,:,t+jnk) = T*squeeze(aK(jnk-1,:,t+jnk-1))';
end
P(:,:,t+1) = T*P(:,:,t)*transpose(T)-T*P(:,mf,t)*transpose(K(:,:,t)) + QQ;
notsteady = ~(max(max(abs(P(:,:,t+1)-P(:,:,t))))<crit);
@ -154,7 +154,7 @@ while t<smpl
a(:,t+1) = T*a(:,t) + K_s*v(:,t);
aK(1,:,t+1) = a(:,t+1);
for jnk=2:nk,
aK(jnk,:,t+jnk) = T*squeeze(aK(jnk-1,:,t+jnk-1));
aK(jnk,:,t+jnk) = T*squeeze(aK(jnk-1,:,t+jnk-1))';
end
end
t = smpl+1;