Quark propagator solution routines
[Fermion action support]


Files

file  asqtad_qprop.h
 Asqtad propagator wrapper.
file  hisq_qprop.h
 Hisq propagator wrapper.

Namespaces

namespace  Chroma::AVPSolver
 Bluegene altivec single-prec solver.
namespace  Chroma::AVPSolver
 Bluegene altivec single-prec solver.
namespace  Chroma::AVPSolver
 Bluegene altivec single-prec solver.
namespace  Chroma::AVPSolver
 Bluegene altivec single-prec solver.
namespace  Chroma::AVPSolver
 Bluegene altivec single-prec solver.
namespace  Chroma::AVPSolver
 Bluegene altivec single-prec solver.
namespace  Chroma::AVPSolverFunctions
 AVP's DWF Solver interface.
namespace  Chroma::AVPSolver
 Bluegene altivec single-prec solver.
namespace  Chroma::AVPSolver
 Bluegene altivec single-prec solver.

Classes

class  Chroma::HtContFrac5DQprop< T, P, Q >
 Propagator of an un-preconditioned Extended-Overlap linear operator. More...
class  Chroma::PrecOvExt5DQprop< T, P, Q >
 Propagator of an un-preconditioned Extended-Overlap linear operator. More...
class  Chroma::ContFrac5DQprop< T, P, Q >
 Propagator of an un-preconditioned Extended-Overlap linear operator. More...
class  Chroma::Ovlap4DQprop
class  Chroma::OvHTCFZ5DQprop< T >
 Propagator of unpreconditioned H_T kernel continued fraction (5D) operator. More...
class  Chroma::OvExt5DQprop< T >
 Propagator of an un-preconditioned Extended-Overlap linear operator. More...
class  Chroma::OvUnprecCF5DQprop< T >
 Propagator of an un-preconditioned Extended-Overlap linear operator. More...
class  Chroma::AsqtadCPSWrapperQprop
 QPROP for wrapping CPS asqtad QPROP. More...
class  Chroma::DWFQprop< T >
 Propagator DWF linear operator. More...
class  Chroma::CGDWFQpropT< SinglePrecSolver, DoublePrecSolver >
 SSE Propagator DWF qpropT. More...
class  Chroma::PrecFermActQprop< T, P, Q >
 Propagator of a generic even-odd preconditioned fermion linear operator. More...
class  Chroma::PrecFermAct5DQprop< T, P, Q >
 Propagator of a generic even-odd preconditioned 5D fermion linear operator. More...
class  Chroma::EvenOddFermActQprop< T, P, Q >
 Propagator of a generic even-odd fermion linear operator. More...
class  Chroma::FermActQprop< T >
 Propagator of a generic non-preconditioned fermion linear operator. More...
class  Chroma::FermAct5DQprop< T >
class  Chroma::MDWFQpropT
 AVP's DWF Solver interface. More...

Functions

void Chroma::dwf_quarkProp4 (LatticePropagator &q_sol, XMLWriter &xml_out, const LatticePropagator &q_src, int t_src, int j_decay, Handle< SystemSolverArray< LatticeFermion > > qpropT, Handle< FermState< LatticeFermion, multi1d< LatticeColorMatrix >, multi1d< LatticeColorMatrix > > > state, const Real &m_q, int &ncg_had)
 Given a complete propagator as a source, this does all the inversions needed.
template<typename T, typename P, typename Q, template< class, class, class > class C>
void Chroma::nef_quarkProp_a (LatticePropagator &q_sol, XMLWriter &xml_out, const LatticePropagator &q_src, int t_src, int j_decay, const C< T, P, Q > &S_f, Handle< FermState< T, P, Q > > state, const GroupXML_t &invParam, int &ncg_had)
 Given a complete propagator as a source, this does all the inversions needed.
void Chroma::nef_quarkProp4 (LatticePropagator &q_sol, XMLWriter &xml_out, const LatticePropagator &q_src, int t_src, int j_decay, const UnprecDWFermActBaseArray< LatticeFermion, multi1d< LatticeColorMatrix >, multi1d< LatticeColorMatrix > > &S_f, Handle< FermState< LatticeFermion, multi1d< LatticeColorMatrix >, multi1d< LatticeColorMatrix > > > state, const GroupXML_t &invParam, int &ncg_had)
 Given a complete propagator as a source, this does all the inversions needed.
void Chroma::nef_quarkProp4 (LatticePropagator &q_sol, XMLWriter &xml_out, const LatticePropagator &q_src, int t_src, int j_decay, const EvenOddPrecDWFermActBaseArray< LatticeFermion, multi1d< LatticeColorMatrix >, multi1d< LatticeColorMatrix > > &S_f, Handle< FermState< LatticeFermion, multi1d< LatticeColorMatrix >, multi1d< LatticeColorMatrix > > > state, const GroupXML_t &invParam, int &ncg_had)
 Given a complete propagator as a source, this does all the inversions needed.
