Clover term. More...

#include <clover_term_bagel_clover.h>

Inheritance diagram for Chroma::BAGELCloverTerm:

Public Types
typedef LatticeFermion	T

typedef multi1d< LatticeColorMatrix >	P

typedef multi1d< LatticeColorMatrix >	Q

Public Member Functions
	BAGELCloverTerm ()
	Empty constructor. Must use create later. More...

	~BAGELCloverTerm ()
	Free the internals. More...

void	create (Handle< FermState< T, P, Q > > fs, const CloverFermActParams &param_, const BAGELCloverTerm &from)
	Create from another. More...

void	create (Handle< FermState< T, P, Q > > fs, const CloverFermActParams &param_)
	Creation routine. More...

void	choles (int cb)
	Computes the inverse of the term on cb using Cholesky. More...

Double	cholesDet (int cb) const
	Computes the inverse of the term on cb using Cholesky. More...

void	apply (LatticeFermion &chi, const LatticeFermion &psi, enum PlusMinus isign, int cb) const

void	applySite (LatticeFermion &chi, const LatticeFermion &psi, enum PlusMinus isign, int site) const

void	triacntr (LatticeColorMatrix &B, int mat, int cb) const
	Calculates Tr_D ( Gamma_mat L ) More...

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

Public Member Functions inherited from Chroma::CloverTermBase< LatticeFermion, LatticeColorMatrix >
virtual	~CloverTermBase ()
	No real need for cleanup here. More...

const Subset &	subset () const
	Subset is all here. More...

void	deriv (multi1d< LatticeColorMatrix > &ds_u, const LatticeFermion &chi, const LatticeFermion &psi, enum PlusMinus isign) const
	Take deriv of D. More...

void	deriv (multi1d< LatticeColorMatrix > &ds_u, const LatticeFermion &chi, const LatticeFermion &psi, enum PlusMinus isign, int cb) const
	Take deriv of D. More...

void	derivMultipole (multi1d< LatticeColorMatrix > &ds_u, const multi1d< LatticeFermion > &chi, const multi1d< LatticeFermion > &psi, enum PlusMinus isign) const
	Take deriv of D. More...

void	derivMultipole (multi1d< LatticeColorMatrix > &ds_u, const multi1d< LatticeFermion > &chi, const multi1d< LatticeFermion > &psi, enum PlusMinus isign, int cb) const
	Take deriv of D. More...

void	derivTrLn (multi1d< LatticeColorMatrix > &ds_u, enum PlusMinus isign, int cb) const
	Take derivative of TrLn D. More...

void	deriv_loops (const int u, const int mu, const int cb, LatticeColorMatrix &ds_u_mu, LatticeColorMatrix &ds_u_nu, const LatticeColorMatrix &Lambda) const

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

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

virtual void	operator() (T &d, const T &psi, enum PlusMinus isign) const
	Apply operator on both checkerboards (entire lattice) More...

virtual void	deriv (P &ds_u, const T &chi, const T &psi, enum PlusMinus isign) const
	Take deriv of D. More...

virtual void	deriv (P &ds_u, const T &chi, const T &psi, enum PlusMinus isign, int cb) const
	Take deriv of D. More...

virtual void	derivMultipole (P &ds_u, const multi1d< T > &chi, const multi1d< T > &psi, enum PlusMinus isign) const
	Return the force for multiple poles. More...

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

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

