portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-04-09 21:50:45 +01:00
Code Releases Activity

Namespace

Browse Source
This commit is contained in:
paboyle 2018-01-14 22:02:09 +00:00
parent 9aa34dc803
commit 8bf78846ee
1 changed files with 12 additions and 14 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

26
lib/qcd/utils/CovariantLaplacian.h
View File

@ -25,13 +25,12 @@ with this program; if not, write to the Free Software Foundation, Inc.,
See the full license in the file "LICENSE" in the top level distribution See the full license in the file "LICENSE" in the top level distribution
directory directory
*************************************************************************************/ *************************************************************************************/
/* END LEGAL */ /* END LEGAL */
#ifndef COVARIANT_LAPLACIAN_H #ifndef COVARIANT_LAPLACIAN_H
#define COVARIANT_LAPLACIAN_H #define COVARIANT_LAPLACIAN_H
namespace Grid { NAMESPACE_BEGIN(Grid);
namespace QCD {
struct LaplacianParams : Serializable { struct LaplacianParams : Serializable {
GRID_SERIALIZABLE_CLASS_MEMBERS(LaplacianParams, GRID_SERIALIZABLE_CLASS_MEMBERS(LaplacianParams,
@ -80,19 +79,19 @@ class LaplacianAdjointField: public Metric<typename Impl::Field> {
MultiShiftFunction PowerHalf; MultiShiftFunction PowerHalf;
MultiShiftFunction PowerInvHalf; MultiShiftFunction PowerInvHalf;
public: public:
INHERIT_GIMPL_TYPES(Impl); INHERIT_GIMPL_TYPES(Impl);
LaplacianAdjointField(GridBase* grid, OperatorFunction<GaugeField>& S, LaplacianParams& p, const RealD k = 1.0) LaplacianAdjointField(GridBase* grid, OperatorFunction<GaugeField>& S, LaplacianParams& p, const RealD k = 1.0)
: U(Nd, grid), Solver(S), param(p), kappa(k){ : U(Nd, grid), Solver(S), param(p), kappa(k){
AlgRemez remez(param.lo,param.hi,param.precision); AlgRemez remez(param.lo,param.hi,param.precision);
std::cout<<GridLogMessage << "Generating degree "<<param.degree<<" for x^(1/2)"<<std::endl; std::cout<<GridLogMessage << "Generating degree "<<param.degree<<" for x^(1/2)"<<std::endl;
remez.generateApprox(param.degree,1,2); remez.generateApprox(param.degree,1,2);
PowerHalf.Init(remez,param.tolerance,false); PowerHalf.Init(remez,param.tolerance,false);
PowerInvHalf.Init(remez,param.tolerance,true); PowerInvHalf.Init(remez,param.tolerance,true);
}; };
void Mdir(const GaugeField&, GaugeField&, int, int){ assert(0);} void Mdir(const GaugeField&, GaugeField&, int, int){ assert(0);}
void Mdiag(const GaugeField&, GaugeField&){ assert(0);} void Mdiag(const GaugeField&, GaugeField&){ assert(0);}
@ -186,12 +185,11 @@ class LaplacianAdjointField: public Metric<typename Impl::Field> {
private: private:
RealD kappa; RealD kappa;
std::vector<GaugeLinkField> U; std::vector<GaugeLinkField> U;
}; };
} NAMESPACE_END(Grid);
}
#endif #endif