subset_assignment.hpp

Go to the documentation of this file.

#ifndef CHIHAYA_SUBSET_ASSIGNMENT_HPP
#define CHIHAYA_SUBSET_ASSIGNMENT_HPP
 
#include "H5Cpp.h"
#include "ritsuko/ritsuko.hpp"
 
#include <stdexcept>
#include <exception>
#include <vector>
#include <algorithm>
 
#include "utils_public.hpp"
#include "utils_list.hpp"
#include "utils_misc.hpp"
#include "utils_subset.hpp"
#include "utils_dimensions.hpp"
 
namespace chihaya {
 
inline ArrayDetails validate_subset_assignment(const H5::Group& group, const ritsuko::Version& version, const Options& options) {
    ArrayDetails seed_details;
    try {
        auto shandle = group.openGroup("seed");
        seed_details = validate(shandle, version, options);
    } catch (...) {
        std::throw_with_nested(std::runtime_error("failed to validate 'seed'"));
    }
 
    ArrayDetails value_details;
    try {
        auto vhandle = group.openGroup("value");
        value_details = validate(vhandle, version, options);
    } catch (...) {
        std::throw_with_nested(std::runtime_error("failed to validate 'value'"));
    }
 
    if (!options.details_only) {
        if ((value_details.type == STRING) != (seed_details.type == STRING)) {
            throw std::runtime_error("both or neither of the 'seed' and 'value' arrays should contain strings");
        }
 
        const auto& seed_dims = seed_details.dimensions;
        if (seed_dims.size() != value_details.dimensions.size()) {
            throw std::runtime_error("'seed' and 'value' arrays should have the same dimensionality");
        }
 
        std::vector<std::pair<std::size_t, std::size_t> > collected;
        try {
            auto ihandle = group.openGroup("index");
            collected = validate_subset_index_list(ihandle, seed_dims, version, options.contiguous_chunk_size);
        } catch (...) {
            std::throw_with_nested(std::runtime_error("failed to validate 'index'"));
        }
 
        auto expected_dims = seed_dims;
        for (auto p : collected) {
            expected_dims[p.first] = p.second;
        }
 
        if (!are_dimensions_equal(expected_dims, value_details.dimensions)) {
            throw std::runtime_error("'value' dimension extents are not consistent with lengths of indices in 'index'");
        }
    }
 
    // Promotion.
    seed_details.type = std::max(seed_details.type, value_details.type);
    return seed_details;
}
 
}
 
#endif