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 "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 */