From a623cbc30d46aa34e147b2f82743b48816b3839f Mon Sep 17 00:00:00 2001 From: Sofie Martins Date: Wed, 11 Sep 2024 10:27:36 +0200 Subject: [PATCH] Removed warnings and double reading of input files --- Disconnected/compute_loops.c | 2 - HMC/hmc_forces.c | 73 ------------------- LibHR/Update/Dphi.c | 2 + ModeNumber/mk_eigvals.c | 1 - ModeNumber/mk_modenumber.c | 3 +- Reweight/reweight.c | 28 +------ Scattering/scatter.c | 2 - Scattering/scattering_lengths.c | 4 +- Scattering/sigma_triangle.c | 4 +- Spectrum/measure_formfactor.c | 2 - Spectrum/measure_spectrum.c | 2 - Spectrum/mk_mesons.c | 3 - Spectrum/mk_mesons_with_z2semwall.c | 60 +-------------- Spectrum/mk_mesons_with_z2semwall_new.c | 3 - Spectrum/trunc_mesons.c | 57 +-------------- StaticPotential/tune_HYP_smearing.c | 62 +--------------- TestProgram/Deflate/check_deflate.c | 19 ----- TestProgram/Propagator/check_tm_correlator.c | 2 - TestProgram/RotatedSF/check_rotatedsf_1.c | 18 +---- TestProgram/RotatedSF/check_rotatedsf_10.c | 29 +------- TestProgram/RotatedSF/check_rotatedsf_2.c | 37 +--------- TestProgram/RotatedSF/check_rotatedsf_3.c | 30 +------- TestProgram/RotatedSF/check_rotatedsf_4.c | 30 +------- TestProgram/RotatedSF/check_rotatedsf_5.c | 31 +------- TestProgram/RotatedSF/check_rotatedsf_6.c | 30 +------- TestProgram/RotatedSF/check_rotatedsf_7.c | 28 +------ TestProgram/RotatedSF/check_rotatedsf_8.c | 29 +------- TestProgram/RotatedSF/check_rotatedsf_9.c | 28 +------ .../Scattering/check_scattering_length_I0.c | 4 - .../Scattering/check_scattering_length_I2.c | 2 - .../Scattering/check_scattering_rhopipi.c | 2 - TestProgram/Scattering/check_sigma_triangle.c | 2 - 32 files changed, 19 insertions(+), 610 deletions(-) diff --git a/Disconnected/compute_loops.c b/Disconnected/compute_loops.c index 65564fe40..ac1316c51 100644 --- a/Disconnected/compute_loops.c +++ b/Disconnected/compute_loops.c @@ -79,9 +79,7 @@ int main(int argc, char *argv[]) { setup_process(&argc, &argv); setup_gauge_fields(); - read_input(glb_var.read, get_input_filename()); read_input(disc_var.read, get_input_filename()); - read_input(rlx_var.read, get_input_filename()); strcpy(list_filename, disc_var.configlist); lprintf("MAIN", 0, "list_filename = %s \n", list_filename, disc_var.configlist); diff --git a/HMC/hmc_forces.c b/HMC/hmc_forces.c index 78088ba20..10d5fdabd 100644 --- a/HMC/hmc_forces.c +++ b/HMC/hmc_forces.c @@ -176,63 +176,9 @@ int main(int argc, char *argv[]) { FILE *list; filename_t fpars; - read_cmdline(argc, argv); - /* setup process communications */ setup_process(&argc, &argv); - /* read global variables file */ - read_input(glb_var.read, input_filename); - - setup_replicas(); - - /* logger setup */ - read_input(logger_var.read, input_filename); - logger_set_input(&logger_var); - if (PID != 0) { - logger_disable(); - } /* disable logger for MPI processes != 0 */ - else { - FILE *stderrp; - sprintf(sbuf, ">>%s", output_filename); - logger_stdout(sbuf); - stderrp = freopen(error_filename, "w", stderr); - error(stderrp == NULL, 1, "main [hmc.c]", "Cannot redirect the stderr"); - } - - print_compiling_info_short(); - lprintf("MAIN", 0, "[RepID: %d][world_size: %d]\n[MPI_ID: %d][MPI_size: %d]\n\n", RID, WORLD_SIZE, MPI_PID, MPI_WORLD_SIZE); - - // lprintf("MAIN",0,"Logger lelvel: %d\n",logger_getlevel(0)); - - /* setup lattice geometry */ - if (geometry_init() == 1) { - finalize_process(); - return 0; - } - geometry_mpi_eo(); - /* test_geometry_mpi_eo(); */ - /* setup random numbers */ - read_input(rlx_var.read, input_filename); - //slower(rlx_var.rlxd_start); //convert start variable to lowercase - if (strcmp(rlx_var.rlxd_start, "continue") == 0 && rlx_var.rlxd_state[0] != '\0') { - /*load saved state*/ - lprintf("MAIN", 0, "Loading rlxd state from file [%s]\n", rlx_var.rlxd_state); - read_ranlxd_state(rlx_var.rlxd_state); - } else { - lprintf("MAIN", 0, "RLXD [%d,%d]\n", rlx_var.rlxd_level, rlx_var.rlxd_seed + MPI_PID); - rlxd_init(rlx_var.rlxd_level, rlx_var.rlxd_seed); - } - -#ifdef GAUGE_SUN - lprintf("MAIN", 0, "Gauge group: SU(%d)\n", NG); -#elif GAUGE_SON - lprintf("MAIN", 0, "Gauge group: SO(%d)\n", NG); -#else - lprintf("MAIN", 0, "Default gauge group: SU(%d)\n", NG); -#endif - lprintf("MAIN", 0, "Fermion representation: " REPR_NAME " [dim=%d]\n", NF); - /* Init Monte Carlo */ init_mc(&flow, input_filename); parse_cnfg_filename(cnfg_filename, &fpars); @@ -241,25 +187,6 @@ int main(int argc, char *argv[]) { error((list = fopen(list_filename, "r")) == NULL, 1, "main [mk_mesons.c]", "Failed to open list file\n"); } - print_compiling_info_short(); - lprintf("MAIN", 0, "PId = %d [world_size: %d]\n\n", PID, WORLD_SIZE); - lprintf("MAIN", 0, "input file [%s]\n", input_filename); - lprintf("MAIN", 0, "output file [%s]\n", output_filename); - if (list_filename[0] != 0) { - lprintf("MAIN", 0, "list file [%s]\n", list_filename); - } else { - lprintf("MAIN", 0, "cnfg file [%s]\n", cnfg_filename); - } - - parse_cnfg_filename(cnfg_filename, &fpars); - - GLB_T = fpars.t; - GLB_X = fpars.x; - GLB_Y = fpars.y; - GLB_Z = fpars.z; - error(fpars.type == UNKNOWN_CNFG, 1, "measure_spectrum.c", "Bad name for a configuration file"); - error(fpars.nc != NG, 1, "measure_spectrum.c", "Bad NG"); - int i = 0; while (1) { double times[num_mon()]; diff --git a/LibHR/Update/Dphi.c b/LibHR/Update/Dphi.c index 1a2ac0f8f..0b8742d97 100644 --- a/LibHR/Update/Dphi.c +++ b/LibHR/Update/Dphi.c @@ -371,7 +371,9 @@ void Dphi_cpu_(spinor_field *restrict out, spinor_field *restrict in) { // this is achieved with comparing the condition to be different than repeat=0,1 _MASTER_FOR(out->type, ix) { +#ifdef _OPENMP register int thread0 = hr_threadId(); +#endif register suNf_spinor *r = _FIELD_AT(out, ix); if (repeat == 0) { _spinor_zero_f(*r); } diff --git a/ModeNumber/mk_eigvals.c b/ModeNumber/mk_eigvals.c index 9ead60572..9921c4682 100644 --- a/ModeNumber/mk_eigvals.c +++ b/ModeNumber/mk_eigvals.c @@ -66,7 +66,6 @@ int main(int argc, char *argv[]) { /* setup process id and communications */ setup_process(&argc, &argv); setup_gauge_fields(); - read_input(glb_var.read, get_input_filename()); read_input(eig_var.read, get_input_filename()); lprintf("MAIN", 0, "list file: [%s]\n", eig_var.configlist); if (strcmp(eig_var.configlist, "") != 0) { diff --git a/ModeNumber/mk_modenumber.c b/ModeNumber/mk_modenumber.c index ac9c937fd..aef02aad0 100644 --- a/ModeNumber/mk_modenumber.c +++ b/ModeNumber/mk_modenumber.c @@ -54,7 +54,7 @@ input_nu nu_var = init_input_nu(nu_var); int main(int argc, char *argv[]) { char tmp[1024]; - FILE *list; + FILE *list = NULL; char *cptr; int neig; double M[1024]; @@ -62,7 +62,6 @@ int main(int argc, char *argv[]) { /* setup process id and communications */ setup_process(&argc, &argv); setup_gauge_fields(); - read_input(glb_var.read, get_input_filename()); read_input(nu_var.read, get_input_filename()); hevamass = nu_var.mass; strcpy(tmp, nu_var.list); diff --git a/Reweight/reweight.c b/Reweight/reweight.c index 0f80066ec..acfeb4544 100644 --- a/Reweight/reweight.c +++ b/Reweight/reweight.c @@ -180,40 +180,14 @@ int main(int argc, char *argv[]) { char sbuf[512]; char rbuf[512]; - // Read commandline - read_cmdline(argc, argv); - // Setup process setup_process(&argc, &argv); + setup_gauge_fields(); logger_setlevel(0, 10); - if (PID != 0) { - logger_disable(); - } else { - sprintf(sbuf, ">>%s", output_filename); - logger_stdout(sbuf); - sprintf(sbuf, "err_%d", PID); - freopen(sbuf, "w", stderr); - } - // Read settings - read_input(glb_var.read, input_filename); - read_input(rlx_var.read, input_filename); read_input(rw_var.read, input_filename); - // Initialize stuff - rlxd_init(rlx_var.rlxd_level, rlx_var.rlxd_seed); - - if (geometry_init() == 1) { - finalize_process(); - return 0; - } - - geometry_mpi_eo(); - - // Allocate gauge field - u_gauge = alloc_suNg_field(&glattice); - u_gauge_f = alloc_suNf_field(&glattice); unit_u(u_gauge); represent_gauge_field(); diff --git a/Scattering/scatter.c b/Scattering/scatter.c index 61152a83e..99795c2e1 100644 --- a/Scattering/scatter.c +++ b/Scattering/scatter.c @@ -77,9 +77,7 @@ int main(int argc, char *argv[]) { setup_gauge_fields(); - read_input(glb_var.read, get_input_filename()); read_input(mes_var.read, get_input_filename()); - read_input(rlx_var.read, get_input_filename()); #ifdef WITH_CLOVER clover_init(mes_var.csw); diff --git a/Scattering/scattering_lengths.c b/Scattering/scattering_lengths.c index 0fe5fca5a..9b494327a 100644 --- a/Scattering/scattering_lengths.c +++ b/Scattering/scattering_lengths.c @@ -70,7 +70,7 @@ typedef struct { int main(int argc, char *argv[]) { int k; - FILE *list; + FILE *list = NULL; int nm; double m[256]; int seq_prop = 0; @@ -79,9 +79,7 @@ int main(int argc, char *argv[]) { setup_gauge_fields(); - read_input(glb_var.read, get_input_filename()); read_input(mes_var.read, get_input_filename()); - read_input(rlx_var.read, get_input_filename()); if (strcmp(mes_var.seq_prop, "true") == 0) { seq_prop = 1; } if (strcmp(mes_var.seq_prop, "none") == 0) { seq_prop = 2; } if (strcmp(mes_var.seq_prop, "false") == 0) { seq_prop = 0; } diff --git a/Scattering/sigma_triangle.c b/Scattering/sigma_triangle.c index 31a5f460b..728b5fdb8 100644 --- a/Scattering/sigma_triangle.c +++ b/Scattering/sigma_triangle.c @@ -72,7 +72,7 @@ typedef struct { int main(int argc, char *argv[]) { int k; - FILE *list; + FILE *list = NULL; int nm; double m[256]; int seq_prop = 0; @@ -81,9 +81,7 @@ int main(int argc, char *argv[]) { setup_gauge_fields(); - read_input(glb_var.read, get_input_filename()); read_input(mes_var.read, get_input_filename()); - read_input(rlx_var.read, get_input_filename()); if (strcmp(mes_var.seq_prop, "true") == 0) { seq_prop = 1; } strcpy(list_filename, mes_var.configlist); diff --git a/Spectrum/measure_formfactor.c b/Spectrum/measure_formfactor.c index 9d65a3d7d..0ffdbaa23 100644 --- a/Spectrum/measure_formfactor.c +++ b/Spectrum/measure_formfactor.c @@ -170,8 +170,6 @@ int main(int argc, char *argv[]) { read_cmdline(argc, argv); setup_process(&argc, &argv); - read_input(glb_var.read, input_filename); - /* logger setup */ /* disable logger for MPI processes != 0 */ logger_setlevel(0, 30); diff --git a/Spectrum/measure_spectrum.c b/Spectrum/measure_spectrum.c index 6c6a9804f..94e6d304b 100644 --- a/Spectrum/measure_spectrum.c +++ b/Spectrum/measure_spectrum.c @@ -135,9 +135,7 @@ int main(int argc, char *argv[]) { setup_gauge_fields(); - read_input(glb_var.read, get_input_filename()); read_input(mes_var.read, get_input_filename()); - read_input(rlx_var.read, get_input_filename()); read_input(ff_var.read, get_input_filename()); strcpy(list_filename, mes_var.configlist); diff --git a/Spectrum/mk_mesons.c b/Spectrum/mk_mesons.c index bc609447d..5bf64aa63 100644 --- a/Spectrum/mk_mesons.c +++ b/Spectrum/mk_mesons.c @@ -220,9 +220,6 @@ int main(int argc, char *argv[]) { read_cmdline(argc, argv); setup_process(&argc, &argv); - read_input(glb_var.read, input_filename); - read_input(rlx_var.read, input_filename); - /* logger setup */ /* disable logger for MPI processes != 0 */ logger_setlevel(0, 30); diff --git a/Spectrum/mk_mesons_with_z2semwall.c b/Spectrum/mk_mesons_with_z2semwall.c index 9ac5d19ae..46d31014b 100644 --- a/Spectrum/mk_mesons_with_z2semwall.c +++ b/Spectrum/mk_mesons_with_z2semwall.c @@ -156,50 +156,10 @@ int main(int argc, char *argv[]) { double m[256]; /* setup process id and communications */ - read_cmdline(argc, argv); setup_process(&argc, &argv); - - read_input(glb_var.read, input_filename); - read_input(rlx_var.read, input_filename); - - /* logger setup */ - /* disable logger for MPI processes != 0 */ - logger_setlevel(0, 30); - if (PID != 0) { logger_disable(); } - if (PID == 0) { - sprintf(tmp, ">%s", output_filename); - logger_stdout(tmp); - sprintf(tmp, "err_%d", PID); - if (!freopen(tmp, "w", stderr)) { lprintf("MAIN", 0, "Error out not open\n"); } - } - - print_compiling_info_short(); - lprintf("MAIN", 0, "PId = %d [world_size: %d]\n\n", PID, WORLD_SIZE); - lprintf("MAIN", 0, "input file [%s]\n", input_filename); - lprintf("MAIN", 0, "output file [%s]\n", output_filename); - if (list_filename != NULL) { - lprintf("MAIN", 0, "list file [%s]\n", list_filename); - } else { - lprintf("MAIN", 0, "cnfg file [%s]\n", cnfg_filename); - } - - /* read & broadcast parameters */ - parse_cnfg_filename(cnfg_filename, &fpars); + setup_gauge_fields(); read_input(mes_var.read, input_filename); - GLB_T = fpars.t; - GLB_X = fpars.x; - GLB_Y = fpars.y; - GLB_Z = fpars.z; - error(fpars.type == UNKNOWN_CNFG, 1, "mk_mesons.c", "Bad name for a configuration file"); - error(fpars.nc != NG, 1, "mk_mesons.c", "Bad NG"); - - lprintf("MAIN", 0, "RLXD [%d,%d]\n", rlx_var.rlxd_level, rlx_var.rlxd_seed); - rlxd_init(rlx_var.rlxd_level, rlx_var.rlxd_seed); - srand(rlx_var.rlxd_seed + PID); - - lprintf("MAIN", 0, "Gauge group: SU(%d)\n", NG); - lprintf("MAIN", 0, "Fermion representation: " REPR_NAME " [dim=%d]\n", NF); nm = 0; if (fpars.type == DYNAMICAL_CNFG) { @@ -216,24 +176,6 @@ int main(int argc, char *argv[]) { } } - /* setup communication geometry */ - if (geometry_init() == 1) { - finalize_process(); - return 0; - } - - /* setup lattice geometry */ - geometry_mpi_eo(); - /* test_geometry_mpi_eo(); */ - - init_BCs(NULL); - - /* alloc global gauge fields */ - u_gauge = alloc_suNg_field(&glattice); -#ifdef ALLOCATE_REPR_GAUGE_FIELD - u_gauge_f = alloc_suNf_field(&glattice); -#endif - lprintf("MAIN", 0, "Inverter precision = %e\n", mes_var.precision); for (k = 0; k < nm; k++) { lprintf("MAIN", 0, "Mass[%d] = %f\n", k, m[k]); diff --git a/Spectrum/mk_mesons_with_z2semwall_new.c b/Spectrum/mk_mesons_with_z2semwall_new.c index 15dd393b0..f6285f9e9 100644 --- a/Spectrum/mk_mesons_with_z2semwall_new.c +++ b/Spectrum/mk_mesons_with_z2semwall_new.c @@ -160,9 +160,6 @@ int main(int argc, char *argv[]) { read_cmdline(argc, argv); setup_process(&argc, &argv); - read_input(glb_var.read, input_filename); - read_input(rlx_var.read, input_filename); - /* logger setup */ /* disable logger for MPI processes != 0 */ logger_setlevel(0, 30); diff --git a/Spectrum/trunc_mesons.c b/Spectrum/trunc_mesons.c index 47104e792..d5ceee11e 100644 --- a/Spectrum/trunc_mesons.c +++ b/Spectrum/trunc_mesons.c @@ -169,49 +169,10 @@ int main(int argc, char *argv[]) { double disc; /* setup process id and communications */ - read_cmdline(argc, argv); setup_process(&argc, &argv); + setup_gauge_fields(); - /* logger setup */ - /* disable logger for MPI processes != 0 */ - logger_setlevel(0, 30); - if (PID != 0) { logger_disable(); } - if (PID == 0) { - sprintf(tmp, ">%s", output_filename); - logger_stdout(tmp); - sprintf(tmp, "err_%d", PID); - freopen(tmp, "w", stderr); - } - - print_compiling_info_short(); - lprintf("MAIN", 0, "PId = %d [world_size: %d]\n\n", PID, WORLD_SIZE); - lprintf("MAIN", 0, "input file [%s]\n", input_filename); - lprintf("MAIN", 0, "output file [%s]\n", output_filename); - if (list_filename != NULL) { - lprintf("MAIN", 0, "list file [%s]\n", list_filename); - } else { - lprintf("MAIN", 0, "cnfg file [%s]\n", cnfg_filename); - } - - /* read & broadcast parameters */ - parse_cnfg_filename(cnfg_filename, &fpars); - - read_input(glb_var.read, input_filename); read_input(ata_qprop_var.read, input_filename); - GLB_T = fpars.t; - GLB_X = fpars.x; - GLB_Y = fpars.y; - GLB_Z = fpars.z; - error(fpars.type == UNKNOWN_CNFG, 1, "mk_hairpins.c", "Bad name for a configuration file"); - error(fpars.nc != NG, 1, "mk_hairpins.c", "Bad NG"); - - lprintf("MAIN", 0, "RLXD [%d,%d]\n", glb_var.rlxd_level, glb_var.rlxd_seed); - rlxd_init(glb_var.rlxd_level, glb_var.rlxd_seed); - srand(glb_var.rlxd_seed + PID); - - lprintf("MAIN", 0, "Gauge group: SU(%d)\n", NG); - lprintf("MAIN", 0, "Fermion representation: " REPR_NAME " [dim=%d]\n", NF); - ata_qprop_var.pars.n_masses = 0; if (fpars.type == DYNAMICAL_CNFG) { ata_qprop_var.pars.n_masses = 1; @@ -227,16 +188,6 @@ int main(int argc, char *argv[]) { } } - /* setup communication geometry */ - if (geometry_init() == 1) { - finalize_process(); - return 0; - } - - /* setup lattice geometry */ - geometry_mpi_eo(); - /* test_geometry_mpi_eo(); */ - /* Print inversion/truncation parameters */ lprintf("MAIN", 0, "Simulation masses:"); for (i = 0; i < ata_qprop_var.pars.n_masses; i++) { @@ -253,12 +204,6 @@ int main(int argc, char *argv[]) { lprintf("MAIN", 0, "Inverter precision: %e\n", ata_qprop_var.pars.inverter_precision); lprintf("MAIN", 0, "Dilution: %d\n", ata_qprop_var.pars.dilution); - /* alloc global gauge fields */ - u_gauge = alloc_suNg_field(&glattice); -#ifndef REPR_FUNDAMENTAL - u_gauge_f = alloc_suNf_field(&glattice); -#endif - ata_qprop_init(&(ata_qprop_var.pars)); haircorr = (double *)malloc(GLB_T * sizeof(double)); diff --git a/StaticPotential/tune_HYP_smearing.c b/StaticPotential/tune_HYP_smearing.c index cd2e3125a..5047f434b 100644 --- a/StaticPotential/tune_HYP_smearing.c +++ b/StaticPotential/tune_HYP_smearing.c @@ -129,74 +129,14 @@ void read_cmdline(int argc, char *argv[]) { int main(int argc, char *argv[]) { int i; - char tmp[256]; FILE *list; - filename_t fpars; double mtp[6859], mtp0, mtp1; double weight[3]; int ibest; /* setup process id and communications */ - read_cmdline(argc, argv); setup_process(&argc, &argv); - - /* logger setup */ - /* disable logger for MPI processes != 0 */ - logger_setlevel(0, 10); - if (PID != 0) { logger_disable(); } - if (PID == 0) { - sprintf(tmp, ">%s", output_filename); - logger_stdout(tmp); - sprintf(tmp, "err_%d", PID); - freopen(tmp, "w", stderr); - } - - print_compiling_info_short(); - lprintf("MAIN", 0, "PId = %d [world_size: %d]\n\n", PID, WORLD_SIZE); - lprintf("MAIN", 0, "input file [%s]\n", input_filename); - lprintf("MAIN", 0, "output file [%s]\n", output_filename); - if (strcmp(list_filename, "") != 0) { - lprintf("MAIN", 0, "list file [%s]\n", list_filename); - } else { - lprintf("MAIN", 0, "cnfg file [%s]\n", cnfg_filename); - } - - /* read & broadcast parameters */ - parse_cnfg_filename(cnfg_filename, &fpars); - - read_input(glb_var.read, input_filename); - GLB_T = fpars.t; - GLB_X = fpars.x; - GLB_Y = fpars.y; - GLB_Z = fpars.z; - error(fpars.type == UNKNOWN_CNFG, 1, "tune_HYP_smearing.c", "Bad name for a configuration file"); - error(fpars.nc != NG, 1, "tune_HYP_smearing.c", "Bad NG"); - - lprintf("MAIN", 0, "Gauge group: SU(%d)\n", NG); - lprintf("MAIN", 0, "Fermion representation: " REPR_NAME " [dim=%d]\n", NF); - - /* setup communication geometry */ - if (geometry_init() == 1) { - finalize_process(); - return 0; - } - - /* setup lattice geometry */ - geometry_mpi_eo(); - /* test_geometry_mpi_eo(); */ - - /* setup random numbers */ - read_input(rlx_var.read, input_filename); - lprintf("MAIN", 0, "RLXD [%d,%d]\n", rlx_var.rlxd_level, rlx_var.rlxd_seed + MPI_PID); - rlxd_init(rlx_var.rlxd_level, rlx_var.rlxd_seed); - - init_BCs(NULL); - - /* alloc global gauge fields */ - u_gauge = alloc_suNg_field(&glattice); -#ifdef ALLOCATE_REPR_GAUGE_FIELD - u_gauge_f = alloc_suNf_field(&glattice); -#endif + setup_gauge_fields(); list = NULL; if (strcmp(list_filename, "") != 0) { diff --git a/TestProgram/Deflate/check_deflate.c b/TestProgram/Deflate/check_deflate.c index 7873b23dd..82d9cf167 100644 --- a/TestProgram/Deflate/check_deflate.c +++ b/TestProgram/Deflate/check_deflate.c @@ -189,28 +189,9 @@ int main(int argc, char *argv[]) { /* setup process id and communications */ logger_map("DEBUG", "debug"); - - read_cmdline(argc, argv); setup_process(&argc, &argv); - setup_gauge_fields(); - - /* setup process id and communications */ - setup_process(&argc, &argv); - - read_input(glb_var.read, input_filename); - - /* read & broadcast parameters */ - parse_cnfg_filename(cnfg_filename, &fpars); - read_input(mes_var.read, input_filename); - GLB_T = fpars.t; - GLB_X = fpars.x; - GLB_Y = fpars.y; - GLB_Z = fpars.z; - error(fpars.type == UNKNOWN_CNFG, 1, "check_propagator.c", "Bad name for a configuration file"); - error(fpars.nc != NG, 1, "check_propagator.c", "Bad NG"); - nm = 0; if (fpars.type == DYNAMICAL_CNFG) { nm = 1; diff --git a/TestProgram/Propagator/check_tm_correlator.c b/TestProgram/Propagator/check_tm_correlator.c index e7f62db37..669a2fd6b 100644 --- a/TestProgram/Propagator/check_tm_correlator.c +++ b/TestProgram/Propagator/check_tm_correlator.c @@ -257,9 +257,7 @@ int main(int argc, char *argv[]) { spinor_field *prop = alloc_spinor_field(4, &glattice); spinor_field *propmu0 = alloc_spinor_field(4, &glattice); - read_input(glb_var.read, get_input_filename()); read_input(cor_var.read, get_input_filename()); - read_input(rlx_var.read, get_input_filename()); m = atof(cor_var.mstring); mu = atof(cor_var.mustring); diff --git a/TestProgram/RotatedSF/check_rotatedsf_1.c b/TestProgram/RotatedSF/check_rotatedsf_1.c index 3dd7a932a..9bdaacca7 100644 --- a/TestProgram/RotatedSF/check_rotatedsf_1.c +++ b/TestProgram/RotatedSF/check_rotatedsf_1.c @@ -72,6 +72,7 @@ int main(int argc, char *argv[]) { double acc = 1.e-26, tau, sig; spinor_field *s0, *s1, *s2, *s3; setup_process(&argc, &argv); + setup_gauge_fields(); BCs_pars_t BCs_pars = { .fermion_twisting_theta = { 0., M_PI / 5., M_PI / 5., M_PI / 5. }, .gauge_boundary_improvement_cs = 1., @@ -97,25 +98,8 @@ int main(int argc, char *argv[]) { freopen(sbuf, "w", stderr); } - lprintf("MAIN", 0, "PId = %d [world_size: %d]\n\n", PID, WORLD_SIZE); - - read_input(glb_var.read, "test_input"); - read_input(rlx_var.read, "test_input"); - - rlxd_init(rlx_var.rlxd_level, rlx_var.rlxd_seed); - - /* setup communication geometry */ - if (geometry_init() == 1) { - finalize_process(); - return 0; - } - - geometry_mpi_eo(); - init_BCs(&BCs_pars); - u_gauge = alloc_suNg_field(&glattice); - u_gauge_f = alloc_suNf_field(&glattice); s0 = alloc_spinor_field(4, &glattice); s1 = s0 + 1; s2 = s1 + 1; diff --git a/TestProgram/RotatedSF/check_rotatedsf_10.c b/TestProgram/RotatedSF/check_rotatedsf_10.c index 292ffae60..ba0ed1c86 100644 --- a/TestProgram/RotatedSF/check_rotatedsf_10.c +++ b/TestProgram/RotatedSF/check_rotatedsf_10.c @@ -9,6 +9,7 @@ rhmc_par _update_par; int main(int argc, char *argv[]) { setup_process(&argc, &argv); + setup_gauge_fields(); BCs_pars_t BCs_pars = { .fermion_twisting_theta = { 0., M_PI / 5., M_PI / 5., M_PI / 5. }, @@ -28,42 +29,14 @@ int main(int argc, char *argv[]) { _update_par.SF_ct = 1; _update_par.SF_zf = 1.0; - char tmp[256]; - logger_setlevel(0, 100); /* log all */ - if (PID != 0) { - logger_disable(); - } else { - sprintf(tmp, ">out_%d", PID); - logger_stdout(tmp); - sprintf(tmp, "err_%d", PID); - freopen(tmp, "w", stderr); - } - - lprintf("MAIN", 0, "PId = %d [world_size: %d]\n\n", PID, WORLD_SIZE); - - read_input(glb_var.read, "test_input"); - - rlxd_init(rlx_var.rlxd_level, rlx_var.rlxd_seed); - #if NG != 3 || NF != 3 #error "Can work only with NC=3 and Nf==3" #endif - /* setup communication geometry */ - if (geometry_init() == 1) { - finalize_process(); - return 0; - } - - geometry_mpi_eo(); - init_BCs(&BCs_pars); lprintf("MAIN", 0, "This test implements a comparison with a working code of Stefan Sint\n"); - u_gauge = alloc_suNg_field(&glattice); - u_gauge_f = alloc_suNf_field(&glattice); - read_gauge_field_nocheck("suNg_field_sint_2.dat"); apply_BCs_on_fundamental_gauge_field(); diff --git a/TestProgram/RotatedSF/check_rotatedsf_2.c b/TestProgram/RotatedSF/check_rotatedsf_2.c index 0e058c7a1..7a25dbf02 100644 --- a/TestProgram/RotatedSF/check_rotatedsf_2.c +++ b/TestProgram/RotatedSF/check_rotatedsf_2.c @@ -63,6 +63,7 @@ static void transform_u(void) { int main(int argc, char *argv[]) { setup_process(&argc, &argv); + setup_gauge_fields(); double acc = 1.e-20; @@ -72,49 +73,13 @@ int main(int argc, char *argv[]) { .chiSF_boundary_improvement_ds = 0.5, .SF_BCs = 0 }; - char tmp[256]; - logger_setlevel(0, 100); /* log all */ - if (PID != 0) { - logger_disable(); - } else { - sprintf(tmp, ">out_%d", PID); - logger_stdout(tmp); - sprintf(tmp, "err_%d", PID); - freopen(tmp, "w", stderr); - } - - lprintf("MAIN", 0, "PId = %d [world_size: %d]\n\n", PID, WORLD_SIZE); - - read_input(glb_var.read, "test_input"); - read_input(rlx_var.read, "test_input"); - - rlxd_init(rlx_var.rlxd_level, rlx_var.rlxd_seed); - - /* setup communication geometry */ - if (geometry_init() == 1) { - finalize_process(); - return 0; - } - - geometry_mpi_eo(); - init_BCs(&BCs_pars); - lprintf("MAIN", 0, "Gauge group: SU(%d)\n", NG); - lprintf("MAIN", 0, "Fermion representation: dim = %d\n", NF); - lprintf("MAIN", 0, "The lattice size is %dx%dx%dx%d\n", T, X, Y, Z); - lprintf("MAIN", 0, "The lattice global size is %dx%dx%dx%d\n", GLB_T, GLB_X, GLB_Y, GLB_Z); - lprintf("MAIN", 0, "The lattice borders are (%d,%d,%d,%d)\n", T_BORDER, X_BORDER, Y_BORDER, Z_BORDER); - lprintf("MAIN", 0, "\n"); - fflush(stdout); - double mass = 0.0; _update_par.SF_zf = 6.; _update_par.SF_ds = 3.; _update_par.SF_sign = 1; - u_gauge = alloc_suNg_field(&glattice); - u_gauge_f = alloc_suNf_field(&glattice); random_u(u_gauge); apply_BCs_on_fundamental_gauge_field(); diff --git a/TestProgram/RotatedSF/check_rotatedsf_3.c b/TestProgram/RotatedSF/check_rotatedsf_3.c index 6a6375a85..6327d89c6 100644 --- a/TestProgram/RotatedSF/check_rotatedsf_3.c +++ b/TestProgram/RotatedSF/check_rotatedsf_3.c @@ -12,6 +12,7 @@ int main(int argc, char *argv[]) { double acc = 1.e-20, tau, sig; spinor_field *s0, *s1, *s2, *s3; setup_process(&argc, &argv); + setup_gauge_fields(); BCs_pars_t BCs_pars = { .fermion_twisting_theta = { 0., M_PI / 5., M_PI / 5., M_PI / 5. }, .gauge_boundary_improvement_cs = 1., @@ -25,36 +26,7 @@ int main(int argc, char *argv[]) { _update_par.SF_ct = 1; _update_par.SF_zf = 1.3; - char tmp[256]; - logger_setlevel(0, 100); /* log all */ - if (PID != 0) { - logger_disable(); - } else { - sprintf(tmp, ">out_%d", PID); - logger_stdout(tmp); - sprintf(tmp, "err_%d", PID); - freopen(tmp, "w", stderr); - } - - lprintf("MAIN", 0, "PId = %d [world_size: %d]\n\n", PID, WORLD_SIZE); - - read_input(glb_var.read, "test_input"); - read_input(rlx_var.read, "test_input"); - - rlxd_init(rlx_var.rlxd_level, rlx_var.rlxd_seed); - - /* setup communication geometry */ - if (geometry_init() == 1) { - finalize_process(); - return 0; - } - - geometry_mpi_eo(); - init_BCs(&BCs_pars); - - u_gauge = alloc_suNg_field(&glattice); - u_gauge_f = alloc_suNf_field(&glattice); s0 = alloc_spinor_field(4, &glattice); s1 = s0 + 1; s2 = s1 + 1; diff --git a/TestProgram/RotatedSF/check_rotatedsf_4.c b/TestProgram/RotatedSF/check_rotatedsf_4.c index a57dd0560..eb696214c 100644 --- a/TestProgram/RotatedSF/check_rotatedsf_4.c +++ b/TestProgram/RotatedSF/check_rotatedsf_4.c @@ -12,6 +12,7 @@ int main(int argc, char *argv[]) { double acc = 1.e-20, tau, sig; spinor_field *s0, *s1, *s2, *s3; setup_process(&argc, &argv); + setup_gauge_fields(); BCs_pars_t BCs_pars = { .fermion_twisting_theta = { 0., M_PI / 5., M_PI / 5., M_PI / 5. }, .gauge_boundary_improvement_cs = 1., @@ -25,37 +26,8 @@ int main(int argc, char *argv[]) { _update_par.SF_ct = 1; _update_par.SF_zf = 1.3; - char tmp[256]; - - logger_setlevel(0, 100); /* log all */ - if (PID != 0) { - logger_disable(); - } else { - sprintf(tmp, ">out_%d", PID); - logger_stdout(tmp); - sprintf(tmp, "err_%d", PID); - freopen(tmp, "w", stderr); - } - - lprintf("MAIN", 0, "PId = %d [world_size: %d]\n\n", PID, WORLD_SIZE); - - read_input(glb_var.read, "test_input"); - read_input(rlx_var.read, "test_input"); - - rlxd_init(rlx_var.rlxd_level, rlx_var.rlxd_seed); - - /* setup communication geometry */ - if (geometry_init() == 1) { - finalize_process(); - return 0; - } - - geometry_mpi_eo(); - init_BCs(&BCs_pars); - u_gauge = alloc_suNg_field(&glattice); - u_gauge_f = alloc_suNf_field(&glattice); s0 = alloc_spinor_field(4, &glattice); s1 = s0 + 1; s2 = s1 + 1; diff --git a/TestProgram/RotatedSF/check_rotatedsf_5.c b/TestProgram/RotatedSF/check_rotatedsf_5.c index 9f419745e..447db4749 100644 --- a/TestProgram/RotatedSF/check_rotatedsf_5.c +++ b/TestProgram/RotatedSF/check_rotatedsf_5.c @@ -9,6 +9,7 @@ rhmc_par _update_par; int main(int argc, char *argv[]) { setup_process(&argc, &argv); + setup_gauge_fields(); char sbuf[350]; BCs_pars_t BCs_pars = { .fermion_twisting_theta = { 0., 0.5, 0.5, 0.5 }, @@ -30,44 +31,14 @@ int main(int argc, char *argv[]) { _update_par.SF_ct = 1; _update_par.SF_zf = 1.3; - logger_setlevel(0, 99); /* log all */ - logger_map("DEBUG", "debug"); - - if (PID != 0) { logger_disable(); } - - if (PID == 0) { - sprintf(sbuf, "out_%d", PID); - logger_stdout(sbuf); - sprintf(sbuf, "err_%d", PID); - freopen(sbuf, "w", stderr); - } - - lprintf("MAIN", 0, "PId = %d [world_size: %d]\n\n", PID, WORLD_SIZE); - - read_input(glb_var.read, "test_input"); - read_input(rlx_var.read, "test_input"); - - rlxd_init(rlx_var.rlxd_level, rlx_var.rlxd_seed); - #if NG != 3 || NF != 3 #error "Can work only with NC=3 and Nf==3" #endif - /* setup communication geometry */ - if (geometry_init() == 1) { - finalize_process(); - return 0; - } - - geometry_mpi_eo(); - init_BCs(&BCs_pars); lprintf("MAIN", 0, "This test implements a comparison with a working code of Stefan Sint\n"); - u_gauge = alloc_suNg_field(&glattice); - u_gauge_f = alloc_suNf_field(&glattice); - lprintf("MAIN", 0, "Reading gauge configuration from file :suNg_field_sint.dat\n"); read_gauge_field_nocheck("suNg_field_sint.dat"); diff --git a/TestProgram/RotatedSF/check_rotatedsf_6.c b/TestProgram/RotatedSF/check_rotatedsf_6.c index 9d32b1032..01804ce45 100644 --- a/TestProgram/RotatedSF/check_rotatedsf_6.c +++ b/TestProgram/RotatedSF/check_rotatedsf_6.c @@ -9,7 +9,7 @@ rhmc_par _update_par; int main(int argc, char *argv[]) { setup_process(&argc, &argv); - char tmp[256]; + setup_gauge_fields(); BCs_pars_t BCs_pars = { .fermion_twisting_theta = { 0., M_PI / 5., M_PI / 5., M_PI / 5. }, @@ -31,42 +31,14 @@ int main(int argc, char *argv[]) { _update_par.SF_ct = 1; _update_par.SF_zf = 1.3; - logger_setlevel(0, 100); /* log all */ - if (PID != 0) { - logger_disable(); - } else { - sprintf(tmp, ">out_%d", PID); - logger_stdout(tmp); - sprintf(tmp, "err_%d", PID); - freopen(tmp, "w", stderr); - } - - lprintf("MAIN", 0, "PId = %d [world_size: %d]\n\n", PID, WORLD_SIZE); - - read_input(glb_var.read, "test_input"); - read_input(rlx_var.read, "test_input"); - - rlxd_init(rlx_var.rlxd_level, rlx_var.rlxd_seed); - #if NG != 3 || NF != 3 #error "Can work only with NC=3 and Nf==3" #endif - /* setup communication geometry */ - if (geometry_init() == 1) { - finalize_process(); - return 0; - } - - geometry_mpi_eo(); - init_BCs(&BCs_pars); lprintf("MAIN", 0, "This test implements a comparison with a working code of Stefan Sint\n"); - u_gauge = alloc_suNg_field(&glattice); - u_gauge_f = alloc_suNf_field(&glattice); - /* unit_u(u_gauge); */ read_gauge_field_nocheck("suNg_field_sint.dat"); apply_BCs_on_fundamental_gauge_field(); diff --git a/TestProgram/RotatedSF/check_rotatedsf_7.c b/TestProgram/RotatedSF/check_rotatedsf_7.c index b5363902f..15634e580 100644 --- a/TestProgram/RotatedSF/check_rotatedsf_7.c +++ b/TestProgram/RotatedSF/check_rotatedsf_7.c @@ -11,31 +11,10 @@ int main(int argc, char *argv[]) { char tmp[256]; setup_process(&argc, &argv); - - logger_setlevel(0, 100); /* log all */ - if (PID != 0) { - logger_disable(); - } else { - sprintf(tmp, ">out_%d", PID); - logger_stdout(tmp); - sprintf(tmp, "err_%d", PID); - freopen(tmp, "w", stderr); - } + setup_gauge_fields(); logger_map("DEBUG", "debug"); - lprintf("MAIN", 0, "PId = %d [world_size: %d]\n\n", PID, WORLD_SIZE); - - read_input(glb_var.read, "test_input"); - - rlxd_init(rlx_var.rlxd_level, rlx_var.rlxd_seed); - - /* setup communication geometry */ - if (geometry_init() == 1) { - finalize_process(); - return 0; - } - BCs_pars_t BCs_pars = { /* .fermion_twisting_theta = {0.,M_PI/5.,M_PI/5.,M_PI/5.}, */ .fermion_twisting_theta = { 0., 0., 0., 0. }, @@ -59,15 +38,10 @@ int main(int argc, char *argv[]) { #error "Can work only with NC=3 and Nf==3" #endif - geometry_mpi_eo(); - init_BCs(&BCs_pars); lprintf("MAIN", 0, "This test implements a comparison with a working code of Stefan Sint\n"); - u_gauge = alloc_suNg_field(&glattice); - u_gauge_f = alloc_suNf_field(&glattice); - unit_u(u_gauge); /* read_gauge_field_nocheck("suNg_field_sint.dat"); */ apply_BCs_on_fundamental_gauge_field(); diff --git a/TestProgram/RotatedSF/check_rotatedsf_8.c b/TestProgram/RotatedSF/check_rotatedsf_8.c index 4d5ea9e0c..c00c44ba4 100644 --- a/TestProgram/RotatedSF/check_rotatedsf_8.c +++ b/TestProgram/RotatedSF/check_rotatedsf_8.c @@ -9,6 +9,7 @@ rhmc_par _update_par; int main(int argc, char *argv[]) { setup_process(&argc, &argv); + setup_gauge_fields(); BCs_pars_t BCs_pars = { /* .fermion_twisting_theta = {0.,M_PI/5.,M_PI/5.,M_PI/5.}, */ @@ -29,42 +30,14 @@ int main(int argc, char *argv[]) { _update_par.SF_ct = 1; _update_par.SF_zf = 1.0; - char tmp[256]; - logger_setlevel(0, 100); /* log all */ - if (PID != 0) { - logger_disable(); - } else { - sprintf(tmp, ">out_%d", PID); - logger_stdout(tmp); - sprintf(tmp, "err_%d", PID); - freopen(tmp, "w", stderr); - } - - lprintf("MAIN", 0, "PId = %d [world_size: %d]\n\n", PID, WORLD_SIZE); - - read_input(glb_var.read, "test_input"); - - rlxd_init(rlx_var.rlxd_level, rlx_var.rlxd_seed); - #if NG != 3 || NF != 3 #error "Can work only with NC=3 and Nf==3" #endif - /* setup communication geometry */ - if (geometry_init() == 1) { - finalize_process(); - return 0; - } - - geometry_mpi_eo(); - init_BCs(&BCs_pars); lprintf("MAIN", 0, "This test implements a comparison with a working code of Stefan Sint\n"); - u_gauge = alloc_suNg_field(&glattice); - u_gauge_f = alloc_suNf_field(&glattice); - unit_u(u_gauge); apply_BCs_on_fundamental_gauge_field(); diff --git a/TestProgram/RotatedSF/check_rotatedsf_9.c b/TestProgram/RotatedSF/check_rotatedsf_9.c index 4dcedc018..97d6fa7df 100644 --- a/TestProgram/RotatedSF/check_rotatedsf_9.c +++ b/TestProgram/RotatedSF/check_rotatedsf_9.c @@ -9,6 +9,7 @@ rhmc_par _update_par; int main(int argc, char *argv[]) { setup_process(&argc, &argv); + setup_gauge_fields(); char sbuf[350]; double a = 1.0; BCs_pars_t BCs_pars = { .fermion_twisting_theta = { 0., a * M_PI / 5., a * M_PI / 5., a * M_PI / 5. }, @@ -29,43 +30,16 @@ int main(int argc, char *argv[]) { _update_par.SF_ct = 1.0; _update_par.SF_zf = 1.0; - logger_setlevel(0, 99); /* log all */ logger_map("DEBUG", "debug"); - if (PID != 0) { logger_disable(); } - - if (PID == 0) { - sprintf(sbuf, ">>out_%d", PID); - logger_stdout(sbuf); - sprintf(sbuf, "err_%d", PID); - freopen(sbuf, "w", stderr); - } - - lprintf("MAIN", 0, "PId = %d [world_size: %d]\n\n", PID, WORLD_SIZE); - - read_input(glb_var.read, "test_input"); - - rlxd_init(rlx_var.rlxd_level, rlx_var.rlxd_seed); - #if NG != 3 || NF != 3 #error "Can work only with NC=3 and Nf==3" #endif - /* setup communication geometry */ - if (geometry_init() == 1) { - finalize_process(); - return 0; - } - - geometry_mpi_eo(); - init_BCs(&BCs_pars); lprintf("MAIN", 0, "This test implements a comparison with a working code of Stefan Sint\n"); - u_gauge = alloc_suNg_field(&glattice); - u_gauge_f = alloc_suNf_field(&glattice); - unit_u(u_gauge); apply_BCs_on_fundamental_gauge_field(); diff --git a/TestProgram/Scattering/check_scattering_length_I0.c b/TestProgram/Scattering/check_scattering_length_I0.c index e4def55e6..4f337acd0 100644 --- a/TestProgram/Scattering/check_scattering_length_I0.c +++ b/TestProgram/Scattering/check_scattering_length_I0.c @@ -459,12 +459,8 @@ int main(int argc, char *argv[]) { error(!(GLB_X == GLB_Y && GLB_X == GLB_Z), 1, "main", "This test works only for GLB_X=GLB_Y=GLB_Z"); setup_process(&argc, &argv); - setup_gauge_fields(); - - read_input(glb_var.read, get_input_filename()); read_input(mes_var.read, get_input_filename()); - read_input(rlx_var.read, get_input_filename()); int numsources = mes_var.nhits; diff --git a/TestProgram/Scattering/check_scattering_length_I2.c b/TestProgram/Scattering/check_scattering_length_I2.c index 3f6dbe4ad..68fb93b42 100644 --- a/TestProgram/Scattering/check_scattering_length_I2.c +++ b/TestProgram/Scattering/check_scattering_length_I2.c @@ -374,9 +374,7 @@ int main(int argc, char *argv[]) { setup_gauge_fields(); - read_input(glb_var.read, get_input_filename()); read_input(mes_var.read, get_input_filename()); - read_input(rlx_var.read, get_input_filename()); int numsources = mes_var.nhits; diff --git a/TestProgram/Scattering/check_scattering_rhopipi.c b/TestProgram/Scattering/check_scattering_rhopipi.c index 05d65e352..39ae3497e 100644 --- a/TestProgram/Scattering/check_scattering_rhopipi.c +++ b/TestProgram/Scattering/check_scattering_rhopipi.c @@ -416,9 +416,7 @@ int main(int argc, char *argv[]) { setup_gauge_fields(); - read_input(glb_var.read, get_input_filename()); read_input(mes_var.read, get_input_filename()); - read_input(rlx_var.read, get_input_filename()); int numsources = mes_var.nhits; diff --git a/TestProgram/Scattering/check_sigma_triangle.c b/TestProgram/Scattering/check_sigma_triangle.c index 0c9a4f6e8..05e8f1e8a 100644 --- a/TestProgram/Scattering/check_sigma_triangle.c +++ b/TestProgram/Scattering/check_sigma_triangle.c @@ -555,9 +555,7 @@ int main(int argc, char *argv[]) { setup_process(&argc, &argv); setup_gauge_fields(); - read_input(glb_var.read, get_input_filename()); read_input(mes_var.read, get_input_filename()); - read_input(rlx_var.read, get_input_filename()); int numsources = mes_var.nhits;