From 4043406781cc66b7820d17d3a8e2e128004e4698 Mon Sep 17 00:00:00 2001
From: michel <michel@ac1d8469-bf42-47a9-8791-bf33cf982152>
Date: Mon, 7 Jan 2008 15:59:43 +0000
Subject: [PATCH] v4 dr1.m: added use of options_.qz_criterium with ordqz.
 Important for models with unit root

git-svn-id: https://www.dynare.org/svn/dynare/dynare_v4@1544 ac1d8469-bf42-47a9-8791-bf33cf982152
---
 matlab/dr1.m | 8 +++++---
 1 file changed, 5 insertions(+), 3 deletions(-)

diff --git a/matlab/dr1.m b/matlab/dr1.m
index 8a25e2d07..d9fe4fb48 100644
--- a/matlab/dr1.m
+++ b/matlab/dr1.m
@@ -272,9 +272,11 @@ if exist('ordqz')
   info1 = 0;
   try
     [ss,tt,q,w]=qz(e,d,'real');
-    [ss,tt,q,w]=ordqz(ss,tt,q,w,'udi');
-    dr.eigval = ordeig(ss,tt);
-    sdim = find(abs(dr.eigval)<1);
+    eigval = ordeig(ss,tt);
+    dr.eigval = eigval;
+    select = abs(eigval) < options_.qz_criterium;
+    [ss,tt,q,w]=ordqz(ss,tt,q,w,select);
+    sdim = sum(select);
   catch
     info(1) = 2;
     info(2) = NaN;