11 multi1d<int> boundary;
12 multi1d<int> rng_seed;
29 XMLReader
reader(filename);
57 push(writer,
"params");
58 push(writer,
"lattice");
82 int main(
int argc,
char **argv)
94 QDPIO::cerr <<
"Caught exception " <<
s << std::endl;
100 Layout::setLattSize(
params.nrow);
105 push(xml_out,
"ritzTest");
110 typedef LatticeFermion
T;
111 typedef multi1d<LatticeColorMatrix>
P;
112 typedef multi1d<LatticeColorMatrix>
Q;
118 multi1d<LatticeColorMatrix>
u(
Nd);
119 XMLReader gauge_file_xml, gauge_xml;
123 MesPlq(xml_out,
"Observables",
u);
137 multi1d<Real> lambda(
params.n_eig);
138 multi1d<Real> check_norm(
params.n_eig);
139 multi1d<LatticeFermion>
psi(
params.n_eig);
144 XMLBufferWriter eig_spec_xml;
163 write(xml_out,
"LowestEv",eig_spec_xml);
168 QDPIO::cout <<
"Look for highest ev" << std::endl;
176 QDPIO::cout <<
"ritz call" << std::endl;
178 XMLBufferWriter eig_spec_xml;
194 write(xml_out,
"HighestEv",eig_spec_xml);
Primary include file for CHROMA in application codes.
Create a simple ferm connection state.
virtual FermState< T, P, Q > * createState(const Q &q) const
Given links (coordinates Q) create the state needed for the linear operators.
Class for counted reference semantics.
Unpreconditioned Wilson fermion action.
virtual DiffLinearOperator< T, P, Q > * lMdagM(Handle< FermState< T, P, Q > > state) const
Produce a linear operator M^dag.M for this action.
void EigSpecRitzCG(const LinearOperator< LatticeFermion > &M, multi1d< Real > &lambda_H, multi1d< LatticeFermion > &psi, int n_eig, int n_renorm, int n_min, int MaxCG, const Real &Rsd_r, const Real &Rsd_a, const Real &zero_cutoff, const bool ProjApsiP, int &n_cg_tot, XMLWriter &xml_out)
Compute low lying eigenvalues of the hermitian H.
void read(XMLReader &xml, const std::string &path, AsqtadFermActParams ¶m)
Read parameters.
void write(XMLWriter &xml, const std::string &path, const AsqtadFermActParams ¶m)
Writer parameters.
void gaugeStartup(XMLReader &gauge_file_xml, XMLReader &gauge_xml, multi1d< LatticeColorMatrix > &u, Cfg_t &cfg)
Initialize the gauge fields.
multi1d< LatticeColorMatrix > P
Handle< FermBC< LatticeStaggeredFermion, multi1d< LatticeColorMatrix >, multi1d< LatticeColorMatrix > > > reader(XMLReader &xml_in, const std::string &path)
Helper function for the FermionAction readers.
Asqtad Staggered-Dirac operator.
static multi1d< LatticeColorMatrix > u
LinOpSysSolverMGProtoClover::Q Q
push(xml_out,"Condensates")
LinOpSysSolverMGProtoClover::T T
void initialize(int *argc, char ***argv)
Chroma initialisation routine.
void finalize(void)
Chroma finalization routine.
std::string getXMLInputFileName()
Get input file name.
void MesPlq(const multi1d< LatticeColorMatrixF3 > &u, multi2d< Double > &plane_plaq, Double &link)
XMLFileWriter & getXMLOutputInstance()
Get xml output instance.
multi1d< LatticeFermion > s(Ncb)
Gauge configuration structure.
Parameters for running program.
void readParams(const std::string &filename, Param_t ¶ms)
int main(int argc, char **argv)
void dumpParams(XMLWriter &writer, Param_t ¶ms)