Commit Graph

8706 Commits (6bddbca01d0514e4506ba67bfdbd7ee8278a301b)

Author SHA1 Message Date
Sébastien Villemot 6bddbca01d
bytecode_steadystate: turn into a local function of evaluate_steady_state 2023-02-17 22:01:25 -05:00
Sébastien Villemot 1a421d93dd
steady: turn homotopy{1,2,3} into local functions 2023-02-16 14:22:27 -05:00
Sébastien Villemot 2c44923384
No longer turn off output pagination
The “more off” command is incompatible with MATLAB Online.

Moreover, output pagination is disabled by default on MATLAB and on
Octave (since version 4 for the latter), so in most cases that command would do
nothing.

Incidentally, this allows users to enable output pagination should they want
to.

[skip ci]
2023-02-08 11:02:29 -05:00
Sébastien Villemot 65053667cf
New “homotopy_alt_starting_point” option to “perfect_foresight_solver” command 2023-02-07 14:49:00 -05:00
Sébastien Villemot d6370a6d28 Merge branch 'occbin_init_smo' into 'master'
occbin enhancements and bug fixes

See merge request Dynare/dynare!2123
2023-02-07 14:23:14 +00:00
Sébastien Villemot be698d5f98
Perfect foresight + homotopy: add missing newline 2023-02-06 15:14:02 -05:00
Marco Ratto 63a299f64e new error codes introduced and applied for occbin smoother.
- 321 when simulation within occbin smoother fails
- 322 when occbin smoother does not converge.
2023-02-06 14:58:42 +01:00
Marco Ratto f565a0a84a bug fix and avoid missing output arguments and when smoother exits after non convergence of simulation. 2023-02-06 11:39:24 +01:00
Marco Ratto f0aa2fb86f cosmetic change 2023-02-03 14:25:59 +01:00
Marco Ratto b886de92dc trap case when simulation does not converge within smoother 2023-02-03 14:25:32 +01:00
Marco Ratto cdd195576e store binding_indicator_history once it has been set 2023-02-03 14:24:58 +01:00
Marco Ratto 783c237d17 evaluate set_prior only when needed 2023-02-03 14:21:02 +01:00
Marco Ratto 035db8b8c4 bug fixes when using guess regime 2023-02-03 14:20:51 +01:00
Marco Ratto 77e87be48b try with base regime before returning error 2023-02-03 14:20:39 +01:00
Marco Ratto 22e6fb1814 make sure also ALL updated variables in period 1 are reconstructed with smoother_redux. test model changed accordingly. 2023-02-03 14:20:28 +01:00
Marco Ratto 7132cb593e bug fix for init variables to be defined only with occbin PKF 2023-02-03 14:20:16 +01:00
Marco Ratto 50aa20c742 eliminate buggy and useless check for existence of guess regime variable. this may create bugs in running PKF update step! 2023-02-03 14:20:02 +01:00
Marco Ratto e51a30eaf5 Bug fix for the case where no prior is defined for estim params 2023-02-03 14:19:50 +01:00
Marco Ratto a346639020 when filter is NOT diffuse, allow likelihood to handle violation of constraints in period 1. 2023-02-03 14:19:28 +01:00
Marco Ratto 6598d615d1 in some cases, order is not yet set to 1 2023-02-03 14:19:16 +01:00
Marco Ratto a91c36f920 first_period_occbin_update is an index, not a boolean 2023-02-03 14:19:05 +01:00
Marco Ratto 74acb2d1f5 bug fix using index t 2023-02-03 14:18:54 +01:00
Marco Ratto 21c9d59e8c incorporate information about states in period 0 for occbin with smoother_redux 2023-02-03 14:18:32 +01:00
Marco Ratto 0622f01f4e filter out t=1 when storing aalphahat and eetahat 2023-02-03 14:18:21 +01:00
Marco Ratto b90dff1986 retrieve states in period 0 from smoother and deploy this to:
- allow having binding regimes in period 1 with occbin
- further improve recovering of smoothed variables in period 1 under smoother_redux option
2023-02-03 14:18:10 +01:00
Sébastien Villemot 5aca770931
perfect_foresight_with_expectation_errors_{setup,solver}: fix bugs with several exogenous
Closes: #1883
2023-01-24 15:02:07 +01:00
Sébastien Villemot 538f38d269 Merge branch 'occbin_debug' into 'master'
missing_observations_kalman_filter.m: provide debugging information on PKF problems

