\\package:elldata
{
  forell(E,1,9999,
    N = ellconvertname(E[1])[1];
    M = ellglobalred( ellinit(E[2], 1) )[1];
    if (N != M, print(E," bad for N = ",N))
  );
}
ellidentify(ellinit([1,1]))
ellsearch("11a1")
ellsearch("11a")
ellsearch("11b")
ellsearch("11")
ellsearch([11,0,1])
ellsearch([11,0])
ellsearch([11])
ellsearch(11)
ellinit("11a1").j
iferr(ellsearch([0,0,0,0]),E,E)
iferr(ellsearch([0,0,0]),E,E)
iferr(ellsearch("curve"),E,E)
iferr(ellsearch(I),E,E)
c=0;forell(E,1,100,c++); c
c=0;forell(E,1,100,c++,1); c
ellidentify(ellinit([1,1],nfinit('a^2-5)))
{
  K = nfinit(t);
  forell(E, 1, 1000,
    my(e,N,f,v);
    v = [1/60,7131/3600,11/12,5/12];
    e = ellchangecurve(ellinit(E[2],1), v);
    eK= ellchangecurve(ellinit(E[2],K), v);
    N = ellconvertname(E[1])[1];
    f = factor(N)[,1];
    for(i=1, #f,
      my(p = f[i], P,A,B);
      P=idealprimedec(K,p)[1];
      A=elllocalred(e,P.p); A[3]=0;
      B=elllocalred(eK,P); B[3]=0;
      if(A!=B, error(E[2]));
      if(ellap(e,p) != ellap(eK,P),error(E[2]))
    )
  )
}
