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; }
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); }
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; }
void model_parameters::preliminary_calculations(void){ admaster_slave_variable_interface(*this); }