CHROMA
Public Types | Public Member Functions | Protected Member Functions | Private Member Functions | Private Attributes | List of all members

5D continued fraction overlap action (Borici,Wenger, Edwards) More...

#include <eoprec_ht_contfrac5d_fermact_array_w.h>

Inheritance diagram for Chroma::EvenOddPrecHtContFrac5DFermActArray:
Chroma::EvenOddPrecConstDetWilsonTypeFermAct5D< LatticeFermion, multi1d< LatticeColorMatrix >, multi1d< LatticeColorMatrix > > Chroma::EvenOddPrecWilsonTypeFermAct5D< T, P, Q > Chroma::EvenOddPrecWilsonTypeFermAct5D< T, P, Q > Chroma::WilsonTypeFermAct5D< T, P, Q > Chroma::WilsonTypeFermAct5D< T, P, Q > Chroma::WilsonTypeFermAct5D< T, P, Q > Chroma::WilsonTypeFermAct5D< T, P, Q > Chroma::DiffFermAct5D< T, P, Q > Chroma::DiffFermAct5D< T, P, Q > Chroma::DiffFermAct5D< T, P, Q > Chroma::DiffFermAct5D< T, P, Q > Chroma::DiffFermAct5D< T, P, Q > Chroma::DiffFermAct5D< T, P, Q > Chroma::DiffFermAct5D< T, P, Q > Chroma::DiffFermAct5D< T, P, Q >

Public Types

typedef LatticeFermion T
 
typedef multi1d< LatticeColorMatrix > P
 
typedef multi1d< LatticeColorMatrix > Q
 

Public Member Functions

 EvenOddPrecHtContFrac5DFermActArray (Handle< CreateFermState< T, P, Q > > fs_, const EvenOddPrecHtContFrac5DFermActParams &param)
 
 EvenOddPrecHtContFrac5DFermActArray (const EvenOddPrecHtContFrac5DFermActArray &a)
 Copy constructor. More...
 
int size (void) const
 Expected length of array index. More...
 
Real getQuarkMass () const
 Return the quark mass. More...
 
EvenOddPrecConstDetLinearOperatorArray< T, P, Q > * linOp (Handle< FermState< T, P, Q > > state) const
 Produce a linear operator for this action. More...
 
EvenOddPrecConstDetLinearOperatorArray< T, P, Q > * linOpPV (Handle< FermState< T, P, Q > > state) const
 Produce a Pauli-Villars linear operator for this action. More...
 
LinearOperatorArray< T > * hermitianLinOp (Handle< FermState< T, P, Q > > state) const
 produce hermitian version of linOp More...
 
LinearOperator< LatticeFermion > * linOp4D (Handle< FermState< T, P, Q > > state, const Real &m_q, const GroupXML_t &invParam) const
 Produce an unpreconditioned linear operator projecting 5D to 4D (the inverse of qprop below) More...
 
LinearOperator< T > * DeltaLs (Handle< FermState< T, P, Q > > state, const GroupXML_t &invParam) const
 Produce a DeltaLs = 1-epsilon^2(H) operator. More...
 
SystemSolver< LatticeFermion > * qprop (Handle< FermState< T, P, Q > > state, const GroupXML_t &invParam) const
 Compute quark propagator over base type. More...
 
 ~EvenOddPrecHtContFrac5DFermActArray ()
 Destructor is automatic. More...
 
- Public Member Functions inherited from Chroma::EvenOddPrecConstDetWilsonTypeFermAct5D< LatticeFermion, multi1d< LatticeColorMatrix >, multi1d< LatticeColorMatrix > >
virtual ~EvenOddPrecConstDetWilsonTypeFermAct5D ()
 Virtual destructor to help with cleanup;. More...
 
virtual ~EvenOddPrecConstDetWilsonTypeFermAct5D ()
 Virtual destructor to help with cleanup;. More...
 
- Public Member Functions inherited from Chroma::EvenOddPrecWilsonTypeFermAct5D< T, P, Q >
virtual ~EvenOddPrecWilsonTypeFermAct5D ()
 Virtual destructor to help with cleanup;. More...
 
virtual SystemSolverArray< T > * qpropT (Handle< FermState< T, P, Q > > state, const GroupXML_t &invParam) const
 Return quark prop solver, solution of unpreconditioned system. More...
 
virtual ~EvenOddPrecWilsonTypeFermAct5D ()
 Virtual destructor to help with cleanup;. More...
 
