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

unique string is now used

Browse Source
This commit is contained in:
ferben 2019-02-11 17:39:42 +00:00
parent dff7d9261d
commit d889cebc60
2 changed files with 21 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -50,6 +50,9 @@ public:
GRID_SERIALIZABLE_CLASS_MEMBERS(PerambLightPar, GRID_SERIALIZABLE_CLASS_MEMBERS(PerambLightPar,
std::string, eigenPack, std::string, eigenPack,
std::string, PerambFileName, std::string, PerambFileName,
std::string, ConfigFileDir,
std::string, ConfigFileName,
std::string, UniqueIdentifier,
bool, multiFile, bool, multiFile,
int, nvec, int, nvec,
int, Ls, // For makeFiveDimGrid int, Ls, // For makeFiveDimGrid
@ -189,10 +192,17 @@ void TPerambLight<FImpl>::execute(void)
const int Nt_inv{Distil.Nt_inv}; // TODO: PROBABLY BETTER: if (full_tdil) Nt_inv=1; else Nt_inv = TI; const int Nt_inv{Distil.Nt_inv}; // TODO: PROBABLY BETTER: if (full_tdil) Nt_inv=1; else Nt_inv = TI;
const int tsrc{Distil.tsrc}; const int tsrc{Distil.tsrc};
const int Ns{Distil.Ns}; const int Ns{Distil.Ns};
const Real mass{Solver.mass};
const Real M5 {Solver.M5};
const bool full_tdil{TI==Nt}; const bool full_tdil{TI==Nt};
const bool exact_distillation{full_tdil && LI==nvec}; const bool exact_distillation{full_tdil && LI==nvec};
const std::string &ConfigFileDir{par().ConfigFileDir};
const std::string &ConfigFileName{par().ConfigFileName};
const std::string &UniqueIdentifier{par().UniqueIdentifier};
//auto &noise = envGet(std::vector<std::vector<std::vector<SpinVector>>>, par().noise); //auto &noise = envGet(std::vector<std::vector<std::vector<SpinVector>>>, par().noise);
auto &noise = envGet(std::vector<Complex>, getName() + "_noise"); auto &noise = envGet(std::vector<Complex>, getName() + "_noise");
auto &perambulator = envGet(Perambulator<SpinVector COMMA 6 COMMA sizeof(Real)>, auto &perambulator = envGet(Perambulator<SpinVector COMMA 6 COMMA sizeof(Real)>,
@ -229,7 +239,8 @@ void TPerambLight<FImpl>::execute(void)
Umu = where(coor==t,Usft,Umu); Umu = where(coor==t,Usft,Umu);
} }
} else { } else {
std::string fileName( "/home/dp008/dp008/dc-rich6/Scripts/ConfigsDeflQED/ckpoint_lat.3000" ); //std::string fileName( "/home/dp008/dp008/dc-rich6/Scripts/ConfigsDeflQED/ckpoint_lat.3000" );
std::string fileName(ConfigFileDir + ConfigFileName);
std::cout << GridLogMessage << "Loading NERSC configuration from '" << fileName << "'" << std::endl; std::cout << GridLogMessage << "Loading NERSC configuration from '" << fileName << "'" << std::endl;
NerscIO::readConfiguration(Umu, header, fileName); NerscIO::readConfiguration(Umu, header, fileName);
std::cout << GridLogMessage << "reading done." << std::endl; std::cout << GridLogMessage << "reading done." << std::endl;
@ -238,7 +249,8 @@ void TPerambLight<FImpl>::execute(void)
//Create Noises //Create Noises
//std::cout << pszGaugeConfigFile << std::endl; //std::cout << pszGaugeConfigFile << std::endl;
//GridSerialRNG sRNG; sRNG.SeedUniqueString(std::string(pszGaugeConfigFile)); //GridSerialRNG sRNG; sRNG.SeedUniqueString(std::string(pszGaugeConfigFile));
GridSerialRNG sRNG; sRNG.SeedUniqueString("unique_string"); // TODO: Proper unique string. Include quark mass, gauge field? Maybe also nvec, but in a way that more nvec would only add noises, not change all of them??? GridSerialRNG sRNG;
sRNG.SeedUniqueString(ConfigFileName + "_" + std::to_string(mass) + "_" + UniqueIdentifier);
Real rn; Real rn;
for (int inoise=0;inoise<nnoise;inoise++) { for (int inoise=0;inoise<nnoise;inoise++) {
@ -289,8 +301,6 @@ void TPerambLight<FImpl>::execute(void)
const int Ntlocal{grid4d->LocalDimensions()[3]}; const int Ntlocal{grid4d->LocalDimensions()[3]};
const int Ntfirst{grid4d->LocalStarts()[3]}; const int Ntfirst{grid4d->LocalStarts()[3]};
const Real mass{Solver.mass};
const Real M5 {Solver.M5};
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;

6
tests/hadrons/Test_hadrons_distil.cc
View File

@ -89,6 +89,9 @@ void test_Perambulators(Application &application)
MDistil::PerambLight::Par PerambPar; MDistil::PerambLight::Par PerambPar;
PerambPar.eigenPack="LapEvec"; PerambPar.eigenPack="LapEvec";
PerambPar.PerambFileName="peramb.bin"; PerambPar.PerambFileName="peramb.bin";
PerambPar.ConfigFileDir="/home/dp008/dp008/dc-rich6/Scripts/ConfigsDeflQED/";
PerambPar.ConfigFileName="ckpoint_lat.3000";
PerambPar.UniqueIdentifier="full_dilution";
PerambPar.Distil.tsrc = 0; PerambPar.Distil.tsrc = 0;
PerambPar.Distil.nnoise = 1; PerambPar.Distil.nnoise = 1;
PerambPar.Distil.LI=5; PerambPar.Distil.LI=5;
@ -133,6 +136,9 @@ void test_PerambulatorsS(Application &application)
MDistil::PerambLight::Par PerambPar; MDistil::PerambLight::Par PerambPar;
PerambPar.eigenPack="LapEvec"; PerambPar.eigenPack="LapEvec";
PerambPar.PerambFileName="perambS.bin"; PerambPar.PerambFileName="perambS.bin";
PerambPar.ConfigFileDir="/home/dp008/dp008/dc-rich6/Scripts/ConfigsDeflQED/";
PerambPar.ConfigFileName="ckpoint_lat.3000";
PerambPar.UniqueIdentifier="full_dilution";
PerambPar.Distil.tsrc = 0; PerambPar.Distil.tsrc = 0;
PerambPar.Distil.nnoise = 1; PerambPar.Distil.nnoise = 1;
PerambPar.Distil.LI=3; PerambPar.Distil.LI=3;