portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2024-11-10 07:55:35 +00:00
Code Releases Activity

enseble can now be specified in LapEvec

Browse Source
This commit is contained in:
ferben 2019-02-13 13:54:31 +00:00
parent 48ec937c55
commit 9f2ca98dfc
2 changed files with 16 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

17
Hadrons/Modules/MDistil/LapEvec.hpp
View File

@ -84,10 +84,12 @@ struct LanczosParameters: Serializable {
class LapEvecPar: Serializable
{
public:
GRID_SERIALIZABLE_CLASS_MEMBERS(LapEvecPar
,std::string, gauge
GRID_SERIALIZABLE_CLASS_MEMBERS(LapEvecPar,
std::string, gauge,
std::string, ConfigFileDir,
std::string, ConfigFileName,
//,std::string, EigenPackName
,StoutParameters, Stout
StoutParameters, Stout
,ChebyshevParameters, Cheby
,LanczosParameters, Lanczos
//,DistilParameters, Distil
@ -229,10 +231,14 @@ void TLapEvec<GImpl>::execute(void)
//else
//assert(nnoise>1);
const std::string &ConfigFileDir{par().ConfigFileDir};
const std::string &ConfigFileName{par().ConfigFileName};
// Debugging only
//envGetTmp(GaugeField, Umu);
auto &Umu = envGet(GaugeField, par().gauge);
envGetTmp(GaugeField, Umu_smear);
FieldMetaData header;
if((1)) {
const std::vector<int> seeds({1, 2, 3, 4, 5});
GridParallelRNG pRNG4d(gridHD);
@ -258,6 +264,11 @@ void TLapEvec<GImpl>::execute(void)
Umu = where(coor==t,Umu_smear,Umu);
}
// std::cout << "Umu is "<<Umu<<std::endl;
} else {
std::string fileName(ConfigFileDir + ConfigFileName);
std::cout << GridLogMessage << "Loading NERSC configuration from '" << fileName << "'" << std::endl;
NerscIO::readConfiguration(Umu, header, fileName);
std::cout << GridLogMessage << "reading done." << std::endl;
}
// Stout smearing

2
tests/hadrons/Test_hadrons_distil.cc
View File

@ -60,6 +60,8 @@ void test_LapEvec(Application &application)
application.createModule<MGauge::Random>(szGaugeName);
// Now make an instance of the LapEvec object
MDistil::LapEvecPar p;
p.ConfigFileDir="/home/dp008/dp008/dc-rich6/Scripts/ConfigsDeflQED/";
p.ConfigFileName="ckpoint_lat.3000";
p.gauge = szGaugeName;
//p.EigenPackName = "ePack";
//p.Distil.TI = 8;