CHROMA
eoprec_constdet_one_flavor_ratio_rat_conv_monomial5d_w.cc
Go to the documentation of this file.
1 /*! @file
2  * @brief One-flavor collection of even-odd preconditioned 5D ferm monomials
3  */
4 
7 
10 
13 
14 namespace Chroma
15 {
16 
17  namespace EvenOddPrecConstDetOneFlavorWilsonTypeFermRatioRatConvMonomial5DEnv
18  {
19  namespace
20  {
21  //! Callback
23  multi1d<LatticeColorMatrix> >* createMonomial(XMLReader& xml, const std::string& path)
24  {
27  }
28 
29  //! Local registration flag
30  bool registered = false;
31  }
32 
33  const std::string name("ONE_FLAVOR_EOPREC_CONSTDET_FERM_RATIO_RAT_CONV_MONOMIAL5D");
34 
35  //! Register all the factories
36  bool registerAll()
37  {
38  bool success = true;
39  if (! registered)
40  {
42  success &= TheMonomialFactory::Instance().registerObject(name, createMonomial);
43  registered = true;
44  }
45  return success;
46  }
47  } //end namespace EvenOddPrecConstDet OneFlavorWilsonFermRatioRatConvMonomialEnv
48 
49 
50  // Constructor
53  {
54  START_CODE();
55 
56  QDPIO::cout << "Constructor: " << __func__ << std::endl;
57 
60  num_pf = param.num_pf;
61 
62  //*********************************************************************
63  // Fermion action
64  {
65  std::istringstream is(param.numer.fermact.xml);
66  XMLReader fermact_reader(is);
67  QDPIO::cout << "Construct fermion action= " << param.numer.fermact.id << std::endl;
68 
69  WilsonTypeFermAct5D<T,P,Q>* tmp_act =
70  TheWilsonTypeFermAct5DFactory::Instance().createObject(param.numer.fermact.id, fermact_reader, param.numer.fermact.path);
71 
73  dynamic_cast<EvenOddPrecConstDetWilsonTypeFermAct5D<T,P,Q>*>(tmp_act);
74 
75  // Check success of the downcast
76  if( downcast == 0x0 ) {
77  QDPIO::cerr << __func__ << ": unable to downcast FermAct to EvenOddPrecConstDetWilsonTypeFermAct5D" << std::endl;
78  QDP_abort(1);
79  }
80 
81  fermact_num = downcast;
82  }
83 
84  //*********************************************************************
85  // Action rational approx
86  {
87  std::istringstream is(param.numer.action.ratApprox.xml);
88  XMLReader approx_reader(is);
89  QDPIO::cout << "Construct action rational approx= " << param.numer.action.ratApprox.id << std::endl;
90 
92  param.numer.action.ratApprox.id,
93  approx_reader,
94  param.numer.action.ratApprox.path));
95 
96  (*approx)(spfe_num, sipfe_num);
97  }
98 
99  //*********************************************************************
100  // Force rational approx
101  {
102  std::istringstream is(param.numer.force.ratApprox.xml);
103  XMLReader approx_reader(is);
104  QDPIO::cout << "Construct force rational approx= " << param.numer.force.ratApprox.id << std::endl;
105 
107  param.numer.force.ratApprox.id,
108  approx_reader,
109  param.numer.force.ratApprox.path));
110 
111  RemezCoeff_t fipfe_num; // discard
112  (*approx)(fpfe_num, fipfe_num);
113  }
114  //*********************************************************************
115 
116  //*********************************************************************
117  // Fermion action
118  {
119  std::istringstream is(param.denom.fermact.xml);
120  XMLReader fermact_reader(is);
121  QDPIO::cout << "Construct fermion action= " << param.denom.fermact.id << std::endl;
122 
123  WilsonTypeFermAct5D<T,P,Q>* tmp_act =
124  TheWilsonTypeFermAct5DFactory::Instance().createObject(param.denom.fermact.id, fermact_reader, param.denom.fermact.path);
125 
127  dynamic_cast<EvenOddPrecConstDetWilsonTypeFermAct5D<T,P,Q>*>(tmp_act);
128 
129  // Check success of the downcast
130  if( downcast == 0x0 ) {
131  QDPIO::cerr << __func__ << ": unable to downcast FermAct to EvenOddPrecConstDetWilsonTypeFermAct5D" << std::endl;
132  QDP_abort(1);
133  }
134 
135  fermact_den = downcast;
136  }
137  //*********************************************************************
138 
139  QDPIO::cout << "Finished constructing: " << __func__ << std::endl;
140 
141  END_CODE();
142  }
143 
144 } //end namespace Chroma
145 
146 
Even-odd preconditioned Wilson-like fermion actions including derivatives.
Class for counted reference semantics.
Definition: handle.h:33
An abstract monomial class, for inexact algorithms.
Definition: abs_monomial.h:43
static T & Instance()
Definition: singleton.h:432
Wilson-like fermion actions.
Definition: fermact.orig.h:403
One-flavor collection of even-odd preconditioned 5D ferm monomials.
Fermion action factories.
All Wilson-type fermion actions.
Monomial factories.
static bool registered
Local registration flag.
const std::string name("ONE_FLAVOR_EOPREC_CONSTDET_FERM_RATIO_RAT_CONV_MONOMIAL5D")
Monomial< multi1d< LatticeColorMatrix >, multi1d< LatticeColorMatrix > > * createMonomial(XMLReader &xml, const std::string &path)
Callback function for the factory.
bool registerAll()
Register all the factories.
Asqtad Staggered-Dirac operator.
Definition: klein_gord.cc:10
START_CODE()
::std::string string
Definition: gtest.h:1979
Rational approximation aggregator.
Rational approximation factories.
GroupXML_t fermact
Definition: comp_approx.h:53
GroupXML_t fermact
Definition: comp_approx.h:36
TermApprox_t action
Definition: comp_approx.h:37
TermApprox_t force
Definition: comp_approx.h:38
Convenient structure to package Remez coeffs.
Definition: remez_coeff.h:19
GroupXML_t invParam
Definition: comp_approx.h:21
GroupXML_t ratApprox
Definition: comp_approx.h:20