CHROMA
eoprec_parwilson_fermact_w.h
Go to the documentation of this file.
1 // -*- C++ -*-
2 /*! \file
3  * \brief Even-odd preconditioned Wilson fermion action with parity breaking term
4  */
5 
6 #ifndef __prec_parwilson_fermact_w_h__
7 #define __prec_parwilson_fermact_w_h__
8 
10 
11 namespace Chroma
12 {
13  //! Name and registration
14  namespace EvenOddPrecParWilsonFermActEnv
15  {
16  extern const std::string name;
17  bool registerAll();
18  }
19 
20 
21  //! Params for wilson ferm acts
23  {
25  EvenOddPrecParWilsonFermActParams(XMLReader& in, const std::string& path);
26 
27  Real Mass;
28  Real H;
29  };
30 
31 
32  // Reader/writers
33  void read(XMLReader& xml, const std::string& path, EvenOddPrecParWilsonFermActParams& param);
34  void write(XMLWriter& xml, const std::string& path, const EvenOddPrecParWilsonFermActParams& param);
35 
36 
37  //! Even-odd preconditioned Wilson fermion action with parity breaking term
38  /*! \ingroup fermacts
39  *
40  * Even-odd preconditioned wilson fermion action with parity breaking term
41  * Only defined on odd subset.
42  *
43  * The kernel for Wilson fermions with a parity breaking term is
44  *
45  * M = (d+M) + i*H*gamma_5 - (1/2) D'
46  */
48  multi1d<LatticeColorMatrix>, multi1d<LatticeColorMatrix> >
49  {
50  public:
51  // Typedefs to save typing
52  typedef LatticeFermion T;
53  typedef multi1d<LatticeColorMatrix> P;
54  typedef multi1d<LatticeColorMatrix> Q;
55 
56  //! General FermBC
58  const EvenOddPrecParWilsonFermActParams& param_) :
59  fs(fs_), param(param_) {}
60 
61  //! Copy constructor
63  fs(a.fs), param(a.param) {}
64 
65  //! Produce a linear operator for this action
67 
68  //! Produce the gamma_5 hermitin op gamma_5 M
70  {
71  QDP_error_exit("gamma5HermLinOp not implemented yet for this action\n");
72  return 0;
73  }
74 
75  //! Destructor is automatic
77 
78  protected:
79  //! Return the fermion create state for this action
80  const CreateFermState<T,P,Q>& getCreateState() const {return *fs;}
81 
82  //! Partial constructor
84  //! Hide =
86 
87  private:
90  };
91 
92 }
93 
94 
95 #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 Wilson fermion action with parity breaking term.
void operator=(const EvenOddPrecParWilsonFermAct &a)
Hide =.
Handle< CreateFermState< T, P, Q > > fs
const CreateFermState< T, P, Q > & getCreateState() const
Return the fermion create state for this action.
EvenOddPrecConstDetLinearOperator< T, P, Q > * linOp(Handle< FermState< T, P, Q > > state) const
Produce a linear operator for this action.
EvenOddPrecParWilsonFermAct(const EvenOddPrecParWilsonFermAct &a)
Copy constructor.
LinearOperator< T > * hermitianLinOp(Handle< FermState< T, P, Q > > state) const
Produce the gamma_5 hermitin op gamma_5 M.
~EvenOddPrecParWilsonFermAct()
Destructor is automatic.
EvenOddPrecParWilsonFermActParams param
EvenOddPrecParWilsonFermAct(Handle< CreateFermState< T, P, Q > > fs_, const EvenOddPrecParWilsonFermActParams &param_)
General FermBC.
Support class for fermion actions and linear operators.
Definition: state.h:94
Class for counted reference semantics.
Definition: handle.h:33
Linear Operator.
Definition: linearop.h:27
Even-odd const determinant Wilson-like fermact.
void read(XMLReader &xml, const std::string &path, AsqtadFermActParams &param)
Read parameters.
void write(XMLWriter &xml, const std::string &path, const AsqtadFermActParams &param)
Writer parameters.
bool registerAll()
Register all the factories.
Asqtad Staggered-Dirac operator.
Definition: klein_gord.cc:10
QDP_error_exit("too many BiCG iterations", n_count, rsd_sq, cp, c, re_rvr, im_rvr, re_a, im_a, re_b, im_b)
Complex a
Definition: invbicg.cc:95
const WilsonTypeFermAct< multi1d< LatticeFermion > > Handle< const ConnectState > state
Definition: pbg5p_w.cc:28
static QDP_ColorVector * in
::std::string string
Definition: gtest.h:1979