Example #1
0
model_parameters::model_parameters(int sz,int argc,char * argv[]) : 
 model_data(argc,argv) , function_minimizer(sz)
{
  initializationfunction();
  lnK.allocate(-10000.,maxB,"lnK");
  lnMSY.allocate("lnMSY");
  n.allocate(0.1,30.,-2,"n");
  f.allocate("f");
  prior_function_value.allocate("prior_function_value");
  likelihood_function_value.allocate("likelihood_function_value");
  z.allocate("z");
  #ifndef NO_AD_INITIALIZE
  z.initialize();
  #endif
  zMSY.allocate("zMSY");
  #ifndef NO_AD_INITIALIZE
  zMSY.initialize();
  #endif
  meanP.allocate("meanP");
  #ifndef NO_AD_INITIALIZE
  meanP.initialize();
  #endif
  sig.allocate("sig");
  #ifndef NO_AD_INITIALIZE
  sig.initialize();
  #endif
  K.allocate("K");
  #ifndef NO_AD_INITIALIZE
  K.initialize();
  #endif
  MSY.allocate("MSY");
  #ifndef NO_AD_INITIALIZE
  MSY.initialize();
  #endif
  SSQ.allocate("SSQ");
  #ifndef NO_AD_INITIALIZE
  SSQ.initialize();
  #endif
  BMSY.allocate("BMSY");
  #ifndef NO_AD_INITIALIZE
  BMSY.initialize();
  #endif
  FMSY.allocate("FMSY");
  #ifndef NO_AD_INITIALIZE
  FMSY.initialize();
  #endif
  epsilon.allocate(1,nobs-1,"epsilon");
  #ifndef NO_AD_INITIALIZE
    epsilon.initialize();
  #endif
  Ppred.allocate(1,nobs-1,"Ppred");
  #ifndef NO_AD_INITIALIZE
    Ppred.initialize();
  #endif
  Pobs.allocate(1,nobs-1,"Pobs");
  #ifndef NO_AD_INITIALIZE
    Pobs.initialize();
  #endif
  sd_sig.allocate("sd_sig");
}
Example #2
0
model_parameters::model_parameters(int sz,int argc,char * argv[]) : 
 model_data(argc,argv) , function_minimizer(sz)
{
  initializationfunction();
  a1.allocate("a1");
  b1.allocate("b1");
  b2.allocate("b2");
  sdobs.allocate(1.0e-16,1.0e+4,"sdobs");
  nll.allocate("nll");
}
Example #3
0
model_parameters::model_parameters(int sz,int argc,char * argv[]) : 
 model_data(argc,argv) , function_minimizer(sz)
{
  initializationfunction();
  logSL50.allocate(logSL50Min,logSL50Max,1,"logSL50");
  logDelta.allocate(logDeltaMin,logDeltaMax,1,"logDelta");
  logFM.allocate(2,"logFM");
  SDLinf.allocate("SDLinf");
  #ifndef NO_AD_INITIALIZE
  SDLinf.initialize();
  #endif
  GTGLinf.allocate("GTGLinf");
  #ifndef NO_AD_INITIALIZE
  GTGLinf.initialize();
  #endif
  LinfdL.allocate("LinfdL");
  #ifndef NO_AD_INITIALIZE
  LinfdL.initialize();
  #endif
  Delta.allocate("Delta");
  #ifndef NO_AD_INITIALIZE
  Delta.initialize();
  #endif
  EP0.allocate("EP0");
  #ifndef NO_AD_INITIALIZE
  EP0.initialize();
  #endif
  EPf.allocate("EPf");
  #ifndef NO_AD_INITIALIZE
  EPf.initialize();
  #endif
  MatDelta.allocate("MatDelta");
  #ifndef NO_AD_INITIALIZE
  MatDelta.initialize();
  #endif
  DiffLinfs.allocate(1,NGTG,"DiffLinfs");
  #ifndef NO_AD_INITIALIZE
    DiffLinfs.initialize();
  #endif
  RecProbs.allocate(1,NGTG,"RecProbs");
  #ifndef NO_AD_INITIALIZE
    RecProbs.initialize();
  #endif
  PredLenComp.allocate(1,NLenMids,"PredLenComp");
  #ifndef NO_AD_INITIALIZE
    PredLenComp.initialize();
  #endif
  PredUnfishedComp.allocate(1,NLenMids,"PredUnfishedComp");
  #ifndef NO_AD_INITIALIZE
    PredUnfishedComp.initialize();
  #endif
  Vul.allocate(1,NLenMids+1,"Vul");
  #ifndef NO_AD_INITIALIZE
    Vul.initialize();
  #endif
  MkL.allocate(1,NLenMids+1,"MkL");
  #ifndef NO_AD_INITIALIZE
    MkL.initialize();
  #endif
  FkL.allocate(1,NLenMids+1,"FkL");
  #ifndef NO_AD_INITIALIZE
    FkL.initialize();
  #endif
  Mat.allocate(1,NLenMids,"Mat");
  #ifndef NO_AD_INITIALIZE
    Mat.initialize();
  #endif
  Wght.allocate(1,NLenMids,"Wght");
  #ifndef NO_AD_INITIALIZE
    Wght.initialize();
  #endif
  Fec.allocate(1,NLenMids,"Fec");
  #ifndef NO_AD_INITIALIZE
    Fec.initialize();
  #endif
  EP0_gtg.allocate(1,NGTG,"EP0_gtg");
  #ifndef NO_AD_INITIALIZE
    EP0_gtg.initialize();
  #endif
  EPf_gtg.allocate(1,NGTG,"EPf_gtg");
  #ifndef NO_AD_INITIALIZE
    EPf_gtg.initialize();
  #endif
  PUnFished.allocate(1,NLenMids+1,"PUnFished");
  #ifndef NO_AD_INITIALIZE
    PUnFished.initialize();
  #endif
  PFished.allocate(1,NLenMids+1,"PFished");
  #ifndef NO_AD_INITIALIZE
    PFished.initialize();
  #endif
  NUnFished.allocate(1,NLenMids,"NUnFished");
  #ifndef NO_AD_INITIALIZE
    NUnFished.initialize();
  #endif
  NFished.allocate(1,NLenMids,"NFished");
  #ifndef NO_AD_INITIALIZE
    NFished.initialize();
  #endif
  currMkL.allocate(1,NLenMids+1,"currMkL");
  #ifndef NO_AD_INITIALIZE
    currMkL.initialize();
  #endif
  currZkL.allocate(1,NLenMids+1,"currZkL");
  #ifndef NO_AD_INITIALIZE
    currZkL.initialize();
  #endif
  UnfishedMatrix.allocate(1,NGTG,1,NLenMids,"UnfishedMatrix");
  #ifndef NO_AD_INITIALIZE
    UnfishedMatrix.initialize();
  #endif
  FishedMatrix.allocate(1,NGTG,1,NLenMids,"FishedMatrix");
  #ifndef NO_AD_INITIALIZE
    FishedMatrix.initialize();
  #endif
  MKLMat.allocate(1,NGTG,1,NLenMids+1,"MKLMat");
  #ifndef NO_AD_INITIALIZE
    MKLMat.initialize();
  #endif
  ZKLMat.allocate(1,NGTG,1,NLenMids+1,"ZKLMat");
  #ifndef NO_AD_INITIALIZE
    ZKLMat.initialize();
  #endif
  FMpar.allocate("FMpar");
  SPR.allocate("SPR");
  SL50.allocate("SL50");
  SL95.allocate("SL95");
  temp.allocate(1,NLenMids+1,"temp");
  #ifndef NO_AD_INITIALIZE
    temp.initialize();
  #endif
  tempNum.allocate("tempNum");
  #ifndef NO_AD_INITIALIZE
  tempNum.initialize();
  #endif
  obj_fun.allocate("obj_fun");
  prior_function_value.allocate("prior_function_value");
  likelihood_function_value.allocate("likelihood_function_value");
}
Example #4
0
model_parameters::model_parameters(int sz,int argc,char * argv[]) : 
 model_data(argc,argv) , function_minimizer(sz)
{
  initializationfunction();
  no_f.allocate("no_f");
  prior_function_value.allocate("prior_function_value");
  likelihood_function_value.allocate("likelihood_function_value");
  kappa.allocate("kappa");
  #ifndef NO_AD_INITIALIZE
  kappa.initialize();
  #endif
  phie.allocate("phie");
  #ifndef NO_AD_INITIALIZE
  phie.initialize();
  #endif
  So.allocate("So");
  #ifndef NO_AD_INITIALIZE
  So.initialize();
  #endif
  Bo.allocate("Bo");
  #ifndef NO_AD_INITIALIZE
  Bo.initialize();
  #endif
  beta.allocate("beta");
  #ifndef NO_AD_INITIALIZE
  beta.initialize();
  #endif
  tBo.allocate("tBo");
  #ifndef NO_AD_INITIALIZE
  tBo.initialize();
  #endif
  lxo.allocate(sage,nage,"lxo");
  #ifndef NO_AD_INITIALIZE
    lxo.initialize();
  #endif
  za.allocate(sage,nage,"za");
  #ifndef NO_AD_INITIALIZE
    za.initialize();
  #endif
  SB.allocate(1,ntstp,"SB");
  #ifndef NO_AD_INITIALIZE
    SB.initialize();
  #endif
  varPos.allocate(sage,nage,"varPos");
  #ifndef NO_AD_INITIALIZE
    varPos.initialize();
  #endif
  maxPos.allocate(sage,nage,"maxPos");
  #ifndef NO_AD_INITIALIZE
    maxPos.initialize();
  #endif
  NationVulB.allocate(1,ntstp,1,nations,"NationVulB");
  #ifndef NO_AD_INITIALIZE
    NationVulB.initialize();
  #endif
  Nage.allocate(1,ntstp,sage,nage,"Nage");
  #ifndef NO_AD_INITIALIZE
    Nage.initialize();
  #endif
  VulB.allocate(1,ntstp,sage,nage,"VulB");
  #ifndef NO_AD_INITIALIZE
    VulB.initialize();
  #endif
  PosX.allocate(1,ntstp,sage,nage,"PosX");
  #ifndef NO_AD_INITIALIZE
    PosX.initialize();
  #endif
  Effage.allocate(1,ntstp,sage,nage,"Effage");
  #ifndef NO_AD_INITIALIZE
    Effage.initialize();
  #endif
  VBarea.allocate(1,ntstp,sarea,narea,"VBarea");
  #ifndef NO_AD_INITIALIZE
    VBarea.initialize();
  #endif
  Effarea.allocate(1,ntstp,sarea,narea,"Effarea");
  #ifndef NO_AD_INITIALIZE
    Effarea.initialize();
  #endif
  NAreaAge.allocate(1,ntstp,sarea,narea,sage,nage,"NAreaAge");
  #ifndef NO_AD_INITIALIZE
    NAreaAge.initialize();
  #endif
  CatchAreaAge.allocate(1,ntstp,sarea,narea,sage,nage,"CatchAreaAge");
  #ifndef NO_AD_INITIALIZE
    CatchAreaAge.initialize();
  #endif
  CatchNatAge.allocate(1,ntstp,1,nations,sage,nage,"CatchNatAge");
  #ifndef NO_AD_INITIALIZE
    CatchNatAge.initialize();
  #endif
  EffNatAge.allocate(1,nations,1,ntstp,sage-2,nage,"EffNatAge");
  #ifndef NO_AD_INITIALIZE
    EffNatAge.initialize();
  #endif
  obsCatchNatAge.allocate(1,tot_pcat,sage-3,nage,"obsCatchNatAge");
  #ifndef NO_AD_INITIALIZE
    obsCatchNatAge.initialize();
  #endif
}