int main (int *arga,int *argb,int *argc,int *argd,double *OUT, char **HEADER,int *QFAI,double *QFADy,double *QFADx,int *QFADNoORF,int *QFADNoTIME,double *PRIORS)/****/ { GetRNGstate(); struct_data *data= malloc(sizeof(struct_data)); struct_para *para= malloc(sizeof(struct_para)); struct_priors *priors= malloc(sizeof(struct_priors)); struct_MH *MH = malloc(sizeof(struct_MH)); int burn,iters,thin, CAPL; /*testargc(argc);*//********/ burn=*arga; /**********/ iters=*argb; /**************/ thin=*argc; /**************/ CAPL=*argd; /**************/ inzstruct_data(data,QFAI,QFADy,QFADx,QFADNoORF,QFADNoTIME); inzstruct_priors(priors,PRIORS); inzstruct_para(para,data,priors); inzstruct_MH(MH); gibbsandMHloop(burn,1,data,para,priors,MH,CAPL,0,OUT,HEADER); gibbsandMHloop(iters,thin,data,para,priors,MH,CAPL,1,OUT,HEADER); PutRNGstate(); return 0; }
int main (int argc,char *argv[]) { struct_data *data= malloc(sizeof(struct_data)); struct_para *para= malloc(sizeof(struct_para)); struct_priors *priors= malloc(sizeof(struct_priors)); struct_MH *MH = malloc(sizeof(struct_MH)); int burn,iters,thin; long seed; const gsl_rng_type * T; gsl_rng * RNG; testargc(argc); gsl_rng_env_setup (); T = gsl_rng_default; RNG = gsl_rng_alloc (T); seed = time (NULL) * getpid(); gsl_rng_set (RNG, seed); /*seed*/ burn=atoi(argv[1]); /*Burn in*/ iters=atoi(argv[2]); /*iterations*/ thin=atoi(argv[3]); /*thining*/ inzstruct_data(data); inzstruct_para(para,data); inzstruct_priors(priors); inzstruct_MH(MH); gibbsandMHloop(burn,1,RNG,data,para,priors,MH,0); gibbsandMHloop(iters,thin,RNG,data,para,priors,MH,1); gsl_rng_free(RNG); return 0; }