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;
}
