get_1d_length.hpp

Go to the documentation of this file.

#ifndef RITSUKO_HDF5_GET_1D_LENGTH_HPP
#define RITSUKO_HDF5_GET_1D_LENGTH_HPP
 
#include "H5Cpp.h"
#include <stdexcept>
 
namespace ritsuko {
 
namespace hdf5 {
 
inline hsize_t get_1d_length(const H5::DataSpace& space, bool allow_scalar) {
    int ndims = space.getSimpleExtentNdims();
    if (ndims == 0) {
        if (!allow_scalar) {
            throw std::runtime_error("expected a 1-dimensional dataset, got a scalar instead");
        } 
        return 0;
    }
    if (ndims != 1) {
        throw std::runtime_error("expected a 1-dimensional dataset, got " + std::to_string(ndims) + " dimensions instead");
    }
 
    hsize_t dims;
    space.getSimpleExtentDims(&dims);
    return dims;
}
 
inline hsize_t get_1d_length(const H5::DataSet& handle, bool allow_scalar) {
    return get_1d_length(handle.getSpace(), allow_scalar);
}
 
inline hsize_t get_1d_length(const H5::Attribute& handle, bool allow_scalar) {
    return get_1d_length(handle.getSpace(), allow_scalar);
}
 
inline bool is_scalar(const H5::DataSpace& space) {
    return space.getSimpleExtentNdims() == 0;
}
 
inline bool is_scalar(const H5::DataSet& handle) {
    return is_scalar(handle.getSpace());
}
 
inline bool is_scalar(const H5::Attribute& handle) {
    return is_scalar(handle.getSpace());
}
 
}
 
}
 
#endif