4D Even Odd preconditioned NEF domain-wall Dirac operator More...

#include <eoprec_nef_linop_array_w.h>

Inheritance diagram for Chroma::EvenOddPrecNEFDWLinOpArray:

Public Types
typedef LatticeFermion	T

typedef multi1d< LatticeColorMatrix >	P

typedef multi1d< LatticeColorMatrix >	Q

Public Member Functions
	EvenOddPrecNEFDWLinOpArray (Handle< FermState< T, P, Q > > fs, const Real &WilsonMass_, const multi1d< Real > &b5_, const multi1d< Real > &c5_, const Real &m_q, int N5_)

	EvenOddPrecNEFDWLinOpArray (Handle< FermState< T, P, Q > > fs, const Real &WilsonMass_, const Real &b5_, const Real &c5_, const Real &m_q, int N5_)

void	create (Handle< FermState< T, P, Q > > fs, const Real &WilsonMass_, const Real &b5_, const Real &c5_, const Real &m_q_, int N5_)
	Creation routine. More...

	~EvenOddPrecNEFDWLinOpArray ()
	set b5 and c5 given kappa and a5 More...

const FermBC< T, P, Q > &	getFermBC () const
	Return the fermion BC object for this linear operator. More...

int	size () const
	Length of DW flavor index/space. More...

void	evenEvenLinOp (multi1d< LatticeFermion > &chi, const multi1d< LatticeFermion > &psi, enum PlusMinus isign) const
	Apply the even-even block onto a source std::vector. More...

void	evenEvenInvLinOp (multi1d< LatticeFermion > &chi, const multi1d< LatticeFermion > &psi, enum PlusMinus isign) const
	Apply the inverse of the even-even block onto a source std::vector. More...

void	evenOddLinOp (multi1d< LatticeFermion > &chi, const multi1d< LatticeFermion > &psi, enum PlusMinus isign) const
	Apply the the even-odd block onto a source std::vector. More...

void	oddEvenLinOp (multi1d< LatticeFermion > &chi, const multi1d< LatticeFermion > &psi, enum PlusMinus isign) const
	Apply the the odd-even block onto a source std::vector. More...

void	oddOddLinOp (multi1d< LatticeFermion > &chi, const multi1d< LatticeFermion > &psi, enum PlusMinus isign) const
	Apply the the odd-odd block onto a source std::vector. More...

void	oddOddInvLinOp (multi1d< LatticeFermion > &chi, const multi1d< LatticeFermion > &psi, enum PlusMinus isign) const
	Apply the inverse of the odd-odd block onto a source std::vector. More...

void	Dminus (LatticeFermion &chi, const LatticeFermion &psi, enum PlusMinus isign, int s5) const
	Apply the Dminus operator on a lattice fermion. See my notes ;-) More...

Public Member Functions inherited from Chroma::EvenOddPrecConstDetLinearOperatorArray< LatticeFermion, multi1d< LatticeColorMatrix >, multi1d< LatticeColorMatrix > >
virtual	~EvenOddPrecConstDetLinearOperatorArray ()
	Virtual destructor to help with cleanup;. More...

const Subset &	subset () const
	Only defined on the odd lattice. More...

virtual void	deriv (multi1d< LatticeColorMatrix > &ds_u, const multi1d< LatticeFermion > &chi, const multi1d< LatticeFermion > &psi, enum PlusMinus isign) const
	Apply the operator onto a source std::vector. More...

virtual void	derivEvenEvenLinOp (multi1d< LatticeColorMatrix > &ds_u, const multi1d< LatticeFermion > &chi, const multi1d< LatticeFermion > &psi, enum PlusMinus isign) const
	Apply the even-even block onto a source std::vector. More...

virtual void	derivEvenOddLinOp (multi1d< LatticeColorMatrix > &ds_u, const multi1d< LatticeFermion > &chi, const multi1d< LatticeFermion > &psi, enum PlusMinus isign) const
	Apply the the even-odd block onto a source std::vector. More...

virtual void	derivOddEvenLinOp (multi1d< LatticeColorMatrix > &ds_u, const multi1d< LatticeFermion > &chi, const multi1d< LatticeFermion > &psi, enum PlusMinus isign) const
	Apply the the odd-even block onto a source std::vector. More...

virtual void	derivOddOddLinOp (multi1d< LatticeColorMatrix > &ds_u, const multi1d< LatticeFermion > &chi, const multi1d< LatticeFermion > &psi, enum PlusMinus isign) const
	Apply the the odd-odd block onto a source std::vector. More...

Public Member Functions inherited from Chroma::EvenOddPrecLinearOperatorArray< T, P, Q >
virtual	~EvenOddPrecLinearOperatorArray ()
	Virtual destructor to help with cleanup;. More...

