John Cremona on Wed, 16 May 2012 14:46:40 +0200 |
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
Re: elliptic curves with prescribed periods |
On 16 May 2012 12:54, Karim Belabas <Karim.Belabas@math.u-bordeaux1.fr> wrote: > > Funny. The whole point of elleisnum(..., flag = 1) is to activate > the following code snippet: > > else if (k==4 && flag) y = gdivgs(y, 12); > else if (k==6 && flag) y = gdivgs(y,-216); > > which your code then gets rid of. Final shortcut: > > x = 0.0017237131288014931634375894696418139876; > y = 0.0088075252675047227730811538038700492673; > om = [2*x,x+y*I]; > ellminimalmodel(ellinit([-27*round(elleisnum(om,4)),54*round(elleisnum(om,6))])) > > :-) > Funny indeed! And I was looking at the code for elleisnum() to see why (in this example) my own C++ code was not getting enough precision to find c4c6 while libpari was doing so easily. John > Cheers, > > K.B. > -- > Karim Belabas, IMB (UMR 5251) Tel: (+33) (0)5 40 00 26 17 > Universite Bordeaux 1 Fax: (+33) (0)5 40 00 69 50 > 351, cours de la Liberation http://www.math.u-bordeaux1.fr/~belabas/ > F-33405 Talence (France) http://pari.math.u-bordeaux1.fr/ [PARI/GP] > `