Bill Allombert on Sat, 09 Aug 2014 22:53:10 +0200


[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]

Re: New ellL1 implementation by Pascal Molin


On Wed, Jun 25, 2014 at 02:24:09PM +0200, Bill Allombert wrote:
> On Sun, Jun 22, 2014 at 12:21:25PM +0200, Bill Allombert wrote:
> > Dear PARI developers,
> > 
> > I have merged Pascal Molin code for computing ellL1 and ellanalyticrank.
> > The new code is much faster than the old one for derivatives.
> > For the value of L itself, it is slightly slower, but only because the
> > older code used unproven (and sometimes wrong) bounds for the summation.
> > 
> > This makes ellheegner faster when the height of the point is large.
> 
> I have commited a second patch that adds an alternative summation strategy
> suggested by Pascal. ellL1 will choose the strategy that requires the less
> amount of memory.

I have commited a third patch that implement Brent and Kung summation for
ellheegner too, when the naive height of the point is > 1900*log(2).
This is significantly faster than the previous summation strategy.

The bound 1900 is not very sharp. It is not very easy to tune this.

Cheers,
Bill.