EvenOddPreconditioned Extended-Overlap (N&N) linear operator. More...

#include <eoprec_ovext_linop_array_w.h>

Inheritance diagram for Chroma::EvenOddPrecOvExtLinOpArray:

Public Types
typedef LatticeFermion	T

typedef multi1d< LatticeColorMatrix >	P

typedef multi1d< LatticeColorMatrix >	Q

Public Member Functions
	EvenOddPrecOvExtLinOpArray (Handle< FermState< T, P, Q > > fs, const int Npoles_, const Real &coeffP_, const multi1d< Real > &resP_, const multi1d< Real > &rootQ_, const multi1d< Real > &beta_, const Real &OverMass_, const Real &Mass_, const Real &b5_, const Real &c5_)
	Full constructor. More...

void	create (Handle< FermState< T, P, Q > > fs, const int Npoles_, const Real &coeffP_, const multi1d< Real > &resP_, const multi1d< Real > &rootQ_, const multi1d< Real > &beta_, const Real &OverMass_, const Real &m_q_, const Real &b5_, const Real &c5_)
	Creation routine. More...

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

	~EvenOddPrecOvExtLinOpArray ()
	Destructor is automatic. More...

const FermBC< T, P, Q > &	getFermBC () const
	Return the fermion BC object for this linear operator. 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	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	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	evenEvenInvLinOp (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	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	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...

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...

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...

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...

unsigned long	evenEvenNFlops (void) const
	Return flops performed by the evenEvenLinOp. More...

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

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

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

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

unsigned long	nFlops () const
	Return flops performed by the operator() 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...

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

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 off diagonal block. More...

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

unsigned long	diagNFlops (void) const
	Return flops performed by the diagonal part. More...

unsigned long	offDiagNFlops (void) const
	Return flops performed by the off diagonal part. More...

unsigned long	diagInvNFlops (void) const
	Return flops performed by the diag inv part. More...

	EvenOddPrecOvExtLinOpArray ()
	Partial constructor. More...

void	operator= (const EvenOddPrecOvExtLinOpArray &)
	Hide =. More...

Private Attributes
int	Npoles

int	N5

Real	A

multi1d< Real >	B

multi1d< Real >	D

multi1d< Real >	C

Real	E

WilsonDslashArray	Dslash

Real	Aprime

Real	Eprime

multi1d< Real >	Bprime

multi1d< Real >	Dprime

multi1d< Real >	Cprime

multi1d< Real >	Atilde

multi1d< Real >	Btilde

multi1d< Real >	Ctilde

multi1d< Real >	D_bd_inv

Real	S

Detailed Description

EvenOddPreconditioned Extended-Overlap (N&N) linear operator.

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

This operator implements hep-lat/0005004

Definition at line 29 of file eoprec_ovext_linop_array_w.h.

Member Typedef Documentation

◆ P

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

Definition at line 35 of file eoprec_ovext_linop_array_w.h.

◆ Q

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

Definition at line 36 of file eoprec_ovext_linop_array_w.h.

◆ T

typedef LatticeFermion Chroma::EvenOddPrecOvExtLinOpArray::T

Definition at line 34 of file eoprec_ovext_linop_array_w.h.

Constructor & Destructor Documentation

◆ EvenOddPrecOvExtLinOpArray() [1/2]

Chroma::EvenOddPrecOvExtLinOpArray::EvenOddPrecOvExtLinOpArray	(	Handle< FermState< T, P, Q > >	fs,
		const int	Npoles_,
		const Real &	coeffP_,
		const multi1d< Real > &	resP_,
		const multi1d< Real > &	rootQ_,
		const multi1d< Real > &	beta_,
		const Real &	OverMass_,
		const Real &	Mass_,
		const Real &	b5_,
		const Real &	c5_
	)

inline

Full constructor.

Definition at line 39 of file eoprec_ovext_linop_array_w.h.

◆ ~EvenOddPrecOvExtLinOpArray()

Chroma::EvenOddPrecOvExtLinOpArray::~EvenOddPrecOvExtLinOpArray ( )

inline

Destructor is automatic.

Definition at line 68 of file eoprec_ovext_linop_array_w.h.

◆ EvenOddPrecOvExtLinOpArray() [2/2]

Chroma::EvenOddPrecOvExtLinOpArray::EvenOddPrecOvExtLinOpArray ( )

inlineprotected

Partial constructor.

Definition at line 269 of file eoprec_ovext_linop_array_w.h.

Member Function Documentation

◆ applyDerivOffDiag()

void Chroma::EvenOddPrecOvExtLinOpArray::applyDerivOffDiag	(	multi1d< LatticeColorMatrix > &	ds_u,
		const multi1d< LatticeFermion > &	chi,
		const multi1d< LatticeFermion > &	psi,
		enum PlusMinus	isign,
		int	cb
	)		const

protected

Apply the even-odd (odd-even) coupling piece of the NEF operator.

Parameters

ds_u	conjugate momenta (Read)
psi	left pseudofermion field (Read)
psi	right pseudofermion field (Read)
isign	Flag ( PLUS \| MINUS ) (Read)
cb	checkerboard ( 0 \| 1 ) (Read)

Definition at line 431 of file eoprec_ovext_linop_array_w.cc.

References Chroma::END_CODE(), and Chroma::START_CODE().

◆ applyDiag()

void Chroma::EvenOddPrecOvExtLinOpArray::applyDiag	(	multi1d< LatticeFermion > &	chi,
		const multi1d< LatticeFermion > &	psi,
		enum PlusMinus	isign,
		const int	cb
	)		const

protected

Apply the even-even (odd-odd) coupling piece of the domain-wall fermion operator.

Parameters

chi	result (Modify)
psi	source (Read)
isign	Flag ( PLUS \| MINUS ) (Read)
cb	checkerboard ( 0 \| 1 ) (Read)
chi	result (Modify)
psi	source (Read)
isign	Flag ( PLUS \| MINUS ) (Read)
cb	checkerboard ( 0 \| 1 ) (Read)

Flopcount = 10*Nc*Ns*(N5-1) + 2*Nc*Ns = 10*Nc*Ns*N5 - 8*Nc*Ns = (10N5 - 8)*Nc*Ns

Definition at line 109 of file eoprec_ovext_linop_array_w.cc.

References A, B, C, Chroma::cb, Chroma::chi(), D, E, Chroma::END_CODE(), Chroma::G5, Chroma::i, Chroma::isign, Chroma::MINUS, N5, Chroma::p, Chroma::PLUS, Chroma::psi, and Chroma::START_CODE().

◆ applyDiagInv()

void Chroma::EvenOddPrecOvExtLinOpArray::applyDiagInv	(	multi1d< LatticeFermion > &	chi,
		const multi1d< LatticeFermion > &	psi,
		enum PlusMinus	isign,
		const int	cb
	)		const

protected

Apply the inverse even-even (odd-odd) coupling piece of the domain-wall fermion operator.

Parameters

chi	result (Modify)
psi	source (Read)
isign	Flag ( PLUS \| MINUS ) (Read)
cb	checkerboard ( 0 \| 1 ) (Read)

Definition at line 266 of file eoprec_ovext_linop_array_w.cc.

References Atilde, Btilde, Chroma::cb, Chroma::chi(), Ctilde, D_bd_inv, Chroma::END_CODE(), ftmp2, Chroma::i, Chroma::isign, Chroma::MINUS, N5, Npoles, Chroma::p, Chroma::PLUS, Chroma::psi, S, and Chroma::START_CODE().

◆ applyOffDiag()

void Chroma::EvenOddPrecOvExtLinOpArray::applyOffDiag	(	multi1d< LatticeFermion > &	chi,
		const multi1d< LatticeFermion > &	psi,
		enum PlusMinus	isign,
		const int	cb
	)		const

protected

Apply the off diagonal block.

Parameters

chi	result (Modify)
psi	source (Read)
isign	Flag ( PLUS \| MINUS ) (Read)
cb	checkerboard ( 0 \| 1 ) (Read)

Definition at line 176 of file eoprec_ovext_linop_array_w.cc.

References Chroma::QDPWilsonDslashArrayOpt::apply(), Aprime, Bprime, Chroma::cb, Chroma::chi(), Cprime, Dprime, Dslash, Chroma::END_CODE(), Eprime, Chroma::G5, Chroma::i, Chroma::isign, Chroma::MINUS, N5, Chroma::p, Chroma::PLUS, Chroma::psi, Chroma::START_CODE(), and Chroma::tmp.

◆ derivEvenEvenLinOp()

void Chroma::EvenOddPrecOvExtLinOpArray::derivEvenEvenLinOp	(	multi1d< LatticeColorMatrix > &	ds_u,
		const multi1d< LatticeFermion > &	chi,
		const multi1d< LatticeFermion > &	psi,
		enum PlusMinus	isign
	)		const

inlinevirtual

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

Reimplemented from Chroma::EvenOddPrecConstDetLinearOperatorArray< LatticeFermion, multi1d< LatticeColorMatrix >, multi1d< LatticeColorMatrix > >.

Definition at line 126 of file eoprec_ovext_linop_array_w.h.

References Nd, and Chroma::zero.

◆ derivEvenOddLinOp()

void Chroma::EvenOddPrecOvExtLinOpArray::derivEvenOddLinOp	(	multi1d< LatticeColorMatrix > &	ds_u,
		const multi1d< LatticeFermion > &	chi,
		const multi1d< LatticeFermion > &	psi,
		enum PlusMinus	isign
	)		const

inlinevirtual

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

Reimplemented from Chroma::EvenOddPrecConstDetLinearOperatorArray< LatticeFermion, multi1d< LatticeColorMatrix >, multi1d< LatticeColorMatrix > >.

Definition at line 146 of file eoprec_ovext_linop_array_w.h.

References chi, Chroma::isign, and psi.

◆ derivOddEvenLinOp()

void Chroma::EvenOddPrecOvExtLinOpArray::derivOddEvenLinOp	(	multi1d< LatticeColorMatrix > &	ds_u,
		const multi1d< LatticeFermion > &	chi,
		const multi1d< LatticeFermion > &	psi,
		enum PlusMinus	isign
	)		const

inlinevirtual

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

Reimplemented from Chroma::EvenOddPrecConstDetLinearOperatorArray< LatticeFermion, multi1d< LatticeColorMatrix >, multi1d< LatticeColorMatrix > >.

Definition at line 154 of file eoprec_ovext_linop_array_w.h.

References chi, Chroma::isign, and psi.

◆ derivOddOddLinOp()

void Chroma::EvenOddPrecOvExtLinOpArray::derivOddOddLinOp	(	multi1d< LatticeColorMatrix > &	ds_u,
		const multi1d< LatticeFermion > &	chi,
		const multi1d< LatticeFermion > &	psi,
		enum PlusMinus	isign
	)		const

inlinevirtual

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

Reimplemented from Chroma::EvenOddPrecConstDetLinearOperatorArray< LatticeFermion, multi1d< LatticeColorMatrix >, multi1d< LatticeColorMatrix > >.

Definition at line 137 of file eoprec_ovext_linop_array_w.h.

References Nd, and Chroma::zero.

◆ diagInvNFlops()

unsigned long Chroma::EvenOddPrecOvExtLinOpArray::diagInvNFlops ( void ) const

inlineprotected

Return flops performed by the diag inv part.

Definition at line 262 of file eoprec_ovext_linop_array_w.h.

References N5.

◆ diagNFlops()

unsigned long Chroma::EvenOddPrecOvExtLinOpArray::diagNFlops ( void ) const

inlineprotected

Return flops performed by the diagonal part.

Definition at line 250 of file eoprec_ovext_linop_array_w.h.

References N5.

◆ evenEvenInvLinOp()

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

inlinevirtual

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

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

Definition at line 108 of file eoprec_ovext_linop_array_w.h.

References chi, Chroma::isign, and psi.

◆ evenEvenInvNFlops()

unsigned long Chroma::EvenOddPrecOvExtLinOpArray::evenEvenInvNFlops ( void ) const

inlinevirtual

Return flops performed by the evenEvenInvLinOp.

Reimplemented from Chroma::EvenOddPrecLinearOperatorArray< T, P, Q >.

Definition at line 183 of file eoprec_ovext_linop_array_w.h.

◆ evenEvenLinOp()

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

inlinevirtual

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

This does not need to be optimized

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

Definition at line 74 of file eoprec_ovext_linop_array_w.h.

References chi, Chroma::isign, and psi.

◆ evenEvenNFlops()

unsigned long Chroma::EvenOddPrecOvExtLinOpArray::evenEvenNFlops ( void ) const

inlinevirtual

Return flops performed by the evenEvenLinOp.

Reimplemented from Chroma::EvenOddPrecLinearOperatorArray< T, P, Q >.

Definition at line 163 of file eoprec_ovext_linop_array_w.h.

◆ evenOddLinOp()

void Chroma::EvenOddPrecOvExtLinOpArray::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 91 of file eoprec_ovext_linop_array_w.h.

References chi, Chroma::isign, and psi.

◆ evenOddNFlops()

unsigned long Chroma::EvenOddPrecOvExtLinOpArray::evenOddNFlops ( void ) const

inlinevirtual

Return flops performed by the evenOddLinOp.

Reimplemented from Chroma::EvenOddPrecLinearOperatorArray< T, P, Q >.

Definition at line 173 of file eoprec_ovext_linop_array_w.h.

◆ getFermBC()

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

inlinevirtual

Return the fermion BC object for this linear operator.

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

Definition at line 71 of file eoprec_ovext_linop_array_w.h.

◆ nFlops()

unsigned long Chroma::EvenOddPrecOvExtLinOpArray::nFlops ( void ) const

inlinevirtual

Return flops performed by the operator()

Reimplemented from Chroma::EvenOddPrecLinearOperatorArray< T, P, Q >.

Definition at line 188 of file eoprec_ovext_linop_array_w.h.

References N5.

◆ oddEvenLinOp()

void Chroma::EvenOddPrecOvExtLinOpArray::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 99 of file eoprec_ovext_linop_array_w.h.

References chi, Chroma::isign, and psi.

◆ oddEvenNFlops()

unsigned long Chroma::EvenOddPrecOvExtLinOpArray::oddEvenNFlops ( void ) const

inlinevirtual

Return flops performed by the oddEvenLinOp.

Reimplemented from Chroma::EvenOddPrecLinearOperatorArray< T, P, Q >.

Definition at line 178 of file eoprec_ovext_linop_array_w.h.

◆ oddOddInvLinOp()

void Chroma::EvenOddPrecOvExtLinOpArray::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 117 of file eoprec_ovext_linop_array_w.h.

References chi, Chroma::isign, and psi.

◆ oddOddLinOp()

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