mortality_factors.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: 2020 W. Traylor <wolfgang.traylor@senckenberg.de>
//
// SPDX-License-Identifier: LGPL-3.0-or-later
 
#ifndef FAUNA_MORTALITY_FACTORS_H
#define FAUNA_MORTALITY_FACTORS_H
 
#include <stdexcept>
 
namespace Fauna {
 
class GetBackgroundMortality {
 public:
 
  GetBackgroundMortality(const double annual_mortality_1st_year,
                         const double annual_mortality);
 
 
  double operator()(const int age_days) const;
 
 private:
  double annual_mortality;
  double annual_mortality_1st_year;
};
 
class GetSimpleLifespanMortality {
 public:
 
  GetSimpleLifespanMortality(const int lifespan_years);
 
 
  double operator()(const int age_days) const {
    if (age_days < 0)
      throw std::invalid_argument(
          "Fauna::GetSimpleLifespanMortality::operator()() "
          "age_days < 0");
    if (age_days >= lifespan_years * 365)
      return 1.0;  // all dead
    else
      return 0.0;  // all survive
  }
 
 private:
  int lifespan_years;
};
 
 
class GetStarvationIlliusOConnor2000 {
 public:
 
  GetStarvationIlliusOConnor2000(const double fat_standard_deviation = .0125,
                                 const bool shift_body_condition = true);
 
 
  double operator()(const double body_condition,
                    double& new_body_condition) const;
 
 private:
  const double fat_standard_deviation;
  bool shift_body_condition;
 
 
  static double cumulative_normal_distribution(double x);
};
 
 
class GetStarvationMortalityThreshold {
 public:
  static const double DEFAULT_MIN_BODYFAT;
 
 
  GetStarvationMortalityThreshold(
      const double min_bodyfat = DEFAULT_MIN_BODYFAT);
 
 
  double operator()(const double bodyfat) const;
 
 private:
  double min_bodyfat;
};
}  // namespace Fauna
 
#endif  // FAUNA_MORTALITY_FACTORS_H