#13003: numerical evaluation of `erf` crashes PARI at large numbers
-------------------------------------------------------------+--------------
Reporter: benjaminfjones | Owner:
burcin
Type: defect | Status:
new
Priority: minor | Milestone:
sage-5.1
Component: symbolics | Resolution:
Keywords: erf pari sd40.5 | Work issues:
Report Upstream: Reported upstream. Little or no feedback. | Reviewers:
Authors: | Merged in:
Dependencies: | Stopgaps:
-------------------------------------------------------------+--------------
Changes (by jdemeyer):
* upstream: N/A => Reported upstream. Little or no feedback.
Old description:
> 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)
> }}}
New description:
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)
}}}
Reported upstream: [http://pari.math.u-bordeaux.fr/cgi-
bin/bugreport.cgi?bug=1320]
--
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/13003#comment:4>
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.
