dense_array.hpp

Go to the documentation of this file.

#ifndef CHIHAYA_DENSE_ARRAY_HPP
#define CHIHAYA_DENSE_ARRAY_HPP
 
#include "H5Cpp.h"
#include "ritsuko/ritsuko.hpp"
#include "ritsuko/hdf5/hdf5.hpp"
 
#include <vector>
#include <cstdint>
 
#include "utils_public.hpp"
#include "utils_type.hpp"
#include "utils_dimnames.hpp"
 
namespace chihaya {
 
namespace dense_array {
 
inline ArrayDetails validate(const H5::Group& handle, const ritsuko::Version& version, [[maybe_unused]] Options& options) {
    ArrayDetails output;
 
    {
        auto dhandle = ritsuko::hdf5::open_dataset(handle, "data");
        auto dspace = dhandle.getSpace();
        auto ndims = dspace.getSimpleExtentNdims();
        if (ndims == 0) {
            throw std::runtime_error("'data' should have non-zero dimensions for a dense array");
        }
 
        std::vector<hsize_t> dims(ndims);
        dspace.getSimpleExtentDims(dims.data());
        output.dimensions.insert(output.dimensions.end(), dims.begin(), dims.end());
 
        try {
            if (version.lt(1, 1, 0)) {
                output.type = internal_type::translate_type_0_0(dhandle.getTypeClass());
                if (internal_type::is_boolean(dhandle)) {
                    output.type = BOOLEAN;
                }
            } else {
                auto type = ritsuko::hdf5::open_and_load_scalar_string_attribute(dhandle, "type");
                output.type = internal_type::translate_type_1_1(type);
                internal_type::check_type_1_1(dhandle, output.type);
            }
 
            if (!options.details_only) {
                internal_misc::validate_missing_placeholder(dhandle, version);
            }
 
            if (dhandle.getTypeClass() == H5T_STRING) {
                ritsuko::hdf5::validate_nd_string_dataset(dhandle, dims, 1000000);
            }
 
        } catch (std::exception& e) {
            throw std::runtime_error("failed to validate 'data'; " + std::string(e.what()));
        }
    }
 
    bool native;
    {
        auto nhandle = ritsuko::hdf5::open_dataset(handle, "native");
        if (!ritsuko::hdf5::is_scalar(nhandle)) {
            throw std::runtime_error("'native' should be a scalar");
        }
 
        if (version.lt(1, 1, 0)) {
            if (nhandle.getTypeClass() != H5T_INTEGER) {
                throw std::runtime_error("'native' should have an integer datatype");
            }
            native = ritsuko::hdf5::load_scalar_numeric_dataset<int>(nhandle);
        } else {
            if (ritsuko::hdf5::exceeds_integer_limit(nhandle, 8, true)) {
                throw std::runtime_error("'native' should have a datatype that fits into an 8-bit signed integer");
            }
            native = ritsuko::hdf5::load_scalar_numeric_dataset<int8_t>(nhandle);
        }
    }
 
    // Do this before the 'native' check.
    if (!options.details_only) {
        if (handle.exists("dimnames")) {
            internal_dimnames::validate(handle, output.dimensions, version);
        }
    }
 
    if (!native) {
        std::reverse(output.dimensions.begin(), output.dimensions.end());
    }
 
    return output;
}
 
}
 
}
 
#endif