int inzstruct_para(struct_para *para,struct_data *data,struct_priors *priors)
{
long size;
size=data->SHIFTlmn;
para->K_lm=malloc(size*sizeof(double));
para->r_lm=malloc(size*sizeof(double));
size=data->L;
para->tau_K_l=malloc(size*sizeof(double));
para->tau_r_l=malloc(size*sizeof(double));
para->K_o_l=malloc(size*sizeof(double));
para->r_o_l=malloc(size*sizeof(double));
para->nu_l=malloc(size*sizeof(double));
fillpara(para,data,priors);
return 0;
}
int inzstruct_para(struct_para *para,struct_data *data)
{
long size;
size=data->L*2;
para->tau_K_cl=malloc(size*sizeof(double));
para->tau_r_cl=malloc(size*sizeof(double));
size=17308*2;/*inputfromfile*/
para->K_clm=malloc(size*sizeof(double));
para->r_clm=malloc(size*sizeof(double));
size=data->L;
para->delta=malloc(size*sizeof(double));
para->gamma=malloc(size*sizeof(double));
para->omega=malloc(size*sizeof(double));
para->K_o_l=malloc(size*sizeof(double));
para->r_o_l=malloc(size*sizeof(double));
para->nu_l=malloc(size*sizeof(double));
size=2;
para->alpha=malloc(size*sizeof(double));
para->beta=malloc(size*sizeof(double));
para->upsilon_c=malloc(size*sizeof(double));
fillpara(para,data);
return 0;
}
/* Justify the entire buffer instead of just a paragraph */
int
fillbuf(int f, int n)
{
LINE *eobline;
REGION region;
if(curbp->b_mode&MDVIEW){ /* don't allow this command if */
return(rdonly()); /* we are in read only mode */
}
else if (fillcol == 0) { /* no fill column set */
mlwrite_utf8("No fill column set", NULL);
return(FALSE);
}
if((lastflag & CFFILL) && (lastflag & CFFLBF)){
/* no use doing a full justify twice */
thisflag |= (CFFLBF | CFFILL);
return(TRUE);
}
/* record the pointer of the last line */
if(gotoeob(FALSE, 1) == FALSE)
return(FALSE);
eobline = curwp->w_dotp; /* last line of buffer */
if(!llength(eobline))
eobline = lback(eobline);
/* and back to the beginning of the buffer */
gotobob(FALSE, 1);
thisflag |= CFFLBF; /* CFFILL also gets set in fillpara */
if(!Pmaster)
sgarbk = TRUE;
curwp->w_flag |= WFMODE;
/*
* clear the kill buffer, that's where we'll store undo
* information, we can't do the fill buffer because
* fillpara relies on its contents
*/
kdelete();
curwp->w_doto = 0;
getregion(®ion, eobline, llength(eobline));
/* Put full message in the kill buffer for undo */
if(!ldelete(region.r_size, kinsert))
return(FALSE);
/* before yank'ing, clear lastflag so we don't just unjustify */
lastflag &= ~(CFFLBF | CFFILL);
/* Now in kill buffer, bring back text to use in fillpara */
yank(FALSE, 1);
gotobob(FALSE, 1);
/* call fillpara until we're at the end of the buffer */
while(curwp->w_dotp != curbp->b_linep)
if(!(fillpara(FALSE, 1)))
return(FALSE);
return(TRUE);
}