template<typename T, typename P, typename Q>
void Chroma::quarkProp_a (LatticeStaggeredPropagator &q_sol, XMLWriter &xml_out, const LatticeStaggeredPropagator &q_src, const StaggeredTypeFermAct< T, P, Q > &S_f, Handle< FermState< T, P, Q > > state, const GroupXML_t &invParam, QuarkSpinType quarkSpinType, int &ncg_had)
 Given a complete propagator as a source, this does all the inversions needed.
void Chroma::quarkProp4 (LatticeStaggeredPropagator &q_sol, XMLWriter &xml_out, const LatticeStaggeredPropagator &q_src, const StaggeredTypeFermAct< LatticeStaggeredFermion, multi1d< LatticeColorMatrix >, multi1d< LatticeColorMatrix > > &S_f, Handle< FermState< LatticeStaggeredFermion, multi1d< LatticeColorMatrix >, multi1d< LatticeColorMatrix > > > state, const GroupXML_t &invParam, QuarkSpinType quarkSpinType, int &ncg_had)
 Given a complete propagator as a source, this does all the inversions needed.
template<typename T>
void Chroma::quarkProp4_a (LatticePropagator &q_sol, XMLWriter &xml_out, const LatticePropagator &q_src, Handle< SystemSolver< T > > qprop, QuarkSpinType quarkSpinType, int &ncg_had)
 Given a complete propagator as a source, this does all the inversions needed.
void Chroma::quarkProp4 (LatticePropagator &q_sol, XMLWriter &xml_out, const LatticePropagator &q_src, Handle< SystemSolver< LF > > qprop, QuarkSpinType quarkSpinType, int &ncg_had)
 Given a complete propagator as a source, this does all the inversions needed.
void Chroma::quarkProp4 (LatticePropagator &q_sol, XMLWriter &xml_out, const LatticePropagator &q_src, Handle< SystemSolver< LatticeFermion > > qprop, QuarkSpinType quarkSpinType, int &ncg_had)
 Given a complete propagator as a source, this does all the inversions needed.
void Chroma::OverlapFermActBase::multiQuarkProp (multi1d< LatticePropagator > &q_sol, XMLWriter &xml_out, const LatticePropagator &q_src, Handle< FermState< T, P, Q > > state, const multi1d< Real > &masses, const GroupXML_t &invParam, const int n_soln, int &ncg_had)
 Given a complete propagator as a source, this does all the inversions needed.

Detailed Description

Routines for computing a quark propagator with various fermion actions

Function Documentation

void Chroma::dwf_quarkProp4 ( LatticePropagator &  q_sol,
XMLWriter &  xml_out,
const LatticePropagator &  q_src,
int  t_src,
int  j_decay,
Handle< SystemSolverArray< LatticeFermion > >  qpropT,
Handle< FermState< LatticeFermion, multi1d< LatticeColorMatrix >, multi1d< LatticeColorMatrix > > >  state,
const Real &  m_q,
int &  ncg_had 
)

Given a complete propagator as a source, this does all the inversions needed.

This routine is actually generic to Domain Wall fermions (Array) fermions

Parameters:
q_sol quark propagator ( Write )
q_src source ( Read )
t_src time slice of source ( Read )
j_decay direction of decay ( Read )
qpropT 5D inverter ( Read )
ncg_had number of CG iterations ( Write )

Definition at line 220 of file dwf_quarkprop4_w.cc.

References Chroma::check_dwf_ward_identity(), Chroma::chi(), Chroma::dwf_conserved_axial_ps_corr(), END_CODE, Chroma::FermToProp(), N5, Chroma::SystemSolverResults_t::n_count, Chroma::PropToFerm(), Chroma::psi(), Chroma::SystemSolverResults_t::resid, START_CODE, Chroma::state, and Chroma::write().

void Chroma::OverlapFermActBase::multiQuarkProp ( multi1d< LatticePropagator > &  q_sol,
XMLWriter &  xml_out,
const LatticePropagator &  q_src,
Handle< FermState< T, P, Q > >  state,
const multi1d< Real > &  masses,
const GroupXML_t &  invParam,
const int  n_soln,
int &  ncg_had 
) [inherited]

Given a complete propagator as a source, this does all the inversions needed.

Parameters:
q_sol quark propagator ( Write )
q_src source ( Read )
invParam inverter parameters ( Read )
ncg_had number of CG iterations ( Write )

Definition at line 114 of file quarkprop4_multi_w.cc.

References Chroma::state.

Referenced by Chroma::InlineMultiPropagator::func().

