********** 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 = 9109334831
Generating log(gamma(a/q)) and g^k%q [decimated in frequency]Start FFT lngamma plan 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
Populate lngamma ... Done ... . Done
Start first FFT transform -- End first FFT transform.
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
Populate akvalues_trasf ... 
Done
Start second FFT transform -- End second FFT transform.
Computing the sum over odd characters. Done
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
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(9109334831) = -0.248739

EK(9109334831)^+ = 12.128187

EK(9109334831) - EK(9109334831)^+ = -12.376927

*** TIMES:
FFT - total computation time: 496.000000 min. 40.000000 sec. 863.208000 msec. 
PLAN FFT - total computation time: 491.000000 min. 19.000000 sec. 582.661000 msec. 
Final comp EKq & EKq^+ - total computation time: 12.000000 min. 44.000000 sec. 230.196000 msec. 
Input data from file - time: 0.000000 min. 0.000000 sec. 0.001000 msec. 
Total elapsed time: 1000.000000 min. 44.000000 sec. 676.066000 msec. 
********** END PROGRAM  ***********