CHROMA
eoprec_wilson_coarse_fine_fermact_w.h
Go to the documentation of this file.
1 // -*- C++ -*-
2 /*! \file
3  * \brief Even-odd preconditioned Wilson fermion action supporting 2+2 anisotropy
4  */
5 
6 #ifndef __eoprec_wilson_coarse_fine_fermact_w_h__
7 #define __eoprec_wilson_coarse_fine_fermact_w_h__
8 
12 
13 namespace Chroma
14 {
15  //! Name and registration
16  /*! \ingroup fermacts */
17  namespace EvenOddPrecWilsonCoarseFineFermActEnv
18  {
19  extern const std::string name;
20  bool registerAll();
21  }
22 
23 
24  //! Even-odd preconditioned WilsonCoarseFine fermion action
25  /*! \ingroup fermacts
26  *
27  * Even-odd preconditioned wilson fermion action.
28  * Only defined on odd subset.
29  */
31  public EvenOddPrecConstDetWilsonTypeFermAct<LatticeFermion,
32  multi1d<LatticeColorMatrix>, multi1d<LatticeColorMatrix> >
33  {
34  public:
35  // Typedefs to save typing
36  typedef LatticeFermion T;
37  typedef multi1d<LatticeColorMatrix> P;
38  typedef multi1d<LatticeColorMatrix> Q;
39 
40  //! General FermBC with Anisotropy
43  cfs(cfs_) {init(p);}
44 
45  //! Copy constructor
47  cfs(a.cfs), param(a.param) {}
48 
49  //! Produce a linear operator for this action
51 
52  //! Produce the gamma_5 hermitian operator H_w
54  {
55  return new lgherm<LatticeFermion>(linOp(state));
56  }
57 
58  //! Return a linear operator solver for this action to solve M*psi=chi
60  const GroupXML_t& invParam) const;
61 
62  //! Return a linear operator solver for this action to solve MdagM*psi=chi
64  const GroupXML_t& invParam) const;
65 
66  //! Destructor is automatic
68 
69  Double getQuarkMass(void) const {
70  return param.Mass;
71  }
72 
73  protected:
74  //! Init params
75  void init(const WilsonCoarseFineFermActParams& param_);
76 
77  //! Return the fermion BC object for this action
78  const CreateFermState<T,P,Q>& getCreateState() const {return *cfs;}
79 
80  //! Partial constructor
82  //! Assignment
84 
85  private:
86  struct Param_t
87  {
88  Real Mass;
89  multi1d<Real> coeffs;
90  };
91 
94  };
95 
96 }
97 
98 #endif
Create a fermion connection state.
Definition: create_state.h:69
Even-odd preconditioned linear operator.
Even-odd preconditioned Wilson-like fermion actions specialised to Wilson Like (gauge independent dia...
Even-odd preconditioned WilsonCoarseFine fermion action.
void init(const WilsonCoarseFineFermActParams &param_)
Init params.
EvenOddPrecConstDetLinearOperator< T, P, Q > * linOp(Handle< FermState< T, P, Q > > state) const
Produce a linear operator for this action.
MdagMSystemSolver< T > * invMdagM(Handle< FermState< T, P, Q > > state, const GroupXML_t &invParam) const
Return a linear operator solver for this action to solve MdagM*psi=chi.
LinearOperator< T > * hermitianLinOp(Handle< FermState< T, P, Q > > state) const
Produce the gamma_5 hermitian operator H_w.
LinOpSystemSolver< T > * invLinOp(Handle< FermState< T, P, Q > > state, const GroupXML_t &invParam) const
Return a linear operator solver for this action to solve M*psi=chi.
EvenOddPrecWilsonCoarseFineFermAct(Handle< CreateFermState< T, P, Q > > cfs_, const WilsonCoarseFineFermActParams &p)
General FermBC with Anisotropy.
EvenOddPrecWilsonCoarseFineFermAct(const EvenOddPrecWilsonCoarseFineFermAct &a)
Copy constructor.
void operator=(const EvenOddPrecWilsonCoarseFineFermAct &a)
Assignment.
const CreateFermState< T, P, Q > & getCreateState() const
Return the fermion BC object for this action.
Support class for fermion actions and linear operators.
Definition: state.h:94
Class for counted reference semantics.
Definition: handle.h:33
SystemSolver disambiguator.
Linear Operator.
Definition: linearop.h:27
SystemSolver disambiguator.
Gamma(5) hermitian linear operator.
Definition: lgherm_w.h:22
Even-odd const determinant Wilson-like fermact.
Asqtad Staggered-Dirac operator.
Definition: klein_gord.cc:10
Complex a
Definition: invbicg.cc:95
const WilsonTypeFermAct< multi1d< LatticeFermion > > Handle< const ConnectState > state
Definition: pbg5p_w.cc:28
FloatingPoint< double > Double
Definition: gtest.h:7351
::std::string string
Definition: gtest.h:1979
Hold group xml and type id.
Wilson coarse-fine 2+2 anisotropic lattice fermion action parameters.