Chroma::EvenOddPrecHtContFrac5DFermActArray

Chroma::EvenOddPrecHtContFrac5DFermActArray Class Reference
[Fermion actions]

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

#include <eoprec_ht_contfrac5d_fermact_array_w.h>

Inheritance diagram for Chroma::EvenOddPrecHtContFrac5DFermActArray:


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.
int	size (void) const
Real	getQuarkMass () const
	Return the quark mass.
EvenOddPrecConstDetLinearOperatorArray< T, P, Q > *	linOp (Handle< FermState< T, P, Q > > state) const
	Produce a linear operator for this action.
EvenOddPrecConstDetLinearOperatorArray< T, P, Q > *	linOpPV (Handle< FermState< T, P, Q > > state) const
	Produce a Pauli-Villars linear operator for this action.
LinearOperatorArray< T > *	hermitianLinOp (Handle< FermState< T, P, Q > > state) const
	produce hermitian version of linOp
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).
LinearOperator< T > *	DeltaLs (Handle< FermState< T, P, Q > > state, const GroupXML_t &invParam) const
	Produce a DeltaLs = 1-epsilon^2(H) operator.
SystemSolver< LatticeFermion > *	qprop (Handle< FermState< T, P, Q > > state, const GroupXML_t &invParam) const
	Propagator of an un-preconditioned Extended-Overlap linear operator.
	~EvenOddPrecHtContFrac5DFermActArray ()
	Destructor is automatic.
Protected Member Functions
const CreateFermState< T, P, Q > &	getCreateState () const
	Return the fermion BC object for this action.
void	init (Real &scale_fac, multi1d< Real > &alpha, multi1d< Real > &beta) const
	Helper in construction.
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 66 of file eoprec_ht_contfrac5d_fermact_array_w.h.

Member Typedef Documentation

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

Definition at line 73 of file eoprec_ht_contfrac5d_fermact_array_w.h.

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

Definition at line 74 of file eoprec_ht_contfrac5d_fermact_array_w.h.

typedef LatticeFermion Chroma::EvenOddPrecHtContFrac5DFermActArray::T

Definition at line 72 of file eoprec_ht_contfrac5d_fermact_array_w.h.

Constructor & Destructor Documentation

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

Definition at line 156 of file eoprec_ht_contfrac5d_fermact_array_w.cc.

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

Chroma::EvenOddPrecHtContFrac5DFermActArray::EvenOddPrecHtContFrac5DFermActArray ( const EvenOddPrecHtContFrac5DFermActArray & a ) [inline]

Copy constructor.

Definition at line 82 of file eoprec_ht_contfrac5d_fermact_array_w.h.

Chroma::EvenOddPrecHtContFrac5DFermActArray::~EvenOddPrecHtContFrac5DFermActArray ( ) [inline]

Destructor is automatic.

Definition at line 128 of file eoprec_ht_contfrac5d_fermact_array_w.h.

Chroma::EvenOddPrecHtContFrac5DFermActArray::EvenOddPrecHtContFrac5DFermActArray ( ) [private]

Member Function Documentation

LinearOperator<T>* Chroma::EvenOddPrecHtContFrac5DFermActArray::DeltaLs	(	Handle< FermState< T, P, Q > >	state,
		const GroupXML_t &	invParam
	)			const `[inline]`

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

Definition at line 116 of file eoprec_ht_contfrac5d_fermact_array_w.h.

References linOp4D(), and Chroma::state.

const CreateFermState<T,P,Q>& Chroma::EvenOddPrecHtContFrac5DFermActArray::getCreateState ( ) const [inline, protected]

Return the fermion BC object for this action.

Definition at line 132 of file eoprec_ht_contfrac5d_fermact_array_w.h.

References fs.

Real Chroma::EvenOddPrecHtContFrac5DFermActArray::getQuarkMass ( ) const [inline]

Return the quark mass.

Definition at line 88 of file eoprec_ht_contfrac5d_fermact_array_w.h.

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

Referenced by qprop().

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

produce hermitian version of linOp

Definition at line 97 of file eoprec_ht_contfrac5d_fermact_array_w.h.

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

Helper in construction.

Definition at line 180 of file eoprec_ht_contfrac5d_fermact_array_w.cc.

References Chroma::EvenOddPrecHtContFrac5DFermActParams::approximation_type, Chroma::EvenOddPrecHtContFrac5DFermActParams::ApproxMax, Chroma::EvenOddPrecHtContFrac5DFermActParams::ApproxMin, Chroma::COEFF_TYPE_TANH_UNSCALED, Chroma::COEFF_TYPE_ZOLOTAREV, Chroma::ExternalFieldEnv::epsilon(), higham(), isLastZeroP, Chroma::EvenOddPrecHtContFrac5DFermActParams::Mass, N5, Chroma::EvenOddPrecHtContFrac5DFermActParams::OverMass, params, Chroma::EvenOddPrecHtContFrac5DFermActParams::RatPolyDeg, and zolotarev().

Referenced by linOp(), and linOpPV().

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

Produce a linear operator for this action.

Parameters:

state

gauge field (Read)

Definition at line 320 of file eoprec_ht_contfrac5d_fermact_array_w.cc.

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

Referenced by qprop().

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

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

Definition at line 106 of file eoprec_ht_contfrac5d_fermact_array_w.h.

Referenced by DeltaLs().

Produce a Pauli-Villars linear operator for this action.

Parameters:

state

gauge field (Read)

Definition at line 350 of file eoprec_ht_contfrac5d_fermact_array_w.cc.

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

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