Skip to content
Snippets Groups Projects
Commit a070f5cb authored by Marco Ratto's avatar Marco Ratto
Browse files

Fixed initialization of matrices for posterior irfs.

parent 4480f5b4
Branches
No related tags found
No related merge requests found
...@@ -172,6 +172,7 @@ localVars.IRUN = IRUN; ...@@ -172,6 +172,7 @@ localVars.IRUN = IRUN;
localVars.irun = irun; localVars.irun = irun;
localVars.irun2=irun2; localVars.irun2=irun2;
localVars.nosaddle=nosaddle; localVars.nosaddle=nosaddle;
localVars.npar = npar;
localVars.type=type; localVars.type=type;
if strcmpi(type,'posterior') if strcmpi(type,'posterior')
......
...@@ -53,20 +53,12 @@ IRUN = myinputs.IRUN; ...@@ -53,20 +53,12 @@ IRUN = myinputs.IRUN;
irun =myinputs.irun; irun =myinputs.irun;
irun2=myinputs.irun2; irun2=myinputs.irun2;
nosaddle=myinputs.nosaddle; nosaddle=myinputs.nosaddle;
npar=myinputs.npar;
type=myinputs.type; type=myinputs.type;
if ~strcmpi(type,'prior'), if ~strcmpi(type,'prior'),
x=myinputs.x; x=myinputs.x;
end end
if options_.dsge_var
gend=myinputs.gend;
nvobs=myinputs.nvobs;
NumberOfParametersPerEquation = myinputs.NumberOfParametersPerEquation;
NumberOfLags = myinputs.NumberOfLags;
NumberOfLagsTimesNvobs = myinputs.NumberOfLagsTimesNvobs;
Companion_matrix = myinputs.Companion_matrix;
end
nvar=myinputs.nvar; nvar=myinputs.nvar;
IndxVariables=myinputs.IndxVariables; IndxVariables=myinputs.IndxVariables;
MAX_nirfs_dsgevar=myinputs.MAX_nirfs_dsgevar; MAX_nirfs_dsgevar=myinputs.MAX_nirfs_dsgevar;
...@@ -77,6 +69,17 @@ NumberOfIRFfiles_dsge=myinputs.NumberOfIRFfiles_dsge; ...@@ -77,6 +69,17 @@ NumberOfIRFfiles_dsge=myinputs.NumberOfIRFfiles_dsge;
NumberOfIRFfiles_dsgevar=myinputs.NumberOfIRFfiles_dsgevar; NumberOfIRFfiles_dsgevar=myinputs.NumberOfIRFfiles_dsgevar;
ifil2=myinputs.ifil2; ifil2=myinputs.ifil2;
if options_.dsge_var
gend=myinputs.gend;
nvobs=myinputs.nvobs;
NumberOfParametersPerEquation = myinputs.NumberOfParametersPerEquation;
NumberOfLags = myinputs.NumberOfLags;
NumberOfLagsTimesNvobs = myinputs.NumberOfLagsTimesNvobs;
Companion_matrix = myinputs.Companion_matrix;
stock_irf_bvardsge = zeros(options_.irf,nvobs,M_.exo_nbr,MAX_nirfs_dsgevar);
end
if whoiam if whoiam
Parallel=myinputs.Parallel; Parallel=myinputs.Parallel;
end end
...@@ -129,7 +132,8 @@ if whoiam ...@@ -129,7 +132,8 @@ if whoiam
end end
% Parallel 'while' very good!!! % Parallel 'while' very good!!!
stock_param=zeros(MAX_nruns,npar);
stock_irf_dsge=zeros(options_.irf,nvar,M_.exo_nbr,MAX_nirfs_dsge);
while fpar<npar while fpar<npar
fpar = fpar + 1; fpar = fpar + 1;
irun = irun+1; irun = irun+1;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment