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); }