/* Copyright (C) 2020 Alessandro Languasco */ /**************** A. LANGUASCO ******************** ************* COMPUTATION OF min_chi |L'/L| MOD q(PRIME) *******/ /******* Direct use of the pari/gp lfun functions *****/ {min_direct(r1,r2,defaultprecision)=local(minutes, millisec, seconds, rows, n,m,L,evenmin,oddmin,q, G, chi, elaptimecomp,absderlog); default(realprecision,defaultprecision); print("************ A. LANGUASCO *************"); print("********* COMPUTATION OF THE MIN |L'/L|**********"); print("******* FOR q IN ONE INTERVAL ********"); print("********* direct use of the lfun command **********"); r1=nextprime(r1); r2=precprime(r2); if (r1>r2, print("error: no odd primes in this interval. END PROGRAM");return); if (r1 < 3, r1=3); print("Interval=[",r1,",",r2,"]"); gettime(); forprime(q=r1,r2, oddmin=10^(10000); evenmin=10^(10000); G= znstar(q, 1); for(m=2, q-1, chi=znconreychar(G,m); L=lfuncreate(Mod(m,q)); \\ defines L as the Dirichlet L-function associated with a non-trivial Dirichlet character (Conrey notation) absderlog=abs(lfun(L,1,1)/lfun(L,1)); \\ computes the absolute value logarithmic derivative of L at 1 if (zncharisodd(G,chi), if(absderlog