virtual SystemSolverArray< T > * qpropT (Handle< FermState< T, P, Q > > state, const GroupXML_t &invParam) const
 Return quark prop solver, solution of unpreconditioned system. More...
 
SystemSolverArray< LF > * qpropT (Handle< FermState< LF, LCM, LCM > > state, const GroupXML_t &invParam) const
 Propagator of a generic even-odd preconditioned fermion linear operator. More...
 
- Public Member Functions inherited from Chroma::WilsonTypeFermAct5D< T, P, Q >
virtual ~WilsonTypeFermAct5D ()
 Virtual destructor to help with cleanup;. More...
 
virtual DiffLinearOperatorArray< T, P, Q > * lMdagM (Handle< FermState< T, P, Q > > state) const
 Produce a linear operator M^dag.M for this action. More...
 
virtual LinOpSystemSolverArray< T > * invLinOp (Handle< FermState< T, P, Q > > state, const GroupXML_t &invParam) const
 Return a linear operator solver for this action to solve M*psi=chi. More...
 
virtual MdagMSystemSolverArray< T > * invMdagM (Handle< FermState< T, P, Q > > state, const GroupXML_t &invParam) const
 Return a linear operator solver for this action to solve MdagM*psi=chi. More...
 
virtual LinOpSystemSolverArray< T > * invLinOpPV (Handle< FermState< T, P, Q > > state, const GroupXML_t &invParam) const
 Return a linear operator solver for this action to solve PV*psi=chi. More...
 
virtual MdagMSystemSolverArray< T > * invMdagMPV (Handle< FermState< T, P, Q > > state, const GroupXML_t &invParam) const
 Return a linear operator solver for this action to solve PV^dag*PV*psi=chi. More...
 
virtual MdagMMultiSystemSolverArray< T > * mInvMdagM (Handle< FermState< T, P, Q > > state, const GroupXML_t &invParam) const
 Return a multi-shift linear operator solver for this action to solve (MdagM+shift)*psi=chi. More...
 
virtual MdagMMultiSystemSolverArray< T > * mInvMdagMPV (Handle< FermState< T, P, Q > > state, const GroupXML_t &invParam) const
 Return a multi-shift linear operator solver for this action to solve (PV^dag*PV+shift)*psi=chi. More...
 
virtual void quarkProp (typename PropTypeTraits< T >::Type_t &q_sol, XMLWriter &xml_out, const typename PropTypeTraits< T >::Type_t &q_src, Handle< FermState< T, P, Q > > state, const GroupXML_t &invParam, QuarkSpinType quarkSpinType, int &ncg_had) const
 Given a complete propagator as a source, this does all the inversions needed. More...
 
virtual ~WilsonTypeFermAct5D ()
 Virtual destructor to help with cleanup;. More...
 
virtual DiffLinearOperatorArray< T, P, Q > * lMdagM (Handle< FermState< T, P, Q > > state) const
 Produce a linear operator M^dag.M for this action. More...
 
virtual LinOpSystemSolverArray< T > * invLinOp (Handle< FermState< T, P, Q > > state, const GroupXML_t &invParam) const
 Return a linear operator solver for this action to solve M*psi=chi. More...
 
virtual MdagMSystemSolverArray< T > * invMdagM (Handle< FermState< T, P, Q > > state, const GroupXML_t &invParam) const
 Return a linear operator solver for this action to solve MdagM*psi=chi. More...
 
virtual LinOpSystemSolverArray< T > * invLinOpPV (Handle< FermState< T, P, Q > > state, const GroupXML_t &invParam) const
 Return a linear operator solver for this action to solve PV*psi=chi. More...
 
virtual MdagMSystemSolverArray< T > * invMdagMPV (Handle< FermState< T, P, Q > > state, const GroupXML_t &invParam) const
 Return a linear operator solver for this action to solve PV^dag*PV*psi=chi. More...
 
virtual MdagMMultiSystemSolverArray< T > * mInvMdagM (Handle< FermState< T, P, Q > > state, const GroupXML_t &invParam) const
 Return a multi-shift linear operator solver for this action to solve (MdagM+shift)*psi=chi. More...
 
virtual MdagMMultiSystemSolverAccumulateArray< T > * mInvMdagMAcc (Handle< FermState< T, P, Q > > state, const GroupXML_t &invParam) const
 Return a multi-shift linear operator solver for this action to solve (MdagM+shift)*psi=chi. More...
 
