From 4f1121057cf15538ba5615b76c634d1b7a3fd348 Mon Sep 17 00:00:00 2001
From: Marco Ratto <marco.ratto@ec.europa.eu>
Date: Fri, 3 Feb 2023 14:24:58 +0100
Subject: [PATCH] store binding_indicator_history once it has been set
(cherry picked from commit cdd195576ef0019efe4a161ebbc540cc95058b0a)
---
matlab/+occbin/solve_one_constraint.m | 5 ++---
1 file changed, 2 insertions(+), 3 deletions(-)
diff --git a/matlab/+occbin/solve_one_constraint.m b/matlab/+occbin/solve_one_constraint.m
index a55ebd0e06..2f7834df9d 100644
--- a/matlab/+occbin/solve_one_constraint.m
+++ b/matlab/+occbin/solve_one_constraint.m
@@ -152,9 +152,7 @@ for shock_period = 1:n_shocks_periods
if length(binding_indicator)<(nperiods_0 + 1)
binding_indicator=[binding_indicator; false(nperiods_0 + 1-length(binding_indicator),1)];
end
-
- binding_indicator_history{iter}=binding_indicator;
-
+
if iter==1 && guess_history_it
regime = regime_history_guess(shock_period).regime;
regime_start = regime_history_guess(shock_period).regimestart;
@@ -164,6 +162,7 @@ for shock_period = 1:n_shocks_periods
end
nperiods_0 = size(binding_indicator,1)-1; %if history is present, update may be required
end
+ binding_indicator_history{iter}=binding_indicator;
% analyze when each regime starts based on current guess
[regime, regime_start, error_code_period]=occbin.map_regime(binding_indicator,opts_simul_.debug);
regime_history(shock_period).regime = regime;
--
GitLab