virtual void	deriv (P &ds_u, const T &chi, const 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::LinearOperator< T >
virtual	~LinearOperator ()
	Virtual destructor to help with cleanup;. More...

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

Protected Member Functions
void	makeClov (const multi1d< LatticeColorMatrix > &f, const Real &diag_mass)
	Create the clover term on cb. More...

void	ldagdlinv (LatticeReal &tr_log_diag, int cb)
	Invert the clover term on cb using LDL^\dagger decomp. More...

void	chlclovms (LatticeReal &log_diag, int cb)
	Invert the clover term on cb using Cholesky decomp. More...

const multi1d< LatticeColorMatrix > &	getU () const
	Get the u field. More...

Real	getCloverCoeff (int mu, int nu) const
	Calculates Tr_D ( Gamma_mat L ) More...

Private Attributes
Handle< FermBC< T, P, Q > >	fbc

multi1d< LatticeColorMatrix >	u

CloverFermActParams	param

LatticeReal	tr_log_diag_

multi1d< bool >	choles_done

PrimitiveClovDiag *	tri_diag

PrimitiveClovOffDiag *	tri_off_diag

Detailed Description

Clover term.

Definition at line 44 of file clover_term_bagel_clover.h.

Member Typedef Documentation

◆ P

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

Definition at line 49 of file clover_term_bagel_clover.h.

◆ Q

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

Definition at line 50 of file clover_term_bagel_clover.h.

◆ T

typedef LatticeFermion Chroma::BAGELCloverTerm::T

Definition at line 48 of file clover_term_bagel_clover.h.

Constructor & Destructor Documentation

◆ BAGELCloverTerm()

Chroma::BAGELCloverTerm::BAGELCloverTerm ( )

Empty constructor. Must use create later.

Definition at line 79 of file clover_term_bagel_clover.cc.

References Chroma::DEFAULT, tri_diag, and tri_off_diag.

◆ ~BAGELCloverTerm()

Chroma::BAGELCloverTerm::~BAGELCloverTerm ( )

Free the internals.

Definition at line 124 of file clover_term_bagel_clover.cc.

References tri_diag, and tri_off_diag.

Member Function Documentation

◆ apply()

void Chroma::BAGELCloverTerm::apply	(	LatticeFermion &	chi,
		const LatticeFermion &	psi,
		enum PlusMinus	isign,
		int	cb
	)		const

virtual

Apply a dslash

Performs the operation

chi <- (L + D + L^dag) . psi

where L is a lower triangular matrix D is the real diagonal. (stored together in type TRIANG)

Arguments:

Parameters

chi	result (Write)
psi	source (Read)
isign	D'^dag or D' ( MINUS \| PLUS ) resp. (Read)
cb	Checkerboard of OUTPUT std::vector (Read)

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

Definition at line 970 of file clover_term_bagel_clover.cc.

References Chroma::cb, Chroma::chi(), Chroma::END_CODE(), getFermBC(), n, Chroma::psi, Chroma::QDP_error_exit(), Chroma::START_CODE(), tri_diag, and tri_off_diag.

Referenced by create().

◆ applySite()

void Chroma::BAGELCloverTerm::applySite	(	LatticeFermion &	chi,
		const LatticeFermion &	psi,
		enum PlusMinus	isign,
		int	site
	)		const

virtual

Apply a dslash

Performs the operation

chi <- (L + D + L^dag) . psi

where L is a lower triangular matrix D is the real diagonal. (stored together in type TRIANG)

Arguments:

Parameters

chi	result (Write)
psi	source (Read)
isign	D'^dag or D' ( MINUS \| PLUS ) resp. (Read)
cb	Checkerboard of OUTPUT std::vector (Read)

Implements Chroma::CloverTermBase< LatticeFermion, LatticeColorMatrix >.

Definition at line 1043 of file clover_term_bagel_clover.cc.

References Chroma::chi(), Chroma::END_CODE(), n, Chroma::psi, Chroma::QDP_error_exit(), Chroma::START_CODE(), tri_diag, and tri_off_diag.

Referenced by create().

◆ chlclovms()

void Chroma::BAGELCloverTerm::chlclovms	(	LatticeReal &	tr_log_diag,
		int	cb
	)

protected

Invert the clover term on cb using Cholesky decomp.

CHLCLOVMS - Cholesky decompose the clover mass term and uses it to compute lower(A^-1) = lower((L.L^dag)^-1) Adapted from Golub and Van Loan, Matrix Computations, 2nd, Sec 4.2.4

Arguments:

Parameters