languasco@OptiPlex-3050:~/Desktop/EK/main_and_precomputations$ chmod 775 compila-optiplex-precomp.sh 
languasco@OptiPlex-3050:~/Desktop/EK/main_and_precomputations$ ./compila-optiplex-precomp.sh 
languasco@OptiPlex-3050:~/Desktop/EK/main_and_precomputations$ ./EKmain.out
Give me an odd prime = 1000003
Computing the primitive root for q = 1000003 and saving on file
Primitive root generation time: 0 min, 0 sec, 0 millisec
******************************************
Now you can call the precomputation programs
******************************************

languasco@OptiPlex-3050:~/Desktop/EK/main_and_precomputations$ chmod 775 lancia.sh 
 
languasco@OptiPlex-3050:~/Desktop/EK/main_and_precomputations$ ./lancia.sh ./m_preclngamma.out 0 1000001 1000002
PRG=./m_preclngamma.out, BEGIN=0, END=1000001, STEP=1000002
x = y = Precomputation lngamma-values for q = 1000003 and saving on file
Precomputation time (I/O included): 0 min, 28 sec, 332 millisec

languasco@OptiPlex-3050:~/Desktop/EK/main_and_precomputations$ ./lancia.sh ./m_precpsi.out 0 1000001 1000002
PRG=./m_precpsi.out, BEGIN=0, END=1000001, STEP=1000002
x = y = Precomputation psi-values for q = 1000003 and saving on file
Precomputation time (I/O included): 0 min, 16 sec, 870 millisec

languasco@OptiPlex-3050:~/Desktop/EK/main_and_precomputations$ ./lancia.sh ./m_precS.out 0 1000001 100000
PRG=./m_precS.out, BEGIN=0, END=1000001, STEP=100000
x = y = Precomputation S-values for q = 1000003 and saving on file
Precomputation time (I/O included): 13 min, 17 sec, 890 millisec
x = y = Precomputation S-values for q = 1000003 and saving on file
Precomputation time (I/O included): 13 min, 23 sec, 417 millisec
x = y = Precomputation S-values for q = 1000003 and saving on file
Precomputation time (I/O included): 13 min, 21 sec, 736 millisec
x = y = Precomputation S-values for q = 1000003 and saving on file
Precomputation time (I/O included): 13 min, 24 sec, 123 millisec
x = y = Precomputation S-values for q = 1000003 and saving on file
Precomputation time (I/O included): 13 min, 20 sec, 308 millisec
x = y = Precomputation S-values for q = 1000003 and saving on file
Precomputation time (I/O included): 13 min, 22 sec, 409 millisec
x = y = Precomputation S-values for q = 1000003 and saving on file

connection dropped
languasco@OptiPlex-3050:~/Desktop/EK/main_and_precomputations$ ./lancia.sh ./m_precS.out 800000 1000001 100000
PRG=./m_precS.out, BEGIN=800000, END=1000001, STEP=100000
x = y = Precomputation S-values for q = 1000003 and saving on file
Precomputation time (I/O included): 13 min, 18 sec, 333 millisec
x = y = Precomputation S-values for q = 1000003 and saving on file
Precomputation time (I/O included): 13 min, 0 sec, 133 millisec
x = y = Precomputation S-values for q = 1000003 and saving on file
Precomputation time (I/O included): 0 min, 0 sec, 16 millisec


languasco@OptiPlex-3050:~/Desktop/EK/main_and_precomputations$ ./lancia.sh ./m_precT.out 0 1000001 100000
PRG=./m_precT.out, BEGIN=0, END=1000001, STEP=100000
x = y = Precomputation T-values for q = 1000003 and saving on file
Precomputation time (I/O included): 11 min, 5 sec, 443 millisec
x = y = Precomputation T-values for q = 1000003 and saving on file
Precomputation time (I/O included): 11 min, 5 sec, 550 millisec
x = y = Precomputation T-values for q = 1000003 and saving on file
Precomputation time (I/O included): 11 min, 4 sec, 132 millisec
x = y = Precomputation T-values for q = 1000003 and saving on file
Precomputation time (I/O included): 11 min, 4 sec, 942 millisec
x = y = Precomputation T-values for q = 1000003 and saving on file
Precomputation time (I/O included): 11 min, 4 sec, 207 millisec
x = y = Precomputation T-values for q = 1000003 and saving on file
Precomputation time (I/O included): 11 min, 4 sec, 269 millisec
x = y = Precomputation T-values for q = 1000003 and saving on file
Precomputation time (I/O included): 11 min, 3 sec, 993 millisec
x = y = Precomputation T-values for q = 1000003 and saving on file

