forked from NOAA-FIMS/MAS
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Fecundity.hpp
112 lines (90 loc) · 3.59 KB
/
Fecundity.hpp
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
/*
* File: Fecundity.hpp
*
* Author: Matthew Supernaw
* National Oceanic and Atmospheric Administration
* National Marine Fisheries Service
* Sustainable Fisheries Division
* St. Petersburg, FL, 33701
*
* Created on September 16, 2016, 12:34 PM
*
* This File is part of the NOAA, National Marine Fisheries Service
* Metapopulation Assessment System project.
*
* This software is a "United States Government Work" under the terms of the
* United States Copyright Act. It was written as part of the author's official
* duties as a United States Government employee and thus cannot be copyrighted.
* This software is freely available to the public for use. The National Oceanic
* And Atmospheric Administration and the U.S. Government have not placed any
* restriction on its use or reproduction. Although all reasonable efforts have
* been taken to ensure the accuracy and reliability of the software and data,
* the National Oceanic And Atmospheric Administration and the U.S. Government
* do not and cannot warrant the performance or results that may be obtained by
* using this software or data. The National Oceanic And Atmospheric
* Administration and the U.S. Government disclaim all warranties, express or
* implied, including warranties of performance, merchantability or fitness
* for any particular purpose.
*
* Please cite the author(s) in any work or product based on this material.
*
*/
#ifndef MAS_FECUNDITY_HPP
#define MAS_FECUNDITY_HPP
#include "Common.hpp"
namespace mas {
template<typename REAL_T>
struct FecundityBase : mas::ModelObject<REAL_T> {
typedef typename VariableTrait<REAL_T>::variable variable;
virtual const variable Evaluate(const int& sex, const variable& age) = 0;
virtual const std::string Name() {
return "FecundityBase";
}
};
template<typename REAL_T>
struct LogisticFec : FecundityBase<REAL_T> {
typedef typename VariableTrait<REAL_T>::variable variable;
variable a50; //age of 50% maturity
variable s; //the rate of increase in maturity at a50
/**
* Age based logistic maturity
*
* @param sex
* @param age
* @return fraction_mature
*/
virtual const variable Evaluate(const int& sex, const variable& age) {
return static_cast<REAL_T> (1.0) / (static_cast<REAL_T> (1.0) + std::exp(-1.0*(age - a50)/s));
}
virtual const std::string Name() {
return "Logistic Fecundity";
}
};
template<typename REAL_T>
struct DoubleLogisticFec : FecundityBase<REAL_T> {
typedef typename VariableTrait<REAL_T>::variable variable;
variable alpha_asc; //ascending alpha
variable beta_asc; // ascending beta
variable alpha_desc; // descending alpha
variable beta_desc; // descending beta
/**
* Age based double logistic maturity
*
* @param sex
* @param age
* @return fraction_mature
*/
virtual const variable Evaluate(const int& sex, const variable& age) {
return (static_cast<REAL_T> (1.0) /
(static_cast<REAL_T> (1.0) +
std::exp(-beta_asc * (age - alpha_asc)))) *
(static_cast<REAL_T> (1.0) - (static_cast<REAL_T> (1.0) /
(static_cast<REAL_T> (1.0) +
std::exp(-beta_desc * (age - alpha_desc)))));
}
virtual const std::string Name() {
return "Double Logistic Fecundity";
}
};
}
#endif /* MAS_FECUNDITY_HPP */