Even-odd preconditioned 5D linear operator. More...

#include <eoprec_logdet_linop.h>

Inheritance diagram for Chroma::EvenOddPrecLogDetLinearOperatorArray< T, P, Q >:

Public Member Functions
virtual	~EvenOddPrecLogDetLinearOperatorArray ()
	Virtual destructor to help with cleanup;. More...

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

virtual int	size () const =0
	Get the szie expected of arrays. More...

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

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

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

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

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

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

virtual void	derivLogDetEvenEvenLinOp (P &ds_u, enum PlusMinus isign) const
	Get the force from the EvenEven Trace Log. More...

virtual Double	logDetEvenEvenLinOp (void) const =0
	Get the log det of the even even part. More...

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

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

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

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

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

virtual void	oddOddLinOp (multi1d< T > &chi, const multi1d< T > &psi, enum PlusMinus isign) const =0
	Apply the the odd-odd block onto a source std::vector. 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...

Detailed Description

template<typename T, typename P, typename Q>
class Chroma::EvenOddPrecLogDetLinearOperatorArray< T, P, Q >

Even-odd preconditioned 5D linear operator.

Support for even-odd preconditioned linear operators Given a matrix M written in block form:

 [      A             D        ]
 [       E,E           E,O     ]
 [                             ]
 [      D             A        ]
 [       O,E           O,O     ]

The preconditioning consists of using the triangular matrices

 [      1              0        ]
 [       E,E            E,O     ]

L = [ ] [ D A^(-1) 1 ] [ O,E E,E O,O ]

and

 [      1            A^{-1}  D       ]
 [       E,E          EE      E,O    ]

U = [ ] [ 0 1 ] [ O,E O,O ]

The preconditioned matrix is formed from

M' = L^-1 * M * U^-1

where

      [      1              0        ]
      [       E,E            E,O     ]

L^(-1) = [ ] [ - D A^(-1) 1 ] [ O,E E,E O,O ]

and

      [      1            - A^(-1) D       ]
      [       E,E            E,E    E,O    ]

U^(-1) = [ ] [ 0 1 ] [ O,E O,O ]

Resulting in a new M

 [      A                    0                      ]
 [       E,E                  E,O                   ]
 [                                                  ]
 [      0                A     -  D    A^(-1)  D    ]
 [       O,E              O,O      O,E   E,E    E,O ]

This class is used to implement the resulting linear operator

 ~
 M  =  A(o,o) - D(o,e) . A^-1(e,e) . D(e,o)

where A^{-1}_{ee} does depend on the gauge fields but we can exactly simulate the determinant without pseudofermion fields since we know det A_{ee} and can write it in the action as exp( log det A ) = exp( Tr Ln A )

Since we can explicitly evaluate Tr Ln A for the action, we can also evaluate the force contribution

d/dt ( Tr Ln A ) = Tr ( (d/dt)A A^{-1} )

and d/dt ( A^{-1} ) = A^{-1} [ (d/dt)A ] A^{-1}

hence we have functions lnDetEvenEven() and derivEvenEven() and derivEvenEvenInv() — this needs work, may not even need more than lnDetEvenEven

Definition at line 395 of file eoprec_logdet_linop.h.

Constructor & Destructor Documentation

◆ ~EvenOddPrecLogDetLinearOperatorArray()

template<typename T , typename P , typename Q >

virtual Chroma::EvenOddPrecLogDetLinearOperatorArray< T, P, Q >::~EvenOddPrecLogDetLinearOperatorArray ( )

inlinevirtual

Virtual destructor to help with cleanup;.

Definition at line 399 of file eoprec_logdet_linop.h.

Member Function Documentation

◆ deriv()

template<typename T , typename P , typename Q >

virtual void Chroma::EvenOddPrecLogDetLinearOperatorArray< T, P, Q >::deriv	(	P &	ds_u,
		const multi1d< T > &	chi,
		const multi1d< T > &	psi,
		enum PlusMinus	isign
	)		const

inlinevirtual

