********** 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 = 1806830951
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(1806830951) = 0.880396

EK(1806830951)^+ = 11.973128

EK(1806830951) - EK(1806830951)^+ = -11.092733

*** TIMES:
FFT - total computation time: 44.000000 min. 58.000000 sec. 358.816000 msec. 
PLAN FFT - total computation time: 1.000000 min. 41.000000 sec. 292.183000 msec. 
Final comp EKq & EKq^+ - total computation time: 1.000000 min. 7.000000 sec. 454.726000 msec. 
Input data from file - time: 0.000000 min. 0.000000 sec. 0.002000 msec. 
Total elapsed time: 47.000000 min. 47.000000 sec. 105.727000 msec. 
********** END PROGRAM  ***********