From 6e555a1539fc7b97a0ae533ee1c931840638b6f3 Mon Sep 17 00:00:00 2001 From: Michel Juillard Date: Sat, 2 Jun 2012 14:14:24 +0200 Subject: [PATCH] removed specialized code for solving purely forward models: it is true that it is possible to derive the solution function by a simple inversion, but this doesn't check that all eigenvalues are larger than one. It is better to use the standard algorithm with QZ decomposition --- matlab/stochastic_solvers.m | 6 ------ 1 file changed, 6 deletions(-) diff --git a/matlab/stochastic_solvers.m b/matlab/stochastic_solvers.m index 044cc2232..f97c4a6fa 100644 --- a/matlab/stochastic_solvers.m +++ b/matlab/stochastic_solvers.m @@ -156,12 +156,6 @@ if M_.maximum_endo_lead == 0 error(['2nd and 3rd order approximation not implemented for purely ' ... 'backward models']) end -elseif M_.maximum_endo_lag == 0 - % purely forward model - dr.ghx = []; - dr.ghu = -b\jacobia_(:,nz+1:end); - dr.eigval = []; - dr.rank = 0; elseif options_.risky_steadystate [dr,info] = dyn_risky_steadystate_solver(oo_.steady_state,M_,dr, ... options_,oo_);