Last login: Tue Jan 29 13:33:02 on ttys009 air:~ languasc$ cd /Users/languasc/Documents/LANGUASCO/matematica/lavori/\[55\]Euler-kronecker-constants/programmi-finali/main_and_precomputations/compiled-mac air:compiled-mac languasc$ ./EKmain.exe 9854964401 Give me an odd prime = 9854964401 Computing the primitive root for q = 9854964401 and saving on file Primitive root generation time: 0 min, 0 sec, 0 millisec ****************************************** Now you can call the precomputation programs ****************************************** ********** 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 = 9854964401 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(9854964401) = -0.096465 EK(9854964401)^+ = 12.807752 EK(9854964401) - EK(9854964401)^+ = -12.904217 *** TIMES: FFT - total computation time: 304.000000 min. 2.000000 sec. 736.437000 msec. PLAN FFT - total computation time: 13.000000 min. 17.000000 sec. 659.174000 msec. Final comp EKq & EKq^+ - total computation time: 8.000000 min. 42.000000 sec. 548.136000 msec. Input data from file - time: 0.000000 min. 0.000000 sec. 0.002000 msec. Total elapsed time: 326.000000 min. 2.000000 sec. 943.749000 msec. ********** END PROGRAM ***********