From 81ba952414f658fb2ece4e99d45d4f4100ab6e4f Mon Sep 17 00:00:00 2001
From: Marco Ratto <marco.ratto@ec.europa.eu>
Date: Fri, 18 Nov 2022 10:52:37 +0100
Subject: [PATCH] big fix in diffuse filter initialization: incorrect mapping
of observables when removing zero columns of T
---
matlab/compute_Pinf_Pstar.m | 1 +
1 file changed, 1 insertion(+)
diff --git a/matlab/compute_Pinf_Pstar.m b/matlab/compute_Pinf_Pstar.m
index aa56b39e69..150acac187 100644
--- a/matlab/compute_Pinf_Pstar.m
+++ b/matlab/compute_Pinf_Pstar.m
@@ -169,6 +169,7 @@ if np0
QTinf(1:np0,np0+1:np0+nk) = STinf0;
QTinf([indx0(:); indx(:)],:) = QTinf;
STinf1 = [zeros(np0+np,np0) [STinf0; eye(nk); zeros(np-nk,nk)] zeros(np0+np,np-nk)];
+ mf = ismember([indx0(:); indx(:)],mf);
for k = 1:nk
if norm(QTinf(mf,:)*ST([indx0(:); indx(:)],k+np0)) < 1e-8
Pinf(k+np0,k+np0) = 0;
--
GitLab