virtual void	operator() (multi1d< T > &chi, const multi1d< T > &psi, enum PlusMinus isign) const
	Apply the operator onto a source std::vector. More...

virtual void	unprecLinOp (multi1d< T > &chi, const multi1d< T > &psi, enum PlusMinus isign) const
	Apply the UNPRECONDITIONED operator onto a source std::vector. More...

virtual void	derivUnprecLinOp (P &ds_u, const multi1d< T > &chi, const multi1d< T > &psi, enum PlusMinus isign) const

virtual unsigned long	evenEvenNFlops () const
	Apply the derivative of the UNPRECONDITIONED operator onto a source std::vector. More...

virtual unsigned long	evenOddNFlops () const
	Return flops performed by the evenOddLinOp. More...

virtual unsigned long	oddEvenNFlops () const
	Return flops performed by the oddEvenLinOp. More...

virtual unsigned long	oddOddNFlops () const
	Return flops performed by the oddOddLinOp. More...

virtual unsigned long	evenEvenInvNFlops () const
	Return flops performed by the evenEvenInvLinOp. More...

virtual unsigned long	nFlops () const
	Return flops performed by the operator() More...

Public Member Functions inherited from Chroma::DiffLinearOperatorArray< T, P, Q >
virtual	~DiffLinearOperatorArray ()
	Virtual destructor to help with cleanup;. More...

virtual void	deriv (P &ds_u, const multi1d< T > &chi, const multi1d< T > &psi, enum PlusMinus isign, const Real &epsilon) const
	Apply the derivative of the operator onto a source std::vector to some precision. More...

Public Member Functions inherited from Chroma::LinearOperatorArray< T >
virtual	~LinearOperatorArray ()
	Virtual destructor to help with cleanup;. More...

virtual void	operator() (multi1d< T > &chi, const multi1d< T > &psi, enum PlusMinus isign, Real epsilon) const
	Apply the operator onto a source std::vector to some precision. More...

Protected Member Functions
void	applyDiag (multi1d< LatticeFermion > &chi, const multi1d< LatticeFermion > &psi, enum PlusMinus isign, const int cb) const
	Apply the even-even (odd-odd) coupling piece of the domain-wall fermion operator. More...

void	applyDiagInv (multi1d< LatticeFermion > &chi, const multi1d< LatticeFermion > &psi, enum PlusMinus isign, const int cb) const
	Apply the inverse even-even (odd-odd) coupling piece of the domain-wall fermion operator. More...

void	applyOffDiag (multi1d< LatticeFermion > &chi, const multi1d< LatticeFermion > &psi, enum PlusMinus isign, const int cb) const
	Apply the even-odd (odd-even) coupling piece of the NEF operator. More...

	EvenOddPrecNEFDWLinOpArray ()
	Partial constructor. More...

void	operator= (const EvenOddPrecNEFDWLinOpArray &)
	Partial constructor. More...

Private Attributes
Real	WilsonMass

Real	c5

Real	b5

Real	m_q

int	N5

Real	c5TwoKappa

Real	c5InvTwoKappa

Real	b5TwoKappa

Real	b5InvTwoKappa

Real	TwoKappa

Real	Kappa

Real	invDfactor

WilsonDslashArray	D

Detailed Description

4D Even Odd preconditioned NEF domain-wall Dirac operator

This routine is specific to Wilson fermions!

Definition at line 22 of file eoprec_nef_linop_array_w.h.

Member Typedef Documentation

◆ P

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

Definition at line 28 of file eoprec_nef_linop_array_w.h.

◆ Q

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

Definition at line 29 of file eoprec_nef_linop_array_w.h.

◆ T

typedef LatticeFermion Chroma::EvenOddPrecNEFDWLinOpArray::T

Definition at line 27 of file eoprec_nef_linop_array_w.h.

Constructor & Destructor Documentation

◆ EvenOddPrecNEFDWLinOpArray() [1/3]

Chroma::EvenOddPrecNEFDWLinOpArray::EvenOddPrecNEFDWLinOpArray	(	Handle< FermState< T, P, Q > >	fs,
		const Real &	WilsonMass_,
		const multi1d< Real > &	b5_,
		const multi1d< Real > &	c5_,
		const Real &	m_q,
		int	N5_
	)

inline

Definition at line 32 of file eoprec_nef_linop_array_w.h.

References Chroma::QDP_error_exit().

◆ EvenOddPrecNEFDWLinOpArray() [2/3]

Chroma::EvenOddPrecNEFDWLinOpArray::EvenOddPrecNEFDWLinOpArray	(	Handle< FermState< T, P, Q > >	fs,
		const Real &	WilsonMass_,
		const Real &	b5_,
		const Real &	c5_,
		const Real &	m_q,
		int	N5_
	)

