Re: bug in ellpointtoz?

[Bernadette.PERRIN-RIOU@math.u-psud.fr ( Bernadette PERRIN-RIOU )]

Perhaps you can take a more recent version


                  GP/PARI CALCULATOR Version 2.2.6 (development CHANGES-1.705)
                              Power PC (PPC kernel) 32-bit version
                      compiled: Mar 12 2003, gcc-3.1 20020420 (prerelease)
                       (readline not compiled in, extended help available)

                                Copyright (C) 2003 The PARI Group

PARI/GP is free software, covered by the GNU General Public License, and comes
WITHOUT ANY 
WARRANTY WHATSOEVER.

Type ? for help, \q to quit.
Type ?12 for how to get moral (and possibly technical) support.

   realprecision = 28 significant digits
   seriesprecision = 16 significant terms
   format = g0.28

parisize = 10000000, primelimit = 500000
(23:23) gp > e=ellinit([0,1,1,-2,0])
%1 = [0, 1, 1, -2, 0, 4, -4, 1, -3, 112, -856, 389, 1404928/389,
[0.9048929600637101729532368497, 0.1354092402217534448906525937,
-2.040302200285463617843889443]~, 2.490212560855055075321357792,
1.971737701551648204422407698*I, -0.5487442089818677019469720750,
-1.696068988382761177923670949*I, 4.910045991115390176268332642]
(23:24) gp > p=[-2,0]
%2 = [-2, 0]
(23:24) gp > ellisoncurve(e,p)
%3 = 1
(23:24) gp > z=ellpointtoz(e,p)
%4 = 0.07973603912227631670323956179 + 0.9858688507758241022112038492*I
(23:24) gp > ellztopoint(e,z)
%5 = [-2.000000000000000000000000000 + 5.04870979 E-29*I, -1.325286321 E-28 +
2.903008131 E-28*I]
(23:24) gp > 

Bernadette