bevholtDa removed, cannot solve

flr · Dec 2, 2023 · 72f09f1 · 72f09f1
1 parent 87a74d5
commit 72f09f1
Show file tree

Hide file tree

Showing 2 changed files with 1 addition and 12 deletions.
diff --git a/DESCRIPTION b/DESCRIPTION
@@ -1,6 +1,6 @@
 Package: FLBRP
 Title: Reference Points for Fisheries Management
-Version: 2.5.9.9020
+Version: 2.5.9.9021
 Authors@R: c(
   person("Iago", "Mosqueira", email = "[email protected]", role = "cre"),
   person("Laurence T.", "Kell", email = "[email protected]", role = "aut"),

diff --git a/src/FLBRP.cpp b/src/FLBRP.cpp
@@ -435,12 +435,6 @@ double FLBRP::Recruits(double FMult, int iUnit, int iIter)
          recruits = sr_params(1,1,iUnit,1,1,iIter)*ssb/(ssb+sr_params(2,1,iUnit,1,1,iIter));
       break;
 
-      // rec = a / (1 + (b / srp) ^ d)
-      case FLRConst_BevHoltDa:
-         ssb      = spr*sr_params(1,1,iUnit,1,1,iIter)-sr_params(2,1,iUnit,1,1,iIter);
-         recruits = sr_params(1,1,iUnit,1,1,iIter) / (1.0 + pow(sr_params(2,1,iUnit,1,1,iIter) / ssb, sr_params(3,1,iUnit,1,1,iIter)));
-      break;
-
       case FLRConst_Ricker:
          ssb      = log(spr*sr_params(1,1,iUnit,1,1,iIter))/sr_params(2,1,iUnit,1,1,iIter);
          recruits = sr_params(1,1,iUnit,1,1,iIter)*ssb*exp(-sr_params(2,1,iUnit,1,1,iIter)*ssb);
@@ -1880,11 +1874,6 @@ adouble FLBRP::Recruits(adouble FMult, int iUnit, int iIter)
          ssb      = spr*sr_params(1,1,iUnit,1,1,iIter)-sr_params(2,1,iUnit,1,1,iIter);
          recruits = sr_params(1,1,iUnit,1,1,iIter)*ssb/(ssb+sr_params(2,1,iUnit,1,1,iIter));
       break;
-      // rec = a / (1 + (b / srp) ^ d)
-      case FLRConst_BevHoltDa:
-         ssb      = spr*sr_params(1,1,iUnit,1,1,iIter)-sr_params(2,1,iUnit,1,1,iIter);
-         recruits = sr_params(1,1,iUnit,1,1,iIter) / (1.0 + pow(sr_params(2,1,iUnit,1,1,iIter) / ssb, sr_params(3,1,iUnit,1,1,iIter)));
-      break;
 
       case FLRConst_Ricker:
          ssb      = adtl::log(spr*sr_params(1,1,iUnit,1,1,iIter))/sr_params(2,1,iUnit,1,1,iIter);