#9969: Update extension code for mpmath-0.17
---------------------------------+------------------------------------------
Reporter: fredrik.johansson | Owner:
Type: enhancement | Status: needs_review
Priority: major | Milestone: sage-4.7
Component: packages | Keywords:
Author: fredrik.johansson | Upstream: N/A
Reviewer: | Merged:
Work_issues: |
---------------------------------+------------------------------------------
Comment(by fbissey):
I got a number of failing mpmath.doctests on my x86 box. A formatting
issue really. Does any of the code interacts with ecl or maxima? The
install I tested it on is based on 4.6.2.alpha4 with ecl-11.1.1 (#10766)
and maxima-5.23.2 (#10773) and a couple of other patch that should be
orthogonal (ppl #10039 and associated patches).
{{{
almosteq 0.004
findroot
**********************************************************************
File "/home/francois/Work/sandbox/sage-4.6.2.alpha4/local/lib/python2.6
/site-packages/mpmath/calculus/optimization.py", line 835, in NoName
Failed example:
findroot(f, -10, solver='anewton', verbose=True)
Expected:
x: -9.88888888888888888889
error: 0.111111111111111111111
converging slowly
x: -9.77890011223344556678
error: 0.10998877665544332211
converging slowly
x: -9.67002233332199662166
error: 0.108877778911448945119
converging slowly
accelerating convergence
x: -9.5622443299551077669
error: 0.107778003366888854764
converging slowly
x: 0.99999999999999999214
error: 10.562244329955107759
x: 1.0
error: 7.8598304758094664213e-18
1.0
Got:
('x: ', -9.88888888888888888889)
('error:', 0.111111111111111111111)
converging slowly
('x: ', -9.77890011223344556678)
('error:', 0.10998877665544332211)
converging slowly
('x: ', -9.67002233332199662166)
('error:', 0.108877778911448945119)
converging slowly
accelerating convergence
('x: ', -9.5622443299551077669)
('error:', 0.107778003366888854764)
converging slowly
('x: ', 0.99999999999999999214)
('error:', 10.562244329955107759)
('x: ', 1.0)
('error:', 7.8598304758094664213e-18)
1.0
0.132
__name__ 0.001
}}}
and
{{{
superfac 0.05
secondzeta
**********************************************************************
File "/home/francois/Work/sandbox/sage-4.6.2.alpha4/local/lib/python2.6
/site-packages/mpmath/functions/zeta.py", line 983, in NoName
Failed example:
secondzeta(0.5 + 40j, error=True, verbose=True)
Expected:
main term = (-30190318549.138656312556 - 13964804384.624622876523j)
computed using 19 zeros of zeta
prime term = (132717176.89212754625045 + 188980555.17563978290601j)
computed using 9 values of the von Mangoldt function
exponential term = (542447428666.07179812536 +
362434922978.80192435203j)
singular term = (512124392939.98154322355 + 348281138038.65531023921j)
((0.059471043 + 0.3463514534j), 1.455191523e-11)
Got:
('main term =', (-30190318549.138656312556 -
13964804384.624622876523j))
(' computed using', 19, 'zeros of zeta')
('prime term =', (132717176.89212754625045 +
188980555.17563978290601j))
(' computed using', 9, 'values of the von Mangoldt function')
('exponential term =', (542447428666.07179812536 +
362434922978.80192435203j))
('singular term =', (512124392939.98154322355 +
348281138038.65531023921j))
((0.059471043 + 0.3463514534j), 1.455191523e-11)
**********************************************************************
File "/home/francois/Work/sandbox/sage-4.6.2.alpha4/local/lib/python2.6
/site-packages/mpmath/functions/zeta.py", line 992, in NoName
Failed example:
secondzeta(0.5 + 40j, a=0.04, error=True, verbose=True)
Expected:
main term = (-151962888.19606243907725 - 217930683.90210294051982j)
computed using 9 zeros of zeta
prime term = (2476659342.3038722372461 + 28711581821.921627163136j)
computed using 37 values of the von Mangoldt function
exponential term = (178506047114.7838188264 +
819674143244.45677330576j)
singular term = (175877424884.22441310708 + 790744630738.28669174871j)
((0.059471043 + 0.3463514534j), 1.455191523e-11)
Got:
('main term =', (-151962888.19606243907725 -
217930683.90210294051982j))
(' computed using', 9, 'zeros of zeta')
('prime term =', (2476659342.3038722372461 +
28711581821.921627163136j))
(' computed using', 37, 'values of the von Mangoldt function')
('exponential term =', (178506047114.7838188264 +
819674143244.45677330576j))
('singular term =', (175877424884.22441310708 +
790744630738.28669174871j))
((0.059471043 + 0.3463514534j), 1.455191523e-11)
34.812
inf 0.002
}}}
I am preparing a run on amd64 with a fresh 4.6.2.rc0 without any patches.
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/9969#comment:23>
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.
