From b2e3a481e364c472f6daf19c2c6b5c102becc64f Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?S=C3=A9bastien=20Villemot?= Date: Thu, 20 Dec 2018 15:00:50 +0100 Subject: [PATCH] Fix for histval with auxiliary lead variables corresponding to non-linear terms The M_.params structure was not correctly passed to dynamic_set_auxiliary_series.m, leading to crashes if a parameter was needed to compute an auxiliary variable (typically an auxiliayr lead variable corresponding to a non-linear term in stochastic setup). --- src/NumericalInitialization.cc | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/NumericalInitialization.cc b/src/NumericalInitialization.cc index 4d1dc8b5..8665e48f 100644 --- a/src/NumericalInitialization.cc +++ b/src/NumericalInitialization.cc @@ -378,7 +378,7 @@ HistValStatement::writeOutput(ostream &output, const string &basename, bool mini } output << "if exist(['+' M_.fname '/dynamic_set_auxiliary_series'])" << endl - << " eval(['M_.histval_dseries = ' M_.fname '.dynamic_set_auxiliary_series(M_.histval_dseries, []);']);" << endl + << " eval(['M_.histval_dseries = ' M_.fname '.dynamic_set_auxiliary_series(M_.histval_dseries, M_.params);']);" << endl << "end" << endl << "M_.endo_histval = M_.histval_dseries{M_.endo_names{:}}(dates(sprintf('%dY', 1-M_.maximum_lag)):dates('0Y')).data';" << endl << "M_.endo_histval(isnan(M_.endo_histval)) = 0;" << endl; // Ensure that lead aux variables do not have a NaN