I found the problem. In evalf.py line 626
re, im, re_acc, im_acc = evalf(func, prec+15, {'subs':{x:t}})
should be
re, im, re_acc, im_acc = evalf(func, mp.prec, {'subs':{x:t}})
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.
It's good to pay attention to performance regressions; quite often
there are real bugs hiding beneath...
Fredrik
--~--~---------~--~----~------------~-------~--~----~
You received this message because you are subscribed to the Google Groups
"sympy" 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/sympy?hl=en
-~----------~----~----~----~------~----~------~--~---
