Karim Belabas on Mon, 29 Nov 2004 13:35:35 +0100


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

Re: extreme elllseries() behavior


* Igor Schein [2003-08-20 10:18]:
> On Sat, Aug 02, 2003 at 04:00:01PM +0200, Karim BELABAS wrote:
> > On Fri, 1 Aug 2003, Igor Schein wrote:
> > > ? forstep(k=23,24,.1,print1(k" ");print(elllseries(e,k)))
> > > 23 0.9999999999999715783271077903
> > > 23.10000000000000000000000000 0.9999999999999752574949227059
> > > 23.20000000000000000000000000 0.9999999999999784603968424610
> > > 23.30000000000000000000000000 0.9999999999999812486647844467
> > > 23.40000000000000000000000000 0
> > > 23.50000000000000000000000000
> > >   ***   object too big, length can't fit in a codeword
> > 
> > And what is 'e' here ?
> 
> sorry, ellinit([0,0,1,0,10])
> 
> Using cut-off point !=1 allows to work around it, but the higher s is,
> the further A needs to be from 1, slowing things down considerably, as
> documented.  In the actual example, error, though too general, at least
> prevents the disaster, while value of 0 at s=23.4 is unacceptable IMO.

I believe the above is now fixed in CVS [ and elllseries should be able to
handle complex arguments ]

Cheers,

    Karim.
-- 
Karim Belabas                     Tel: (+33) (0)1 69 15 57 48
Dep. de Mathematiques, Bat. 425   Fax: (+33) (0)1 69 15 60 19
Universite Paris-Sud              http://www.math.u-psud.fr/~belabas/
F-91405 Orsay (France)            http://pari.math.u-bordeaux.fr/  [PARI/GP]