Apply the operator onto a source std::vector.

User should make sure deriv routines do a resize

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

Definition at line 412 of file eoprec_logdet_linop.h.

References chi, Chroma::isign, Chroma::MINUS, Chroma::PLUS, psi, tmp2, and tmp3.

◆ derivEvenEvenLinOp()

template<typename T , typename P , typename Q >

virtual void Chroma::EvenOddPrecLogDetLinearOperatorArray< T, P, Q >::derivEvenEvenLinOp	(	P &	ds_u,
		const multi1d< T > &	chi,
		const multi1d< T > &	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 459 of file eoprec_logdet_linop.h.

◆ derivEvenOddLinOp()

template<typename T , typename P , typename Q >

virtual void Chroma::EvenOddPrecLogDetLinearOperatorArray< T, P, Q >::derivEvenOddLinOp	(	P &	ds_u,
		const multi1d< T > &	chi,
		const multi1d< T > &	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 467 of file eoprec_logdet_linop.h.

◆ derivLogDetEvenEvenLinOp()

template<typename T , typename P , typename Q >

virtual void Chroma::EvenOddPrecLogDetLinearOperatorArray< T, P, Q >::derivLogDetEvenEvenLinOp	(	P &	ds_u,
		enum PlusMinus	isign
	)		const

inlinevirtual

Get the force from the EvenEven Trace Log.

Definition at line 491 of file eoprec_logdet_linop.h.

◆ derivOddEvenLinOp()

template<typename T , typename P , typename Q >

virtual void Chroma::EvenOddPrecLogDetLinearOperatorArray< T, P, Q >::derivOddEvenLinOp	(	P &	ds_u,
		const multi1d< T > &	chi,
		const multi1d< T > &	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 475 of file eoprec_logdet_linop.h.

◆ derivOddOddLinOp()

template<typename T , typename P , typename Q >

virtual void Chroma::EvenOddPrecLogDetLinearOperatorArray< T, P, Q >::derivOddOddLinOp	(	P &	ds_u,
		const multi1d< T > &	chi,
		const multi1d< T > &	psi,
		enum PlusMinus	isign
	)		const

inlinevirtual

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

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

Definition at line 483 of file eoprec_logdet_linop.h.

◆ getFermBC()

template<typename T , typename P , typename Q >

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

pure virtual

Return the fermion BC object for this linear operator.

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

◆ logDetEvenEvenLinOp()

template<typename T , typename P , typename Q >

virtual Double Chroma::EvenOddPrecLogDetLinearOperatorArray< T, P, Q >::logDetEvenEvenLinOp ( void ) const

pure virtual

Get the log det of the even even part.

◆ size()

template<typename T , typename P , typename Q >

virtual int Chroma::EvenOddPrecLogDetLinearOperatorArray< T, P, Q >::size ( ) const

pure virtual

Get the szie expected of arrays.

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

◆ subset()

template<typename T , typename P , typename Q >

const Subset& Chroma::EvenOddPrecLogDetLinearOperatorArray< T, P, Q >::subset ( void ) const

inlinevirtual

Only defined on the odd lattice.

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

Definition at line 402 of file eoprec_logdet_linop.h.

The documentation for this class was generated from the following file:

eoprec_logdet_linop.h

Public Member Functions

Detailed Description

template<typename T, typename P, typename Q> class Chroma::EvenOddPrecLogDetLinearOperatorArray< T, P, Q >

Constructor & Destructor Documentation

◆ ~EvenOddPrecLogDetLinearOperatorArray()

Member Function Documentation

◆ deriv()

◆ derivEvenEvenLinOp()

◆ derivEvenOddLinOp()

◆ derivLogDetEvenEvenLinOp()

◆ derivOddEvenLinOp()

◆ derivOddOddLinOp()

◆ getFermBC()

◆ logDetEvenEvenLinOp()

◆ size()

◆ subset()

template<typename T, typename P, typename Q>
class Chroma::EvenOddPrecLogDetLinearOperatorArray< T, P, Q >