details: http://hg.sympy.org/sympy/rev/18e3ec461fa2
changeset: 1820:18e3ec461fa2
user: Ondrej Certik <[EMAIL PROTECTED]>
date: Sun Oct 19 16:15:19 2008 +0200
description:
Fixes an old bug in sympy/core/evalf.py:do_integral().
prec+15 has to be changed to mp.prec. It needs to be mp.prec for oscillatory
quadrature because mpmath significantly increases the precision when quadosc is
called.
This has been broken in SymPy all along. It only worked in the
previous tests by accident. The change of quadrature algorithm in
mpmath 0.10 exposed the problem.
Signed-off-by: Ondrej Certik <[EMAIL PROTECTED]>
Signed-off-by: Fredrik Johansson <[EMAIL PROTECTED]>
diffs (12 lines):
diff -r 3f7eab4594d2 -r 18e3ec461fa2 sympy/core/evalf.py
--- a/sympy/core/evalf.py Sun Oct 19 16:15:19 2008 +0200
+++ b/sympy/core/evalf.py Sun Oct 19 16:15:19 2008 +0200
@@ -623,7 +623,7 @@
max_imag_term = [MINUS_INF]
def f(t):
- re, im, re_acc, im_acc = evalf(func, prec+15, {'subs':{x:t}})
+ re, im, re_acc, im_acc = evalf(func, mp.prec, {'subs':{x:t}})
have_part[0] = re or have_part[0]
have_part[1] = im or have_part[1]
--~--~---------~--~----~------------~-------~--~----~
You received this message because you are subscribed to the Google Groups
"sympy-commits" group.
To post to this group, send email to sympy-commits@googlegroups.com
To unsubscribe from this group, send email to [EMAIL PROTECTED]
For more options, visit this group at
http://groups.google.com/group/sympy-commits?hl=en
-~----------~----~----~----~------~----~------~--~---
