size_grid.hpp

Go to the documentation of this file.

#ifndef FIMS_POPULATION_DYNAMICS_SIZE_GRID_HPP
#define FIMS_POPULATION_DYNAMICS_SIZE_GRID_HPP
 
#include <cstddef>
 
#include "../../common/fims_vector.hpp"
 
namespace fims_popdy {

struct SizeGrid {
  std::size_t n_bins = 0;       
  fims::Vector<double> edges;   
  fims::Vector<double> centers; 

  bool HasStrictlyIncreasingEdges() const {
    if (n_bins == 0) {
      return edges.size() == 0;
    }
 
    if (edges.size() != n_bins + 1) {
      return false;
    }
 
    for (std::size_t i = 1; i < edges.size(); ++i) {
      if (!(edges[i] > edges[i - 1])) {
        return false;
      }
    }
 
    return true;
  }

  bool CentersFallWithinBins() const {
    if (n_bins == 0) {
      return centers.size() == 0;
    }
 
    if (centers.size() != n_bins || edges.size() != n_bins + 1) {
      return false;
    }
 
    for (std::size_t i = 0; i < n_bins; ++i) {
      if (!(centers[i] > edges[i] && centers[i] < edges[i + 1])) {
        return false;
      }
    }
 
    return true;
  }

  bool IsConsistent() const {
    if (n_bins == 0) {
      return edges.size() == 0 && centers.size() == 0;
    }
 
    return HasStrictlyIncreasingEdges() && CentersFallWithinBins();
  }
};
 
}  // namespace fims_popdy
#endif /* FIMS_POPULATION_DYNAMICS_SIZE_GRID_HPP */