if(g_eps_sq_force1 < 0) g_eps_sq_force1 = g_eps_sq_force; if(g_eps_sq_force2 < 0) g_eps_sq_force2 = g_eps_sq_force; if(g_eps_sq_force3 < 0) g_eps_sq_force3 = g_eps_sq_force; if(g_eps_sq_acc1 < 0) g_eps_sq_acc1 = g_eps_sq_acc; if(g_eps_sq_acc2 < 0) g_eps_sq_acc2 = g_eps_sq_acc; if(g_eps_sq_acc3 < 0) g_eps_sq_acc3 = g_eps_sq_acc;
int read_input(char * conf_file){  /********************************************   * Setting default values!   ********************************************/  T_global = _default_T_global;  L = _default_L;  LX = _default_LX;  LY = _default_LY;  LZ = _default_LZ;  N_PROC_X = _default_N_PROC_X;  g_kappa = _default_g_kappa;  g_mu = _default_g_mu;  g_mu1 = _default_g_mu1;  g_mu2 = _default_g_mu2;  g_mu3 = _default_g_mu3;  g_beta = _default_g_beta;  g_c_sw = _default_g_c_sw;  q_off = _default_q_off;  q_off2 = _default_q_off2;  dtau = _default_dtau;  tau = _default_tau;  Nsteps = _default_Nsteps;  nsmall = _default_nsmall;  integtyp = _default_integtyp;  random_seed = _default_random_seed;  matrix_element_flag = _default_matrix_element_flag;  solver_flag = _default_solver_flag;  operator_flag = _default_operator_flag;  startoption = _default_startoption;  Ntherm = _default_Ntherm;  Nmeas = _default_Nmeas;  Nskip = _default_Nskip;  save_config_flag = _default_save_config_flag;  save_prop_flag = _default_save_prop_flag;  save_prop_g2_flag = _default_save_prop_g2_flag;  write_cp_flag = _default_write_cp_flag;  cp_interval = _default_cp_interval;  nstore = _default_nstore;  strcpy(rlxd_input_filename, _default_rlxd_input_filename);  strcpy(gauge_input_filename, _default_gauge_input_filename);  g_stdio_proc = _default_g_stdio_proc;  index_start = _default_index_start;  index_end = _default_index_end;  first_prop_flag = _default_first_prop_flag;  ITER_MAX_CG = _default_ITER_MAX_CG;  ITER_MAX_BCG = _default_ITER_MAX_BCG;  X0 = _default_X0;  max_solver_iterations = _default_max_solver_iterations;  solver_precision = _default_solver_precision;  mass_number = _default_mass_number;  g_rgi_C1 = _default_g_rgi_C1;  read_source_flag= _default_read_source_flag;  strcpy(source_input_filename, _default_source_filename);  g_eps_sq_force = _default_g_eps_sq_force;  g_eps_sq_acc = _default_g_eps_sq_acc;  g_relative_precision_flag = _default_g_relative_precision_flag;  return_check_flag = _default_return_check_flag;  return_check_interval = _default_return_check_interval;  /********************************************/  if ((yyin = fopen(conf_file, "rt")) == NULL){    return(2);  }  yyout = fopen("/dev/null", "w");  parse_config();    if(LX == 0) {    LX = L;  }  if(LY == 0) {    LY = L;  }  if(LZ == 0) {    LZ = L;  }  g_rgi_C0 = 1. - 8.*g_rgi_C1;  g_ka_csw_8 = g_kappa*g_c_sw/8.;    fclose(yyout);  fclose(yyin);  return(0);}