virtual MdagMMultiSystemSolverArray< T > * mInvMdagMPV (Handle< FermState< T, P, Q > > state, const GroupXML_t &invParam) const
 Return a multi-shift linear operator solver for this action to solve (MdagM+shift)*psi=chi. More...
 
virtual MdagMMultiSystemSolverAccumulateArray< T > * mInvMdagMPVAcc (Handle< FermState< T, P, Q > > state, const GroupXML_t &invParam) const
 Return a multi-shift linear operator solver for this action to solve (PV^dag*PV+shift)*psi=chi. More...
 
virtual void quarkProp (typename PropTypeTraits< T >::Type_t &q_sol, XMLWriter &xml_out, const typename PropTypeTraits< T >::Type_t &q_src, Handle< FermState< T, P, Q > > state, const GroupXML_t &invParam, QuarkSpinType quarkSpinType, int &ncg_had) const
 Given a complete propagator as a source, this does all the inversions needed. More...
 
void quarkProp (LatticePropagator &q_sol, XMLWriter &xml_out, const LatticePropagator &q_src, Handle< FermState< LF, LCM, LCM > > state, const GroupXML_t &invParam, QuarkSpinType quarkSpinType, int &ncg_had) const
 Given a complete propagator as a source, this does all the inversions needed. More...
 
LinOpSystemSolverArray< LF > * invLinOp (Handle< FermState< LF, LCM, LCM > > state, const GroupXML_t &invParam) const
 
MdagMSystemSolverArray< LF > * invMdagM (Handle< FermState< LF, LCM, LCM > > state, const GroupXML_t &invParam) const
 Return a linear operator solver for this action to solve MdagM*psi=chi. More...
 
LinOpSystemSolverArray< LF > * invLinOpPV (Handle< FermState< LF, LCM, LCM > > state, const GroupXML_t &invParam) const
 
MdagMSystemSolverArray< LF > * invMdagMPV (Handle< FermState< LF, LCM, LCM > > state, const GroupXML_t &invParam) const
 Return a linear operator solver for this action to solve PV^dag*PV*psi=chi. More...
 
MdagMMultiSystemSolverArray< LF > * mInvMdagM (Handle< FermState< LF, LCM, LCM > > state, const GroupXML_t &invParam) const
 Return a linear operator solver for this action to solve (MdagM+shift_i)*psi_i = chi. More...
 
MdagMMultiSystemSolverAccumulateArray< LF > * mInvMdagMAcc (Handle< FermState< LF, LCM, LCM > > state, const GroupXML_t &invParam) const
 Return a linear operator solver for this action to solve (MdagM+shift_i)*psi_i = chi. More...
 
MdagMMultiSystemSolverArray< LF > * mInvMdagMPV (Handle< FermState< LF, LCM, LCM > > state, const GroupXML_t &invParam) const
 Return a linear operator solver for this action to solve (MdagM+shift_i)*psi_i = chi. More...
 
MdagMMultiSystemSolverAccumulateArray< LF > * mInvMdagMPVAcc (Handle< FermState< LF, LCM, LCM > > state, const GroupXML_t &invParam) const
 Return a linear operator solver for this action to solve (MdagM+shift_i)*psi_i = chi. More...
 
- Public Member Functions inherited from Chroma::DiffFermAct5D< T, P, Q >
virtual ~DiffFermAct5D ()
 Virtual destructor to help with cleanup;. More...
 
virtual ~DiffFermAct5D ()
 Virtual destructor to help with cleanup;. More...
 
- Public Member Functions inherited from Chroma::FermAct5D< T, P, Q >
virtual ~FermAct5D ()
 Virtual destructor to help with cleanup;. More...
 
virtual ~FermAct5D ()
 Virtual destructor to help with cleanup;. More...
 
SystemSolverArray< LatticeFermion > * qpropT (Handle< FermState< LatticeFermion, multi1d< LatticeColorMatrix >, multi1d< LatticeColorMatrix > > > state, const GroupXML_t &invParam) const
 
- Public Member Functions inherited from Chroma::FermionAction< T, P, Q >
virtual ~FermionAction ()
 Virtual destructor to help with cleanup;. More...
 
virtual FermState< T, P, Q > * createState (const Q &q) const
 Given links (coordinates Q) create the state needed for the linear operators. More...
 
virtual FermState< T, P, Q > * createState (const Q &q, XMLReader &reader, const std::string &path) const
 Given links (coordinates Q) create a state with additional info held by the XMLReader. More...
 
