diff options
author | 2016-04-09 20:07:36 +0100 | |
---|---|---|
committer | 2016-04-09 20:07:36 +0100 | |
commit | de057ebe541d5a6c1297ea94a89dcaf35582d44e (patch) | |
tree | d43459882cc5c1df686afddb69f1aeb77b2277f3 /Eigen | |
parent | a350c25a396aa4fdef4878d165bb3dbaedf0a4bb (diff) |
Added nans to zeta function.
Diffstat (limited to 'Eigen')
-rw-r--r-- | Eigen/src/Core/SpecialFunctions.h | 5 |
1 files 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 . |