示例#1
0
void model_parameters::preliminary_calculations(void)
{

#if defined(USE_ADPVM)

  admaster_slave_variable_interface(*this);

#endif
  // get the data out of the columns of the data matrix
  landings=column(data,1);
  biomass=column(data,2);
  
  // calculate observed surplus production (Pobs)
  for(int i=1;i<=nobs-1;i++)
  {
  Pobs(i)=biomass(i+1)-biomass(i)+landings(i);			
  }
  
  // set starting values
  lnMSY=log(0.6*max(Pobs));
  lnK=log(max(biomass));
  n=2;
  
  
}
示例#2
0
void model_parameters::preliminary_calculations(void)
{

#if defined(USE_ADPVM)

  admaster_slave_variable_interface(*this);

#endif
	incidence_functions();
	initialization();
	move_grow_die();
	clean_catage();
	output_true();
	output_dat();
	output_pin();
	exit(1);
}
示例#3
0
void model_parameters::preliminary_calculations(void)
{

  admaster_slave_variable_interface(*this);
 int X;
 double pi = 3.14159265358979323844;
 MatDelta = L95 - L50;
 Mat = 1.0/(1+mfexp(-log(19)*(LenMids-L50)/MatDelta));
 Fec = elem_prod(Mat, pow(LenMids, FecB));
 Fec = Fec/max(Fec);
 SDLinf = CVLinf * Linf;
 LinfdL = ((Linf + MaxSD * SDLinf) - (Linf - MaxSD * SDLinf))/(NGTG-1);
 for (X=0;X<NGTG;X++) {
   DiffLinfs(X+1) = (Linf - MaxSD * SDLinf) + X * LinfdL;  	 
 }
 
 RecProbs = 1/(sqrt(2*pi*SDLinf*SDLinf)) * mfexp(-(elem_prod((DiffLinfs-Linf),(DiffLinfs-Linf)))/(2*SDLinf*SDLinf));
 RecProbs = RecProbs/sum(RecProbs);
 
 cout << "damnit" << endl;
}
示例#4
0
void model_parameters::preliminary_calculations(void){
  admaster_slave_variable_interface(*this);
}