Jackson Vanover
Jackson Vanover
SciPy, GSL, and jmat all support evaluating the limit of the besselk function as it approaches infinity, returning zero. However, neither the fp nor the mp contexts for mpmath do...
mpmath version: 1.1.0 When the third argument is a NaN, the gegenbauer function returns numerical results without complaining; perhaps propagating the NaN or throwing an exception is more appropriate? ```...
When evaluating the limit approaching infinity of the hermite polynomial of degree zero, SciPy, GSL, and even the fp version of mpmath function agree that the output should be 1....
Though mpmath appears to support evaluating limits approaching positive infinity, matching the results returned by SciPy and GSL, and jmat, it diverges from those three libraries when evaluating the limit...
When the second argument is negative infinity, `ff` returns a zero regardless of the first argument. When the second argument is a positive infinity, `ff` returns a NaN which seems...
The mp and fp versions of the psi function differ greatly for certain inputs. Here are a few examples that I found of their divergence (starting with one example that...
mpmath version: 1.1.0 For certain inputs, e.g. `hyperu(1, 2, inf)`, the implementations for mpmath, scipy, and gsl all agree on the output (0 in this case). However, when the first...
mpmath version: 1.1.0 I'm no expert on this function, but it appears that the expected behavior for `laguerre(1, a, b)` is to return `1 + a - b`. SciPy agrees...
mpmath version: 1.1.0 `mp.rgamma(inf)` evaluates to 0, but `fp.rgamma(inf)` throws the exception `cannot convert float infinity to integer`. I mention this as a potential issue because I've observed that in...
mpmath version: 1.1.0 When the last argument to hyp2f2 is a NaN, hyp2f2 returns 1. If a NaN is fed to any other argument, a ValueError is thrown. For what...