double_logistic.hpp

Go to the documentation of this file.

 
#ifndef POPULATION_DYNAMICS_SELECTIVITY_DOUBLE_LOGISTIC_HPP
#define POPULATION_DYNAMICS_SELECTIVITY_DOUBLE_LOGISTIC_HPP
 
// #include "../../../interface/interface.hpp"
#include "../../../common/fims_math.hpp"
#include "../../../common/fims_vector.hpp"
#include "selectivity_base.hpp"
 
namespace fims_popdy {
 
template <typename Type>
struct DoubleLogisticSelectivity : public SelectivityBase<Type> {
  fims::Vector<Type> inflection_point_asc; 
  fims::Vector<Type> slope_asc; 
  fims::Vector<Type> inflection_point_desc; 
  fims::Vector<Type> slope_desc; 
  DoubleLogisticSelectivity() : SelectivityBase<Type>() {}
 
  virtual ~DoubleLogisticSelectivity() {}
 
  virtual const Type evaluate(const Type& x) {
    return fims_math::double_logistic<Type>(
        inflection_point_asc[0], slope_asc[0], inflection_point_desc[0],
        slope_desc[0], x);
  }
 
  virtual const Type evaluate(const Type& x, size_t pos) {
    return fims_math::double_logistic<Type>(
        inflection_point_asc.get_force_scalar(pos),
        slope_asc.get_force_scalar(pos),
        inflection_point_desc.get_force_scalar(pos),
        slope_desc.get_force_scalar(pos), x);
  }
 
  virtual void create_report_vectors(
      std::map<std::string, fims::Vector<fims::Vector<Type>>>& report_vectors) {
    report_vectors["inflection_point_asc"].emplace_back(
        inflection_point_asc.to_tmb());
    report_vectors["slope_asc"].emplace_back(slope_asc.to_tmb());
    report_vectors["inflection_point_desc"].emplace_back(
        inflection_point_desc.to_tmb());
    report_vectors["slope_desc"].emplace_back(slope_desc.to_tmb());
  }
 
  virtual void get_report_vector_count(
      std::map<std::string, size_t>& report_vector_count) {
    report_vector_count["inflection_point_asc"] += 1;
    report_vector_count["slope_asc"] += 1;
    report_vector_count["inflection_point_desc"] += 1;
    report_vector_count["slope_desc"] += 1;
  }
};
 
}  // namespace fims_popdy
 
#endif /* POPULATION_DYNAMICS_SELECTIVITY_DOUBLE_LOGISTIC_HPP */