From c8568a7fb8ae62f21e02dc43d65adf923d652efe Mon Sep 17 00:00:00 2001
From: Marco Ratto <marco.ratto@jrc.ec.europa.eu>
Date: Fri, 12 Feb 2010 19:30:27 +0100
Subject: [PATCH] fixed bug: transpose after squeeze
---
matlab/DiffuseKalmanSmoother1.m | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/matlab/DiffuseKalmanSmoother1.m b/matlab/DiffuseKalmanSmoother1.m
index ef64185df0..def1b080cf 100644
--- a/matlab/DiffuseKalmanSmoother1.m
+++ b/matlab/DiffuseKalmanSmoother1.m
@@ -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;
--
GitLab