From beee1d2f822e197b15a21c53394cf812c0f26eef Mon Sep 17 00:00:00 2001
From: ferhat <ferhat.mihoubi@univ-evry.fr>
Date: Thu, 11 Sep 2014 18:19:31 +0200
Subject: [PATCH] Correct the test on the number of eigenvalues outside the
 unit circle: do the test on each block instead of at the global level

---
 matlab/dr_block.m | 6 +++---
 1 file changed, 3 insertions(+), 3 deletions(-)

diff --git a/matlab/dr_block.m b/matlab/dr_block.m
index 1878d572c..c729d6b00 100644
--- a/matlab/dr_block.m
+++ b/matlab/dr_block.m
@@ -463,17 +463,17 @@ for i = 1:Size;
             %First order approximation
             if task ~= 1
                 if nba ~= nyf
-                    sorted_roots = sort(abs(dr.eigval));
                     if isfield(options_,'indeterminacy_continuity')
                         if options_.indeterminacy_msv == 1
-                            [ss,tt,w,q] = qz(e',d');
+                            [ss,tt,w,q] = qz(E',D');
                             [ss,tt,w,junk] = reorder(ss,tt,w,q);
                             ss = ss';
                             tt = tt';
                             w  = w';
-                            %nba = nyf;
+                            nba = nyf;
                         end
                     else
+					    sorted_roots = sort(abs(data(i).eigval));
                         if nba > nyf
                             temp = sorted_roots(nd-nba+1:nd-nyf)-1-options_.qz_criterium;
                             info(1) = 3;