Exemple #1
0
int main(int argc, char *argv[])
{ 

  Meas.TaskList.TaskList_len = 1;
  Meas.TaskList.TaskList_val = &Task;

  Meas.GaugeStem = "";
  Meas.RNGStem = "";

  Task.OutputFilestem = "";
  Task.ArgFilename = "arg_stem";
  ws_out.cg = "";
  ws_out.cg2 = "filename";
  ws_out.pbp = "filename";
  ws_out.mid_point = "filename";
  ws_out.nucleon = "filename";
  ws_out.nucleon_prime = "filename";
  ws_out.delta_x = "filename";
  ws_out.delta_y = "filename";
  ws_out.delta_z = "filename";
  ws_out.delta_t = "filename";
  evo_arg.ensemble_id = "id";
  evo_arg.ensemble_label = "label";
  evo_arg.creator = "creator";
  evo_arg.gauge_file_stem = "ckpoint_lat";
  evo_arg.rng_file_stem = "ckpoint_rng";
  evo_arg.plaquette_stem = "plaq";
  evo_arg.evo_stem = "hmd";
  evo_arg.work_directory = "";

  ws_arg.CgArgFile = "CgFile";

  rat_arg.resize(4);
  rat_arg.resize(0, 10, 16);
  rat_arg.resize(1, 10, 16);
  rat_arg.resize(2, 9, 14);
  rat_arg.resize(3, 6, 9);
  bsn_arg.resize(1);
  frm_arg.resize(1);
  split_arg.resize(2);
  quo_arg.resize(1);
  rat_quo_arg.resize(1);
  rat_quo_arg.resize(0, 9, 15, 6, 9);
  float_array.resize(4);
  param_array.resize(3);
  param_array.params.params_val[0].name="mass";
  param_array.params.params_val[1].name="epsilon";
  param_array.params.params_val[2].name="u0";
  stag_quark_src.type = S_QUARK_POINT;
  for(int i = 0; i < 4; i++)
    {
      stag_quark_src.origin[i] = 0;
      stag_quark_src.end[i] = 0;
    }
  stag_quark_src.dir = HDM_T;
  stag_quark_arg.qid = 0;
  stag_quark_arg.cg.mass = 0.1;
  stag_quark_arg.cg.max_num_iter = 5000;
  stag_quark_arg.cg.stop_rsd = 1e-10;
  stag_quark_arg.cg.true_rsd = 1e-10;
  stag_quark_arg.cg.RitzMatOper = MATPCDAG_MATPC;
  stag_quark_arg.cg.Inverter = CG;
  stag_quark_arg.cg.bicgstab_n = 1;
  stag_quark_arg.src = stag_quark_src;
  stag_quark_arg.sln = LOCAL;
  stag_meson_arg.qid0 = 0;
  stag_meson_arg.qid1 = 1;
  stag_meson_arg.dir = HDM_T;
  stag_meson_arg.meson_buf = 0;
  stag_mom_meson_arg.qid0 = 0;
  stag_mom_meson_arg.qid1 = 1;
  stag_mom_meson_arg.dir = HDM_T;
  stag_mom_meson_arg.no_of_momenta = 4;
  stag_mom_meson_arg.meson_buf = 0;
  stag_nucleon_arg.qid0 = 0;
  stag_nucleon_arg.qid1 = 1;
  stag_nucleon_arg.qid2 = 2;
  stag_nucleon_arg.dir = HDM_T;
  stag_nucleon_arg.nucleon_buf = 0;
  stag_non_local_arg.qid0 = 0;
  stag_non_local_arg.qid1 = 1;
  stag_non_local_arg.qid2 = 2;
  stag_non_local_arg.dir = HDM_T;
  stag_non_local_arg.nlocal_buf = 0;
  for(int i = 0; i < 8; i++)
    nlstag_meson_arg.qid0[i] = i;
  nlstag_meson_arg.dir = HDM_T;
  nlstag_meson_arg.nlocal_buf = 0;
  nuc3pt_arg.cname="AlgNuc3pt";
  nuc3pt_arg.ensemble_label="AlgNuc3pt";

  mdwf_arg.M5 = 1.8;
 // mdwf_arg.cg_arg_p = &cg_arg;

  cg_arg.Encode("cg_arg.vml","cg_arg");
  mdwf_arg.Encode("mdwf_arg.vml","mdwf_arg");
  hmd_arg.Encode("hmd_arg.vml","hmd_arg");
  evo_arg.Encode("evo_arg.vml","evo_arg");
  do_arg.Encode("do_arg.vml","do_arg");
  ws_arg.Encode("w_spect_arg.vml","w_spect_arg");
  ws_out.Encode("w_spect_output.vml","w_spect_output");
  Meas.Encode("meas_arg.vml","meas_arg");
  pa.Encode("pot_arg.vml","pot_arg");
  pbp.Encode("pbp_arg.vml","pbp_arg");
  fga.Encode("fix_gauge_arg.vml","fix_gauge_arg");
  hmc_arg.Encode("hmc_arg.vml","hmc_arg");
  quo_arg.Encode("quo_arg.vml","quo_arg");
  frm_arg.Encode("frm_arg.vml","frm_arg");
  bsn_arg.Encode("bsn_arg.vml","bsn_arg");
  rat_arg.Encode("rat_arg.vml","rat_arg");
  rat_quo_arg.Encode("rat_quo_arg.vml","rat_quo_arg");
  split_arg.Encode("split_arg.vml","split_arg");
  gauge_arg.Encode("gauge_arg.vml","gauge_arg");
  ab_arg.Encode("ab_arg.vml","ab_arg");
  eig_arg.Encode("eig_arg.vml","eig_arg");
  float_array.Encode("float_array.vml","float_array");
  param_array.Encode("param_array.vml","param_array");
  stag_quark_src.Encode("stag_quark_src.vml","stag_quark_src");
  stag_quark_arg.Encode("stag_quark_arg.vml","stag_quark_arg");
  stag_meson_arg.Encode("stag_meson_arg.vml","stag_meson_arg");
  stag_mom_meson_arg.Encode("stag_mom_meson_arg.vml","stag_mom_meson_arg");
  stag_nucleon_arg.Encode("stag_nucleon_arg.vml","stag_nucleon_arg");
  stag_non_local_arg.Encode("stag_non_local_arg.vml","stag_non_local_arg");
  nlstag_meson_arg.Encode("nlstag_meson_arg.vml","nlstag_meson_arg");
  qpropw_arg.Encode("qpropw_arg.vml","qpropw_arg");
  eigcg_arg.Encode("eigcg_arg.vml","eigcg_arg");
  nuc3pt_arg.Encode("nuc3pt_arg.vml","nuc3pt_arg");

  vml_markup_type(VML_XML);
  
  cg_arg.Encode("cg_arg.xml","cg_arg");
  hmd_arg.Encode("hmd_arg.xml","hmd_arg");
  evo_arg.Encode("evo_arg.xml","evo_arg");
  do_arg.Encode("do_arg.xml","do_arg");
  ws_arg.Encode("w_spect_arg.xml","w_spect_arg");
  ws_out.Encode("w_spect_output.xml","w_spect_output");
  hmc_arg.Encode("hmc_arg.xml","hmc_arg");
  quo_arg.Encode("quo_arg.xml","quo_arg");
  frm_arg.Encode("frm_arg.xml","frm_arg");
  bsn_arg.Encode("bsn_arg.xml","bsn_arg");
  rat_arg.Encode("rat_arg.xml","rat_arg");
  rat_quo_arg.Encode("rat_quo_arg.xml","rat_quo_arg");
  split_arg.Encode("split_arg.xml","split_arg");
  gauge_arg.Encode("gauge_arg.xml","gauge_arg");
  ab_arg.Encode("ab_arg.xml","ab_arg");
  eig_arg.Encode("eig_arg.xml","eig_arg");
  stag_quark_src.Encode("stag_quark_src.xml","stag_quark_src");
  stag_quark_arg.Encode("stag_quark_arg.xml","stag_quark_arg");
  stag_meson_arg.Encode("stag_meson_arg.xml","stag_meson_arg");
  stag_mom_meson_arg.Encode("stag_mom_meson_arg.xml","stag_mom_meson_arg");
  stag_nucleon_arg.Encode("stag_nucleon_arg.xml","stag_nucleon_arg");
  stag_non_local_arg.Encode("stag_non_local_arg.xml","stag_non_local_arg");
  nlstag_meson_arg.Encode("nlstag_meson_arg.xml","nlstag_meson_arg");
  qpropw_arg.Encode("qpropw_arg.xml","qpropw_arg");
  eigcg_arg.Encode("eigcg_arg.xml","eigcg_arg");
  nuc3pt_arg.Encode("nuc3pt_arg.xml","nuc3pt_arg");

  return(0);
}