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

checksummed HDF5 IO

Browse Source
This commit is contained in:
Antonin Portelli 2018-11-27 17:43:19 +00:00
parent c7f33ca2a8
commit a4c55406ed
2 changed files with 5 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
Grid/serialisation/Hdf5IO.h
View File

@ -123,9 +123,11 @@ namespace Grid
if (flatx.size() > dataSetThres_) if (flatx.size() > dataSetThres_)
{ {
H5NS::DataSet dataSet; H5NS::DataSet dataSet;
H5NS::DSetCreatPropList plist;
dataSet = group_.createDataSet(s, Hdf5Type<Element>::type(), dataSpace); plist.setFletcher32();
dataSet = group_.createDataSet(s, Hdf5Type<Element>::type(), dataSpace, plist);
dataSet.write(flatx.data(), Hdf5Type<Element>::type()); dataSet.write(flatx.data(), Hdf5Type<Element>::type());
} }
else else

3
Hadrons/A2AMatrix.hpp
View File

@ -442,6 +442,7 @@ void A2AMatrixIo<T>::initFile(const MetadataType &d, const unsigned int chunkSiz
push(reader, dataname_); push(reader, dataname_);
auto &group = reader.getGroup(); auto &group = reader.getGroup();
plist.setChunk(chunk.size(), chunk.data()); plist.setChunk(chunk.size(), chunk.data());
plist.setFletcher32();
dataset = group.createDataSet(HADRONS_A2AM_NAME, Hdf5Type<T>::type(), dataspace, plist); dataset = group.createDataSet(HADRONS_A2AM_NAME, Hdf5Type<T>::type(), dataspace, plist);
#else #else
HADRONS_ERROR(Implementation, "all-to-all matrix I/O needs HDF5 library"); HADRONS_ERROR(Implementation, "all-to-all matrix I/O needs HDF5 library");
@ -502,14 +503,12 @@ void A2AMatrixIo<T>::load(Vec<VecT> &v, double *tRead)
H5NS::DataSet dataset; H5NS::DataSet dataset;
H5NS::DataSpace dataspace; H5NS::DataSpace dataspace;
H5NS::CompType datatype; H5NS::CompType datatype;
H5NS::DSetCreatPropList plist;
push(reader, dataname_); push(reader, dataname_);
auto &group = reader.getGroup(); auto &group = reader.getGroup();
dataset = group.openDataSet(HADRONS_A2AM_NAME); dataset = group.openDataSet(HADRONS_A2AM_NAME);
datatype = dataset.getCompType(); datatype = dataset.getCompType();
dataspace = dataset.getSpace(); dataspace = dataset.getSpace();
plist = dataset.getCreatePlist();
hdim.resize(dataspace.getSimpleExtentNdims()); hdim.resize(dataspace.getSimpleExtentNdims());
dataspace.getSimpleExtentDims(hdim.data()); dataspace.getSimpleExtentDims(hdim.data());
if ((nt_*ni_*nj_ != 0) and if ((nt_*ni_*nj_ != 0) and