portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-04-25 13:15:55 +01:00
Code Releases Activity

Namespace, indent

Browse Source
This commit is contained in:
paboyle 2018-01-14 23:21:40 +00:00
parent 8178a17b88
commit 017dcd69a6
1 changed files with 36 additions and 39 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

75
lib/qcd/action/fermion/MobiusFermion.h
View File

@ -1,4 +1,4 @@
/************************************************************************************* /*************************************************************************************
Grid physics library, www.github.com/paboyle/Grid Grid physics library, www.github.com/paboyle/Grid
@ -24,57 +24,54 @@ Author: Peter Boyle <paboyle@ph.ed.ac.uk>
51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA. 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
See the full license in the file "LICENSE" in the top level distribution directory See the full license in the file "LICENSE" in the top level distribution directory
*************************************************************************************/ *************************************************************************************/
/* END LEGAL */ /* END LEGAL */
#ifndef GRID_QCD_MOBIUS_FERMION_H #ifndef GRID_QCD_MOBIUS_FERMION_H
#define GRID_QCD_MOBIUS_FERMION_H #define GRID_QCD_MOBIUS_FERMION_H
#include <Grid/qcd/action/fermion/FermionCore.h> #include <Grid/qcd/action/fermion/FermionCore.h>
namespace Grid { NAMESPACE_BEGIN(Grid);
namespace QCD { template<class Impl>
class MobiusFermion : public CayleyFermion5D<Impl>
{
public:
INHERIT_IMPL_TYPES(Impl);
public:
template<class Impl> virtual void Instantiatable(void) {};
class MobiusFermion : public CayleyFermion5D<Impl> // Constructors
{ MobiusFermion(GaugeField &_Umu,
public: GridCartesian &FiveDimGrid,
INHERIT_IMPL_TYPES(Impl); GridRedBlackCartesian &FiveDimRedBlackGrid,
public: GridCartesian &FourDimGrid,
GridRedBlackCartesian &FourDimRedBlackGrid,
RealD _mass,RealD _M5,
RealD b, RealD c,const ImplParams &p= ImplParams()) :
virtual void Instantiatable(void) {}; CayleyFermion5D<Impl>(_Umu,
// Constructors FiveDimGrid,
MobiusFermion(GaugeField &_Umu, FiveDimRedBlackGrid,
GridCartesian &FiveDimGrid, FourDimGrid,
GridRedBlackCartesian &FiveDimRedBlackGrid, FourDimRedBlackGrid,_mass,_M5,p)
GridCartesian &FourDimGrid,
GridRedBlackCartesian &FourDimRedBlackGrid,
RealD _mass,RealD _M5,
RealD b, RealD c,const ImplParams &p= ImplParams()) :
CayleyFermion5D<Impl>(_Umu, {
FiveDimGrid, RealD eps = 1.0;
FiveDimRedBlackGrid,
FourDimGrid,
FourDimRedBlackGrid,_mass,_M5,p)
{ std::cout<<GridLogMessage << "MobiusFermion (b="<<b<<",c="<<c<<") with Ls= "<<this->Ls<<" Tanh approx"<<std::endl;
RealD eps = 1.0; Approx::zolotarev_data *zdata = Approx::higham(eps,this->Ls);// eps is ignored for higham
assert(zdata->n==this->Ls);
std::cout<<GridLogMessage << "MobiusFermion (b="<<b<<",c="<<c<<") with Ls= "<<this->Ls<<" Tanh approx"<<std::endl; // Call base setter
Approx::zolotarev_data *zdata = Approx::higham(eps,this->Ls);// eps is ignored for higham this->SetCoefficientsTanh(zdata,b,c);
assert(zdata->n==this->Ls);
// Call base setter Approx::zolotarev_free(zdata);
this->SetCoefficientsTanh(zdata,b,c);
Approx::zolotarev_free(zdata);
}
};
} }
}
};
NAMESPACE_END(Grid);
#endif #endif