connection dropped
languasco@OptiPlex-3050:~/Desktop/EK/main_and_precomputations$ ./lancia.sh ./m_precT.out 900000 1000001 100000
PRG=./m_precT.out, BEGIN=900000, END=1000001, STEP=100000
x = y = Precomputation T-values for q = 1000003 and saving on file
Precomputation time (I/O included): 11 min, 4 sec, 533 millisec
x = y = Precomputation T-values for q = 1000003 and saving on file
Precomputation time (I/O included): 0 min, 0 sec, 9 millisec

languasco@OptiPlex-3050:~/Desktop/EK/main_and_precomputations$ chmod 775 unisci.sh 
languasco@OptiPlex-3050:~/Desktop/EK/main_and_precomputations$ ./unisci.sh 



languasco@OptiPlex-3050:~/Desktop/EK/for_fftw$ chmod 775 compila-optiplex-fft.sh 
languasco@OptiPlex-3050:~/Desktop/EK/for_fftw$ chmod 775 finale-fft.sh 
languasco@OptiPlex-3050:~/Desktop/EK/for_fftw$ ./compila-optiplex-fft.sh 
languasco@OptiPlex-3050:~/Desktop/EK/for_fftw$ ./finale-fft.sh 
********** A. LANGUASCO  *********** 
Computation of the Euler-Kronecker constant with the T function. [DOUBLE PRECISION]
Acquiring q from file: q = 1000003
Acquiring precomputed values for T(a/q). Done
Acquiring precomputed values for psi(a/q). Done

EK(1000003) = 17.379970

FFT - computation time (I/O excluded): 
0.000000 min. 0.000000 sec. 80.067000 millisec. 
********** END PROGRAM  *********** 
********** A. LANGUASCO  *********** 
Computation of the Euler-Kronecker constant with the T function. [LONG DOUBLE PRECISION]
Acquiring q from file: q = 1000003
Acquiring precomputed values for T(a/q). Done
Acquiring precomputed values for psi(a/q). Done

EK(1000003) = 17.379970

FFT - computation time (I/O excluded): 
0.000000 min. 0.000000 sec. 374.260000 millisec. 
********** END PROGRAM  *********** 
********** A. LANGUASCO  *********** 
Computation of the Euler-Kronecker constant with the S function. [DOUBLE PRECISION]
Acquiring q from file: q = 1000003
Acquiring precomputed values for S(a/q). Done
Acquiring precomputed values for log(gamma(a/q)). Done
Building the sequence g^k%q. Done

EK(1000003)^+ = 15.298449
EK(1000003) = 17.379970

FFT - computation time (I/O excluded): 
0.000000 min. 0.000000 sec. 122.185000 millisec. 
********** END PROGRAM  *********** 
********** A. LANGUASCO  *********** 
Computation of the Euler-Kronecker constant with the S function. [LONG DOUBLE PRECISION]
Acquiring q from file: q = 1000003
Acquiring precomputed values for S(a/q). Done
Acquiring precomputed values for log(gamma(a/q)). Done
Building the sequence g^k%q. Done

EK(1000003)^+ = 15.298449
EK(1000003) = 17.379970

FFT - computation time (I/O excluded): 
0.000000 min. 0.000000 sec. 558.329000 millisec. 
********** END PROGRAM  *********** 
********** A. LANGUASCO and L. RIGHI *********** 
Computation of the Euler-Kronecker constant with the T function. [QUADRUPLE PRECISION]
Acquiring q from file: q = 1000003
Acquiring precomputed values for T(a/q). Done
Acquiring precomputed values for psi(a/q). Done

EK(1000003) = +17.37997023831120957268                      

FFT - computation time (I/O excluded): 
0.000000 min. 12.000000 sec. 266.218000 millisec. 
********** END PROGRAM  *********** 

********** A. LANGUASCO and L. RIGHI *********** 
Computation of the Euler-Kronecker constant with the S function. [QUADRUPLE PRECISION]
Acquiring q from file: q = 1000003
Acquiring precomputed values for S(a/q). Done
Acquiring precomputed values for log(gamma(a/q)). Done
Building the sequence g^k%q. Done

EK(1000003)^+ = +15.29844947429554917174                      

EK(1000003) = +17.37997023827717073547                      

FFT - computation time (I/O excluded): 
0.000000 min. 18.000000 sec. 461.072000 millisec. 
********** END PROGRAM  *********** 