Paul Masson

Hankel functions incorrect for large imaginary parts and standard precision

1

For example: ``` >>> hankel1(1,1+30j) mpc(real='0.0', imag='0.0') >>> hankel2(1,1-30j) mpc(real='0.0', imag='0.0') ``` whereas Mathematica gives ``` In[1]:= HankelH1[1, 1. + 30 I] Out[1]= -7.25495*10^-15 - 1.17346*10^-14 I In[2]:= HankelH2[1, 1....

Streamline Carlson algorithms

1

@fredrik-johansson in reading Carlson's 1995 paper, I noticed ways in which your algorithm for evaluating the RJ integral can be streamlined to speed up the computation. The other cleanup comes...

Minor error in lambertw at -1/e for branches 0 and -1

1

With version 1.1.0: >>> lambertw(-exp(-1),0) mpc(real='-1.0', imag='8.2200797148366183e-9') >>> lambertw(-exp(-1),-1) mpc(real='-1.0', imag='-8.2200797148366183e-9') These should be exactly -1. Perhaps best handled as a special case.

From direct summation, with more terms than necessary: ``` >>> from mpmath import * >>> x=50+5j >>> 1+1/2**x+1/3**x+1/4**x+1/5**x (0.9999999999999991+2.8288195730055187e-16j) >>> zeta(x) mpc(real='0.99999999999999911', imag='2.8288190741413715e-16') ``` The imaginary part returned by Mathematica...

Imaginary parts of error functions are incorrect

1

The error function is odd in its argument, so that comparing values of the complementary error function on either side of the complex plane should merely have a change of...

elliprc has incorrect imaginary part for negative y-axis

3

Comparing values above and below the negative y-axis, >>> elliprc(-1,-2+.001j,pv=False) mpc(real='0.00014269906542846742', imag='-0.78539811887314348') >>> elliprc(-1,-2,pv=False) mpc(real='0.0', imag='0.78539816339744828') >>> elliprc(-1,-2-.001j,pv=False) mpc(real='0.0014280962796194981', imag='2.3561933566202637') the imaginary part on the axis has the wrong sign....

Numerical difference in elliptic pi from Mathematica

6

In mpmath, `ellippi(2,-20,.5)` returns `mpc(real='3.4158110893487281', imag='23.579391735044833')` In Mathematica 10, `EllipticPi[2, -20, .5]` returns `-0.346725 + 1.8138 I`. Wolfram|Alpha returns the same real part but omits the imaginary part. Mathematica apparently...

The command `npm run build` fails due to removed Babel Stage presets, even after attempting update with `npx babel-upgrade`. - Platform / Device: Windows 10 - Node 10 - Babel...

Unusable on Oculus Go

2

When the website is viewed with an Oculus Go, the user is in the middle of the main keyboard. The controller is visible to the far right but can’t access...