Re: your mail

[Bill Allombert <Bill.Allombert@math.u-bordeaux.fr>]

On Fri, Feb 17, 2017 at 11:40:35PM +0330, Benyamin Gholami wrote:
> Hi
>  i wanted to compute mestre- nagao sum with pari as sieving method for
> finding high rank fibrations from an elliptic surface and asked my question
> and get this codes:
> 
 S(E, N) =
 { my (s = 0.0);
   forprime(p = 2, N, my(a = ellap(E,p)); s += (2-a)/(p+1-a));
   return (s);
 }
> 
> for(t=-10,10, E=ellinit([0,t,0,t,1]); if(E==[],,print(t,": ",S(E,100))))
> 
> (for example in interval integers between -10 and 10 for y^2=x^3+t*x^2+t*x+1
> )
> it is completely working but i have two questions:
> first: when i enlarge the interval pari although compute and print all
> cases but ater that i cant see all of them . how can fix this problem?

Write the output to a file using write ?
Or use
\l somefile

Or put it in a list:
{
  my(L=List());
  for(t=-10,10, E=ellinit([0,t,0,t,1]); 
    if(E==[],,
      listput(L,[t,S(E,100)])));
  L
}

> second:  i really need to use rational intervals because integer numbers
> induce large coefficients and mwrank will be useless. i need this
> computations for$ t \in a/b$ witch where a and b are restricted and then
> print only curves that have large summations , for example in t=a/b witch
> 1<a,b,<1000 and just print curves with S(E_t,100)>3 with associated t
> value. can any body tell me the codes?

do a double loop.
{
  for(a=1, 1000,
    for(b=1, 1000, if(gcd(a,b)==1,
      my(t=a/b);
      my(Et=ellinit([0,t,0,t,1]));
      if(Et!=[],
        my(E=ellminimalmodel(Et));
        my(s=S(E,100));
        if(s>3, print(t,": ",s);
  )))))
}

Cheers,
Bill.