Cohen Henri on Wed, 22 Oct 2025 23:51:41 +0200


[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]

Re: Question on finding a Riemann Zeta function zero for high values of s


Sorry, should have replied to list. In integrand_h0, replace line 309

p1=gmul(expIxy(pmd,gsqr(zn),prec), by

p1 = gmul(pmd, mulcxI(gsqr(zn))); p1r = greal(p1);
 if (gcmpgs(mpabs(p1r), 2*prec2nbits(prec)) > 0) p1 = p1r;
 p1 = gmul(gexp(p1, prec),
 

and declare p1r as a GEN. This is not good, but works here. It is to avoid an

overflow in gsincos when taking the exponential of A+I*B where B is huge and

A is huge negative or positive. The condition I wrote is ridiculous, but I don't know what

to put else.

Henri