typedef float pooh; #ifndef VMLH #include "conf.h" #include "precision.h" #endif /*!\file \brief Magic numbers. */ /*------------------------------------------------------------------*/ /*! The directions in the lattice*/ /*------------------------------------------------------------------*/ enum DirType {DIR_X, DIR_Y, DIR_Z, DIR_T, DIR_S}; // chiral projectors 1+-g5 enum ChiralProj{PL=-6,PR=-7}; enum PrecType { PREC_SINGLE, PREC_DOUBLE }; /*------------------------------------------------------------------*/ /*! The types of fermion action*/ /* ~~F_CLASS_WILSON_TM added for twisted mass Wilson fermions */ /*------------------------------------------------------------------*/ enum FclassType { F_CLASS_NONE, F_CLASS_STAG, F_CLASS_WILSON, F_CLASS_CLOVER, F_CLASS_DWF, F_CLASS_ASQTAD, F_CLASS_P4, F_CLASS_HISQ, F_CLASS_WILSON_TM, F_CLASS_MDWF, F_CLASS_BFM, // F_CLASS_BFM_TYPE2, F_CLASS_MOBIUS, F_CLASS_DWF4D, F_CLASS_DWF4D_PAIR, F_CLASS_ZMOBIUS, F_CLASS_NAIVE, F_CLASS_GRID, F_CLASS_GRID_GPARITY_MOBIUS, F_CLASS_GRID_MOBIUS, F_CLASS_GRID_ZMOBIUS, F_CLASS_GRID_GPARITY_WILSON_TM, F_CLASS_GRID_WILSON_TM }; /*------------------------------------------------------------------*/ /*! The types of gauge action.*/ /*------------------------------------------------------------------*/ enum GclassType { G_CLASS_NONE, G_CLASS_WILSON, G_CLASS_POWER_PLAQ, G_CLASS_IMPR_RECT, G_CLASS_POWER_RECT, G_CLASS_IMPR_OLSYM, G_CLASS_TADPOLE_RECT }; /*------------------------------------------------------------------*/ /*! The lattice storage orders.*/ /*------------------------------------------------------------------*/ enum StrOrdType { CANONICAL = 0, /*!< Canonical storage order: %Lattice sites are ordered so that the site with cartesian coordinates (x, y, z, t) on a local lattice of dimensions (N_x, N_y, N_z, N_t) has the index n = x + N_x y + N_xN_y z + N_xN_yN_z t i.e the x coordinate runs fastest, then y, then z and t slowest. The gauge link direction index is 0, 1, 2 or 3 for direction X, Y, Z and T respectively. This order is expected by most Lattice functions. */ STAG = 1, /*!< This order is used by actions with staggered fermions. Gauge links are stored in the CANONICAL site order, but they are multiplied by the staggered phases and the hermitian conjugate of each link is stored. Fermions are stored in an odd-even order, where lattice sites are ordered by parity. The parity of a site with cartesian coordinates (x, y, z, t) is even(odd) if x+y+z+t is even(odd) Even parity sites are numbered before odd parity sites. On a local lattice of dimensions (N_x, N_y, N_z, N_t) has the index n = [t + N_t x + N_t N_x y + N_t N_x N_y z + N_t N_x N_y N_z ((x+y+z+t)%2)]/2 i.e the t coordinate runs fastest, then x, then y and z slowest. */ WILSON = 2, /*!< This is used with actions with wilson-type fermions. %Lattice sites for both gauge and fermion fields are ordered by parity. The parity of a site with cartesian coordinates (x, y, z, t) is even(odd) if x+y+z+t is even(odd) For gauge fields, even parity sites are numbered before odd parity sites. On a local lattice of dimensions (N_x, N_y, N_z, N_t) each site has the index n = [x-x%2 + N_x y + N_x N_y z + N_x N_y N_z t + N_x N_y N_z N_t ((x+y+z+t)%2)]/2 For fermion fields, odd parity sites are numbered before even parity sites. Each site has the index n = [x-x%2 + N_x y + N_x N_y z + N_x N_y N_z t + N_x N_y N_z N_t ((x+y+z+t+1)%2)]/2 i.e in both cases, the x coordinate runs fastest, then y, then z and t slowest. */ G_WILSON_HB = 3, /*!< Storage order for the %Wilson gauge action heat bath: Site ordering is CANONICAL but the link direction indices are 0, 1, 2 or 3 for direction T, X, Y and Z respectively. This order is expected by all QncWilsonHb functions.*/ STAG_BLOCK = 4, DWF_5D_EOPREC = WILSON, /*!< DWF with 5 dimensional even/odd precondition, it's actually same as WILSON. */ DWF_4D_EOPREC = 5, /*!< DWF with 4 dimensional even/odd precondition, for odd-odd preconditioning. This should work for Mobius as well. */ DWF_4D_EOPREC_EE = 6, /*!< DWF with 4 dimensional even/odd precondition, for even-even preconditioning. This should work for Mobius as well. */ S_INNER = 7 // s innermost, for Force term, bfm ordering, etc. }; /*------------------------------------------------------------------*/ /*! The storage order conversion flags.*/ /*------------------------------------------------------------------*/ enum CnvFrmType { CNV_FRM_NO = 0, /*!< Do not convert fermion field */ CNV_FRM_YES = 1 /*!< Convert fermion field assuming that it is defined on the whole lattice \e i.e. one spinor per site and not just on the sites of a single checkerboard parity.*/ }; /*------------------------------------------------------------------*/ /*! The dimensionality of the fermion field.*/ /*------------------------------------------------------------------*/ enum FermionFieldDimension { FOUR_D, FIVE_D }; /*------------------------------------------------------------------*/ /*! The kinds of preservation */ /*------------------------------------------------------------------*/ enum PreserveType {PRESERVE_NO = 0, /*!< Do not preserve. */ PRESERVE_YES = 1 /*!< Preserve. */ }; /*------------------------------------------------------------------*/ /*! The kinds of starting configurations.*/ /**/ /* START_CONF_ORD Ordered start. After the configuration is set*/ /* the GJP.StartConfKind() is set to START_CONF_MEM.*/ /**/ /* START_CONF_DISORD Disordered start. After the configuration*/ /* is set the GJP.StartConfKind() is set to START_CONF_MEM.*/ /**/ /* START_CONF_FILE Read from file. After the configuration is set*/ /* the GJP.StartConfKind() is set to START_CONF_MEM.*/ /**/ /* START_CONF_LOAD Memory is not allocated for the gauge field.*/ /* Instead the gauge_field pointer is set to GJP.StartConfLoadAddr().*/ /* After the configuration is set the GJP.StartConfKind() is set */ /* to START_CONF_MEM.*/ /**/ /* START_CONF_MEM Memory is allocated for the gauge field with */ /* pmalloc. When the program begins execution it checcks if the*/ /* pmalloc address is the same as the GJP.StartConfLoadAddr(). */ /* If not it exits with an error. Subsequent calls to the */ /* constructor with START_CONF_MEM do nothing except to set the */ /* gauge field to canonical order.*/ /**/ /* For implementation details see the Lattice constructor in */ /* util/lattice/lattice_base/lattice_base.C*/ /**/ /*------------------------------------------------------------------*/ enum StartConfType { START_CONF_ORD = 0, /*!< Ordered start */ START_CONF_DISORD = 1, /*!< Disordered start */ START_CONF_FILE = 2, /*!< Read from file */ START_CONF_LOAD = 3, /*!< The gauge_field pointer is set to the address specified in the DoArg structure but no memory is allocated for it.*/ START_CONF_MEM = 4 /*!< The gauge_field pointer is set to the address specified in the DoArg structure and memory is allocated for it. */ }; /*------------------------------------------------------------------*/ /*! The possible kinds of initial RNG seeds.*/ /*------------------------------------------------------------------*/ enum StartSeedType { START_SEED_FIXED = 0, /*!< Seed has a fixed value but it is different on each node */ START_SEED_FIXED_UNIFORM = 1, /*!< Seed has a fixed value and is the same on all nodes */ START_SEED = 2, /*!< Seed has a time-dependent value and it is different on each node */ START_SEED_UNIFORM = 3, /*!< Seed has a time-dependent value and is the same on all nodes */ START_SEED_INPUT = 4, /*!< Seed value comes from the DoArg structure and is different on all nodes*/ START_SEED_INPUT_UNIFORM = 5, /*!< Seed value comes from the DoArg structure and is the same on all nodes*/ START_SEED_INPUT_NODE = 6, START_SEED_FILE = 7 /*!< Seed structure comes from the DoArg structure and is different on all nodes. This should be used when in a farming mode. */ }; /*------------------------------------------------------------------*/ /* The possible types of checker boards (chequer boards?)*/ /*------------------------------------------------------------------*/ enum ChkbType { CHKB_EVEN =0, /* Even checkerboard*/ CHKB_ODD =1}; /* Odd checkerboard*/ /*------------------------------------------------------------------*/ /*! Hermitian conjugate flag*/ /*------------------------------------------------------------------*/ enum DagType { DAG_NO = 0, /*!< No hermitian conjugate */ DAG_YES = 1 /*!< Hermitian conjugate */ }; /*------------------------------------------------------------------*/ /*! The kinds of boundary conditions.*/ /*------------------------------------------------------------------*/ enum BndCndType { BND_CND_PRD, /*!< Periodic */ BND_CND_APRD, /*!< Antiperiodic */ BND_CND_TWISTED, /*!< Twisted BCs */ BND_CND_GPARITY, /*!< G-parity */ BND_CND_GPARITY_TWISTED /*!< G-parity + Twisted BCs*/ }; /*------------------------------------------------------------------*/ /*! The possible kinds of gauge fixing*/ /*------------------------------------------------------------------*/ enum FixGaugeType{ FIX_GAUGE_NONE = -2, /*!< No gauge fixing */ FIX_GAUGE_LANDAU = -1, /*!< Fixing to Landau gauge */ FIX_GAUGE_COULOMB_X = 0, /*!< Fixing to Coulomb gauge */ FIX_GAUGE_COULOMB_Y = 1, /*!< Fixing to Coulomb gauge */ FIX_GAUGE_COULOMB_Z = 2, /*!< Fixing to Coulomb gauge */ FIX_GAUGE_COULOMB_T = 3 /*!< Fixing to Coulomb gauge */ }; /*------------------------------------------------------------------*/ /*! The possible kinds of spin projections*/ /*------------------------------------------------------------------*/ enum SprojType { SPROJ_XM = 0, /*!< (1 - gamma_0) projection */ SPROJ_YM = 1, /*!< (1 - gamma_1) projection */ SPROJ_ZM = 2, /*!< (1 - gamma_2) projection */ SPROJ_TM = 3, /*!< (1 - gamma_3) projection */ SPROJ_XP = 4, /*!< (1 + gamma_0) projection */ SPROJ_YP = 5, /*!< (1 + gamma_1) projection */ SPROJ_ZP = 6, /*!<(1 + gamma_2) projection */ SPROJ_TP = 7 /*!< (1 + gamma_3) projection */ }; /*------------------------------------------------------------------*/ /*! The possible kinds of Sigma projections*/ /*------------------------------------------------------------------*/ enum SigmaprojType { SIGMAPROJ_XY = 0, /*!< Sigma_{0,1} projection */ SIGMAPROJ_XZ = 1, /*!< Sigma_{0,2} projection */ SIGMAPROJ_XT = 2, /*!< Sigma_{0,3} projection */ SIGMAPROJ_YZ = 3, /*!< Sigma_{1,2} projection */ SIGMAPROJ_YT = 4, /*!< Sigma_{1,3} projection */ SIGMAPROJ_YX = 5, /*!< Sigma_{1,0} projection */ SIGMAPROJ_ZT = 6, /*!< Sigma_{2,3} projection */ SIGMAPROJ_ZX = 7, /*!< Sigma_{2,0} projection */ SIGMAPROJ_ZY = 8, /*!< Sigma_{2,1} projection */ SIGMAPROJ_TX = 9, /*!< Sigma_{3,0} projection */ SIGMAPROJ_TY =10, /*!< Sigma_{3,1} projection */ SIGMAPROJ_TZ =11 /*!< Sigma_{3,2} projection */ }; /*------------------------------------------------------------------*/ /*! Operators, in terms of the fermion matrix \e M, of which the eigenvalues/vectors can be found.*/ /*------------------------------------------------------------------*/ enum RitzMatType { NONE, /*!< No eigenvalues requested */ MAT_HERM, /*!< The hermitian matrix on the full lattice */ MATPC_HERM, /*!< The preconditioned hermitian matrix on a single parity. For Wilson case, it's the gamma_5 D_oo. For DWF, it's the gamma_5 R D4_eo, where D4_ee is the 4dim even/odd preconditioned matrix */ MATPCDAG_MATPC, /*!< The preconditioned \f$M^\dagger M\f$ on a single parity */ NEG_MATPCDAG_MATPC, /*!< The preconditioned \f$-M^\dagger M\f$ on a single parity */ MATDAG_MAT, /*!< \f$M^\dagger M\f$ on the full lattice */ NEG_MATDAG_MAT, /*!< \f$-M^\dagger M\f$ on the full lattice*/ MATDAG_MAT_NORM, /*!< \f$cM^\dagger M\f$ on the full lattice (normalised)*/ NEG_MATDAG_MAT_NORM, /*!< \f$-cM^\dagger M\f$ on the full lattice (normalised)*/ MATPCDAG_MATPC_SHIFT, /*!< The preconditioned \f$M^\dagger M\f$ on a single parity, Shift for eigen spectrum will be applied, currently DWF only. When the original (unshifted) opertor is \f$M^\dagger M = H^2, H = \Gamma_5 M$\f, the shifted operator will be \f (H-\mu)(H-\mu)\f, here \f$M$\f is the 4d-even/odd precondition operator. For Wilson types, it's easy to extend, I don't know for staggered.*/ RitzMatType_LAST /* sentinel */ }; /*------------------------------------------------------------------*/ /*! The types of rational approximation for RHMC*/ /*------------------------------------------------------------------*/ enum RatApproxType { CONSTANT, /*!< CONSTANT - approximation is held constant*/ DYNAMIC /*!< DYNAMIC - approximation recalculated at the end of the trajectory*/ }; /*------------------------------------------------------------------*/ /*! The types multishift solve to perform*/ /*------------------------------------------------------------------*/ enum MultiShiftSolveType { SINGLE, /*!< SINGLE - solutions are summed to a common vector*/ MULTI, /*!< MULTI - solutions are summed to separate vectors*/ GENERAL /*!< GENERAL - non-zero initial guess (R2 aglorithm) */ }; enum MassRenormaliseDir { RENORM_BACKWARDS = 0, RENORM_FORWARDS = 1 }; /*------------------------------------------------------------------*/ /*! The field type (valid for dwf rhmc)*/ /*------------------------------------------------------------------*/ enum FieldType { FERMION, /*< FERMION - field is of fermion type*/ BOSON /*< BOSON - field is of boson type*/ }; /*------------------------------------------------------------------*/ /*! The type rhmc being done*/ /*------------------------------------------------------------------*/ enum RatType { RATIONAL_STANDARD, /*< Standard rhmc, single timescale and single approx */ RATIONAL_QUOTIENT, /*< Single timescale with double fermion/boson approximation */ RATIONAL_SPLIT /* Multiple timescale with single approx */ }; enum WbaryonFold {BARYON_FOLD, BARYON_RAW, BARYON_PAST}; enum SourceKind { POINT_W = 0, WALL_W, BOX_W, JACOBI_W, /*added by Thomas and Xiaodong*/ MAX_NUM_SINK, /*- N/A yet */ Z2, /* added by Meifeng Lin */ COMPLEX_Z2, KURAMASHI }; enum SinkKind { W_POINT, W_WALL, W_BOX }; enum MomentumKind {MOM_000 = 0, /* average over permutaions*/ MOM_001, MOM_002, MOM_011, MOM_022, MOM_111, MOM_222, MAX_NUM_MOMENTA }; /*-----------------------------------------------------------------------------*/ /* Related to extended mesons constructed from derivative operators*/ /*-----------------------------------------------------------------------------*/ /*Note: The actual derivative direction depends on propagation direction*/ /* See the end of w_quark.C for detail. */ /* DEV1DEV2 means DEV2 first DEV1 second !!*/ /*The order is very important, changes will affect WspectExtendedMesons*/ /*DEVIDEVI not used now*/ enum DEVOperatorKind{ /*Group 0: Individual Operators*/ UNIT=0, DEV1, DEV2, DEV3, DEV1DEV2, DEV2DEV1, DEV2DEV3, DEV3DEV2, DEV1DEV3, DEV3DEV1, DEV1DEV1, DEV2DEV2, DEV3DEV3, DEV_OP_NUM, /*Naming convetion: F--first deriv S--second deriv(followed by ANTISYM/SYM/DIAG)*/ /*Group 1: Sum over polarizations*/ SUM_F, /*sum of first derivatives*/ SUM_S_ANTISYM, /*sum of second derivatives(antisymmtric combination: DIDJ-DJDI)*/ SUM_S_SYM, /*sum of second derivatives(symmtric combination: DIDJ+DJDI)*/ SUM_S_DIAG, /*sum of second derivatives(diagonal combination: DIDI)*/ /*Group 2: */ SUM_F_S_ANTISYM, /*SUM_F+SUM_S_ANTISYM*/ /*SUM_S_SYM */ /*SUM_S_DIAG*/ /*Group 3:*/ /*SUM_F_S_ANTISYM*/ SUM_S_SYM_DIAG, /*Group 4:*/ SUM_UNIT_F_S_ANTISYM, /*SUM_S_SYM_DIAG*/ END_SUM_OP, /*Related to operators with local B/E fields*/ /*group 1*/ BEGIN_BE_OP, FB1_OP=0,FB2_OP,FB3_OP, FE1_OP,FE2_OP,FE3_OP, FUNIT_OP, /*for mixing*/ /*SUM operators*/ /*group 2*/ SUM_MAGN_OP, SUM_ELEC_OP, /*group 3*/ SUM_MAGN_ELEC_OP, /**/ END_BE_OP }; enum WMesonOpKind{ /*each name starts with MO(Meson Operator)*/ /*Normal mesons*/ /* MO_a0, MO_a0_prime, MO_a1_x, MO_a1_y, MO_a1_z, MO_b1_x, MO_b1_y, MO_b1_z, MO_rho_x, MO_rho_y, MO_rho_z, MO_rho_prime_x, MO_rho_prime_y, MO_rho_prime_z, MO_pion, MO_pion_prime, */ /* extended mesons*/ MO_a0xP_x, MO_a0xP_y, MO_a0xP_z, MO_pionxP_x, MO_pionxP_y, MO_pionxP_z, MO_a0_primexP_x, MO_a0_primexP_y, MO_a0_primexP_z, MO_rhoxP_A1, MO_rhoxP_T1_x, MO_rhoxP_T1_y, MO_rhoxP_T1_z, MO_rhoxP_T2_x, MO_rhoxP_T2_y, MO_rhoxP_T2_z, MO_a1xP_A1, MO_a1xP_T2_x, MO_a1xP_T2_y, MO_a1xP_T2_z, MO_a1xP_E_1,MO_a1xP_E_2, MO_b1xP_T1_x, MO_b1xP_T1_y, MO_b1xP_T1_z, MO_b1xD_A2, MO_b1xD_T1_x, MO_b1xD_T1_y, MO_b1xD_T1_z, MO_b1xD_T2_x, MO_b1xD_T2_y, MO_b1xD_T2_z, MO_b1xD_E_1, MO_b1xD_E_2, MO_a0_primexD_x, MO_a0_primexD_y, MO_a0_primexD_z, MO_rhoxB_T1_x, MO_rhoxB_T1_y, MO_rhoxB_T1_z, MO_rhoxB_T2_x, MO_rhoxB_T2_y, MO_rhoxB_T2_z, MO_a1xB_A1, MO_a1xB_T1_x, MO_a1xB_T1_y, MO_a1xB_T1_z, MO_a1xB_T2_x, MO_a1xB_T2_y, MO_a1xB_T2_z, MO_a1xD_A2, MO_a1xD_T1_x, MO_a1xD_T1_y, MO_a1xD_T1_z, MO_a1xD_T2_x, MO_a1xD_T2_y, MO_a1xD_T2_z, MO_a1xD_E_1, MO_a1xD_E_2, MO_rhoxD_A2, MO_rhoxD_T1_x,MO_rhoxD_T1_y,MO_rhoxD_T1_z, MO_rhoxD_T2_x,MO_rhoxD_T2_y,MO_rhoxD_T2_z, /*added*/ MO_pionxB_T1_x, MO_pionxB_T1_y, MO_pionxB_T1_z, MO_pionxD_T2_x, MO_pionxD_T2_y, MO_pionxD_T2_z, NUM_WMESON_OP_KIND }; /*--------------------------------*/ /*enum WMesonState (with polarization)*/ /*--------------------------------*/ enum WMesonState{ /*Normal Mesons*/ /* MS_a0, MS_a0_prime, MS_a1_x, MS_a1_y, MS_a1_z, MS_b1_x, MS_b1_y, MS_b1_z, MS_rho_x, MS_rho_y, MS_rho_z, MS_rho_prime_x, MS_rho_prime_y, MS_rho_prime_z, MS_pion, MS_pion_prime, */ /*ExtendedMesons*/ MS_a0xP_x, MS_a0xP_y, MS_a0xP_z, MS_pionxP_x, MS_pionxP_y, MS_pionxP_z, MS_a0_primexP_x, MS_a0_primexP_y, MS_a0_primexP_z, MS_rhoxP_A1_1, MS_rhoxP_T1_x, MS_rhoxP_T1_y, MS_rhoxP_T1_z, MS_rhoxP_T2_x, MS_rhoxP_T2_y, MS_rhoxP_T2_z, MS_a1xP_A1_1, MS_a1xP_T2_x, MS_a1xP_T2_y, MS_a1xP_T2_z, MS_a1xP_E_1, MS_a1xP_E_2, MS_b1xP_T1_x, MS_b1xP_T1_y, MS_b1xP_T1_z, MS_b1xD_A2_1, MS_b1xD_T1_x, MS_b1xD_T1_y, MS_b1xD_T1_z, MS_b1xD_T2_x, MS_b1xD_T2_y, MS_b1xD_T2_z, MS_b1xD_E_1, MS_b1xD_E_2, MS_a0_primexD_x, MS_a0_primexD_y, MS_a0_primexD_z, MS_rhoxB_T1_x, MS_rhoxB_T1_y, MS_rhoxB_T1_z, MS_rhoxB_T2_x, MS_rhoxB_T2_y, MS_rhoxB_T2_z, MS_a1xB_A1_1, MS_a1xB_T1_x, MS_a1xB_T1_y, MS_a1xB_T1_z, MS_a1xB_T2_x, MS_a1xB_T2_y, MS_a1xB_T2_z, MS_a1xD_A2_1, MS_a1xD_T1_x,MS_a1xD_T1_y,MS_a1xD_T1_z, MS_a1xD_T2_x,MS_a1xD_T2_y,MS_a1xD_T2_z, MS_a1xD_E_1, MS_a1xD_E_2, MS_rhoxD_A2_1, MS_rhoxD_T1_x,MS_rhoxD_T1_y,MS_rhoxD_T1_z, MS_rhoxD_T2_x,MS_rhoxD_T2_y,MS_rhoxD_T2_z, MS_pionxB_T1_x,MS_pionxB_T1_y,MS_pionxB_T1_z, MS_pionxD_T2_x, MS_pionxD_T2_y, MS_pionxD_T2_z, /* Mixing states, syntax: source_mix_sink_polarisation*/ /* add desired combinations from mixing_terms.h here*/ NUM_WMESON_STATE /*should equal number of states in WGinfo.h?*/ }; /*the average of all polarizations, written out to file*/ /* the first EXTMESONS enum's have to be in the same order as in WspectOutput*/ #define MAX_FUZZING_C_NUM 10 enum WMesonOutputName{ /*actual outputs(average over polarizarion)*/ /* a0, a0_prime, a1, b1, rho, rho_prime, pion, pion_prime, */ /*Extended Meson*/ a0xP, pionxP, a0_primexP, rhoxP_A1, rhoxP_T1, rhoxP_T2, a1xP_A1, a1xP_T2, a1xP_E, b1xP_T1, b1xD_A2, b1xD_T1, b1xD_T2, b1xD_E, a0_primexD, rhoxB_T1, rhoxB_T2, a1xB_A1, a1xB_T1, a1xB_T2, a1xD_A2, a1xD_T1, a1xD_T2, a1xD_E, rhoxD_A2, rhoxD_T1, rhoxD_T2, pionxB_T1, pionxD_T2, /* Mixing states*/ /* add desired combinations from mixing_terms.h here*/ NUM_WMESON_OUTPUT /* should be here when done with all 15 mesons*/ /*all polarizations*/ }; enum WMesonCategory{ NORMALMESON, EXT_FIRSTDEV_MESON, EXT_SECONDDEV_SYM_MESON, EXT_SECONDDEV_ANTISYM_MESON, EXT_SECONDDEV_DIAG_MESON, MIXING }; /*------------------------------------------------------------------------*/ /* Related to extended mesons constructed from local fields*/ /*------------------------------------------------------------------------*/ /*output index(average over all polarisations)*/ enum WExtMesonBEOutputName{ BE_pionxB=0, BE_rhoxB_T1, NUM_WEXTMESON_BE_OUTPUT }; /*state IDs(with polarisation)*/ enum WExtMesonBEState{ BE_MS_pionxB_x=0, BE_MS_pionxB_y, BE_MS_pionxB_z, BE_MS_rhoxB_T1_x, BE_MS_rhoxB_T1_y, BE_MS_rhoxB_T1_z, NUM_WEXTMESON_BE_STATES }; /*meson operator*/ enum WExtMesonBEOp{ BE_MO_pionxB_x=0,BE_MO_pionxB_y, BE_MO_pionxB_z, BE_MO_rhoxB_T1_x, BE_MO_rhoxB_T1_y, BE_MO_rhoxB_T1_z, NUM_WEXTMESON_BE_OPS }; /*extmeson category*/ enum WExtMesonBECategory{ ELEC_HYBRID_BE=0, MAG_HYBRID_BE, MIXING_BE }; enum FieldTensorId{ /*used in state table*/ /*group 1*/ FB1=0,FB2,FB3, FE1,FE2,FE3, NUM_FLDS, FUNIT, /*for mixing*/ /*SUM operators*/ /*group 2*/ SUM_MAGN, SUM_ELEC, /*group 3*/ SUM_MAGN_ELEC, /**/ NUM_FLD_OPS }; /*! How to obtain the masses at which the observable is measured.*/ enum PatternType { LIN=0, /*!< Masses are elements of an arithmetic progression. */ ARRAY, /*!< Masses are defined in a list. */ LOG, /*!< Masses are elements of a geometric progression. */ FLOW /*!< Used for the spetral flow */ }; /*! Different types of integrators and actions. */ enum IntegratorType { INT_LEAP=0, INT_OMELYAN, INT_CAMPOSTRINI, INT_OMELYAN_44, INT_OMELYAN_45, INT_FORCE_GRAD_PQPQP, INT_FORCE_GRAD_QPQPQ, INT_FORCE_GRAD_QPQPQPQ, INT_FORCE_GRAD_PQPQPQPQP, INT_SUM, INT_MOM, INT_GAUGE, INT_FERMION, INT_BOSON, INT_QUOTIENT, INT_RATIONAL, INT_RATIONAL_SPLIT, INT_RATIONAL_QUOTIENT, INT_EOFA }; /*! Define the level of the integrator (used to determine whether a CSM.comment is necesssary. */ enum IntegratorLevel { EMBEDDED_INTEGRATOR=0, TOP_LEVEL_INTEGRATOR }; /* Moved these from hmd_arg.x to allow use in other classes. */ enum ReunitarizeType { REUNITARIZE_NO = 0, REUNITARIZE_YES = 1 }; /* Perform a reproducibility test? */ enum ReproduceTest { REPRODUCE_NO = 0, REPRODUCE_YES = 1 }; /* Perform a reproducibility test? */ enum TestReproduceTest { TEST_REPRODUCE_NO = 0, TEST_REPRODUCE_YES = 1 }; /* Perform a reversibility test? */ enum ReverseTest { REVERSE_NO = 0, REVERSE_YES = 1 }; /* Perform a Metropolis accept/reject? */ enum MetropolisType { METROPOLIS_NO = 0, METROPOLIS_YES = 1 }; /* Measure the magnitude of the force? */ enum ForceMeasure { FORCE_MEASURE_NO = 0, FORCE_MEASURE_YES = 1 }; /* Measure the eigenvalue bounds? */ enum EigenMeasure { EIGEN_MEASURE_NO = 0, EIGEN_MEASURE_YES = 1 }; enum RhmcPolesAction { RHMC_POLES_CALC = 0, RHMC_POLES_READ = 1, RHMC_POLES_CALC_WRITE = 2 }; enum HmdLimits { MAX_HMD_MASSES=200 , /* The maximum number of dynamical masses.*/ MAX_RAT_DEGREE=30 /* The maximum degree of the rational approximation.*/ }; /* Which type of solver to use? */ enum InverterType { CG = 0, /* Conjugate Gradients. */ BICGSTAB = 1, /* BiCGstab(n). */ EIGCG = 2, /* EigCG */ LOWMODEAPPROX = 3, /* Low Mode Approximation */ CG_LOWMODE_DEFL = 4, /* CG accelerated using low-mode deflation */ CG_FIXED_ITER = 5, /* CG accelerated using low-mode deflation */ HDCG = 6, /* HDCG */ FAKE = 7 /* For testing */ }; /* CG_FIXED_ITER = 6, */ /* Which type of approximation to use? */ enum RationalApproxType { RATIONAL_APPROX_POWER = 0, RATIONAL_APPROX_QUOTIENT = 1, RATIONAL_APPROX_ZERO_POLE = 2 }; /* Are set automatically (staggered only) or not? */ enum RationalBoundsType { RATIONAL_BOUNDS_AUTOMATIC = 0, RATIONAL_BOUNDS_MANUAL = 1 }; /* for AlgStaticB & QPropS by YA */ enum StaticBActionLinkSmearType { SB_ALS_NONE, SB_ALS_APE, SB_ALS_APE_NO_PROJ, SB_ALS_APE_OLEG, SB_ALS_HYP_HK, SB_ALS_HYP_L, SB_ALS_HYP_2, SB_ALS_STOUT }; /* for AlgStaticB & QPropW by YA */ // to apply links in Gaussian smearing for src and sink enum GaussianKernelLinkSmearType { GKLS_NONE = 0, GKLS_APE = 1, GKLS_STOUT = 2 }; /*! IO type for the quark propagator calculations in AlgNuc3pt. */ enum CalcQpropType { READ_QPROP, // load quark propagator NOIO_QPROP, // calculate quark propagator without IO WRITE_QPROP // calculate quark propagator and write it to disk }; enum CalcSeqType { READ_SEQ, // load sequential propagator NOIO_SEQ, // calculate sequential propagator for each quark propagator without IO WRITE_SEQ, // calculate sequential propagator and write it to disk MULT_SEQ, // calculate multiple sequential propagators without IO READ_MULT_SEQ, // read a sequential propagators calculated using multiple-sink method WRITE_MULT_SEQ // calculate multiple sequential propagators and write to disk } ; /* we have to add a prefix (here I choose it to be BFM_) to all names to avoid naming collision. */ /* keep this list in sync with bfm releases. */ enum BfmSolverType { BFM_DWF, BFM_DWFrb4d, BFM_WilsonFermion, BFM_WilsonTM, BFM_WilsonNN, BFM_HwPartFracZolo, BFM_HwContFracZolo, BFM_HwPartFracTanh, BFM_HwContFracTanh, BFM_HwCayleyZolo, BFM_HtCayleyZolo, BFM_HwCayleyTanh, BFM_HmCayleyTanh, BFM_HtCayleyTanh, BFM_DWFTransfer, BFM_DWFTransferInv, BFM_HtContFracTanh, BFM_HtContFracZolo };