age_specific.hpp

Go to the documentation of this file.

 
#ifndef POPULATION_DYNAMICS_AGE_SPECIFIC_SELECTIVITY_HPP
#define POPULATION_DYNAMICS_AGE_SPECIFIC_SELECTIVITY_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 AgeSpecificSelectivity : public SelectivityBase<Type> {
  fims::Vector<Type> logit_sel_at_age;
 
  size_t n_ages;
 
  size_t min_age;
 
  AgeSpecificSelectivity() : SelectivityBase<Type>() {}
 
  virtual ~AgeSpecificSelectivity() {}
 
  virtual const Type evaluate(const Type &x) {
    Type a = static_cast<Type>(0.0);
    Type b = static_cast<Type>(1.0);
    double x_temp = fims_math::Value(x);
    size_t x_final = static_cast<size_t>(x_temp);
    size_t i_age = x_final - min_age;
    return fims_math::inv_logit<Type>(a, b, logit_sel_at_age[i_age]);
  }
 
  virtual const Type evaluate(const Type &x, size_t pos) {
    Type a = static_cast<Type>(0.0);
    Type b = static_cast<Type>(1.0);
    double x_temp = fims_math::Value(x);
    size_t x_final = static_cast<size_t>(x_temp);
    size_t i_age_year = pos * n_ages + x_final - min_age;
    return fims_math::inv_logit<Type>(
        a, b, logit_sel_at_age.get_force_scalar_wrap(i_age_year));
  }
};
 
}  // namespace fims_popdy
 
#endif /* POPULATION_DYNAMICS_AGE_SPECIFIC_SELECTIVITY_HPP */