logistic.hpp

Go to the documentation of this file.

 
#ifndef POPULATION_DYNAMICS_SELECTIVITY_LOGISTIC_HPP
#define POPULATION_DYNAMICS_SELECTIVITY_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 LogisticSelectivity : public SelectivityBase<Type> {
  fims::Vector<Type>
      inflection_point;     
  fims::Vector<Type> slope; 
  LogisticSelectivity() : SelectivityBase<Type>() {}
 
  virtual ~LogisticSelectivity() {}
 
  virtual const Type evaluate(const Type &x) {
    return fims_math::logistic<Type>(inflection_point[0], slope[0], x);
  }
 
  virtual const Type evaluate(const Type &x, size_t pos) {
    return fims_math::logistic<Type>(inflection_point.get_force_scalar(pos),
                                     slope.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"].emplace_back(inflection_point);
    report_vectors["slope"].emplace_back(slope);
  }
 
  virtual void get_report_vector_count(
      std::map<std::string, size_t> &report_vector_count) {
    report_vector_count["inflection_point"] += 1;
    report_vector_count["slope"] += 1;
  }
};
 
}  // namespace fims_popdy
 
#endif /* POPULATION_DYNAMICS_SELECTIVITY_LOGISTIC_HPP */