hdf5_utils.hpp

Go to the documentation of this file.

 
#ifndef READ_FIELDS_HPP
#define READ_FIELDS_HPP
 
#include <cassert>
#include <memory>
#include <string>
#include <cstdlib>
 
#include "H5Cpp.h"
 
namespace hdf5Utils {
 
template<typename TH5, typename Trequ>
static Trequ
read_scalar_attr (H5::Group &header, const std::string &name)
{// {{{
    TH5 out;
    auto attr = header.openAttribute(name);
    assert (attr.getDataType().getSize() == sizeof(TH5));
    attr.read(attr.getDataType(), &out);
    return (Trequ)(out);
}// }}}
 
template<typename TH5, typename Trequ>
static Trequ
read_vector_attr (H5::Group &header, const std::string &name, size_t idx)
{// {{{
    auto attr = header.openAttribute(name);
    auto aspace = attr.getSpace();
    hsize_t dim_lenghts[16];
    auto Ndims = aspace.getSimpleExtentDims(dim_lenghts);
    assert(Ndims == 1);
    assert(dim_lenghts[0] > idx);
    assert(attr.getDataType().getSize() == sizeof(TH5));
    TH5 out[dim_lenghts[0]];
    attr.read(attr.getDataType(), out);
    return (Trequ)(out[idx]);
}// }}}
 
template<typename TH5, typename Trequ>
static size_t
read_vector_attr (H5::Group &header, const std::string &name, Trequ *out)
{// {{{
    auto attr = header.openAttribute(name);
    auto aspace = attr.getSpace();
    hsize_t dim_lenghts[16];
    auto Ndims = aspace.getSimpleExtentDims(dim_lenghts);
    assert(Ndims == 1);
    assert(attr.getDataType().getSize() == sizeof(TH5));
    
    if constexpr (std::is_same_v<TH5, Trequ>)
        attr.read(attr.getDataType(), out);
    else
    {
        TH5 out1[dim_lenghts[0]];
        attr.read(attr.getDataType(), out1);
        for (size_t ii=0; ii != dim_lenghts[0]; ++ii)
            out[ii] = (Trequ)out1[ii];
    }
 
    return dim_lenghts[0];
}// }}}
 
} // namespace hdf5Utils
 
#endif // READ_FIELDS_HPP