void Chroma::nef_quarkProp4 ( LatticePropagator &  q_sol,
XMLWriter &  xml_out,
const LatticePropagator &  q_src,
int  t_src,
int  j_decay,
const EvenOddPrecDWFermActBaseArray< LatticeFermion, multi1d< LatticeColorMatrix >, multi1d< LatticeColorMatrix > > &  S_f,
Handle< FermState< LatticeFermion, multi1d< LatticeColorMatrix >, multi1d< LatticeColorMatrix > > >  state,
const GroupXML_t &  invParam,
int &  ncg_had 
)

Given a complete propagator as a source, this does all the inversions needed.

This routine is actually generic to Domain Wall fermions (Array) fermions

Parameters:
q_sol quark propagator ( Write )
q_src source ( Read )
t_src time slice of source ( Read )
j_decay direction of decay ( Read )
invParam inverter params ( Read )
ncg_had number of CG iterations ( Write )

Definition at line 349 of file nef_quarkprop4_w.cc.

References Chroma::S_f, and Chroma::state.

void Chroma::nef_quarkProp4 ( LatticePropagator &  q_sol,
XMLWriter &  xml_out,
const LatticePropagator &  q_src,
int  t_src,
int  j_decay,
const UnprecDWFermActBaseArray< LatticeFermion, multi1d< LatticeColorMatrix >, multi1d< LatticeColorMatrix > > &  S_f,
Handle< FermState< LatticeFermion, multi1d< LatticeColorMatrix >, multi1d< LatticeColorMatrix > > >  state,
const GroupXML_t &  invParam,
int &  ncg_had 
)

Given a complete propagator as a source, this does all the inversions needed.

This routine is actually generic to Domain Wall fermions (Array) fermions

Parameters:
q_sol quark propagator ( Write )
q_src source ( Read )
t_src time slice of source ( Read )
j_decay direction of decay ( Read )
invParam inverter params ( Read )
ncg_had number of CG iterations ( Write )

Definition at line 312 of file nef_quarkprop4_w.cc.

References Chroma::S_f, and Chroma::state.

Referenced by Chroma::UnprecZoloNEFFermActArray::quarkProp(), Chroma::UnprecNEFFermActArray::quarkProp(), Chroma::EvenOddPrecZoloNEFFermActArray::quarkProp(), Chroma::EvenOddPrecNEFFermActArray::quarkProp(), and Chroma::EvenOddPrecKNOFermActArray::quarkProp().

template<typename T, typename P, typename Q, template< class, class, class > class C>
void Chroma::nef_quarkProp_a ( LatticePropagator &  q_sol,
XMLWriter &  xml_out,
const LatticePropagator &  q_src,
int  t_src,
int  j_decay,
const C< T, P, Q > &  S_f,
Handle< FermState< T, P, Q > >  state,
const GroupXML_t &  invParam,
int &  ncg_had 
)

Given a complete propagator as a source, this does all the inversions needed.

This routine is actually generic to Domain Wall fermions (Array) fermions

Parameters:
q_sol quark propagator ( Write )
q_src source ( Read )
t_src time slice of source ( Read )
j_decay direction of decay ( Read )
invParam inverter params ( Read )
ncg_had number of CG iterations ( Write )

constuct the conserved axial current correlator nef_conserved_axial_ps_corr(cfield,state->getLinks(),prop5d,j_decay);

push(xml_out, "DWF_ConservedAxial"); write(xml_out, "mesprop", mesprop); pop(xml_out);

