stepan
/
dynare
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

posterior moments: fix bugs that cause crashes 

- Delete non-existing output argument
- prevent info from being overwritten
trustregion
Johannes Pfeifer 2022-03-14 14:50:49 +01:00
parent 882091cc97
commit 28536314c0
2 changed files with 6 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
matlab/compute_decision_rules.m
View File

@ -1,5 +1,5 @@
function [dr,info,M_,oo_] =compute_decision_rules(M_,options_,oo_) function [dr,info,M_,oo_] =compute_decision_rules(M_,options_,oo_)
% function [dr,info,M_,options_,oo_] =compute_decision_rules(M_,options_,oo_) % function [dr,info,M_,oo_] =compute_decision_rules(M_,options_,oo_)
% INPUTS % INPUTS
% - M_ [structure] Matlab's structure describing the model (M_). % - M_ [structure] Matlab's structure describing the model (M_).
% - options_ [structure] Matlab's structure describing the current options (options_). % - options_ [structure] Matlab's structure describing the current options (options_).

10
matlab/selec_posterior_draws.m
View File

@ -20,7 +20,7 @@ function SampleAddress = selec_posterior_draws(SampleSize,drsize)
% None. % None.
% %
% Copyright (C) 2006-2017 Dynare Team % Copyright (C) 2006-2022 Dynare Team
% %
% This file is part of Dynare. % This file is part of Dynare.
% %
@ -113,9 +113,9 @@ if info
load([BaseName '_mh' num2str(mhfile) '_blck' num2str(mhblck) '.mat'],'x2') load([BaseName '_mh' num2str(mhfile) '_blck' num2str(mhblck) '.mat'],'x2')
end end
pdraws(i,1) = {x2(SampleAddress(i,4),:)}; pdraws(i,1) = {x2(SampleAddress(i,4),:)};
if info-1 if info==2
set_parameters(pdraws{i,1}); set_parameters(pdraws{i,1});
[dr,info,M_,oo_] =compute_decision_rules(M_,options_,oo_); [dr,~,M_,oo_] =compute_decision_rules(M_,options_,oo_);
pdraws(i,2) = { dr }; pdraws(i,2) = { dr };
end end
old_mhfile = mhfile; old_mhfile = mhfile;
@ -140,9 +140,9 @@ if info
load([BaseName '_mh' num2str(mhfile) '_blck' num2str(mhblck) '.mat'],'x2') load([BaseName '_mh' num2str(mhfile) '_blck' num2str(mhblck) '.mat'],'x2')
end end
pdraws(linee,1) = {x2(SampleAddress(i,4),:)}; pdraws(linee,1) = {x2(SampleAddress(i,4),:)};
if info-1 if info==2
set_parameters(pdraws{linee,1}); set_parameters(pdraws{linee,1});
[dr,info,M_,options_,oo_] = compute_decision_rules(M_,options_,oo_); [dr,~,M_,oo_] = compute_decision_rules(M_,options_,oo_);
pdraws(linee,2) = { dr }; pdraws(linee,2) = { dr };
end end
old_mhfile = mhfile; old_mhfile = mhfile;