sage: pari(3*i).erfc()
-1.76710569338983 E-16 - 1629.86732385786*I
sage: mpmath.erfc(3*i)
mpc(real='1.0', imag='-1629.9946226015657')
sage: 1-pari(3*i).erfc()
1.00000000000000 + 1629.86732385786*I
sage: mpmath.erf(3*i)
mpc(real='0.0', imag='1629.9946226015657')
Now, in 5.1.beta6's gp I get
$ Downloads/sage-5.1.beta6/sage -gp
GP/PARI CALCULATOR Version 2.5.1 (development git-5c5e253)
i386 running darwin (x86-64/GMP-5.0.2 kernel) 64-bit version
compiled: Jun 25 2012, gcc-4.6.3 (GCC)
(readline v6.2 enabled, extended help enabled)
Copyright (C) 2000-2011 The PARI Group
PARI/GP is free software, covered by the GNU General Public License, and
comes
WITHOUT ANY WARRANTY WHATSOEVER.
Type ? for help, \q to quit.
Type ?12 for how to get moral (and possibly technical) support.
parisize = 8000000, primelimit = 500509
? erfc(3*I)
%1 = 0.E-35 - 1629.5516567497094550267455532288372861*I
so apparently this is still an issue. See #11948, which exposed this
(didn't exactly *cause* it).
This would be fixed by #13050 or #1173. Or is it possible that this is not
a Pari bug, but in our using it this way? Note that Wolfram Alpha also
gives the mpmath answers. I'd be grateful for any feedback.
- kcrisman
--
To post to this group, send an email to [email protected]
To unsubscribe from this group, send an email to
[email protected]
For more options, visit this group at http://groups.google.com/group/sage-devel
URL: http://www.sagemath.org
