********** 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 = 2830676081
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(2830676081) = 1.254528

EK(2830676081)^+ = 12.438044

EK(2830676081) - EK(2830676081)^+ = -11.183516

*** TIMES:
FFT - total computation time: 82.000000 min. 13.000000 sec. 168.957000 msec. 
PLAN FFT - total computation time: 4.000000 min. 46.000000 sec. 847.154000 msec. 
Final comp EKq & EKq^+ - total computation time: 2.000000 min. 58.000000 sec. 925.570000 msec. 
Input data from file - time: 0.000000 min. 0.000000 sec. 0.003000 msec. 
Total elapsed time: 89.000000 min. 58.000000 sec. 941.684000 msec. 
********** END PROGRAM  ***********