Skip to content
Snippets Groups Projects
Commit aeb1b92e authored by Johannes Pfeifer's avatar Johannes Pfeifer
Browse files

:bug: make_ex_.m: correct matrix dimensions with datafile and exo_nbr==0

homotopy expects a matrix conformable with the initval_series
parent 840d938f
Branches
Tags
1 merge request!2292:bug: make_ex_.m: correct matrix dimensions with datafile and exo_nbr==0
...@@ -49,22 +49,26 @@ if isempty(oo_.initval_series) ...@@ -49,22 +49,26 @@ if isempty(oo_.initval_series)
error('histval and endval cannot be used simultaneously') error('histval and endval cannot be used simultaneously')
end end
end end
elseif M_.exo_nbr > 0 else
x = oo_.initval_series{M_.exo_names{:}}.data; if M_.exo_nbr > 0
oo_.exo_simul = x(M_.orig_maximum_lag-M_.maximum_lag+1:M_.orig_maximum_lag + options_.periods + M_.maximum_lead,:); x = oo_.initval_series{M_.exo_names{:}}.data;
if ~isempty(M_.exo_histval) oo_.exo_simul = x(M_.orig_maximum_lag-M_.maximum_lag+1:M_.orig_maximum_lag + options_.periods + M_.maximum_lead,:);
oo_.exo_simul(1:M_.maximum_lag, :) ... if ~isempty(M_.exo_histval)
= M_.exo_histval(:, 1:M_.maximum_lag)'; oo_.exo_simul(1:M_.maximum_lag, :) ...
= M_.exo_histval(:, 1:M_.maximum_lag)';
end
else
oo_.exo_simul=zeros(M_.maximum_lag + options_.periods + M_.maximum_lead,M_.exo_nbr);
end end
elseif M_.exo_det_nbr > 0 if M_.exo_det_nbr > 0
x_det = oo_.initval_series{M_.exo_det_names{:}}.data; x_det = oo_.initval_series{M_.exo_det_names{:}}.data;
oo_.exo_det_simul = x_det(M_.orig_maximum_lag-M_.maximum_lag+1:M_.orig_maximum_lag + options_.periods + M_.maximum_lead,:); oo_.exo_det_simul = x_det(M_.orig_maximum_lag-M_.maximum_lag+1:M_.orig_maximum_lag + options_.periods + M_.maximum_lead,:);
if ~isempty(M_.exo_det_histval) if ~isempty(M_.exo_det_histval)
oo_.exo_det_simul(1:M_.maximum_lag, :) ... oo_.exo_det_simul(1:M_.maximum_lag, :) ...
= M_.exo_det_histval(:, 1:M_.maximum_lag)'; = M_.exo_det_histval(:, 1:M_.maximum_lag)';
end
end end
end end
% Initialize oo_.exo_det_simul % Initialize oo_.exo_det_simul
if M_.exo_det_nbr > 0 if M_.exo_det_nbr > 0
if isempty(M_.exo_det_histval) if isempty(M_.exo_det_histval)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment