********** A. LANGUASCO and L. RIGHI *********** Computation of the Euler-Kronecker constant with the S function. [LONG DOUBLE PRECISION - GURU64] Acquiring q from program argument: q = 212634221 Start FFT plans preparation Mapping in memory (RW) values array (first step: logGamma values to be computed) Map of lnGamma file/array done Create FFT lnGamma plan... Calling fftwl_plan_guru64_dft_r2c() for lnGamma plan... FFT lnGamma plan created Mapping in memory (RW) values array (second step: akValues values to be computed) Map of akValues file/array done Create FFT akValues plan... Calling fftwl_plan_guru64_dft() for akValues plan... FFT akValues plan created Mapping in memory (RW) values array (third step: precomputed Svalues) Map of Svalues file/array done Create FFT Svalues plan... Calling fftwl_plan_guru64_dft_r2c() for Svalues plan... FFT Svalues plan created . Done Generating log(gamma(a/q)) and g^k%q [decimated in frequency]. Done Start first FFT transform -- End first FFT transform. Start second FFT transform -- End second FFT transform. Computing the sum over odd characters. Done Acquiring precomputed values for S(a/q) [decimated in frequency]. Done Start third FFT transform -- End third FFT transform. Computing the sum over even characters. Done *** RESULTS: EK(212634221) = 1.435141 EK(212634221)^+ = 11.883540 EK(212634221) - EK(212634221)^+ = -10.448398 *** TIMES: FFT - total computation time: 3.000000 min. 48.000000 sec. 606.787000 msec. PLAN FFT - total computation time: 0.000000 min. 32.000000 sec. 503.017000 msec. Final comp EKq & EKq^+ - total computation time: 0.000000 min. 6.000000 sec. 124.798000 msec. Input data from file - time: 0.000000 min. 0.000000 sec. 0.001000 msec. Total elapsed time: 4.000000 min. 27.000000 sec. 234.603000 msec. ********** END PROGRAM ***********