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

Hadrons: multi-file eigenpack in separate directory

Browse Source
This commit is contained in:
Antonin Portelli 2018-05-07 17:52:54 +01:00
parent 2f92721249
commit 27a4d4c951
1 changed files with 19 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

22
extras/Hadrons/EigenPack.hpp
View File

@ -107,9 +107,15 @@ protected:
std::string evecFilename(const std::string stem, const int vec, const int traj) std::string evecFilename(const std::string stem, const int vec, const int traj)
{ {
std::string t = (traj < 0) ? "" : ("." + std::to_string(traj)); std::string t = (traj < 0) ? "" : ("." + std::to_string(traj));
std::string v = (traj < 0) ? "" : (".v" + std::to_string(vec));
return stem + "_evec" + t + v + ".bin"; if (vec == -1)
{
return stem + t + ".bin";
}
else
{
return stem + t + "/v" + std::to_string(vec) + ".bin";
}
} }
template <typename T> template <typename T>
@ -163,8 +169,13 @@ protected:
const std::vector<double> &eval, const unsigned int size) const std::vector<double> &eval, const unsigned int size)
{ {
ScidacWriter binWriter(evec[0]._grid->IsBoss()); ScidacWriter binWriter(evec[0]._grid->IsBoss());
XmlWriter xmlWriter("", "eigenPackPar"); XmlWriter xmlWriter("", "eigenPackPar");
std::string d = dirname(filename);
if (mkdir(d))
{
HADRONS_ERROR(Io, "cannot create directory '" + d + "'");\
}
xmlWriter.pushXmlString(record.operatorXml); xmlWriter.pushXmlString(record.operatorXml);
xmlWriter.pushXmlString(record.solverXml); xmlWriter.pushXmlString(record.solverXml);
binWriter.open(filename); binWriter.open(filename);
@ -188,7 +199,12 @@ protected:
ScidacWriter binWriter(evec._grid->IsBoss()); ScidacWriter binWriter(evec._grid->IsBoss());
XmlWriter xmlWriter("", "eigenPackPar"); XmlWriter xmlWriter("", "eigenPackPar");
VecRecord vecRecord; VecRecord vecRecord;
std::string d = dirname(filename);
if (mkdir(d))
{
HADRONS_ERROR(Io, "cannot create directory '" + d + "'");\
}
xmlWriter.pushXmlString(record.operatorXml); xmlWriter.pushXmlString(record.operatorXml);
xmlWriter.pushXmlString(record.solverXml); xmlWriter.pushXmlString(record.solverXml);
binWriter.open(filename); binWriter.open(filename);