From b8194f78452bb1f5a0f43b76d358d25fd3560ef2 Mon Sep 17 00:00:00 2001 From: Michel Juillard Date: Thu, 23 Jul 2015 14:25:18 +0200 Subject: [PATCH] histval: replaced initialization with NaN by zeros --- matlab/perfect-foresight-models/make_y_.m | 2 +- preprocessor/NumericalInitialization.cc | 4 ++-- 2 files changed, 3 insertions(+), 3 deletions(-) diff --git a/matlab/perfect-foresight-models/make_y_.m b/matlab/perfect-foresight-models/make_y_.m index 7611c9e97..1085b9030 100644 --- a/matlab/perfect-foresight-models/make_y_.m +++ b/matlab/perfect-foresight-models/make_y_.m @@ -53,6 +53,6 @@ else end % the first NaNs take care of the case where there are lags > 1 on % exogenous variables - oo_.endo_simul = [NaN(M_.endo_nbr,M_.maximum_lag-1) M_.endo_histval ... + oo_.endo_simul = [M_.endo_histval ... oo_.steady_state*ones(1,options_.periods+M_.maximum_lead)]; end diff --git a/preprocessor/NumericalInitialization.cc b/preprocessor/NumericalInitialization.cc index 85055905f..f8491d8b0 100644 --- a/preprocessor/NumericalInitialization.cc +++ b/preprocessor/NumericalInitialization.cc @@ -274,7 +274,7 @@ HistValStatement::writeOutput(ostream &output, const string &basename, bool mini output << "%" << endl << "% HISTVAL instructions" << endl << "%" << endl - << "M_.endo_histval = zeros(M_.endo_nbr,M_.maximum_endo_lag);" << endl; + << "M_.endo_histval = zeros(M_.endo_nbr,M_.maximum_lag);" << endl; for (hist_values_t::const_iterator it = hist_values.begin(); it != hist_values.end(); it++) @@ -310,7 +310,7 @@ HistValStatement::writeOutput(ostream &output, const string &basename, bool mini int tsid = symbol_table.getTypeSpecificID(symb_id) + 1; if (type == eEndogenous) - output << "M_.endo_histval( " << tsid << ", M_.maximum_endo_lag + " << lag << ") = "; + output << "M_.endo_histval( " << tsid << ", M_.maximum_lag + " << lag << ") = "; else if (type == eExogenous) output << "oo_.exo_simul( M_.maximum_lag + " << lag << ", " << tsid << " ) = "; else if (type != eExogenousDet)