combine.hpp

Go to the documentation of this file.

#ifndef CHIHAYA_COMBINE_HPP
#define CHIHAYA_COMBINE_HPP
 
#include "H5Cpp.h"
#include "ritsuko/ritsuko.hpp"
#include "ritsuko/hdf5/hdf5.hpp"
 
#include <stdexcept>
#include <vector>
#include <string>
#include <cstdint>
 
#include "utils_list.hpp"
#include "utils_misc.hpp"
 
namespace chihaya {
 
inline ArrayDetails validate(const H5::Group&, const ritsuko::Version&, Options&);
namespace combine {
 
inline ArrayDetails validate(const H5::Group& handle, const ritsuko::Version& version, Options& options) {
    uint64_t along = internal_misc::load_along(handle, version);
 
    auto shandle = ritsuko::hdf5::open_group(handle, "seeds");
    internal_list::ListDetails list_params;
    try {
        list_params = internal_list::validate(shandle, version);
    } catch (std::exception& e) {
        throw std::runtime_error(std::string("failed to load 'seeds' list; ") + e.what());
    }
    if (list_params.present.size() != list_params.length) {
        throw std::runtime_error("missing elements in the 'seeds' list");
    }
 
    std::vector<size_t> dimensions;
    ArrayType type = BOOLEAN;
    {
        bool first = true;
        size_t num_strings = 0;
 
        for (const auto& p : list_params.present) {
            auto current = ritsuko::hdf5::open_group(shandle, p.second.c_str());
 
            ArrayDetails cur_seed;
            try {
                cur_seed = ::chihaya::validate(current, version, options);
            } catch (std::exception& e) {
                throw std::runtime_error("failed to validate 'seeds/" + p.second + "'; " + std::string(e.what()));
            }
 
            if (first) {
                type = cur_seed.type;
                dimensions = cur_seed.dimensions;
                if (static_cast<size_t>(along) >= dimensions.size()) {
                    throw std::runtime_error("'along' should be less than the seed dimensionality");
                }
                first = false;
            } else {
                if (type < cur_seed.type) {
                    type = cur_seed.type;
                }
                if (dimensions.size() != cur_seed.dimensions.size()) {
                    throw std::runtime_error("dimensionality mismatch between seeds");
                }
                for (size_t d = 0; d < dimensions.size(); ++d) {
                    if (d == static_cast<size_t>(along)) {
                        dimensions[d] += cur_seed.dimensions[d];
                    } else if (dimensions[d] != cur_seed.dimensions[d]) {
                        throw std::runtime_error("inconsistent dimension extents between seeds");
                    }
                }
            }
 
            num_strings += (cur_seed.type == STRING);
        }
 
        if (num_strings != 0 && num_strings != list_params.length) {
            throw std::runtime_error("either none or all of the arrays to be combined should contain strings");
        }
    }
 
    return ArrayDetails(type, std::move(dimensions));
}
 
}
 
}
 
#endif