From de057ebe541d5a6c1297ea94a89dcaf35582d44e Mon Sep 17 00:00:00 2001
From: Till Hoffmann <tillahoffmann@gmail.com>
Date: Sat, 9 Apr 2016 20:07:36 +0100
Subject: [PATCH] Added nans to zeta function.

---
 Eigen/src/Core/SpecialFunctions.h | 5 +++--
 1 file changed, 3 insertions(+), 2 deletions(-)

diff --git a/Eigen/src/Core/SpecialFunctions.h b/Eigen/src/Core/SpecialFunctions.h
index 2a0a6ff15..954972cdd 100644
--- a/Eigen/src/Core/SpecialFunctions.h
+++ b/Eigen/src/Core/SpecialFunctions.h
@@ -881,13 +881,14 @@ struct zeta_impl {
         const Scalar maxnum = NumTraits<Scalar>::infinity();
         const Scalar zero = 0.0, half = 0.5, one = 1.0;
         const Scalar machep = igamma_helper<Scalar>::machep();
+        const Scalar nan = NumTraits<Scalar>::quiet_NaN();
         
         if( x == one )
             return maxnum;
         
         if( x < one )
         {
-            return zero;
+            return nan;
         }
         
         if( q <= zero )
@@ -899,7 +900,7 @@ struct zeta_impl {
             p = x;
             r = numext::floor(p);
             if (p != r)
-                return zero;
+                return nan;
         }
         
         /* Permit negative q but continue sum until n+q > +9 .