See merge request Dynare/dynare!2122
2023-01-23 10:39:17 +00:00
Johannes Pfeifer 4307711c60 missing_observations_kalman_filter.m: provide debugging information on PKF problems 2023-01-19 11:57:35 +01:00
Sébastien Villemot cfaf5a685d
Preprocessor: actually obey the “mfs” option for the static model
As a consequence, fix steady state computation when mfs > 0.
2023-01-17 19:08:32 +01:00
Sébastien Villemot 8aabdaee9c
Consolidate functions for solving steady state with block and/or bytecode 2023-01-17 15:25:47 +01:00
Sébastien Villemot 7b6058509f
Drop unused function 2023-01-17 14:48:36 +01:00
Sébastien Villemot 005769210a
Drop unused function 2023-01-17 14:47:52 +01:00
Sébastien Villemot 5342476196
Drop unused function 2023-01-17 14:46:50 +01:00
Sébastien Villemot 023afc7ebb
Drop unused function 2023-01-17 14:46:46 +01:00
Sébastien Villemot 1c5c1639e1
Drop unused function 2023-01-16 14:38:09 +01:00
Sébastien Villemot e5f100b152
Drop unused function 2023-01-16 14:38:06 +01:00
Sébastien Villemot c187d0ead9
Drop unused function
It should probably have been removed in 9c6e219990.
2023-01-16 14:38:03 +01:00
Sébastien Villemot 8f1a4cb363
Drop block-decomposed first-order perturbation solution and kalman filter 2023-01-13 16:57:49 +01:00
Stéphane Adjemian (Ryûk) 91dd917d62
Cosmetic/Efficiency changes.
- Use bsxfun for centering data if possible,
 - Factorise LU decomposition,
 - Remove useless operations during the presampling step.
2023-01-12 17:45:51 +01:00
Sébastien Villemot e250067959
Use sparse representation for block-decomposed perfect foresight and steady state computation
Ref. #1859
2023-01-12 16:28:02 +01:00
Sébastien Villemot 3c55aa57e1
Use sparse representation for evaluating the static model
We also take advantage of the fact that the non-block version is always
available next to the block one, so when we are only interested in the residual
as a whole, we simplify by using the non-block version.
2023-01-10 16:21:44 +01:00
Sébastien Villemot b343f0231d
Ensure that generated sparse MEX files are unloaded before the preprocessor is run
This workaround already existed for the legacy representation.
2023-01-10 16:07:29 +01:00
Johannes Pfeifer 0b0a4ef386 dsge_likelihood.m: flag analytical_derivation as incompatible with endogenous_prior 2023-01-09 17:57:18 +01:00
Sébastien Villemot f84753025d
Bytecode: the block and non-block versions now coexist in parallel
As a consequence, a new “block_decomposed” option of the bytecode MEX has been
introduced to explicitly select the block-decomposed version.

Note that we do not always use the “block_decomposed” option even when the
“block” option has been passed to the user, in situations where the block
decomposition brings nothing (e.g. when evaluating the residuals of the whole
model).
2023-01-09 15:20:46 +01:00
Stéphane Adjemian (Ryûk) b7693c3273
Add routine for conditional likelihood (first order). 2023-01-06 14:57:45 +01:00
Sébastien Villemot f93a7f7192
Merge branch 'disp_aux_4' of git.dynare.org:JohannesPfeifer/dynare
Ref. !2118
2023-01-04 10:59:21 +01:00
Johannes Pfeifer 45450651a2 model_info.m: fix display of block decomposition 2023-01-03 16:17:43 +01:00
Sébastien Villemot ead9023cec Merge branch 'steady_nan' into 'master'
steady state computation: trap cases where parameters are NaN

See merge request Dynare/dynare!2119
2023-01-03 11:24:26 +00:00
Stéphane Adjemian (Guts) 5b72a3041c
Use new unit tests interface. 2023-01-02 22:58:20 +01:00
Johannes Pfeifer 64d61c3268 steady state computation: trap cases where parameters are NaN 2022-12-28 11:48:23 +01:00