From eb6f17e78397ed2dba8e5f79262c2514a925f1db Mon Sep 17 00:00:00 2001
From: Marco Ratto <marco.ratto@ec.europa.eu>
Date: Wed, 19 Oct 2022 18:07:06 +0200
Subject: [PATCH] check if initial states are empty, and set them to zero.

---
 matlab/+occbin/mkdata.m                  | 2 +-
 matlab/+occbin/mkdatap_anticipated_dyn.m | 2 +-
 2 files changed, 2 insertions(+), 2 deletions(-)

diff --git a/matlab/+occbin/mkdata.m b/matlab/+occbin/mkdata.m
index 303ff1ef6c..40af37863d 100644
--- a/matlab/+occbin/mkdata.m
+++ b/matlab/+occbin/mkdata.m
@@ -30,7 +30,7 @@ function [data_mat]=mkdata(n_periods,dr_A,dr_B,endo_names,exo_names,wish_list,sh
 % given decision rule
 neqs = size(endo_names,1);
 
-if nargin<9
+if nargin<9 || isempty(var_init)
     var_init = zeros(neqs,1);
 end
 
diff --git a/matlab/+occbin/mkdatap_anticipated_dyn.m b/matlab/+occbin/mkdatap_anticipated_dyn.m
index b361c62f41..f680851925 100644
--- a/matlab/+occbin/mkdatap_anticipated_dyn.m
+++ b/matlab/+occbin/mkdatap_anticipated_dyn.m
@@ -48,7 +48,7 @@ T = DM.decrulea;
 CONST = zeros(n_vars,1);
 R = DM.decruleb;
 
-if nargin<7
+if nargin<7 || isempty(init)
     init=zeros(n_vars,1);
 end
 
-- 
GitLab