#13003: numerical evaluation of `erf` crashes PARI at large numbers
------------------------------+---------------------------------------------
Reporter: benjaminfjones | Owner: tbd
Type: defect | Status: new
Priority: minor | Milestone: sage-5.1
Component: PLEASE CHANGE | Keywords: erf pari sd40.5
Work issues: | Report Upstream: N/A
Reviewers: | Authors:
Merged in: | Dependencies:
Stopgaps: |
------------------------------+---------------------------------------------
In sage-5.0:
{{{
sage: for i in range(10):
....: print erf(45*10**i).n()
....:
1.00000000000000
1.00000000000000
1.00000000000000
-infinity
1.00000000000000
1.00000000000000
-1.89139086613300e232327576
1.00000000000000
ERROR: An unexpected error occurred while tokenizing input
The following traceback may be corrupted or invalid
The error message is: ('EOF in multi-line statement', (9099, 0))
---------------------------------------------------------------------------
PariError Traceback (most recent call
last)
/Users/jonesbe/sage/sage-5.0/<ipython console> in <module>()
/Users/jonesbe/sage/sage-5.0/local/lib/python2.7/site-
packages/sage/symbolic/expression.so in
sage.symbolic.expression.Expression._numerical_approx
(sage/symbolic/expression.cpp:18303)()
/Users/jonesbe/sage/sage-5.0/local/lib/python2.7/site-
packages/sage/symbolic/expression.so in
sage.symbolic.expression.Expression._convert
(sage/symbolic/expression.cpp:5394)()
/Users/jonesbe/sage/sage-5.0/local/lib/python2.7/site-
packages/sage/functions/other.pyc in _evalf_(self, x, parent)
221 except AttributeError: # not a Sage parent
222 prec = 0
--> 223 return parent(1) - parent(pari(x).erfc(prec))
224
225 def _derivative_(self, x, diff_param=None):
/Users/jonesbe/sage/sage-5.0/local/lib/python2.7/site-
packages/sage/libs/pari/gen.so in sage.libs.pari.gen._pari_trap
(sage/libs/pari/gen.c:49843)()
PariError: (15)
}}}
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/13003>
Sage <http://www.sagemath.org>
Sage: Creating a Viable Open Source Alternative to Magma, Maple, Mathematica,
and MATLAB
--
You received this message because you are subscribed to the Google Groups
"sage-trac" group.
To post to this group, send email to [email protected].
To unsubscribe from this group, send email to
[email protected].
For more options, visit this group at
http://groups.google.com/group/sage-trac?hl=en.
