Differences between revisions 2 and 8 (spanning 6 versions)
 ⇤ ← Revision 2 as of 2008-11-05 23:03:53 → Size: 753 Editor: was Comment: ← Revision 8 as of 2008-11-10 03:50:07 → ⇥ Size: 1767 Editor: SouravSenGupta Comment: Deletions are marked like this. Additions are marked like this. Line 2: Line 2: TODO:    1. Modularize the code (break into the following classes and sort dependencies)        * LSeries -------- Is already in there (have to free it from the other 3 modules)        * GammaFactor ------- Done and documented        * IncompleteGeneralizedGammaFactor ------ To be done        * InverseMellinGammaFactor ------- To be done    2. Stuff still computed using gp:        * Delta polynomials in _recursions_at_infinity (search for comment below)        * _without_gp (gamma_series) has this line                sinser = sage_eval(rs(gp_eval('Vec(sin(Pi*(%s)))'%(z0+x))))        * init_Ginf: still uses pari (see below)        * Ginf: still uses pari to evaluate continued fraction    3. Doctest everything, making sure it all works 100% and fix issues with coercion, complex inputs, etc., as they are systematically uncovered.    4. Possibly maybe change from digits to bits prec.    5. Optimize.----------------------------------------------------- Line 20: Line 42: Dokchiter's Paper: attachment:dokchiter.png Dokchiter's Paper: attachment:dokchitser.pdf

# Dokchitser Project for Sage Days 11

TODO:

1. Modularize the code (break into the following classes and sort dependencies)
• LSeries

Is already in there (have to free it from the other 3 modules)

Done and documented

To be done

To be done

1. Stuff still computed using gp:
• Delta polynomials in _recursions_at_infinity (search for comment below)
• _without_gp (gamma_series) has this line
• sinser = sage_eval(rs(gp_eval('Vec(sin(Pi*(%s)))'%(z0+x))))
• init_Ginf: still uses pari (see below)
• Ginf: still uses pari to evaluate continued fraction
2. Doctest everything, making sure it all works 100% and fix issues with coercion, complex inputs, etc., as they are systematically uncovered.
3. Possibly maybe change from digits to bits prec.
4. Optimize.

From Jen:

Here's the version (closest to Dokchitser's original pari code) that still uses continued fraction approximation:

(needs gamma_series.py to run:

The version with Pade approximation (l5.py) has a negligible speedup but only really works for low precision. I'm not sure if Pade gives us a means of computing bounds (I think Mike Rubinstein said that continued fractions won't). Also, l4.py doesn't work for imaginary inputs yet - some coercion with SymbolicRing that I didn't try.

Dokchiter's Paper: attachment:dokchitser.pdf

days11/projects/dokchitser (last edited 2008-11-14 13:42:12 by anonymous)