portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-04-10 06:00:45 +01:00
Code Releases Activity

LapEvec test works on --grid 4.4.4.8

Browse Source
This commit is contained in:
Michael Marshall 2019-01-24 16:30:13 +00:00
parent cf85f0388d
commit dfb7fb1d9f
1 changed files with 9 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -323,6 +323,8 @@ void TLapEvec<FImpl>::execute(void)
envGetTmp(LatticeColourVector, src);
const int Ntlocal{gridHD->LocalDimensions()[Tdir]};
const int Ntfirst{gridHD->LocalStarts()[Tdir]};
const char DefaultOperatorXml[] = "<OPERATOR>Michael</OPERATOR>";
const char DefaultsolverXml[] = "<SOLVER>Felix</SOLVER>";
for(int t=Ntfirst;bReturnValue && t<Ntfirst+Ntlocal;t++){
std::cout << GridLogMessage << "------------------------------------------------------------" << std::endl;
std::cout << GridLogMessage << " Compute eigenpack, Timeslice = " << t << std::endl;
@ -372,8 +374,8 @@ void TLapEvec<FImpl>::execute(void)
// Write the eigenvectors and eigenvalues to disk
//std::cout << GridLogMessage << "Writing eigenvalues/vectors to " << pszEigenPack << std::endl;
eig[t].record.operatorXml="<OPERATOR>Michael</OPERATOR>";
eig[t].record.solverXml="<SOLVER>Felix</SOLVER>";
eig[t].record.operatorXml = DefaultOperatorXml;
eig[t].record.solverXml = DefaultsolverXml;
eig[t].write(sEigenPackName,false,t);
//std::cout << GridLogMessage << "Written eigenvectors" << std::endl;
}
@ -383,6 +385,11 @@ void TLapEvec<FImpl>::execute(void)
}
}
// Now write out the 4d eigenvectors
eig4d.record.operatorXml = DefaultOperatorXml;
eig4d.record.solverXml = DefaultsolverXml;
eig4d.write(sEigenPackName,false);
// Close the local debugging log file
if( ll ) {
*ll << " Returning " << bReturnValue << std::endl;