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

fixing bugs in evaluation of steady state file introduced in last 

related commit
time-shift
Michel Juillard 2012-04-21 11:39:21 +02:00
parent 6131e597f3
commit d1120da658
2 changed files with 6 additions and 13 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
matlab/evaluate_steady_state.m
View File

@ -55,7 +55,11 @@ function [ys,params,info] = evaluate_steady_state(ys_init,M,options,oo,steadysta
[ys,params] = dyn_ramsey_static(ys_init,M,options,oo);
elseif steadystate_flag
% explicit steady state file
[ys,params1,check] = evaluate_steady_state_file(ys_init,exo_ss,M,options);
[ys,params1,info] = evaluate_steady_state_file(ys_init,exo_ss,M, ...
options);
if info(1)
return;
end
elseif (options.bytecode == 0 && options.block == 0)
if options.linear == 0
% non linear model
@ -105,15 +109,3 @@ function [ys,params,info] = evaluate_steady_state(ys_init,M,options,oo,steadysta
return
end
if options.steadystate_flag && updated_params_flag && ~isreal(params)
info(1) = 23;
info(2) = sum(imag(params).^2);
return
end
if options.steadystate_flag && updated_params_flag && ~isempty(find(isnan(params)))
info(1) = 24;
info(2) = NaN;
return
end

1
matlab/evaluate_steady_state_file.m
View File

@ -79,6 +79,7 @@ function [ys,params,info] = evaluate_steady_state_file(ys_init,exo_ss,M,options)
% adding values for auxiliary variables
if length(M.aux_vars) > 0
h_set_auxiliary_variables = str2func([M.fname '_set_auxiliary_variables']);
ys = h_set_auxiliary_variables(ys,exo_ss,params);
end