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

Hadrons: more verbose eigenpack convert

Browse Source
This commit is contained in:
Antonin Portelli 2018-10-02 12:24:45 +01:00
parent aadd9f4468
commit 4ab8cfbe2a
1 changed files with 17 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
Hadrons/Utilities/EigenPackCast.cc
View File

@ -26,8 +26,14 @@ void convert(const std::string outFilename, const std::string inFilename,
gIn.reset(SpaceTimeGrid::makeFourDimGrid(dim, simdIn, GridDefaultMpi())); gIn.reset(SpaceTimeGrid::makeFourDimGrid(dim, simdIn, GridDefaultMpi()));
FOut bufOut(gOut.get()); FOut bufOut(gOut.get());
FIn bufIn(gIn.get()); FIn bufIn(gIn.get()), testIn(gIn.get());
LOG(Message) << "==== EIGENPACK CONVERSION" << std::endl;
LOG(Message) << "In path : " << inFilename << std::endl;
LOG(Message) << "In type : " << typeName<FIn>() << std::endl;
LOG(Message) << "Out path : " << outFilename << std::endl;
LOG(Message) << "Out type : " << typeName<FOut>() << std::endl;
LOG(Message) << "Multifile: " << (multiFile ? "yes" : "no") << std::endl;
if (multiFile) if (multiFile)
{ {
for(unsigned int k = 0; k < size; ++k) for(unsigned int k = 0; k < size; ++k)
@ -39,6 +45,9 @@ void convert(const std::string outFilename, const std::string inFilename,
std::string outV = outFilename + "/v" + std::to_string(k) + ".bin"; std::string outV = outFilename + "/v" + std::to_string(k) + ".bin";
std::string inV = inFilename + "/v" + std::to_string(k) + ".bin"; std::string inV = inFilename + "/v" + std::to_string(k) + ".bin";
LOG(Message) << "==== Converting vector " << k << std::endl;
LOG(Message) << "In : " << inV << std::endl;
LOG(Message) << "Out: " << outV << std::endl;
makeFileDir(outV, gOut.get()); makeFileDir(outV, gOut.get());
binWriter.open(outV); binWriter.open(outV);
binReader.open(inV); binReader.open(inV);
@ -46,6 +55,9 @@ void convert(const std::string outFilename, const std::string inFilename,
EPOut::writeHeader(binWriter, record); EPOut::writeHeader(binWriter, record);
EPIn::readElement(bufIn, vecRecord, binReader); EPIn::readElement(bufIn, vecRecord, binReader);
precisionChange(bufOut, bufIn); precisionChange(bufOut, bufIn);
precisionChange(testIn, bufOut);
testIn -= bufIn;
LOG(Message) << "Diff norm^2: " << norm2(testIn) << std::endl;
EPOut::writeElement(binWriter, bufOut, vecRecord); EPOut::writeElement(binWriter, bufOut, vecRecord);
binWriter.close(); binWriter.close();
binReader.close(); binReader.close();
@ -66,8 +78,12 @@ void convert(const std::string outFilename, const std::string inFilename,
{ {
VecRecord vecRecord; VecRecord vecRecord;
LOG(Message) << "==== Converting vector " << k << std::endl;
EPIn::readElement(bufIn, vecRecord, binReader); EPIn::readElement(bufIn, vecRecord, binReader);
precisionChange(bufOut, bufIn); precisionChange(bufOut, bufIn);
precisionChange(testIn, bufOut);
testIn -= bufIn;
LOG(Message) << "Diff norm^2: " << norm2(testIn) << std::endl;
EPOut::writeElement(binWriter, bufOut, vecRecord); EPOut::writeElement(binWriter, bufOut, vecRecord);
} }
binWriter.close(); binWriter.close();