-------------------------------------------------------------------------------------------------------------------- Computation of the Euler-Kronecker constant with the S function. [LONG DOUBLE PRECISION - GURU64] Acquiring q from program argument: q = 2488788101 Start FFT plans preparationMapping in memory (RW) values array (first step: logGamma values to be computed) . 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(2488788101) = 0.424880 EK(2488788101)^+ = 12.248837 EK(2488788101) - EK(2488788101)^+ = -11.823958 *** TIMES: FFT - total computation time: 87.000000 min. 22.000000 sec. 75.642000 msec. PLAN FFT - total computation time: 12.000000 min. 42.000000 sec. 61.142000 msec. Final comp EKq & EKq^+ - total computation time: 3.000000 min. 4.000000 sec. 135.783000 msec. Input data from file - time: 0.000000 min. 0.000000 sec. 0.003000 msec. Total elapsed time: 103.000000 min. 8.000000 sec. 272.570000 msec. ********** END PROGRAM *********** --------------------------------------------------------------------------------------------------------------------