From b3f36ab7f53ac0ca9904b2251e67a0d6616a41c2 Mon Sep 17 00:00:00 2001
From: Marco Ratto <marco.ratto@ec.europa.eu>
Date: Thu, 20 Oct 2022 15:14:09 +0200
Subject: [PATCH] impede prior_trunc=0 when doing sensitivity
(cherry picked from commit 1a46c1dc897076ac6f2dec170f847db59d7c772d)
---
matlab/dynare_sensitivity.m | 6 ++++++
1 file changed, 6 insertions(+)
diff --git a/matlab/dynare_sensitivity.m b/matlab/dynare_sensitivity.m
index 9bdb5b61d4..398beaa615 100644
--- a/matlab/dynare_sensitivity.m
+++ b/matlab/dynare_sensitivity.m
@@ -89,6 +89,8 @@ if options_.order~=1
options_.order = 1;
end
+original_prior_trunc = options_.prior_trunc;
+
if ~isempty(options_gsa.datafile) || isempty(bayestopt_) || options_gsa.rmse
if isempty(options_gsa.datafile) && options_gsa.rmse
disp('The data file and all relevant estimation options ')
@@ -134,6 +136,9 @@ else
if isempty(options_.qz_criterium)
options_.qz_criterium = 1+1e-6;
end
+ if options_.prior_trunc==0
+ options_.prior_trunc=1.e-10;
+ end
end
if M_.exo_nbr==0
@@ -430,6 +435,7 @@ if options_gsa.rmse
filt_mc_(OutputDirectoryName,options_gsa,dataset_,dataset_info);
end
options_.opt_gsa = options_gsa;
+options_.prior_trunc=original_prior_trunc;
if options_gsa.glue
--
GitLab