virtual const FermBC< T, P, Q > & getFermBC () const
 Return the fermion BC object for this action. More...
 
virtual void quarkProp (typename PropTypeTraits< T >::Type_t &q_sol, XMLWriter &xml_out, const typename PropTypeTraits< T >::Type_t &q_src, int t_src, int j_decay, Handle< FermState< T, P, Q > > state, const GroupXML_t &invParam, QuarkSpinType quarkSpinType, bool obsvP, int &ncg_had) const
 Given a complete propagator as a source, this does all the inversions needed. More...
 
virtual ~FermionAction ()
 Virtual destructor to help with cleanup;. More...
 
virtual FermState< T, P, Q > * createState (const Q &q) const
 Given links (coordinates Q) create the state needed for the linear operators. More...
 
virtual FermState< T, P, Q > * createState (const Q &q, XMLReader &reader, const std::string &path) const
 Given links (coordinates Q) create a state with additional info held by the XMLReader. More...
 
virtual const FermBC< T, P, Q > & getFermBC () const
 Return the fermion BC object for this action. More...
 
virtual void quarkProp (typename PropTypeTraits< T >::Type_t &q_sol, XMLWriter &xml_out, const typename PropTypeTraits< T >::Type_t &q_src, int t_src, int j_decay, Handle< FermState< T, P, Q > > state, const GroupXML_t &invParam, QuarkSpinType quarkSpinType, bool obsvP, int &ncg_had) const
 Given a complete propagator as a source, this does all the inversions needed. More...
 

Protected Member Functions

const CreateFermState< T, P, Q > & getCreateState () const
 Return the fermion BC object for this action. More...
 
void init (Real &scale_fac, multi1d< Real > &alpha, multi1d< Real > &beta) const
 Helper in construction. More...
 

Private Member Functions

 EvenOddPrecHtContFrac5DFermActArray ()
 

Private Attributes

Handle< CreateFermState< T, P, Q > > fs
 
EvenOddPrecHtContFrac5DFermActParams params
 
int N5
 
bool isLastZeroP
 

Detailed Description

5D continued fraction overlap action (Borici,Wenger, Edwards)

This operator applies the extended version of the hermitian overlap operator Chi = ((1+Mass)/(1-Mass)*gamma_5 + B) . Psi where B is the continued fraction of the zolotarev approx. to eps(H(m))

Definition at line 65 of file eoprec_ht_contfrac5d_fermact_array_w.h.

Member Typedef Documentation

◆ P

typedef multi1d<LatticeColorMatrix> Chroma::EvenOddPrecHtContFrac5DFermActArray::P

Definition at line 72 of file eoprec_ht_contfrac5d_fermact_array_w.h.

◆ Q

typedef multi1d<LatticeColorMatrix> Chroma::EvenOddPrecHtContFrac5DFermActArray::Q

Definition at line 73 of file eoprec_ht_contfrac5d_fermact_array_w.h.

◆ T

Definition at line 71 of file eoprec_ht_contfrac5d_fermact_array_w.h.

Constructor & Destructor Documentation

◆ EvenOddPrecHtContFrac5DFermActArray() [1/3]

Chroma::EvenOddPrecHtContFrac5DFermActArray::EvenOddPrecHtContFrac5DFermActArray ( Handle< CreateFermState< T, P, Q > >  fs_,
const EvenOddPrecHtContFrac5DFermActParams param 
)

◆ EvenOddPrecHtContFrac5DFermActArray() [2/3]

Chroma::EvenOddPrecHtContFrac5DFermActArray::EvenOddPrecHtContFrac5DFermActArray ( const EvenOddPrecHtContFrac5DFermActArray a)
inline

Copy constructor.

Definition at line 81 of file eoprec_ht_contfrac5d_fermact_array_w.h.

◆ ~EvenOddPrecHtContFrac5DFermActArray()

Chroma::EvenOddPrecHtContFrac5DFermActArray::~EvenOddPrecHtContFrac5DFermActArray ( )
inline

Destructor is automatic.

Definition at line 127 of file eoprec_ht_contfrac5d_fermact_array_w.h.

◆ EvenOddPrecHtContFrac5DFermActArray() [3/3]

Chroma::EvenOddPrecHtContFrac5DFermActArray::EvenOddPrecHtContFrac5DFermActArray ( )
private

Member Function Documentation

◆ DeltaLs()

