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

moved hard-coded parameters to module input

Browse Source
This commit is contained in:
ferben 2019-01-23 11:26:07 +00:00
parent b8afa7314c
commit 4cc2ebc9e4
1 changed files with 31 additions and 15 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

46
Hadrons/Modules/MDistil/PerambLight.hpp
View File

@ -26,7 +26,21 @@ public:
GRID_SERIALIZABLE_CLASS_MEMBERS(PerambLightPar, GRID_SERIALIZABLE_CLASS_MEMBERS(PerambLightPar,
std::string, noise, std::string, noise,
std::string, eigenPack, std::string, eigenPack,
bool, multiFile); bool, multiFile,
int, tsrc,
int, nnoise,
int, LI,
int, SI,
int, TI,
int, nvec,
int, Ns,
int, Nt,
int, Nt_inv,
Real, mass,
Real, M5,
int, Ls,
double, CGPrecision,
int, MaxIterations);
}; };
template <typename FImpl> template <typename FImpl>
@ -155,26 +169,28 @@ void TPerambLight<FImpl>::execute(void)
int Ntlocal = grid4d->LocalDimensions()[3]; int Ntlocal = grid4d->LocalDimensions()[3];
int Ntfirst = grid4d->LocalStarts()[3]; int Ntfirst = grid4d->LocalStarts()[3];
int tsrc=0;
int nnoise=1; int tsrc=par().tsrc;
int LI=6; int nnoise=par().nnoise;
int Ns=4; int LI=par().LI;
int Nt_inv=1; int Ns=par().Ns;
int Nt=64; int Nt_inv=par().Nt_inv;
int TI=64; int Nt=par().Nt;
int nvec=6; int TI=par().TI;
bool full_tdil=true; int nvec=par().nvec;
bool full_tdil=(TI==Nt);
Real mass=0.005; // TODO Infile Real mass=par().mass; // TODO Infile
Real M5 =1.8; // TODO Infile Real M5 =par().M5; // TODO Infile
std::cout << "init RBG " << std::endl; std::cout << "init RBG " << std::endl;
GridRedBlackCartesian RBGrid(grid4d); GridRedBlackCartesian RBGrid(grid4d);
std::cout << "init RBG done" << std::endl; std::cout << "init RBG done" << std::endl;
int Ls=16; int Ls=par().Ls;
double CGPrecision = 10e-8; double CGPrecision = par().CGPrecision;
int MaxIterations = 10000; int MaxIterations = par().MaxIterations;
GridCartesian * FGrid = SpaceTimeGrid::makeFiveDimGrid(Ls,grid4d); GridCartesian * FGrid = SpaceTimeGrid::makeFiveDimGrid(Ls,grid4d);
GridRedBlackCartesian * FrbGrid = SpaceTimeGrid::makeFiveDimRedBlackGrid(Ls,grid4d); GridRedBlackCartesian * FrbGrid = SpaceTimeGrid::makeFiveDimRedBlackGrid(Ls,grid4d);