On Monday, July 2, 2012 10:03:40 AM UTC-4, kcrisman wrote: > > 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). >
Jeff Denny has also found some other issues on the imaginary axis. I've opened http://trac.sagemath.org/sage_trac/ticket/13193 for this, as 'critical'. -- To post to this group, send an email to sage-devel@googlegroups.com To unsubscribe from this group, send an email to sage-devel+unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-devel URL: http://www.sagemath.org
