From 894ed518c83364e29131216be6b02efd2ae423ff Mon Sep 17 00:00:00 2001
From: Michel Juillard <michel.juillard@ens.fr>
Date: Fri, 14 May 2010 15:03:15 +0200
Subject: [PATCH] mjdgges.c: eliminates bogus imaginary NaNs in eigenvalues

---
 mex/sources/mjdgges/mjdgges.c | 5 ++++-
 1 file changed, 4 insertions(+), 1 deletion(-)

diff --git a/mex/sources/mjdgges/mjdgges.c b/mex/sources/mjdgges/mjdgges.c
index ecdec0546..8569cb20e 100644
--- a/mex/sources/mjdgges/mjdgges.c
+++ b/mex/sources/mjdgges/mjdgges.c
@@ -61,7 +61,10 @@ mjdgges(double *a, double *b, double *z, double *n, double *sdim, double *eval_r
   for (per = eval_r; per <= &eval_r[i_n-1]; ++per)
     {
       *per = *par / *pb;
-      *pei = *pai / *pb;
+      if (*pai == 0.0 && *pb == 0.0)
+	*pei = 0.0;
+      else
+	*pei = *pai / *pb;
       ++par;
       ++pai;
       ++pb;