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);
  }
};
 
}  // namespace fims_popdy
 
#endif /* POPULATION_DYNAMICS_SELECTIVITY_DOUBLE_LOGISTIC_HPP */