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

#include <eoprec_ovlap_contfrac5d_fermact_array_w.h>

Inheritance diagram for Chroma::EvenOddPrecOvlapContFrac5DFermActArray:

Public Types
typedef LatticeFermion	T

typedef multi1d< LatticeColorMatrix >	P

typedef multi1d< LatticeColorMatrix >	Q

Public Member Functions
	EvenOddPrecOvlapContFrac5DFermActArray (Handle< CreateFermState< T, P, Q > > cfs_, const EvenOddPrecOvlapContFrac5DFermActParams &param)
	Construct the action out of a parameter structure. More...

	EvenOddPrecOvlapContFrac5DFermActArray (const EvenOddPrecOvlapContFrac5DFermActArray &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 a hermitian version of this operator More...

LinearOperator< T > *	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...

	~EvenOddPrecOvlapContFrac5DFermActArray ()
	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^dagPVpsi=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^dagPV+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^dagPVpsi=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^dagPV+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^dagPVpsi=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
	EvenOddPrecOvlapContFrac5DFermActArray ()

Private Attributes
Handle< CreateFermState< T, P, Q > >	cfs

EvenOddPrecOvlapContFrac5DFermActParams	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 63 of file eoprec_ovlap_contfrac5d_fermact_array_w.h.

Member Typedef Documentation

◆ P

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

Definition at line 69 of file eoprec_ovlap_contfrac5d_fermact_array_w.h.

◆ Q

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

Definition at line 70 of file eoprec_ovlap_contfrac5d_fermact_array_w.h.

◆ T

typedef LatticeFermion Chroma::EvenOddPrecOvlapContFrac5DFermActArray::T

Definition at line 68 of file eoprec_ovlap_contfrac5d_fermact_array_w.h.

Constructor & Destructor Documentation

◆ EvenOddPrecOvlapContFrac5DFermActArray() [1/3]

Chroma::EvenOddPrecOvlapContFrac5DFermActArray::EvenOddPrecOvlapContFrac5DFermActArray	(	Handle< CreateFermState< T, P, Q > >	cfs_,
		const EvenOddPrecOvlapContFrac5DFermActParams &	param
	)

Construct the action out of a parameter structure.

Definition at line 138 of file eoprec_ovlap_contfrac5d_fermact_array_w.cc.

References isLastZeroP, N5, params, and Chroma::EvenOddPrecOvlapContFrac5DFermActParams::RatPolyDeg.

◆ EvenOddPrecOvlapContFrac5DFermActArray() [2/3]

Chroma::EvenOddPrecOvlapContFrac5DFermActArray::EvenOddPrecOvlapContFrac5DFermActArray ( const EvenOddPrecOvlapContFrac5DFermActArray & a )

inline

Copy constructor.

Definition at line 78 of file eoprec_ovlap_contfrac5d_fermact_array_w.h.

◆ ~EvenOddPrecOvlapContFrac5DFermActArray()

Chroma::EvenOddPrecOvlapContFrac5DFermActArray::~EvenOddPrecOvlapContFrac5DFermActArray ( )

inline

Destructor is automatic.

Definition at line 122 of file eoprec_ovlap_contfrac5d_fermact_array_w.h.

◆ EvenOddPrecOvlapContFrac5DFermActArray() [3/3]

Chroma::EvenOddPrecOvlapContFrac5DFermActArray::EvenOddPrecOvlapContFrac5DFermActArray ( )

private

Member Function Documentation

◆ DeltaLs()

LinearOperator<T>* Chroma::EvenOddPrecOvlapContFrac5DFermActArray::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 110 of file eoprec_ovlap_contfrac5d_fermact_array_w.h.

References linOp4D(), and Chroma::state.

◆ getCreateState()

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

inlineprotectedvirtual

Return the fermion BC object for this action.

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

Definition at line 127 of file eoprec_ovlap_contfrac5d_fermact_array_w.h.

References cfs.

◆ getQuarkMass()

Real Chroma::EvenOddPrecOvlapContFrac5DFermActArray::getQuarkMass ( ) const

inlinevirtual

Return the quark mass.

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

Definition at line 84 of file eoprec_ovlap_contfrac5d_fermact_array_w.h.

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

Referenced by qprop().

◆ hermitianLinOp()

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

inlinevirtual

produce a hermitian version of this operator

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

Definition at line 94 of file eoprec_ovlap_contfrac5d_fermact_array_w.h.

References linOp(), and Chroma::state.

◆ init()

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

protected

◆ linOp()

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

virtual

Produce a linear operator for this action.

Parameters

state gauge field (Read)

Implements Chroma::EvenOddPrecConstDetWilsonTypeFermAct5D< LatticeFermion, multi1d< LatticeColorMatrix >, multi1d< LatticeColorMatrix > >.

Definition at line 303 of file eoprec_ovlap_contfrac5d_fermact_array_w.cc.

References Chroma::StagPhases::alpha(), Chroma::StagPhases::beta(), init(), isLastZeroP, Chroma::EvenOddPrecOvlapContFrac5DFermActParams::Mass, N5, Chroma::EvenOddPrecOvlapContFrac5DFermActParams::OverMass, params, and Chroma::START_CODE().

Referenced by hermitianLinOp(), main(), and qprop().

◆ linOp4D()

LinearOperator<T>* Chroma::EvenOddPrecOvlapContFrac5DFermActArray::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 100 of file eoprec_ovlap_contfrac5d_fermact_array_w.h.

Referenced by DeltaLs().

◆ linOpPV()

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

virtual

Produce a Pauli-Villars linear operator for this action.

Parameters

state gauge field (Read)

Implements Chroma::EvenOddPrecConstDetWilsonTypeFermAct5D< LatticeFermion, multi1d< LatticeColorMatrix >, multi1d< LatticeColorMatrix > >.

Definition at line 331 of file eoprec_ovlap_contfrac5d_fermact_array_w.cc.

References Chroma::StagPhases::alpha(), Chroma::StagPhases::beta(), init(), isLastZeroP, Chroma::EvenOddPrecOvlapContFrac5DFermActParams::Mass, N5, Chroma::EvenOddPrecOvlapContFrac5DFermActParams::OverMass, and params.

Referenced by main().

◆ qprop()

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