extended_path algorithm

[Barbara73]

Hi there,
I have a model with an occasionally binding constraint and I am using the extended_path algorithm (1600 periods, just one shock with sd 0.002) Unfortunately I get the error message reported below. Any clues on how I can proceed? Any hint and/or suggestion will be very appreciated.
Barbara

Warning: Matrix is singular to working precision.
> In solve_perfect_foresight_model at 77
In extended_path at 213
In binding at 345
In dynare at 120
Warning: Matrix is singular to working precision.
> In solve_perfect_foresight_model at 77
In homotopic_steps at 63
In extended_path at 344
In binding at 345
In dynare at 120
Warning: Matrix is singular to working precision.
> In solve_perfect_foresight_model at 77
In homotopic_steps at 63
In extended_path at 344
In binding at 345
In dynare at 120
Warning: Matrix is singular to working precision.
> In solve_perfect_foresight_model at 77
In homotopic_steps at 63
In extended_path at 344
In binding at 345
In dynare at 120
Warning: Matrix is singular to working precision.
> In solve_perfect_foresight_model at 77
In homotopic_steps at 63
In extended_path at 344
In binding at 345
In dynare at 120
Warning: Matrix is singular to working precision.
> In solve_perfect_foresight_model at 77
In homotopic_steps at 63
In extended_path at 344
In binding at 345
In dynare at 120
Warning: Matrix is singular to working precision.
> In solve_perfect_foresight_model at 77
In homotopic_steps at 63
In extended_path at 344
In binding at 345
In dynare at 120
Warning: Matrix is singular to working precision.
> In solve_perfect_foresight_model at 77
In homotopic_steps at 136
In extended_path at 344
In binding at 345
In dynare at 120
??? Error using ==> homotopic_steps at 151
extended_path::homotopy: Oups! I did my best, but I
am not able to simulate this model...

Error in ==> extended_path at 344
[INFO,tmp] =
homotopic_steps(endo_simul,exo_simul_1,.5,.01,pfm1);

Error in ==> binding at 345
oo_.endo_simul = [ oo_.steady_state,
extended_path([], 1600) ];

Error in ==> dynare at 120
evalin('base',fname) ;

[SébastienVillemot]

Your model seems to be singular in some way. Did you first try to run a "simple" deterministic simulation with the "simul" command?

Also, providing your MOD file would help.