From a07b0ae6141835145d68daf881fd43d65fab592c Mon Sep 17 00:00:00 2001
From: Jacob Smith <jacob.smith@atl.frb.org>
Date: Thu, 16 Jun 2011 10:28:42 -0400
Subject: [PATCH] Bug Fix, when using parameter uncertainty on multiple regimes
---
mex/sources/ms-sbvar/mex_ms_forecast.cc | 1 +
mex/sources/ms-sbvar/mex_ms_irf.cc | 1 +
mex/sources/ms-sbvar/mex_ms_variance_decomposition.cc | 1 +
3 files changed, 3 insertions(+)
diff --git a/mex/sources/ms-sbvar/mex_ms_forecast.cc b/mex/sources/ms-sbvar/mex_ms_forecast.cc
index 678be6b958..08d5eabbe9 100644
--- a/mex/sources/ms-sbvar/mex_ms_forecast.cc
+++ b/mex/sources/ms-sbvar/mex_ms_forecast.cc
@@ -115,6 +115,7 @@ mexFunction(int nlhs, mxArray *plhs[],
for (s = 0; s < nstates; s++)
{
printf("Constructing percentiles for forecast - regime %d\n", s);
+ if (options->simulation_file) rewind(options->simulation_file);
if (histogram = forecast_percentile_regime(options->shocks, options->simulation_file, options->thin, s, T, options->horizon, model, type))
{
for (k = 0; k < options->num_percentiles; k++)
diff --git a/mex/sources/ms-sbvar/mex_ms_irf.cc b/mex/sources/ms-sbvar/mex_ms_irf.cc
index ae607dd570..fe01315628 100644
--- a/mex/sources/ms-sbvar/mex_ms_irf.cc
+++ b/mex/sources/ms-sbvar/mex_ms_irf.cc
@@ -114,6 +114,7 @@ mexFunction(int nlhs, mxArray *plhs[],
for (s = 0; s < nstates; s++)
{
printf("Constructing percentiles for impulse responses - regime %d\n", s);
+ if (options->simulation_file) rewind(options->simulation_file);
if (histogram = impulse_response_percentile_regime(options->simulation_file, options->thin, s, options->horizon, model, type))
{
for (k = 0; k < options->num_percentiles; k++)
diff --git a/mex/sources/ms-sbvar/mex_ms_variance_decomposition.cc b/mex/sources/ms-sbvar/mex_ms_variance_decomposition.cc
index fdf0df0382..9d8e813a61 100644
--- a/mex/sources/ms-sbvar/mex_ms_variance_decomposition.cc
+++ b/mex/sources/ms-sbvar/mex_ms_variance_decomposition.cc
@@ -130,6 +130,7 @@ mexFunction(int nlhs, mxArray *plhs[],
}
else
{
+ if (options->simulation_file) rewind(options->simulation_file);
if (histogram = variance_decomposition_percentiles_regime(options->simulation_file, options->thin, s, options->horizon, model, type))
{
tvd = CreateMatrix(options->horizon, nvars*nvars);
--
GitLab