verify.hpp

Go to the documentation of this file.

#ifndef SINGLEPP_LOADERS_VERIFY_HPP
#define SINGLEPP_LOADERS_VERIFY_HPP
 
#include "rankings.hpp"
#include "singlepp/singlepp.hpp"
 
#include <stdexcept>
#include <vector>
 
namespace singlepp_loaders {
 
template<typename Value_, typename Index_, typename Rank_, typename Label_>
void verify(const RankMatrix<Value_, Index_, Rank_>& rankings, const std::vector<Label_>& labels, const singlepp::Markers<Index_>& markers) {
    size_t NC = rankings.ncol();
    if (labels.size() != NC) {
        throw std::runtime_error("number of columns of 'rankings' does not equal length of 'labels'");
    }
 
    size_t nlabels = 0;
    if (NC) {
        nlabels = static_cast<size_t>(*std::max_element(labels.begin(), labels.end())) + 1;
    }
    if (markers.size() != nlabels) {
        throw std::runtime_error("length of 'markers' does not equal the number of labels in 'labels'");
    }
 
    Index_ NR = rankings.nrow();
    for (const auto& mm : markers) {
        if (mm.size() != nlabels) {
            throw std::runtime_error("number of inner vectors in 'markers' does not equal the number of labels in 'labels'");
        }
        for (const auto& m : mm) {
            for (auto x : m) {
                if (x < 0 || x >= NR) {
                    throw std::runtime_error("gene indices in 'markers' should be non-negative and less than the number of rows in 'rankings'");
                }
            }
        }
    }
}
 
}
 
#endif