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

evaluate_steady_state.m: do not require second solver call if debug is requested 

All output has already been computed before
kalman-mex
Johannes Pfeifer 2023-08-27 11:25:29 +02:00
parent 757b6e10ec
commit dafe3fbdb2
1 changed files with 2 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
matlab/evaluate_steady_state.m
View File

@ -293,22 +293,19 @@ elseif ~options.bytecode && ~options.block
options.mcppath.lb = lb; options.mcppath.lb = lb;
options.mcppath.ub = ub; options.mcppath.ub = ub;
end end
[ys,check,fvec] = dynare_solve(@static_mcp_problem,... [ys,check,fvec, fjac, errorcode] = dynare_solve(@static_mcp_problem,...
ys_init,... ys_init,...
options.steady.maxit, options.solve_tolf, options.solve_tolx, ... options.steady.maxit, options.solve_tolf, options.solve_tolx, ...
options, exo_ss, params,... options, exo_ss, params,...
M.endo_nbr, static_resid, static_g1, ... M.endo_nbr, static_resid, static_g1, ...
M.static_g1_sparse_rowval, M.static_g1_sparse_colval, M.static_g1_sparse_colptr, eq_index); M.static_g1_sparse_rowval, M.static_g1_sparse_colval, M.static_g1_sparse_colptr, eq_index);
else else
[ys, check] = dynare_solve(@static_problem, ys_init, ... [ys, check, fvec, fjac, errorcode] = dynare_solve(@static_problem, ys_init, ...
options.steady.maxit, options.solve_tolf, options.solve_tolx, ... options.steady.maxit, options.solve_tolf, options.solve_tolx, ...
options, exo_ss, params, M.endo_nbr, static_resid, static_g1, ... options, exo_ss, params, M.endo_nbr, static_resid, static_g1, ...
M.static_g1_sparse_rowval, M.static_g1_sparse_colval, M.static_g1_sparse_colptr); M.static_g1_sparse_rowval, M.static_g1_sparse_colval, M.static_g1_sparse_colptr);
end end
if check && options.debug if check && options.debug
[ys, check, fvec, fjac, errorcode] = dynare_solve(@static_problem, ys_init, ...
options.steady.maxit, options.solve_tolf, options.solve_tolx, ...
options, exo_ss, params, M.endo_nbr, static_resid, static_g1, M.static_g1_sparse_rowval, M.static_g1_sparse_colval, M.static_g1_sparse_colptr);
dprintf('Nonlinear solver routine returned errorcode=%i.', errorcode) dprintf('Nonlinear solver routine returned errorcode=%i.', errorcode)
skipline() skipline()
[infrow,infcol]=find(isinf(fjac) | isnan(fjac)); [infrow,infcol]=find(isinf(fjac) | isnan(fjac));