The local axial corruent pseudoscalar correlator int d(1<<j_decay); cfield = trace( adj(q_sol)*Gamma(d)*q_sol ) ; corr = sumMulti(cfield, trick.getSet()) ; for(int t(0);t<length; t++){ int t_eff( (t - t_src + length) % length ) ; mesprop[t_eff] = -real(corr[t]) ; // sign fix } push(xml_out, "DWF_LocalAxial"); write(xml_out, "mesprop", mesprop); pop(xml_out);

check_nef_ward_identity(state->getLinks(),prop5d,q_src, q_sol,q_mp,S_f.quark_mass(), j_decay);

Definition at line 114 of file nef_quarkprop4_w.cc.

References Chroma::chi(), END_CODE, Chroma::FermToProp(), Chroma::SftMom::getSet(), N5, Chroma::SystemSolverResults_t::n_count, Chroma::SftMom::numSubsets(), Chroma::PLUS, Chroma::PropToFerm(), Chroma::psi(), Chroma::SystemSolverResults_t::resid, START_CODE, Chroma::state, and Chroma::write().

void Chroma::quarkProp4 ( LatticePropagator &  q_sol,
XMLWriter &  xml_out,
const LatticePropagator &  q_src,
Handle< SystemSolver< LatticeFermion > >  qprop,
QuarkSpinType  quarkSpinType,
int &  ncg_had 
)

Given a complete propagator as a source, this does all the inversions needed.

This routine is actually generic to all Wilson-like fermions

Parameters:
q_sol quark propagator ( Write )
q_src source ( Read )
qprop inverter ( Read )
ncg_had number of CG iterations ( Write )

void Chroma::quarkProp4 ( LatticePropagator &  q_sol,
XMLWriter &  xml_out,
const LatticePropagator &  q_src,
Handle< SystemSolver< LF > >  qprop,
QuarkSpinType  quarkSpinType,
int &  ncg_had 
)

Given a complete propagator as a source, this does all the inversions needed.

This routine is actually generic to all Wilson-like fermions

Parameters:
q_sol quark propagator ( Write )
q_src source ( Read )
invParam inverter parameters ( Read )
ncg_had number of CG iterations ( Write )

Definition at line 206 of file quarkprop4_w.cc.

void Chroma::quarkProp4 ( LatticeStaggeredPropagator &  q_sol,
XMLWriter &  xml_out,
const LatticeStaggeredPropagator &  q_src,
const StaggeredTypeFermAct< LatticeStaggeredFermion, multi1d< LatticeColorMatrix >, multi1d< LatticeColorMatrix > > &  S_f,
Handle< FermState< LatticeStaggeredFermion, multi1d< LatticeColorMatrix >, multi1d< LatticeColorMatrix > > >  state,
const GroupXML_t &  invParam,
QuarkSpinType  quarkSpinType,
int &  ncg_had 
)

Given a complete propagator as a source, this does all the inversions needed.

This routine is actually generic to all Wilson-like fermions

Parameters:
q_sol quark propagator ( Write )
q_src source ( Read )
invParam inverter parameters ( Read )
ncg_had number of CG iterations ( Write )

Definition at line 125 of file quarkprop4_s.cc.

References Chroma::quarkProp_a(), Chroma::S_f, and Chroma::state.

Referenced by main(), Chroma::UnprecZoloNEFFermActArray::quarkProp(), Chroma::UnprecNEFFermActArray::quarkProp(), Chroma::EvenOddPrecZoloNEFFermActArray::quarkProp(), Chroma::EvenOddPrecNEFFermActArray::quarkProp(), and Chroma::EvenOddPrecKNOFermActArray::quarkProp().

template<typename T>
void Chroma::quarkProp4_a ( LatticePropagator &  q_sol,
XMLWriter &  xml_out,
const LatticePropagator &  q_src,
Handle< SystemSolver< T > >  qprop,
QuarkSpinType  quarkSpinType,
int &  ncg_had 
)

Given a complete propagator as a source, this does all the inversions needed.

This routine is actually generic to all Wilson-like fermions

Parameters:
q_sol quark propagator ( Write )
q_src source ( Read )
RsdCG CG (or MR) residual used here ( Read )
MaxCG maximum number of CG iterations ( Read )
ncg_had number of CG iterations ( Write )

Definition at line 33 of file quarkprop4_w.cc.

References Chroma::chi(), END_CODE, Chroma::FermToProp(), Chroma::SystemSolverResults_t::n_count, Chroma::PropToFerm(), Chroma::psi(), Chroma::QUARK_SPIN_TYPE_FULL, Chroma::QUARK_SPIN_TYPE_LOWER, Chroma::QUARK_SPIN_TYPE_UPPER, Chroma::SystemSolverResults_t::resid, START_CODE, and Chroma::write().

template<typename T, typename P, typename Q>
void Chroma::quarkProp_a ( LatticeStaggeredPropagator &  q_sol,
XMLWriter &  xml_out,
const LatticeStaggeredPropagator &  q_src,
const StaggeredTypeFermAct< T, P, Q > &  S_f,
Handle< FermState< T, P, Q > >  state,
const GroupXML_t &  invParam,
QuarkSpinType  quarkSpinType,
int &  ncg_had 
)

Given a complete propagator as a source, this does all the inversions needed.

This routine is actually generic to all Wilson-like fermions

Parameters:
q_sol quark propagator ( Write )
q_src source ( Read )
invParam inverter parameters ( Read )
RsdCG CG (or MR) residual used here ( Read )
MaxCG maximum number of CG iterations ( Read )
ncg_had number of CG iterations ( Write )

Definition at line 35 of file quarkprop4_s.cc.

References Chroma::chi(), END_CODE, Chroma::FermToProp(), Chroma::SystemSolverResults_t::n_count, Chroma::PropToFerm(), Chroma::psi(), Chroma::SystemSolverResults_t::resid, Chroma::S_f, START_CODE, Chroma::state, and Chroma::write().

Referenced by Chroma::quarkProp4().

Generated on Sun Nov 22 04:38:37 2009 for CHROMA by  doxygen 1.4.7