inline

Full constructor Set b5 = 1.0 and c5=0.0 to get Shamir DWF with a5=1. Set b5 = 1.0 and c5=1.0 to get Borichi DWF.

Definition at line 44 of file eoprec_nef_linop_array_w.h.

References create(), and m_q.

◆ ~EvenOddPrecNEFDWLinOpArray()

Chroma::EvenOddPrecNEFDWLinOpArray::~EvenOddPrecNEFDWLinOpArray ( )

inline

set b5 and c5 given kappa and a5

Destructor is automatic

Definition at line 58 of file eoprec_nef_linop_array_w.h.

◆ EvenOddPrecNEFDWLinOpArray() [3/3]

Chroma::EvenOddPrecNEFDWLinOpArray::EvenOddPrecNEFDWLinOpArray ( )

inlineprotected

Partial constructor.

Definition at line 168 of file eoprec_nef_linop_array_w.h.

Member Function Documentation

◆ Dminus()

void Chroma::EvenOddPrecNEFDWLinOpArray::Dminus	(	LatticeFermion &	chi,
		const LatticeFermion &	psi,
		enum PlusMinus	isign,
		int	s5
	)		const

virtual

Apply the Dminus operator on a lattice fermion. See my notes ;-)

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

Definition at line 421 of file eoprec_nef_linop_array_w.cc.

References chi, Chroma::isign, and psi.

◆ evenEvenInvLinOp()

void Chroma::EvenOddPrecNEFDWLinOpArray::evenEvenInvLinOp	(	multi1d< LatticeFermion > &	chi,
		const multi1d< LatticeFermion > &	psi,
		enum PlusMinus	isign
	)		const

inlinevirtual

Apply the inverse of the even-even block onto a source std::vector.

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

Definition at line 77 of file eoprec_nef_linop_array_w.h.

References applyDiagInv(), Chroma::chi(), Chroma::isign, and Chroma::psi.

◆ evenEvenLinOp()

void Chroma::EvenOddPrecNEFDWLinOpArray::evenEvenLinOp	(	multi1d< LatticeFermion > &	chi,
		const multi1d< LatticeFermion > &	psi,
		enum PlusMinus	isign
	)		const

inlinevirtual

Apply the even-even block onto a source std::vector.

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

Definition at line 68 of file eoprec_nef_linop_array_w.h.

References applyDiag(), Chroma::chi(), Chroma::isign, and Chroma::psi.

◆ evenOddLinOp()

void Chroma::EvenOddPrecNEFDWLinOpArray::evenOddLinOp	(	multi1d< LatticeFermion > &	chi,
		const multi1d< LatticeFermion > &	psi,
		enum PlusMinus	isign
	)		const

inlinevirtual

Apply the the even-odd block onto a source std::vector.

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

Definition at line 85 of file eoprec_nef_linop_array_w.h.

References applyOffDiag(), Chroma::chi(), Chroma::isign, and Chroma::psi.

◆ getFermBC()

const FermBC<T,P,Q>& Chroma::EvenOddPrecNEFDWLinOpArray::getFermBC ( ) const

inlinevirtual

Return the fermion BC object for this linear operator.

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

Definition at line 61 of file eoprec_nef_linop_array_w.h.

References D, and Chroma::QDPWilsonDslashArrayOpt::getFermBC().

◆ oddEvenLinOp()

void Chroma::EvenOddPrecNEFDWLinOpArray::oddEvenLinOp	(	multi1d< LatticeFermion > &	chi,
		const multi1d< LatticeFermion > &	psi,
		enum PlusMinus	isign
	)		const

inlinevirtual

Apply the the odd-even block onto a source std::vector.

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

Definition at line 93 of file eoprec_nef_linop_array_w.h.

References applyOffDiag(), Chroma::chi(), Chroma::isign, and Chroma::psi.

◆ oddOddInvLinOp()

void Chroma::EvenOddPrecNEFDWLinOpArray::oddOddInvLinOp	(	multi1d< LatticeFermion > &	chi,
		const multi1d< LatticeFermion > &	psi,
		enum PlusMinus	isign
	)		const

inline

Apply the inverse of the odd-odd block onto a source std::vector.

Definition at line 111 of file eoprec_nef_linop_array_w.h.

References applyDiagInv(), Chroma::chi(), Chroma::isign, and Chroma::psi.

◆ oddOddLinOp()

void Chroma::EvenOddPrecNEFDWLinOpArray::oddOddLinOp	(	multi1d< LatticeFermion > &	chi,
		const multi1d< LatticeFermion > &	psi,
		enum PlusMinus	isign
	)		const