comparison tools/x11play.c @ 0:5242703e91d3 tip

Initial checkin for AIM92 aimR8.2 (last updated May 1997).

+/****************************************************************************
+*
+*   x11play    -Animate successive plots of n points.
+		The points are read as binary shorts.
+*               Like x11plot, but plots successive blocks of n points
+*               on top of each other, rubbing out the previous plot, so
+*               as to form an "animated" sequence of plots.
+*
+*              -Help:  x11play -H
+*              -Usage, EG:
+*                      fft -n128 -o120 file | x11play -n64 -y0 -c5000
+*                      Quit program with left mouse button.
+*             -Step-mode
+*                      With a -S flag, single-step through frames using
+*                      the mouse as follows:
+*
+*                      left -   quit
+*                      centre - freeze animation (successive centre
+*                               buttons give single successive frames).
+*                      right  - continue animation
+*
+*                      Note: animation is slower in step-mode because
+*                      XSynchronise has to be on.
+****************************************************************************/
+#include <stdio.h>
+#include <math.h>
+#include "x11coord.h"
+
+/****************************************************************************
+*   Defaults
+****************************************************************************/
+#define UNSET        (-9999)
+
+#define DEFAMPL        1000
+#define DEFSCALE        256
+
+float   argXd0=UNSET,   argYd0=UNSET;
+float   argXd1=UNSET,   argYd1=UNSET;
+
+/****************************************************************************
+*   Input-file arguments
+****************************************************************************/
+int     n =      DEFSCALE;      /* num points to read from file */
+int     s =      0;             /* start position in file */
+short   summary= 0;             /* Flag for method of block summary in downsample */
+short   z =      1;             /* Flag for drawing zero line */
+short   stepmode=0;             /* Flag for single-step mode */
+short   overplot=0;             /* Flag for overplot mode */
+short   stepcount=0;
+char    DEFNAME[]="stdin";
+char    *name=DEFNAME;
+short   *data;
+XPoint  *points;
+
+FILE    *fp, *fopen();
+
+/****************************************************************************
+*       main
+*   The window is plotted when an expose event is (automatically) generated.
+*   The plotting routines are called from within the event monitor.
+****************************************************************************/
+main(argc, argv)
+int   argc ;
+char *argv[] ;
+{
+    Window  w;
+    int     npoints;
+
+    while (--argc > 0 && **++argv == '-')
+	switch (*++*argv) {
+	    case 'n':   n       = atoi(++*argv); break;
+	    case 's':   s       = atoi(++*argv); break;
+	    case 'x':   argXd0  = atof(++*argv); break;
+	    case 'y':   argYd0  = atof(++*argv); break;
+	    case 'c':   argXd1  = atof(++*argv); break;
+	    case 'a':   argYd1  = atof(++*argv); break;
+	    case 'b':   summary = atoi(++*argv); break;
+	    case 'z':   z       = 0;             break;
+	    case 'S':   stepmode++;              break;
+	    case 'O':   overplot++;              break;
+	    case 'F':   theFontString = ++*argv; break;
+	    case 'H':
+	    default:    help();
+	}
+
+    if (argXd0==UNSET)  argXd0 = s;
+    if (argXd1==UNSET)  argXd1 = n;
+    else                argXd1 = argXd1-argXd0; /* ie, a range */
+    if (argYd0==UNSET) {
+	if (argYd1==UNSET)  argYd1 = 2*DEFAMPL;
+	else                argYd1 = 2*argYd1;
+	argYd0 = argYd1/2;
+    }
+    else {
+	if (argYd1==UNSET)  argYd1 = DEFAMPL+argYd0;
+	else                argYd1 = argYd1+argYd0;
+    }
+    argYd0 = (-argYd0);
+
+    if ((data = (short*) malloc(n * sizeof(short))) == NULL)
+	error("malloc out of space\n");
+    if ((points = (XPoint*)malloc(n * sizeof(XPoint))) == NULL)
+	error("malloc out of space\n");
+
+    if (argc) {
+	/* Open data file */
+	name = *argv;
+	if ((fp = fopen(name, "r")) == NULL) {
+	    fprintf(stderr,"can't open %s\n", *argv);
+	    exit(1);
+	}
+    }
+    else fp = stdin;
+
+    /* Seek start of data in file */
+    if (s > 0) seekstart(s,fp);
+
+    set_window_parameters(PLOTXORG,PLOTYORG, PLOTWIDTH,PLOTHEIGHT);
+    set_bordered_box(10);
+    set_data_parameters(0,0, argXd0,argYd0, argXd1,argYd1);
+    xopen();
+    xgc( 1, theFontString ) ;
+    w = xcreate(name, ButtonPressMask | ExposureMask);
+    init_xevent(w);
+}
+
+
+/****************************************************************************
+*   Seek s points from current position in stream.
+*   This general seek works even when fp is stdin, unlike fseek.
+****************************************************************************/
+seekstart(s,fp)
+int   s;
+FILE  *fp;
+{
+    int   i;
+    short p;
+
+    for (i=0 ; i<s && fread(&p, sizeof(short), 1, fp) ; i++)
+	;
+    if (i<s) error("seek overshot end-of-file\n");
+}
+
+
+/****************************************************************************
+*       X11 Event Monitor.
+*   The window is plotted when an expose event is (automatically) generated.
+*   This routine catches the initial expose event to first draw the window.
+*   After this, the xevent_monitor is called from upsample, which is called
+*   from draw_plot.
+****************************************************************************/
+init_xevent(w)  /* to catch initial expose event and draw first plot */
+Window w;
+{
+    XEvent event;
+
+    for ( ; ; ) {
+	XNextEvent(theDisplay,&event);
+	switch (event.type) {
+	    case ButtonPress:
+		switch(event.xbutton.button) {
+		    case Button1:           /* Left   */
+			XDestroyWindow(theDisplay,w);
+			exit(0);
+		    case Button2:           /* Middle */
+		    case Button3: break;    /* Right */
+		    default:
+			fprintf(stderr,"button %d not used\n", event.xbutton.button);
+		}
+		break;
+	    case Expose:
+		if (event.xexpose.count == 0) {
+		    draw_calibrated_box(w);
+		    draw_plot(w);
+		}
+		break;
+	    default:
+		fprintf(stderr,"event type %d not found\n", event.type);
+	}
+    }
+}
+
+
+xevent_monitor(w)
+Window w;
+{
+    XEvent event;
+    static short  wait=1;
+
+    if (!stepmode) wait=0;
+    do {
+	/* get next event provided one is on queue; don't wait otherwise */
+	if (wait || XEventsQueued(theDisplay,QueuedAlready) > 0) {
+	    XNextEvent(theDisplay,&event);
+	    switch (event.type) {
+		case ButtonPress:
+		    switch(event.xbutton.button) {
+			case Button1:               /* Left   */
+			    XDestroyWindow(theDisplay,w);
+			    exit(0);
+			case Button2: /* Middle */
+			    if (stepmode) {
+			       if (wait) single_shot(w,fp);
+			       else wait=1;
+			    }
+			    break;
+			case Button3: if (stepmode) wait=0; break;/* Right */
+			default:
+			    fprintf(stderr,"button %d not used\n", event.xbutton.button);
+		    }
+		    break;
+		default:
+		    fprintf(stderr,"event type %d not found\n", event.type);
+	    }
+	}
+    } while (wait);
+}
+
+
+/****************************************************************************
+*   Draw a calibrated box in a window and plot data inside it.
+****************************************************************************/
+draw_plot(w)
+Window  w;
+{
+    int i, npoints;
+
+    if (z) draw_Xaxis(w,0);     /* Draw zero line */
+    if (!overplot) XSetFunction(theDisplay,theGC,GXinvert);
+    if (stepmode) XSynchronize(theDisplay,True);
+    if (n <= pwidth()) npoints = upsample(w,fp);
+    else  fprintf(stderr,"successive downsampling not implemented yet!\n");
+    /*    npoints = downsample(w,fp);   */
+}
+
+
+/****************************************************************************
+*       Up and Down sampling.
+****************************************************************************/
+char stepcountstr[128];
+
+upsample (w,fp)
+Window  w;
+FILE   *fp;
+{
+    int    i, x, j=0;
+    int    lastx,lasty, X,Y;
+    short  y, notfirst=0;
+    float  res, incr, carry;
+
+    res = (float)pwidth()/(n-1);                /* resolution of plot */
+    incr = res;                                 /* x-value increment */
+    x = 0;                                      /* initial x-value */
+    while (fread(&y, sizeof(short), 1, fp)) {
+	if (stepmode) {
+	    if (stepcount > 0) supertopline(w,stepcountstr);
+	    sprintf(stepcountstr,"%d", ++stepcount);
+	    supertopline(w,stepcountstr);
+	}
+	lastx=pleft(); lasty=py(y);
+	carry = res - (int)res;
+	incr = res + carry;
+	x = (int)incr;
+	for (i=1 ; i<n  && fread(&y, sizeof(short), 1, fp) ; i++) {
+	    X = x+pleft();
+	    Y = py(y);
+	    /* Plot only inside y-limits of box */
+	    if (Y > pbottom()) Y = pbottom();
+	    if (Y < ptop())    Y = ptop();
+
+	    /* rub-out a segment of last plot */
+	    if (notfirst)
+		XDrawLine(theDisplay, w, theGC, points[i-1].x, points[i-1].y, points[i].x, points[i].y);
+	    else  notfirst=1;
+	    /* draw a new segment for current plot */
+	    XDrawLine(theDisplay, w, theGC, lastx, lasty, X, Y);
+
+	    points[i-1].x = lastx;
+	    points[i-1].y = lasty;
+	    lastx = X;
+	    lasty = Y;
+
+	    carry = incr - (int)incr;
+	    incr = res + carry;
+	    x += (int)incr;
+	}
+	points[i-1].x = lastx;
+	points[i-1].y = lasty;
+/*        if (XEventsQueued(theDisplay,QueuedAlready) > 0)      */
+	    xevent_monitor(w);
+    }
+}
+
+
+single_shot (w,fp)
+Window  w;
+FILE   *fp;
+{
+    int    i, x, j=0;
+    int    lastx,lasty, X,Y;
+    short  y, notfirst=0;
+    float  res, incr, carry;
+
+    res = (float)pwidth()/(n-1);                /* resolution of plot */
+    incr = res;                                 /* x-value increment */
+    x = 0;                                      /* initial x-value */
+    if (fread(&y, sizeof(short), 1, fp)) {
+	if (stepmode) {
+	    if (stepcount > 0) supertopline(w,stepcountstr);
+	    sprintf(stepcountstr,"%d", ++stepcount);
+	    supertopline(w,stepcountstr);
+	}
+	lastx=pleft(); lasty=py(y);
+	carry = res - (int)res;
+	incr = res + carry;
+	x = (int)incr;
+	for (i=1 ; i<n  && fread(&y, sizeof(short), 1, fp) ; i++) {
+	    X = x+pleft();
+	    Y = py(y);
+	    /* Plot only inside y-limits of box */
+	    if (Y > pbottom()) Y = pbottom();
+	    if (Y < ptop())    Y = ptop();
+
+	    /* rub-out a segment of last plot */
+	    if (notfirst)
+		XDrawLine(theDisplay, w, theGC, points[i-1].x, points[i-1].y, points[i].x, points[i].y);
+	    else  notfirst=1;
+	    /* draw a new segment for current plot */
+	    XDrawLine(theDisplay, w, theGC, lastx, lasty, X, Y);
+
+	    points[i-1].x = lastx;
+	    points[i-1].y = lasty;
+	    lastx = X;
+	    lasty = Y;
+
+	    carry = incr - (int)incr;
+	    incr = res + carry;
+	    x += (int)incr;
+	}
+	points[i-1].x = lastx;
+	points[i-1].y = lasty;
+    }
+}
+
+
+#define BLOCK   256
+float block;            /* for averaging in downsample */
+
+
+downsample (w,fp)
+Window  w;
+FILE   *fp;
+{
+    int   i, x, j = 0;
+    short y;
+    float rate;
+    float res;
+    float incr;
+    float carry;
+
+    rate = (float)n/pwidth();
+    if ((block = rate) >= BLOCK)
+	error ("excessive downsampling; re-define BLOCK\n");
+
+    res = (float)pwidth()/(n-1);                /* resolution of plot */
+    incr = res;                                 /* x-value increment */
+    x = 0;                                      /* initial x-value */
+    for (i=0 ; i<n ; i++) {
+	if (incr >= 1.0) {
+	    if (sample(&y, i, rate)) {
+		points[j].x = x+pleft();
+		points[j].y = py(y);
+
+		/* Plot only inside y-limits of box */
+		if (points[j].y > pbottom()) points[j].y = pbottom();
+		if (points[j].y < ptop())    points[j].y = ptop();
+		j++;
+	    }
+	}
+	carry = incr - (int)incr;
+	incr = res + carry;
+	x += (int)incr;
+    }
+    return j;
+}
+
+sample (y, i, rate)
+short  *y;
+int     i;
+float   rate;
+{
+    short average();
+    float carry;
+
+    if (block <= n-i) {
+	*y = average ((int)block, &data[i]);
+	carry = block - (int)block;
+	block = rate + carry;
+	return 1;
+    }
+    else
+	return 0;
+}
+
+/*
+short
+average (block, p)
+int block;
+short p[];
+{
+    int   i;
+    float sum = 0;
+
+    for (i = 0 ; i < block ; i++)
+	sum += p[i];
+    return (short)(sum / block);
+}
+*/
+
+short
+average (block, p)
+int block;
+short p[];
+{
+    int   i;
+    float sum = 0;
+    short MIN = 32000, MAX = 0;
+
+    for (i = 0 ; i < block ; i++) {
+	sum += p[i];
+	if (p[i] < MIN) MIN = p[i];
+	if (p[i] > MAX) MAX = p[i];
+    }
+    if (summary==0) {                   /* summary of block is largest extremal value */
+	if (MAX >= (-MIN)) return MAX;     /* summary of block is MAX value */
+	else               return MIN;     /* summary of block is MIN value */
+    }
+    if (summary == 1) return (p[0]);    /* summary of block is sample value */
+    return (short)(sum / block);        /* summary of block is mean value */
+}
+
+
+/****************************************************************************
+*       Miscellaneous
+****************************************************************************/
+error(s)
+{
+    fprintf(stderr,"%s",s);
+    exit(1);
+}
+
+
+help ()
+{
+    printf ("\nx11play: animated plot of sets of n points, read as binary shorts.\n");
+    printf ("Usage: a) x11play  [options]  filename\n");
+    printf ("       b) cat filename  |  x11play  [options]\n");
+    printf ("    where filename is an input stream of 16-bit binary numbers (shorts)\n");
+    printf ("Options:\n");
+    printf ("-n [int]      = Number of points plotted (default=%d).\n", DEFSCALE);
+    printf ("-s [int]      = Offset points from start of file (default=0).\n");
+    printf ("-a [float]    = Amplitude between zero-line and top of box (default=%d).\n", DEFAMPL);
+    printf ("-y [int]      = Amplitude between zero-line and bottom of box (default=a).\n");
+    printf ("-c [float]    = Calibration at right end of x-axis (default=n+s).\n");
+    printf ("-x [int]      = Calibration at left end of x-axis (default=s).\n");
+    printf ("-z            = Suppress zero-line.\n");
+    printf ("-b [0,1,2]    = Method of block summary for down-sampling.\n");
+    printf ("                0=MAX, 1=SAMPLE, 2=MEAN, (default=MAX).\n");
+    printf ("-S            = Step-mode, (centre-button=single-step, right-button=play).\n");
+    printf ("-O            = Over-plot successive plots of n points.\n");
+    printf ("-F [font]     = Font name (default=%s)\n", theFontString);
+    exit (0);
+}
+
+