static void write_xvgr_graphs(const char *file, int ngraphs, int nsetspergraph,
const char *title, const char *subtitle,
const char *xlabel, const char **ylabel,
int n, real *x, real **y, real ***sy,
real scale_x, gmx_bool bZero, gmx_bool bSplit,
const output_env_t oenv)
{
FILE *out;
int g,s,i;
real min,max,xsp,ysp;
out=ffopen(file,"w");
if (output_env_get_xvg_format(oenv) == exvgXMGRACE) {
fprintf(out,"@ autoscale onread none\n");
}
for(g=0; g<ngraphs; g++) {
if (y) {
min=y[g][0];
max=y[g][0];
for(i=0; i<n; i++) {
if (y[g][i]<min) min=y[g][i];
if (y[g][i]>max) max=y[g][i];
}
} else {
min=sy[g][0][0];
max=sy[g][0][0];
for(s=0; s<nsetspergraph; s++)
for(i=0; i<n; i++) {
if (sy[g][s][i]<min) min=sy[g][s][i];
if (sy[g][s][i]>max) max=sy[g][s][i];
}
}
if (bZero)
min=0;
else
min=min-0.1*(max-min);
max=max+0.1*(max-min);
xsp=tick_spacing((x[n-1]-x[0])*scale_x,4);
ysp=tick_spacing(max-min,3);
if (output_env_get_print_xvgr_codes(oenv)) {
fprintf(out,"@ with g%d\n@ g%d on\n",g,g);
if (g==0) {
fprintf(out,"@ title \"%s\"\n",title);
if (subtitle)
fprintf(out,"@ subtitle \"%s\"\n",subtitle);
}
if (g==ngraphs-1)
fprintf(out,"@ xaxis label \"%s\"\n",xlabel);
else
fprintf(out,"@ xaxis ticklabel off\n");
if (n > 1) {
fprintf(out,"@ world xmin %g\n",x[0]*scale_x);
fprintf(out,"@ world xmax %g\n",x[n-1]*scale_x);
fprintf(out,"@ world ymin %g\n",min);
fprintf(out,"@ world ymax %g\n",max);
}
fprintf(out,"@ view xmin 0.15\n");
fprintf(out,"@ view xmax 0.85\n");
fprintf(out,"@ view ymin %g\n",0.15+(ngraphs-1-g)*0.7/ngraphs);
fprintf(out,"@ view ymax %g\n",0.15+(ngraphs-g)*0.7/ngraphs);
fprintf(out,"@ yaxis label \"%s\"\n",ylabel[g]);
fprintf(out,"@ xaxis tick major %g\n",xsp);
fprintf(out,"@ xaxis tick minor %g\n",xsp/2);
fprintf(out,"@ xaxis ticklabel start type spec\n");
fprintf(out,"@ xaxis ticklabel start %g\n",ceil(min/xsp)*xsp);
fprintf(out,"@ yaxis tick major %g\n",ysp);
fprintf(out,"@ yaxis tick minor %g\n",ysp/2);
fprintf(out,"@ yaxis ticklabel start type spec\n");
fprintf(out,"@ yaxis ticklabel start %g\n",ceil(min/ysp)*ysp);
if ((min<0) && (max>0)) {
fprintf(out,"@ zeroxaxis bar on\n");
fprintf(out,"@ zeroxaxis bar linestyle 3\n");
}
}
for(s=0; s<nsetspergraph; s++) {
for(i=0; i<n; i++) {
if ( bSplit && i>0 && abs(x[i])<1e-5 )
{
if (output_env_get_print_xvgr_codes(oenv))
fprintf(out,"&\n");
}
fprintf(out,"%10.4f %10.5f\n",
x[i]*scale_x,y ? y[g][i] : sy[g][s][i]);
}
if (output_env_get_print_xvgr_codes(oenv))
fprintf(out,"&\n");
}
}
ffclose(out);
}
void ps_mat(char *outf,int nmat,t_matrix mat[],t_matrix mat2[],
bool bFrame,bool bDiag,bool bFirstDiag,
bool bTitle,bool bTitleOnce,bool bYonce,int elegend,
real size,real boxx,real boxy,char *m2p,char *m2pout,
int mapoffset)
{
char *libm2p;
char buf[256],*legend;
t_psdata out;
t_psrec psrec,*psr;
int W,H;
int i,j,x,y,col,leg=0;
real x0,y0,xx;
real w,h,dw,dh;
int nmap1=0,nmap2=0,leg_nmap;
t_mapping *map1=NULL,*map2=NULL,*leg_map;
bool bMap1,bNextMap1,bDiscrete;
libm2p = m2p ? strdup(libfn(m2p)) : m2p;
get_params(libm2p,m2pout,&psrec);
psr=&psrec;
if (psr->X.major <= 0 )
tick_spacing((mat[0].flags & MAT_SPATIAL_X) ? mat[0].nx + 1 : mat[0].nx,
mat[0].axis_x, psr->X.offset, 'X',
&(psr->X.major), &(psr->X.minor) );
if (psr->X.minor <= 0 )
psr->X.minor = psr->X.major / 2;
if (psr->Y.major <= 0)
tick_spacing((mat[0].flags & MAT_SPATIAL_Y) ? mat[0].ny + 1 : mat[0].ny,
mat[0].axis_y, psr->Y.offset, 'Y',
&(psr->Y.major), &(psr->Y.minor) );
if (psr->Y.minor <= 0)
psr->Y.minor = psr->Y.major / 2;
if (boxx>0) {
psr->xboxsize=boxx;
psr->yboxsize=boxx;
}
if (boxy>0)
psr->yboxsize=boxy;
if (psr->xboxsize==0) {
psr->xboxsize = size/mat[0].nx;
printf("Set the x-size of the box to %.3f\n",psr->xboxsize);
}
if (psr->yboxsize==0) {
psr->yboxsize = size/mat[0].nx;
printf("Set the y-size of the box to %.3f\n",psr->yboxsize);
}
nmap1=0;
for (i=0; (i<nmat); i++)
if (mat[i].nmap>nmap1) {
nmap1=mat[i].nmap;
map1=mat[i].map;
leg=i+1;
}
if (leg!=1)
printf("Selected legend of matrix # %d for display\n",leg);
if (mat2) {
nmap2=0;
for (i=0; (i<nmat); i++)
if (mat2[i].nmap>nmap2) {
nmap2=mat2[i].nmap;
map2=mat2[i].map;
leg=i+1;
}
if (leg!=1)
printf("Selected legend of matrix # %d for second display\n",leg);
}
if ( (mat[0].legend[0]==0) && psr->legend )
strcpy(mat[0].legend, psr->leglabel);
bTitle = bTitle && mat[nmat-1].title[0];
bTitleOnce = bTitleOnce && mat[nmat-1].title[0];
psr->bTitle = bTitle;
psr->bTitleOnce = bTitleOnce;
psr->bYonce = bYonce;
/* Set up size of box for nice colors */
box_dim(nmat,mat,mat2,psr,elegend,bFrame,&w,&h,&dw,&dh);
/* Set up bounding box */
W=w+dw;
H=h+dh;
/* Start box at */
x0=dw;
y0=dh;
x = W+psr->xoffs;
y = H+psr->yoffs;
if (bFrame) {
x += 5*DDD;
y += 4*DDD;
}
out=ps_open(outf,0,0,x,y);
ps_linewidth(out,psr->linewidth);
ps_init_rgb_box(out,psr->xboxsize,psr->yboxsize);
ps_init_rgb_nbox(out,psr->xboxsize,psr->yboxsize);
ps_translate(out,psr->xoffs,psr->yoffs);
if (bFrame) {
ps_comment(out,"Here starts the BOX drawing");
draw_boxes(out,x0,y0,w,nmat,mat,psr);
}
for(i=0; (i<nmat); i++) {
if (bTitle || (bTitleOnce && i==nmat-1) ) {
/* Print title, if any */
ps_rgb(out,BLACK);
ps_strfont(out,psr->titfont,psr->titfontsize);
if (!mat2 || (strcmp(mat[i].title,mat2[i].title) == 0))
strcpy(buf,mat[i].title);
else
sprintf(buf,"%s / %s",mat[i].title,mat2[i].title);
ps_ctext(out,x0+w/2,y0+box_height(&(mat[i]),psr)+psr->titfontsize,
buf,eXCenter);
}
sprintf(buf,"Here starts the filling of box #%d",i);
ps_comment(out,buf);
for(x=0; (x<mat[i].nx); x++) {
int nexty;
int nextcol;
xx=x0+x*psr->xboxsize;
ps_moveto(out,xx,y0);
y=0;
bMap1 = (!mat2 || (x<y || (x==y && bFirstDiag)));
if ((bDiag) || (x!=y))
col = mat[i].matrix[x][y];
else
col = -1;
for(nexty=1; (nexty<=mat[i].ny); nexty++) {
bNextMap1 = (!mat2 || (x<nexty || (x==nexty && bFirstDiag)));
/* TRUE: upper left -> map1 */
/* FALSE: lower right -> map2 */
if ((nexty==mat[i].ny) || (!bDiag && (x==nexty)))
nextcol = -1;
else
nextcol=mat[i].matrix[x][nexty];
if ( (nexty==mat[i].ny) || (col!=nextcol) || (bMap1!=bNextMap1) ) {
if (col >= 0)
if (bMap1)
ps_rgb_nbox(out,&(mat[i].map[col].rgb),nexty-y);
else
ps_rgb_nbox(out,&(mat2[i].map[col].rgb),nexty-y);
else
ps_moverel(out,0,psr->yboxsize);
y=nexty;
bMap1=bNextMap1;
col=nextcol;
}
}
}
y0+=box_height(&(mat[i]),psr)+box_dh(psr)+box_dh_top(IS_ONCE,psr);
}
if (psr->X.lineatzero || psr->Y.lineatzero) {
/* reset y0 for first box */
y0=dh;
ps_comment(out,"Here starts the zero lines drawing");
draw_zerolines(out,x0,y0,w,nmat,mat,psr);
}
if (elegend!=elNone) {
ps_comment(out,"Now it's legend time!");
ps_linewidth(out,psr->linewidth);
if ( mat2==NULL || elegend!=elSecond ) {
bDiscrete = mat[0].bDiscrete;
legend = mat[0].legend;
leg_nmap = nmap1;
leg_map = map1;
} else {
bDiscrete = mat2[0].bDiscrete;
legend = mat2[0].legend;
leg_nmap = nmap2;
leg_map = map2;
}
if (bDiscrete)
leg_discrete(out,psr->legfontsize,DDD,legend,
psr->legfontsize,psr->legfont,leg_nmap,leg_map);
else {
if ( elegend!=elBoth )
leg_continuous(out,x0+w/2,w/2,DDD,legend,
psr->legfontsize,psr->legfont,leg_nmap,leg_map,
mapoffset);
else
leg_bicontinuous(out,x0+w/2,w,DDD,mat[0].legend,mat2[0].legend,
psr->legfontsize,psr->legfont,nmap1,map1,nmap2,map2);
}
ps_comment(out,"Were there, dude");
}
ps_close(out);
}
