From ffebf6b262039a2682ad2cd2f1bbbaaab192dcc1 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?St=C3=A9phane=20Adjemian=20=28Charybdis=29?=
 <stephane.adjemian@univ-lemans.fr>
Date: Wed, 15 Feb 2017 12:41:56 +0100
Subject: [PATCH] Check positive steady state for loglinear option...

1e-9 is too far from zero. Set 1e-16 as a threshold value for the positivity test.
---
 matlab/resol.m | 6 ++++--
 1 file changed, 4 insertions(+), 2 deletions(-)

diff --git a/matlab/resol.m b/matlab/resol.m
index 2ad99cf2c..ec390ba0f 100644
--- a/matlab/resol.m
+++ b/matlab/resol.m
@@ -109,9 +109,11 @@ if info(1)
 end
 
 if options.loglinear
+    threshold = 1e-16;
     % Find variables with non positive steady state. Skip auxiliary
     % variables for lagges/leaded exogenous variables
-    idx = find(dr.ys(get_all_variables_but_lagged_leaded_exogenous(M)) < 1e-9);
+    idx = find(dr.ys(get_all_variables_but_lagged_leaded_exogenous(M))<threshold);
+    idx
     if length(idx)
         if options.debug
             variables_with_non_positive_steady_state = M.endo_names(idx,:);
@@ -132,7 +134,7 @@ if options.loglinear
             end
         end
         info(1)=26;
-        info(2)=sum(dr.ys(dr.ys<1e-9).^2);
+        info(2)=sum(dr.ys(dr.ys<threshold).^2);
         return
     end
 end