Dynare_Forum.zip (806.2 KB)

Hi, I am having troubles solving the perfect foresight path for small shocks in this model. For some reason, the lmmcp algorithm works, but the standard solver does not.

To replicate the results that lmmcp works , run the matfile RUN_0_Replicate_Error.m

To replicate the result that the standard solver does not work, go to .mod file, and exchange line 349 for “perfect_foresight_solver;”. Then RUN_0_Replicate_Error.m

Any idea of why this might be happenning ?

Best,

Pablo

Could it be that the steady state file is missing?

There is a singularity in the stacked Jacobian that throws off the standard solver. For some reason, this does not throw off the lmmcp-solver, most probably due to the least-squares formulation it employs. A similar thing happens in Dynare failed to do the deterministic simulation - #2 by jpfeifer

Thanks for the prompt reply.

Is there a specific reason why this might happen ? The model seems to be stable, and the stochastic simulation works.

How come then the non linear perfect foresight solver finds this singularity and the lmmcp or the stochastic simulation does not ?

I tried simplifying the model further to the point where the model is just a simple Gertler Karadi. Do you have any idea of why this model does not work with standard perfect foresight ?

Best,

Pablo

It works with standard perfect foresight, just not with the standard numerical solver. I don’t have a good intuition for this. Maybe @MichelJuillard has an idea.