LinearOperator<T>* Chroma::EvenOddPrecHtContFrac5DFermActArray::DeltaLs ( Handle< FermState< T, P, Q > >  state,
const GroupXML_t invParam 
) const
inlinevirtual

Produce a DeltaLs = 1-epsilon^2(H) operator.

Implements Chroma::WilsonTypeFermAct5D< T, P, Q >.

Definition at line 115 of file eoprec_ht_contfrac5d_fermact_array_w.h.

References linOp4D(), and Chroma::state.

◆ getCreateState()

const CreateFermState<T,P,Q>& Chroma::EvenOddPrecHtContFrac5DFermActArray::getCreateState ( ) const
inlineprotectedvirtual

Return the fermion BC object for this action.

Implements Chroma::FermionAction< T, P, Q >.

Definition at line 131 of file eoprec_ht_contfrac5d_fermact_array_w.h.

References fs.

◆ getQuarkMass()

Real Chroma::EvenOddPrecHtContFrac5DFermActArray::getQuarkMass ( ) const
inlinevirtual

Return the quark mass.

Implements Chroma::FermAct5D< T, P, Q >.

Definition at line 87 of file eoprec_ht_contfrac5d_fermact_array_w.h.

References Chroma::EvenOddPrecHtContFrac5DFermActParams::Mass, and params.

Referenced by qprop().

◆ hermitianLinOp()

LinearOperatorArray<T>* Chroma::EvenOddPrecHtContFrac5DFermActArray::hermitianLinOp ( Handle< FermState< T, P, Q > >  state) const
inlinevirtual

produce hermitian version of linOp

Implements Chroma::WilsonTypeFermAct5D< T, P, Q >.

Definition at line 96 of file eoprec_ht_contfrac5d_fermact_array_w.h.

◆ init()

void Chroma::EvenOddPrecHtContFrac5DFermActArray::init ( Real &  scale_fac,
multi1d< Real > &  alpha,
multi1d< Real > &  beta 
) const
protected

◆ linOp()

EvenOddPrecConstDetLinearOperatorArray< LatticeFermion, multi1d< LatticeColorMatrix >, multi1d< LatticeColorMatrix > > * Chroma::EvenOddPrecHtContFrac5DFermActArray::linOp ( Handle< FermState< T, P, Q > >  state_) const
virtual

◆ linOp4D()

LinearOperator<LatticeFermion>* Chroma::EvenOddPrecHtContFrac5DFermActArray::linOp4D ( Handle< FermState< T, P, Q > >  state,
const Real &  m_q,
const GroupXML_t invParam 
) const
inlinevirtual

Produce an unpreconditioned linear operator projecting 5D to 4D (the inverse of qprop below)

Implements Chroma::WilsonTypeFermAct5D< T, P, Q >.

Definition at line 105 of file eoprec_ht_contfrac5d_fermact_array_w.h.

Referenced by DeltaLs().

◆ linOpPV()

EvenOddPrecConstDetLinearOperatorArray< LatticeFermion, multi1d< LatticeColorMatrix >, multi1d< LatticeColorMatrix > > * Chroma::EvenOddPrecHtContFrac5DFermActArray::linOpPV ( Handle< FermState< T, P, Q > >  state_) const
virtual

◆ qprop()

SystemSolver< LatticeFermion > * Chroma::EvenOddPrecHtContFrac5DFermActArray::qprop ( Handle< FermState< T, P, Q > >  state,
const GroupXML_t invParam 
) const
virtual

◆ size()

int Chroma::EvenOddPrecHtContFrac5DFermActArray::size ( void  ) const
inlinevirtual

Expected length of array index.

Implements Chroma::FermAct5D< T, P, Q >.

Definition at line 84 of file eoprec_ht_contfrac5d_fermact_array_w.h.

References N5.

Member Data Documentation

◆ fs

Handle< CreateFermState<T,P,Q> > Chroma::EvenOddPrecHtContFrac5DFermActArray::fs
private

Definition at line 143 of file eoprec_ht_contfrac5d_fermact_array_w.h.

Referenced by getCreateState().

◆ isLastZeroP

bool Chroma::EvenOddPrecHtContFrac5DFermActArray::isLastZeroP
private

◆ N5

int Chroma::EvenOddPrecHtContFrac5DFermActArray::N5
private

◆ params

EvenOddPrecHtContFrac5DFermActParams Chroma::EvenOddPrecHtContFrac5DFermActArray::params
private

The documentation for this class was generated from the following files: