| 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]