portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-11-03 21:44:33 +00:00
Code Releases Activity

Added trajectory number to output file

Browse Source
This commit is contained in:
Michael Marshall
2019-01-25 13:26:48 +00:00
parent 2568f5b925
commit d15bf4b8e1
2 changed files with 11 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files

1
Hadrons/Modules/MDistil/Distil.hpp
View File

@@ -36,6 +36,7 @@
/******************************************************************************
This potentially belongs in CartesianCommunicator
Turns out I don't actually need this when running inside hadrons
******************************************************************************/
BEGIN_MODULE_NAMESPACE(Grid)

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

@@ -315,6 +315,8 @@ void TLapEvec<FImpl>::execute(void)
////////////////////////////////////////////////////////////////////////
std::string sEigenPackName(getName());
sEigenPackName.append("_");
sEigenPackName.append(std::to_string(vm().getTrajectory()));
bool bReturnValue = true;
auto & eig4d = envGet(DistilEP, getName() );
envGetTmp(std::vector<DistilEP>, eig); // Eigenpack for each timeslice
@@ -369,12 +371,14 @@ void TLapEvec<FImpl>::execute(void)
// Now rotate the eigenvectors into our phase convention
RotateEigen( eig[t].evec );
// Write the eigenvectors and eigenvalues to disk
//std::cout << GridLogMessage << "Writing eigenvalues/vectors to " << pszEigenPack << std::endl;
eig[t].record.operatorXml = DefaultOperatorXml;
eig[t].record.solverXml = DefaultsolverXml;
eig[t].write(sEigenPackName,false,t);
//std::cout << GridLogMessage << "Written eigenvectors" << std::endl;
if((1)) { // Debugging only
// Write the eigenvectors and eigenvalues to disk
//std::cout << GridLogMessage << "Writing eigenvalues/vectors to " << pszEigenPack << std::endl;
eig[t].record.operatorXml = DefaultOperatorXml;
eig[t].record.solverXml = DefaultsolverXml;
eig[t].write(sEigenPackName,false,t);
//std::cout << GridLogMessage << "Written eigenvectors" << std::endl;
}
}
for (int i=0;i<LPar.Nvec;i++){
std::cout << "Inserting Timeslice " << t << " into vector " << i << std::endl;