void SavePlait(PlaitWS *pws, char *fdir){
FILE *fp;
int i;
char filename[BUFSIZE];
ws = pws;
/// print segment positions
for(i=0;i<ws->Opt.idx;i++){
sprintf(filename, "%ssegment.%d",fdir, i);
if(( fp = fopen (filename, "w") ) == NULL )
error("can not open:",filename);
PrintStEd(fp, ws->Opt.s[i]);
fclose(fp);
}//i
/// print regime labels
sprintf(filename, "%ssegment.labels",fdir);
if(( fp = fopen (filename, "w") ) == NULL )
error("can not open:",filename);
for(i=0;i<ws->Opt.idx;i++)
fprintf(fp, "%d\t\t%s\t\t%.0f\t\t%d \n", i, ws->Opt.s[i]->label, ws->Opt.s[i]->costT, ws->Opt.s[i]->model.k);
fclose(fp);
#if(PRINTHMM)
/// print HMM params
for(i=0;i<ws->Opt.idx;i++){
sprintf(filename, "%smodel.%d",fdir, i);
if(( fp = fopen (filename, "w") ) == NULL )
error("can not open:",filename);
PrintHMM(fp, &ws->Opt.s[i]->model);
fclose(fp);
}//i
#endif
#if(VITPATH)
/// print vit path
_computeViterbiPath(&ws->Opt);
#endif
}
void printModelSet(SegBox *s, int n){
int i;
for(i=0;i<n;i++){
PrintHMM(stdout,&s[i].model);
fprintf(stdout,"\n");
}
return;
}
void SaveModel(HMM *hmm, int g, char *fn){
FILE *fp;
int i;
char filename[BUFSIZE];
for(i=0;i<g;i++) {
sprintf(filename, "%s.model.0.%d",fn, i);
if(( fp = fopen (filename, "w") ) == NULL )
error("can not open:",fn);
PrintHMM(fp, &(hmm[i]));
fclose(fp);
}// i
}
