Bill Allombert on Sat, 20 Dec 2008 20:31:40 +0100


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

Re: Functions elllderiv and ellanalyticrank


On Sat, Dec 20, 2008 at 07:21:22PM +0100, Karim Belabas wrote:
> * Bill Allombert [2008-12-20 16:39]:
> > On Sat, Dec 20, 2008 at 10:16:23AM +0100, Karim Belabas wrote:
> > > * Bill Allombert [2008-12-19 23:20]:
> > > > Also the DEBUG message should be improved to be more useful.
> > > 
> > > Done to some extent, I think. More specific suggestions (svn11475 and above) ?
> > 
> > I'd like timing of each ellL1 computation done by ellanalyticrank 
> > separately.
> 
> Done.
> 
> > With \p38, ellanalytic rank fails on a large number (about 1/3 so far)
> > of rank 5 curve in Womack database (it reports rank 1 instead of 5).  
> > 
> > An example: [0, 0, 0, -5187, 176830]
> > ? \p38
> >    realprecision = 38 significant digits
> > ? ellanalyticrank(ellinit([0, 0, 0, -5187, 176830]))
> > %102 = [1, 1278372.8304848296806080373284485378614]
> [...]
> > Clearly the value L'(1) is worng (independently of the rank of the curve).
> 
> Fixed since svn-11480 (the an cache should have been initialized to 0).

Now, ellanalytic report rank 5 for all the rank 5 curves in Womack
database. Good!

One comment on ellanalyticrank():
the instances of ellL1_i could reuse the vecsmall el->ap, growing it if
necessary.

Cheers,
Bill.