mjdgges.c: eliminates bogus imaginary NaNs in eigenvalues

2010-05-14 15:03:15 +02:00 · 2010-05-14 15:03:15 +02:00 · 894ed518c8
parent 891ad22564
commit 894ed518c8
1 changed files with 4 additions and 1 deletions
--- 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;