default(realprecision,38);
rk(x)=x=ellinit(x);ellanalyticrank(x);
rk([0, -1, 1, -10, -20])
rk([0, 0, 1, -1, 0])
rk([0, 1, 1, -2, 0])
rk([0, 0, 1, -7, 6])
rk([-5187, 176830])

he(x)=x=ellinit(x);ellheegner(x);
he([1, 1, 0, -1297, -18530])
he([0, -1, 1, -33, 93])
he([-157^2,0])
he([0,0,-9/484,0,-27/234256])
getheap()[1]

allocatemem(30000000);
E=ellinit([0,-1437004800,0,458885065605120000,0]);
ellglobalred(E); ellheegner(E)
