v09i062: draw, part 02/02

Brandon S. Allbery - comp.sources.misc allbery at uunet.UU.NET
Wed Dec 13 13:27:16 AEST 1989


Posting-number: Volume 9, Issue 62
Submitted-by: jta at sah.vtt.fi (Juha Takala)
Archive-name: draw_jt/part02

[This is "dr_lib", a library to produce plots in various formats including
plot(3) and HPGL formats, and "draw", a simple front-end for the library.
See the README for more information.  ++bsa]

#! /bin/sh
# This is a shell archive.  Remove anything before this line, then unpack
# it by saving it into a file and typing "sh file".  To overwrite existing
# files, type "sh file -c".  You can also feed this as standard input via
# unshar, or by typing "sh <file", e.g..  If this archive is complete, you
# will see the following message at the end:
#		"End of archive 2 (of 2)."
# Contents:  dr_lib.c draw.c
# Wrapped by allbery at uunet on Tue Dec 12 21:15:58 1989
PATH=/bin:/usr/bin:/usr/ucb ; export PATH
if test -f 'dr_lib.c' -a "${1}" != "-c" ; then 
  echo shar: Will not clobber existing file \"'dr_lib.c'\"
else
echo shar: Extracting \"'dr_lib.c'\" \(31312 characters\)
sed "s/^X//" >'dr_lib.c' <<'END_OF_FILE'
X/*
X
XCopyright (C) 1988, 1989 by Juha Takala, jta at sah.vtt.fi
X
X     This program is free software; you can redistribute it and/or modify
X     it under the terms of the GNU General Public License as published by
X     the Free Software Foundation; version 1.
X
X     This program is distributed in the hope that it will be useful,
X     but WITHOUT ANY WARRANTY; without even the implied warranty of
X     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
X     the file "License" for more details
X
X     You should have received a copy of the GNU General Public License
X     along with this program; if not, write to the Free Software
X     Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
X
X
X */
X
X
X/* Description:
X *  dr_lib.c provides interface between device level and my draw(1).
X *  Supported devices are:
X *  - hpgl output for hp7475 plotter
X *  - gl_plot(3) calls for EGA or matrix printer output
X *  - plot(3) interface
X *
X */
Xstatic char *rcsid =
X  "$Id: dr_lib.c,v 2.4 89/12/08 14:00:34 jta Exp Locker: jta $";
X
X/* define this, if your plotter connection is 'y-cable' (not direct) */
X/* #define YCABLE			/*  */
X
X/* definede this if you want to ask paper size from plotter */
X/* (should not be used, if you are going to use this package as */
X/* part of pipelines) */
X/* #define ASK			/*  */
X
X#include <stdio.h> 
X#include <math.h>
X#ifdef MSDOS
X#include <bios.h>
X#endif
X#include "dr_lib.h"
X
X#define  ESC 0x1b		/* ascii <ESC> */
X
Xstatic int xscale(), yscale();
Xstatic void sends(), invdev();
X
X#ifdef ASK
Xstatic void takes();
X#endif /* ASK */
X
X/*
X * internal global variable declarations
X *
X */
X				/* coordinate transformations */
Xstatic float xuserincr,   yuserincr;
Xstatic float xuserlow,    yuserlow;
Xstatic float xuserhigh,   yuserhigh;
Xstatic float xdevicelow,  ydevicelow;
Xstatic float xdevicehigh, ydevicehigh;
Xstatic float xrelarea,    yrelarea;
Xstatic float xfactor,     yfactor;
Xstatic float shape_factor, size_factor;	/* text 'fatness' and size factors */
Xstatic int lastx=0, lasty=0;	/* device coordinates for last point */
X
X
Xstatic int device = INVDEV;	/* what output device we are working for */
Xstatic int old_linetype = -1;	/* remember last used linetype */
Xstatic int old_pen = -1;	/* ...and pen */
X
Xstatic int errcount, warncount, errfile;
Xstatic char *errfilename = "draw.err";
X
Xstatic FILE *errfilefp;
X
X
X
X/* ********************************************************** */
X/*                                                            */
X/* Following functions are ment to be subroutines to be       */
X/* called from a user supplied main program                   */
X/*                                                            */
X/* ********************************************************** */
X
X/*
X * dr_set_cset()   - select cset
X *
X */
Xvoid dr_set_cset(n)
Xint n;
X{
X    char s[50];
X
X    switch (device)
X      {
X      case HP7475:
X	  if (n < 0 || n == 5 || (n > 9 && n < 30) || n > 39) {
X	      (void)sprintf(s, "Invalid cset code %d, using 0", n);
X	      dr_errmsg(s);
X	      n = 0;
X	  }
X	  (void)sprintf(s, "CS%d;\n", n);
X	  sends(s);
X	  break;
X#ifdef GL_LIB
X      case SCREEN:
X/*	  (void)sprintf(s, "cset not yet for SCREEN");
X	  dr_errmsg(s);		/*  */
X	  break;
X#endif
X#ifdef PLOT
X      case PLOTLIB:
X/*	  (void)sprintf(s, "cset not valid for PLOT");
X	  dr_errmsg(s);		/*  */
X	  break;
X#endif
X      default:
X	  invdev(); break;
X      }
X}
X
X/*
X * dr_set_pen()   - select pen
X *
X */
Xvoid dr_set_pen(n)
Xint n;				/* pen number, 1..6, 0=='put it away' */
X{
X    char s[50];
X
X    switch (device)
X      {
X      case HP7475:
X	  if (n < 0 || n > 6) {
X	      (void)sprintf(s, "Invalid pen code %d, using 1", n);
X	      dr_errmsg(s);
X	      n = 1;
X	  }
X	  if (n != old_pen) {	/* optimize for hpgl... */
X	      old_pen = n;
X	      (void)sprintf(s, "SP%d;", n);
X	      sends(s);
X	  }
X	  break;
X#ifdef GL_LIB
X      case SCREEN:
X/*	  (void)sprintf(s, "dr_set_pen(colour) not yet for SCREEN\n");
X	  dr_errmsg(s);		/*  */
X	  break;
X#endif
X#ifdef PLOT
X      case PLOTLIB:
X/*	  (void)sprintf(s, "dr_set_pen(colour) not valid for PLOT\n");
X	  dr_errmsg(s);		/*  */
X	  break;
X#endif
X      default:
X	  invdev(); break;
X      }
X}
X
X/*
X * dr_goto()  - move pen to point x, y
X *
X */
Xvoid dr_goto(x, y, penup)
Xfloat x, y;			/* where he wants to go */
Xint penup;			/* nonzero == lift pen before moving */
X{
X    char s[60];
X    lastx = xscale(x);
X    lasty = yscale(y);
X    switch (device)
X      {
X      case HP7475:
X	  if (penup) sends ("PU");
X	  else sends ("PA;PD");
X	      
X	  /* go to point, leave pen as told
X	   */
X	  (void)sprintf(s, "%d,%d;\n", lastx, lasty);
X	  sends(s);
X	  break;
X#ifdef GL_LIB
X      case SCREEN:
X	  if (penup)	n_movepen(lastx, lasty);
X	  else		n_draw(lastx, lasty);
X	  break;
X#endif
X#ifdef PLOT
X      case PLOTLIB:
X	  if (penup)	move(lastx, lasty);
X	  else		cont(lastx, lasty);
X	  break;
X#endif
X      default:
X	  invdev(); break;
X      }
X    return;
X}
X
X
X/*
X * dr_draw_line()  - draw line (x1,y1) .. (x2,y2)
X *
X * connect the two data points with line
X */
Xvoid dr_draw_line(x1, y1, x2, y2)
Xfloat x1, y1, x2, y2;		/* data points */
X{
X    char s[60];
X    
X    lastx = xscale(x2);
X    lasty = yscale(y2);
X    switch (device)
X      {
X      case HP7475:
X	  /* lift pen,
X	   * go to first point,
X	   * drop pen there,
X	   */
X	  (void)sprintf(s, "PA;PU%d,%d;PD;", xscale(x1), yscale(y1));
X	  sends(s);
X
X	  /* draw line segment, lift pen
X	   */
X	  (void)sprintf(s, "PD%d,%d;PU;\n", lastx, lasty);
X	  sends(s);
X	  break;
X#ifdef GL_LIB
X      case SCREEN:
X	  n_line(xscale(x1), yscale(y1), lastx, lasty);
X	  break;
X#endif
X#ifdef PLOT
X      case PLOTLIB:
X	  line(xscale(x1), yscale(y1), lastx, lasty);
X	  break;
X#endif
X      default:
X	  invdev(); break;
X      }
X    return;
X}
X
X
X/*
X * dr_draw_array()  - draw y = f(x)
X *
X * connect data points with line
X */
Xvoid dr_draw_array(x, y, n)
Xfloat x[], y[];			/* data points */
Xint n;				/* number of points in x[] and y[] */
X{
X    int i;
X    char s[60];
X    
X    /* if we really have something to draw
X     */
X    if (n < 1)
X      return;
X
X    switch (device)
X      {
X      case HP7475:
X	  /* lift pen,
X	   * go to first point,
X	   * drop pen there,
X	   */
X	  (void)sprintf(s, "PA;PU%d,%d;PD;", xscale(x[0]), yscale(y[0]));
X	  sends(s);
X	  
X	  /* draw line segments
X	   */
X	  for (i=1; i<n; i++) {
X	      if (i % 4 == 0) sends("\n");
X	      (void)sprintf(s, "PD%d,%d;", xscale(x[i]), yscale(y[i]));
X	      sends(s);
X	  }
X	  
X	  /* finally lift pen
X	   */
X	  sends("PU;\n");
X	  break;
X#ifdef GL_LIB
X      case SCREEN:
X	  n_movepen(xscale(x[0]), yscale(y[0]));
X	  for (i=1; i<n; i++)
X	    n_draw(xscale(x[i]), yscale(y[i]));
X	  break;
X#endif
X#ifdef PLOT
X      case PLOTLIB:
X	  move(xscale(x[0]), yscale(y[0]));
X	  for (i=1; i<n; i++)
X	    cont(xscale(x[i]), yscale(y[i]));
X	  break;
X#endif
X      default:
X	  invdev(); break;
X      }
X    lastx = xscale(x[i-1]);
X    lasty = yscale(y[i-1]);
X    return;
X}
X
X
X/*
X * dr_put_mark()  - make a mark at (x,y)
X *
X * put mark at data point
X */
Xvoid dr_put_mark(x, y, mark) 
Xfloat x, y;			/* data point */
Xchar mark;
X{
X    char s[60];
X
X    lastx = xscale(x);
X    lasty = yscale(y);
X    switch (device)
X      {
X      case HP7475:
X	  /* if we really have something to plot
X	   */
X	  
X	  /* lift pen,
X	   * set symbol mode,
X	   */
X	  (void)sprintf(s, "PA;PU;SM%c;", mark);
X	  sends(s);
X	  
X	  /* put symbol at all data points
X	   */
X	  (void)sprintf(s, "PU%d,%d;", lastx, lasty);
X	  sends(s);
X	  
X	  /* reset symbol mode
X	   */
X	  sends("SM;\n");
X	  break;
X#ifdef GL_LIB
X      case SCREEN:
X	  n_movepen(lastx, lasty);
X	  n_grafchar(mark);
X	  break;
X#endif
X#ifdef PLOT
X      case PLOTLIB:
X	  move(lastx, lasty);
X	  *s = mark;
X	  *(s+1) = '\0';
X	  label(s);
X	  break;
X#endif
X      default:
X	  invdev(); break;
X      }
X    return;
X}
X
X
Xvoid dr_set_shape(x)
Xfloat x;
X{
X    if (x > 0.0)
X      shape_factor = x * xrelarea/yrelarea;
X    return;
X}
X
X
Xvoid dr_set_size(x)
Xfloat x;
X{
X    if (x > 0.0)
X      size_factor = x;
X    return;
X}
X
X
X/*
X * dr_put_text()  - put some text
X * plot(3) only has one size and orientation, what can you do!
X *
X */
Xvoid dr_put_text(x, y, xofs, yofs, angle, size, s)
Xfloat x, y;			/* where to put it */
Xfloat xofs, yofs;		/* offset of start position, measured in */
X				/* character size units */
Xfloat angle;			/* radians, 0.0 == horizontal */
Xfloat size;			/* relative size of text */
Xchar *s;			/* the text */
X{
X    char ss[60];
X    int senddir;
X    char dir[30];
X
X    if (size <= 0.0) {
X	sprintf(ss, "invalid size: %f");
X	dr_errmsg(ss);
X	return;
X    }
X    size *= size_factor;
X
X    switch (device)
X      {
X      case HP7475:
X	  senddir = (angle != 0.0);
X	  if (senddir) {
X	      if (angle == M_PI_2)
X		sprintf(dir,"DI0,1");
X	      else if (angle == -M_PI_2)
X		sprintf(dir,"DI0,-1");
X	      else
X		sprintf(dir,"DI1,%f", tan(angle));
X	  }
X
X	  /* lift pen, go to specified position, set direction & size,
X	   */
X
X	  (void)sprintf(ss, "PA;PU%d,%d;%s;SR%4.2f,%4.2f;\n",
X			xscale(x), yscale(y),
X			(senddir ? dir : ""), 
X			size * 0.75 * shape_factor * xrelarea,
X			size * 1.5 * yrelarea);
X	  sends(ss);
X
X	  /* go to offset,
X	   * put text,
X	   * restore 0 degree direction.
X	   */
X	  if ((xofs != 0.0) || (yofs != 0.0)) {
X	      (void)sprintf(ss, "CP%4.2f,%4.2f;\n", xofs, yofs);
X	      sends(ss);
X	  }
X	  sends("LB"); sends(s); sends("\03\n");
X
X	  /* correction for y-offset, if necessary, to leave pen on same */
X	  /* height but after the text */
X	  if (yofs != 0.0) {
X	      (void)sprintf(ss, "CP0,%4.2f;\n", -yofs);
X	      sends(ss);
X	  }
X
X	  if (senddir)
X	    sends("DI1,0;\n");
X/*	  lastx = ....;
X	  lasty = ....;		/* hp7475 does not really need these, */
X				/* it remembers the pen position */
X
X	  break;
X#ifdef GL_LIB
X      case SCREEN:
X	  n_movepen(xscale(x), yscale(y));
X
X#define GL_MODE_255_HAS_BUG
X#ifdef GL_MODE_255_HAS_BUG
X	  /* DT Lewis!!  Wake up! */
X	  /* For some reason gl_lib works strangely with GL_MODE == 255 */
X	  /* (the output is  (almost) microscopic small and weird shape). */
X	  /* We can compensate that strangeness here by modifying size */
X	  /* and sape factors for g_fontctl()...  The values 3.0 and 0.4 */
X	  /* are experimental results (by kekkonen at sah.vtt.fi). */
X	  if (strncmp(getenv("GLMODE"),"256",3)==0)
X	    g_fontctl(3.0*size*yrelarea, 0.4*shape_factor, 1.0, angle, 0.0);
X	  else
X#endif
X#undef GL_MODE_255_HAS_BUG
X	    g_fontctl(size*yrelarea, shape_factor, 1.0, angle, 0.0);
X	  yofs = -0.5 - yofs;
X	  n_ch_ofs(xofs, yofs);
X	  n_grafstr(s);
X	  n_ch_ofs(0.0, -yofs);
X	  g_fontctl(size, 1.0, 1.0, 0.0, 0.0);
X/*	  lastx = ....;
X	  lasty = ....;		/* this does not really need these, */
X				/* it remembers the pen position */
X
X	  break;
X#endif
X#ifdef PLOT
X#define PLOT_CH_XSIZE 70
X#define PLOT_CH_YSIZE 100
X      case PLOTLIB:
X	  {
X	      int i;
X	      float xsize, ysize;
X	      float xbase, ybase;
X	      float xincr, yincr;
X	      float xofs1, yofs1;
X
X	      xsize = PLOT_CH_XSIZE * size * xrelarea;
X	      ysize = PLOT_CH_YSIZE * size * yrelarea;
X
X	      yofs += 0.4;	/* middle --> botom */
X	      xofs1 = (xofs * xsize * cos(angle) * shape_factor -
X		       yofs * ysize * sin(angle));
X	      yofs1 = (xofs * xsize * sin(angle) * shape_factor +
X		       yofs * ysize * cos(angle));
X	      xbase = xscale(x) + xofs1;
X	      ybase = yscale(y) + yofs1;
X
X	      xincr = cos(angle) * xsize * shape_factor;
X	      yincr = sin(angle) * ysize;
X	      
X	      /* This is how we simulate different printing directions */
X	      /* and character sizes: by moving to the position of every */
X	      /* single character one at a time */
X	      *(ss+1) = '\0';
X	      for (i=0; *ss=*s; s++,i) {
X		  move((int) (0.5 + xbase + i * xincr),
X		       (int) (0.5 + ybase + i * yincr));
X		  label(ss);
X		  i++;
X	      }
X
X	      lastx = (int) (0.5 + xbase + i * xincr - xofs1);
X	      lasty = (int) (0.5 + ybase + i * yincr - yofs1);
X	      move(lastx, lasty);
X	  }
X	  break;
X#endif
X      default:
X	  invdev(); break;
X      }
X    return;
X}
X
X/*
X * dr_small_circle()  - draw small circle at current point.
X * This routine needs the help of lastx and lasty in case of plot(3) library.
X *
X */
Xvoid dr_draw_circle()
X{
X    switch (device)
X      {
X      case HP7475:
X	  sends("CI35;\n");	  break;
X#ifdef GL_LIB
X      case SCREEN:
X	  n_ellipse1(100, 100);	  break;
X#endif
X#ifdef PLOT
X      case PLOTLIB:
X	  circle(lastx, lasty, 20);	  break;
X#endif
X      default:
X	  invdev(); break;
X      }
X}
X
X
X/*
X * dr_finish()  - terminate the plotting
X *
X * This routine must be called after all plotting activity is done
X * to restore the Y-cable connection so that terminal gets control
X */
Xint dr_finish(clean_scr)
Xint clean_scr;			/* != 0 if we have clear screen */
X{
X    char s[20];
X
X    switch (device)
X      {
X      case HP7475:
X	  /* deselct pen,
X	   * reset Y-cable connection
X	   *   (<esc>.Z is ignored if not Y-cable connecton )
X	   */
X	  (void)sprintf(s, "SP;%c.Z\n", ESC);
X	  sends(s);
X	  break;
X#ifdef GL_LIB
X      case SCREEN:
X	  if (clean_scr && (strncmp(getenv("GLMODE"), "256", 3) != 0)) {
X#ifdef UNIX
X	      FILE *fp = fopen("/dev/tty", "r");
X	      fgets(s, 2, fp);	/* wait for user response */
X	      fclose(fp);
X#endif /* UNIX */
X#ifdef MSDOS
X	      _bios_keybrd(_KEYBRD_READ); /* wait for keybrd. input forever */
X#endif /* MSDOS */
X	  }
X	  {
X	      int status;
X	      status = g_finish();	/* do clean screen */
X	      if (status != 0)
X		dr_errmsg("dr_lib(): problems with g_finish().");
X	  }
X	  break;
X#endif
X#ifdef PLOT
X      case PLOTLIB:
X	  move(0,0);
X	  closepl();
X	  break;
X#endif
X      default:
X	  break;		/* call to invdev() --> recursion */
X      }
X	  
X    /* check errors, tell if any or remove error log file
X     */
X    if (errcount == 0 && warncount == 0)
X      (void)unlink(errfilename);
X    else {
X	if (errcount)
X	  (void)fprintf(stderr, "%d error%s,",
X			errcount, (errcount == 1 ? "" : "s"));
X	if (warncount)
X	  (void)fprintf(stderr, "%d warning%s,",
X			warncount, (warncount == 1 ? "" : "s"));
X	(void)fprintf (stderr, " see file \"%s\".\n", errfilename);
X    }
X    /* return something useful
X     */
X    return (errcount);
X}
X
X
X/*
X * dr_start  - init the plotting device
X *
X * This routine must be called before any plotting activity is done
X * to set up the Y-cable connction so that plotter starts processing
X * or in case of direct connection, open the device
X */
Xvoid dr_start(xlow, ylow, xhigh, yhigh, dev)
Xint dev;			/* plotting device & paper size */
Xfloat xlow, xhigh;		/* drawing area, relative to.. */
Xfloat ylow, yhigh;		/* ..whole paper, that means: */
X				/* 0.0, 0.5, 0.5, 1.0 will use */
X				/* upper left quarter of paper */
X{
X    int p1x, p1y, p2x, p2y;
X    float t1, t2;
X    char s[200];
X
X    /* open error log file
X     */
X    errcount = 0;
X    warncount = 0;
X    if ((errfilefp = fopen(errfilename, "w")) == NULL) {
X	(void)fprintf(stderr, "\nCan't open %s.  Aborting.\n", errfilename);
X	exit(1);
X    }
X    errfile = fileno(errfilefp);
X    xrelarea = xhigh - xlow;
X    yrelarea = yhigh - ylow;
X    shape_factor = yrelarea/xrelarea; /* back to normal 'fatness' */
X    size_factor = 1.0;		/* actual size normalized to area */
X
X    switch (device = (dev & DEVMASK))
X      {
X      case HP7475:
X	  /* set Y-cable connection,
X	   *   (<esc>.Y is ignored if not Y-cable connecton )
X	   * init the plotter,
X	   * set handshake mode x-on/x-off,
X	   * select pen #1
X	   */
X	  (void)sprintf(s, "\n%c.Y;IN;", ESC); /* Y-connection, plotter on */
X	  sends(s);
X	  (void)sprintf(s, "%c.I80;;17:", ESC);	/* enable x-on/x-off x-on=^q */
X	  sends(s);
X	  (void)sprintf(s, "%c.N;19:\n", ESC);	/* x-off=^s */
X	  sends(s);
X#ifdef ASK
X	  /* ask picture area from plotter
X	   * set up limits (leave some room for axis labels etc.) */
X	  sends("OP;");
X	  takes(s);
X#else /* ASK */
X	   /* we would get one of this */
X	  switch (dev & SIZEMASK)
X	    {
X	    case SIZE_A:
X		(void)sprintf (s, "250,596,10250,7796\n");
X		sends("PS4;\n");
X		break;
X	    case SIZE_A4:
X		(void)sprintf (s, "603,521,10603,7721\n");
X		sends("PS4;\n");
X		break;
X	    case SIZE_B:
X		(void)sprintf (s, "522,259,15722,10259\n");
X		sends("PS3;\n");
X		break;
X	    case SIZE_A3:
X		(void)sprintf (s, "170,602,15370,10602\n");
X		sends("PS3;\n");
X		break;
X	    default:
X		dr_errmsg("Invalid paper size.  Aborting.");
X		(void)dr_finish();
X		exit(1);
X	    }
X#endif /* ASK */
X	  if ( sscanf(s, "%d,%d,%d,%d", &p1x, &p1y, &p2x, &p2y) == 4 ) {
X				/* some room for labels */
X	      t1 = p1x + 0.05*(p2x-p1x); /* y-label */
X	      t2 = p2x - 0.001*(p2x-p1x); /* right margin */
X	      xdevicelow  = t1 + xlow * (t2 - t1);
X	      xdevicehigh = t1 + xhigh * (t2 - t1);
X	      t1 = p1y + 0.05*(p2y-p1y); /* x-label */
X	      t2 = p2y - 0.07*(p2y-p1y); /* name */
X	      ydevicelow  = t1 + ylow * (t2 - t1);
X	      ydevicehigh = t1 + yhigh * (t2 - t1);
X	  }
X	  else {
X	      dr_errmsg("Can't receive plotter dimensions.  Aborting.");
X	      (void) dr_finish();
X	      exit(1);
X	  }
X	  
X	  dr_set_pen(1);		/* select some pen */
X	  break;
X#ifdef GL_LIB
X      case SCREEN:
X	  p1x = 0;	p1y = 32767;
X	  p2x = 32767;	p2y = 0;
X	  t1 = p1x + 0.06*(p2x-p1x); /* room y-axis label */
X	  t2 = p2x - 0.001*(p2x-p1x); /* no room at right margin */
X	  xdevicelow  = t1 + xlow * (t2 - t1);
X	  xdevicehigh = t1 + xhigh * (t2 - t1);
X	  t1 = p1y + 0.10*(p2y-p1y); /* room for x-axis label */
X	  t2 = p2y - 0.05*(p2y-p1y); /* room fo picture's name */
X	  ydevicelow  = t1 + ylow * (t2 - t1);
X	  ydevicehigh = t1 + yhigh * (t2 - t1);
X
X#define ENV_MODE 0		/* from gl.h; it means that g_init() will */
X				/* look operating mode in environment */
X				/* variable GLMODE */
X	  g_init(ENV_MODE);
X#undef ENV_MODE
X	  /* dr_set_pen(1);		/* select some pen */
X	  break;
X#endif
X
X#ifdef PLOT
X      case PLOTLIB:
X				/* Y-cable for Tektornix? */
X	  p1x = 0;	p1y = 0;
X	  p2x = 5450;	p2y = 4095;
X				/* 5% room for labels */
X	  t1 = p1x + 0.07*(p2x-p1x); /* y-label */
X	  t2 = p2x - 0.02*(p2x-p1x);; /* righ margin */
X	  xdevicelow  = t1 + xlow * (t2 - t1);
X	  xdevicehigh = t1 + xhigh * (t2 - t1);
X	  t1 = p1y + 0.10*(p2y-p1y); /* x-label */
X	  t2 = p2y - 0.05*(p2y-p1y); /* name */
X	  ydevicelow  = t1 + ylow * (t2 - t1);
X	  ydevicehigh = t1 + yhigh * (t2 - t1);
X
X	  openpl();
X	  space(0,0,4096,4096);
X	  erase();
X	  break;
X#endif
X      default:
X	  invdev(); break;
X      }
X    
X    /* check coordinate position legality
X     */
X    if (xlow < -0.0001 || xhigh > 1.0001 ||
X	ylow < -0.0001 || yhigh > 1.0001) {
X	dr_errmsg("Illegal x/y low/high position.  Aborting.");
X	(void)dr_finish();
X	exit(1);
X    }
X    return;
X}
X
X
X/* 
X * dr_xgrid  - draw x grid 
X *
X * Draw grid lines parallel to y-axis.  When argument n == 0, makes
X * lines at every tic-mark, n == 1 makes one grid line between every
X * tic mark etc.  To leave grid lines away, just don't call this.
X */
Xvoid dr_xgrid(n)
Xint n;				/* number of grid lines */
X				/* between tic-marks */
X{
X    int i;
X    int x, y1, y2;		/* temps for grid line positions */
X    int nlines;			/* number of grid lines to draw */
X    char s[200];
X
X    /* check interval
X     */
X    if (n < 1) {
X	dr_errmsg("Invalid x-grid interval, using 1");
X	n = 1;
X    }
X
X    /* how many lines we have?
X     */
X    nlines = (int) (n*(xuserhigh-xuserlow)/xuserincr - 0.5);
X
X     switch (device)
X      {
X      case HP7475:
X	  /*  select absolute moving
X	   */
X	  sends("PA;\n");
X#ifdef SLOW
X	  if (slow) {
X	  }
X#endif	  
X	  /* draw desired number of lines, sawing up and down
X	   */
X	  y1 = yscale(yuserlow);
X	  y2 = yscale(yuserhigh);
X	  for (i=1; i<=nlines; i++) {
X				/* go up */
X	      x = xscale(i*xuserincr/n+xuserlow);
X	      (void)sprintf(s, "PU%d,%d;PD%d,%d;", x, y1, x, y2);
X	      sends(s);
X				/* go down */
X	      if (++i <= nlines) {
X		  x = xscale(i*xuserincr/n+xuserlow);
X		  (void)sprintf(s, "PU%d,%d;PD%d,%d;\n", x, y2, x, y1);
X		  sends(s);
X	      }
X	  }
X	  
X	  /* lift pen
X	   */
X	  sends("PU;\n");
X	  break;
X#ifdef GL_LIB
X      case SCREEN:
X	  y1 = yscale(yuserlow);
X	  y2 = yscale(yuserhigh);
X	  for (i=1; i<=nlines; i++) {
X	      x = xscale(i*xuserincr/n+xuserlow);
X	      n_line(x,y1, x,y2);
X	  }
X	  break;
X#endif
X#ifdef PLOT
X      case PLOTLIB:
X	  y1 = yscale(yuserlow);
X	  y2 = yscale(yuserhigh);
X	  for (i=1; i<=nlines; i++) {
X	      x = xscale(i*xuserincr/n+xuserlow);
X	      line(x,y1, x,y2);
X	  }
X	  break;
X#endif
X      default:
X	  invdev(); break;
X      }
X}
X
X
X/* 
X * dr_ygrid  - draw y grid 
X *
X * See dr_xgrid().
X */
Xvoid dr_ygrid(n)
Xint n;
X{
X    int i;
X    int nlines;			/* number of lines */
X    char s[200];
X    int x1, x2, y;
X
X    /* check interval
X     */
X    if (n < 1) {
X	dr_errmsg("Invalid y-grid interval, using 1");
X	n = 1;
X    }
X
X    /* how many lines do we have?
X     */
X    nlines = (int) (n*(yuserhigh-yuserlow)/yuserincr - 0.5);
X
X    switch (device)
X      {
X      case HP7475:
X	  /* select absolute moving
X	   */
X	  sends("PA;\n");
X	  
X	  /* draw desired number of lines, sawing right and left
X	   */
X	  x1 = xscale(xuserlow);
X	  x2 = xscale(xuserhigh);
X	  for (i=1; i<=nlines; i++) {
X				/* go right */
X	      y = yscale(yuserlow+i*yuserincr/n);
X	      (void)sprintf(s, "PU%d,%d;PD%d,%d;", x1, y, x2, y);
X	      sends(s);
X				/* go left */
X	      if (++i <= nlines) {
X		  y = yscale(yuserlow+i*yuserincr/n);
X		  (void)sprintf(s, "PU%d,%d;PD%d,%d;", x2, y, x1, y);
X		  sends(s);
X	      }
X	      sends("\n");
X	  }
X	  
X	  /* lift pen
X	   */
X	  sends("PU;\n");
X	  break;
X#ifdef GL_LIB
X      case SCREEN:
X	  x1 = xscale(xuserlow);
X	  x2 = xscale(xuserhigh);
X	  for (i=1; i<=nlines; i++) {
X	      y = yscale(i*yuserincr/n+yuserlow);
X	      n_line(x1,y, x2,y);
X	  }
X	  break;
X#endif
X#ifdef PLOT
X      case PLOTLIB:
X	  x1 = xscale(xuserlow);
X	  x2 = xscale(xuserhigh);
X	  for (i=1; i<=nlines; i++) {
X	      y = yscale(i*yuserincr/n+yuserlow);
X	      line(x1,y, x2,y);
X	  }
X	  break;
X#endif
X      default:
X	  invdev(); break;
X      }
X}
X
X
X/* 
X * dr_area_def  - establish x and y axis 
X *
X * [xy]low    lower left corner
X * [xy]high   upper right corner
X * [xy]incr   increment of tic-marks
X * no_axes    flag, != 0 means don't draw axes
X *
X * Coordinate values will be printed at every tic-mark, so useful
X * values for incr is such that (high-low)/incr will be some small
X * integer.  With functions dr_xgrid() and dr_ygrid() it is possible to
X * draw grid lines also in between these tic-marcs.
X */
Xvoid dr_area_def(xlow, xincr, xhigh, ylow, yincr, yhigh, no_axes)
Xfloat xlow, xincr, xhigh;
Xfloat ylow, yincr, yhigh;
Xint no_axes;
X{
X    int i;
X    char units[30];		/* place for units */
X    char s[200];
X    float x, y;
X    int ix, iy, nx, ny;
X
X    /* set the globals
X     */
X    xuserlow = xlow;   xuserhigh = xhigh;   xuserincr = xincr;
X    yuserlow = ylow;   yuserhigh = yhigh;   yuserincr = yincr;
X    xfactor = (xdevicehigh - xdevicelow) / (xuserhigh - xuserlow);
X    yfactor = (ydevicehigh - ydevicelow) / (yuserhigh - yuserlow) ;
X
X    /* check legality,
X     * abort on error
X     */
X    if ((xincr == 0.0) ||
X	(xuserlow == xuserhigh) ||
X	((xuserhigh-xuserlow)*xincr < 0.0)) {
X	dr_errmsg("x-axis defined improperly. Aborting.");
X	(void)dr_finish();
X	exit(1);
X    }
X    if ((yincr == 0.0) ||
X	(yuserlow == yuserhigh) ||
X	((yuserhigh-yuserlow)*yincr < 0.0)) {
X	dr_errmsg("y-axis defined improperly. Aborting.");
X	(void)dr_finish();
X	exit(1);
X    }
X    if (no_axes)
X      return;
X    
X    switch (device)
X      {
X      case HP7475:
X
X	  /* draw box
X	   * --------
X	   */
X	  
X	  /* select absolute plotting mode,
X	   * lift pen,
X	   * go to lower left corner
X	   */
X	  (void)sprintf(s,"PA;PU%d,%d;\n",
X			xscale(xlow), yscale(ylow));
X	  sends(s);
X	  
X	  /* draw x-axis with x-tick marks
X	   */
X	  
X	  (void)sprintf(s, "TL%.1f,%.1f;\n", 0.5*xrelarea, 0.5*yrelarea);
X	  sends(s);
X	  for (nx=0, x = xlow;
X	       x < xhigh + xincr/100.0;
X	       x += xincr, nx++) {
X	      (void)sprintf(s, "PD%d,%d;XT;\n", xscale(x), yscale(ylow));
X	      sends(s);
X	  }
X	  (void)sprintf(s, "PD%d,%d;\n", xscale(xhigh), yscale(ylow));
X	  sends(s);
X	  
X	  /* draw y-axis with y-tick marks
X	   */
X	  (void)sprintf(s, "PU%d,%d;\n",
X			xscale(xlow), yscale(ylow));
X	  sends(s);
X	  for (ny=0, y = ylow;
X	       y < yhigh + yincr/100.0;
X	       y += yincr, ny++) {
X	      (void)sprintf(s, "PD%d,%d;YT;\n", xscale(xlow), yscale(y));
X	      sends(s);
X	  }
X	  (void)sprintf(s, "PD%d,%d;\n", xscale(xlow), yscale(yhigh));
X	  sends(s);
X
X	  /* outer borders of the box
X	   */
X	  (void)sprintf(s, "PU%d,%d;PD%d,%d,%d,%d;PU\n",
X			xscale(xlow), yscale(yhigh),
X			xscale(xhigh), yscale(yhigh),
X			xscale(xhigh), yscale(ylow));
X	  sends(s);
X	  break;
X#ifdef GL_LIB
X#define GL_TIC_SIZE 200
X      case SCREEN:
X	  iy = yscale(ylow);
X	  n_movepen(xscale(xlow), iy);
X
X	  for (nx=0, x = xlow;
X	       x < xhigh + xincr/100.0;
X	       x += xincr, nx++) {
X	      ix = xscale(x);
X	      n_draw(ix, iy);
X	      n_line(ix, iy+GL_TIC_SIZE, ix, iy-GL_TIC_SIZE);
X	      n_movepen(ix, iy);
X	  }
X	  n_draw(xscale(xhigh), iy);
X
X	  /* draw y-axis with y-tick marks
X	   */
X	  ix = xscale(xlow);
X	  n_movepen(ix, yscale(ylow));
X	  for (ny=0, y = ylow;
X	       y < yhigh + yincr/100.0;
X	       y += yincr, ny++) {
X	      iy = yscale(y);
X	      n_draw(ix, iy);
X	      n_line(ix+GL_TIC_SIZE, iy, ix-GL_TIC_SIZE, iy);
X	      n_movepen(ix, iy);
X	  }
X	  n_draw(ix, yscale(yhigh));
X	  
X	  /* outer borders of the box
X	   */
X	  n_movepen(xscale(xlow), yscale(yhigh));
X	  n_draw(xscale(xhigh), yscale(yhigh));
X	  n_draw(xscale(xhigh), yscale(ylow));
X	  break;
X#endif
X#ifdef PLOT
X#define PLOT_TIC_SIZE 40
X      case PLOTLIB:
X	  iy = yscale(ylow);
X	  move(xscale(xlow), iy);
X
X	  for (nx=0, x = xlow;
X	       x < xhigh + xincr/100.0;
X	       x += xincr, nx++) {
X	      ix = xscale(x);
X	      cont(ix, iy);
X	      line(ix, iy+PLOT_TIC_SIZE, ix, iy-PLOT_TIC_SIZE);
X	      move(ix, iy);
X	  }
X	  cont(xscale(xhigh), iy);
X
X	  /* draw y-axis with y-tick marks
X	   */
X	  ix = xscale(xlow);
X	  move(ix, yscale(ylow));
X	  for (ny=0, y = ylow;
X	       y < yhigh + yincr/100.0;
X	       y += yincr, ny++) {
X	      iy = yscale(y);
X	      cont(ix, iy);
X	      line(ix+PLOT_TIC_SIZE, iy, ix-PLOT_TIC_SIZE, iy);
X	      move(ix, iy);
X	  }
X	  cont(ix, yscale(yhigh));
X	  
X	  /* outer borders of the box
X	   */
X	  move(xscale(xlow), yscale(yhigh));
X	  cont(xscale(xhigh), yscale(yhigh));
X	  cont(xscale(xhigh), yscale(ylow));
X	  break;
X#endif
X      default:
X	  invdev(); break;
X      }
X
X    
X    /* put y-units
X     */
X    for (i=ny-1; i>=0; i--) {
X 	(void)sprintf(units, "%g", ylow + i * yincr);
X 	dr_put_text(xlow, ylow+i*yincr,
X		    -strlen(units)/2.0, -1.3, -M_PI_2, 1.0, units);
X    }
X
X    /* put x-units
X     */
X    for (i=0; i<nx; i++) {
X	(void)sprintf(units, "%g", xlow + i * xincr);
X	dr_put_text(xlow+i*xincr, ylow,
X		    -strlen(units)/2.0, -1.3, 0.0, 1.0, units);
X    }
X    return;
X}
X
X/*
X * linetypes
X *
X */
X#ifdef GL_LIB
Xstatic long gl_line_types[] = {	/* 76543210765432107654321076543210 */
X    0xf0000000L,		/* ====............................ */
X    0xffff0000L,		/* ================................ */
X    0xffffff00L,		/* ========================........ */
X    0xffffff18L,		/* ========================...==... */
X    0xfffff0f0L,		/* ====================....====.... */
X    0xfffff39cL,		/* ====================..===..===.. */
X    0xffffffffL,		/* ================================ */
X    0x66666666L,		/* ==..==..==..==..==..==..==..==.. */
X    0xf8f8f8f8L,		/* =====...=====...=====...=====... */
X    0xffffff00L,		/* ========================........ */
X    0xfc30fc30L };		/* ======....==....======....==.... */
X#define DEF_GL_TYPE 7		/* this will produce continous line */
X#endif /* GL_PLOT */
X#ifdef PLOT
Xstatic char *plot_line_types[] = {
X    "dotted",			/* 1 */
X    "longdashed",		/* 2 */
X    "shortdashed",		/* 3 */
X    "dotdashed",		/* 4 */
X    "solid",			/* 5 */
X    "solid",			/* 6 */
X    "solid" };			/* 7 */
X#define DEF_PLOT_TYPE 7		/* this will produce continous line */
X#endif /* PLOT */
X
Xint dr_set_ltype(linetype)
Xint linetype;
X{
X    char s[40];
X    int old_val = old_linetype;
X
X    switch (device)
X      {
X      case HP7475:
X	  if (linetype < 0 || linetype > 6) {
X	      linetype = 7;
X	      sends("LT;");
X	  }
X	  else {
X	      (void)sprintf(s, "LT%d;", linetype);
X	      sends(s);
X	  }
X	  break;
X#ifdef GL_LIB
X      case SCREEN:
X	  if (linetype < 1 ||
X	      linetype > (sizeof(gl_line_types) / sizeof(gl_line_types[0]))) {
X	       sprintf (s, "Invalid linetype: %d", linetype);
X	       dr_errmsg(s);
X	       linetype = DEF_GL_TYPE;
X	   }
X	  g_style(gl_line_types[linetype-1]);
X	  break;
X#endif
X#ifdef PLOT	    
X      case PLOTLIB:
X	  if (linetype < 1 ||
X	      linetype > (sizeof(plot_line_types) /
X			  sizeof(plot_line_types[0]))) {
X	       sprintf (s, "Invalid linetype: %d", linetype);
X	       dr_errmsg(s);
X	       linetype = DEF_PLOT_TYPE;
X	   }
X	  linemod(plot_line_types[linetype-1]);
X	  break;
X#endif
X      default:
X	  invdev(); break;
X      }
X    old_linetype = linetype;
X    return old_val;
X}
X
X/*
X * xscale  - scaling in horizontal directon
X *
X * convert user x-coordinate to physical device x-coordinate
X */
Xstatic int xscale(x)
Xfloat x;
X{
X    int temp;
X    char s[60];
X
X    /* calculate it
X     */
X    temp = (int) ((x - xuserlow) * xfactor + xdevicelow + 0.5);
X    /* check, if it is good
X     */
X    if ((temp >= -32768) && (temp <= 32767))
X      return (temp);
X
X    /* when not good make report,
X     * return something
X     */
X    (void)sprintf(s, "x-value out of range: %f", x);
X    dr_warnmsg(s);
X    return (int) (xdevicelow + 0.5);
X}
X
X
X/*
X * yscale  - scaling in vertical directon
X *
X * convert user y-coordinate to physical device y-coordinate
X */
Xstatic int yscale(y)
Xfloat y;
X{
X    int temp;
X    char s[60];
X
X    /* calculate it
X     */
X    temp = (int) ((y - yuserlow) * yfactor + ydevicelow + 0.5);
X
X    /* check, if it is good
X     */
X    if ((temp >= -32768) && (temp <= 32767))
X      return (temp);
X
X    /* when not good make report,
X     * return something
X     */
X    (void)sprintf(s, "y-value out of range: %f", y);
X    dr_warnmsg(s);
X    return (int) (ydevicelow + 0.5);
X}
X
Xstatic void invdev()
X{
X    dr_errmsg("Invalid plotting device.");
X    (void)dr_finish();
X    exit(1);
X}
X
X
X/*
X * dr_errmsg  - send a string to error output
X *
X */
Xvoid dr_errmsg(s)
Xchar *s;
X{
X    char *p = s + strlen(s);
X    if (p != s) p--;
X    while ((p != s) && (*p == '\n')) /* strip the newlines */
X      *(p--) = '\0';
X    fprintf(errfilefp, "dr_lib error: %s\n", s);
X    errcount++;
X    return;
X}
X
X
X/*
X * warnmsg  - send a string to error output
X *
X */
Xvoid dr_warnmsg(s)
Xchar *s;
X{
X    char *p = s + strlen(s);
X    if (p != s) p--;
X    while ((p != s) && (*p == '\n')) *(p--) = '\0';
X    fprintf(errfilefp, "Warning: %s\n", s);
X    warncount++;
X    return;
X}
X
X
X/*
X * sends  - send something to plotting device
X *
X */
Xstatic void sends(s)
Xchar *s;
X{
X    (void)printf("%s", s);
X}
X
X#ifdef ASK
X/*
X * takes  - take a string from plotter (to read plotter coordinates)
X *
X */
X#include <signal.h>
X#include <setjmp.h>
X
Xjmp_buf sjbuf;
Xint (*savAlrm)();		/* save incomming alarm function */
X
Xstatic int timerh()		/* timer interrupt handler */
X{
X    longjmp(sjbuf, 1);
X}
X
Xstatic void takes(s)
Xchar *s;
X{
X#ifdef YCABLE
X    savAlrm = signal(SIGALRM,timerh);
X    alarm(2);			/* give time to wake up */
X    if (setjmp(sjbuf)) {
X	alarm(0);
X	signal(SIGALRM,savAlrm); /* cancel timer */
X	(void)dr_finish();
X	(void)fprintf(stderr, "\n\7Timeout reading stdin.  Aborting.\n");
X	exit(1);
X    }
X    gets(s);
X    alarm(0);
X    signal(SIGALRM,savAlrm);	/* cancel timer */
X#else
X    /* do it yorself */
X#endif
X}
X#endif /* ASK */
END_OF_FILE
if test 31312 -ne `wc -c <'dr_lib.c'`; then
    echo shar: \"'dr_lib.c'\" unpacked with wrong size!
fi
# end of 'dr_lib.c'
fi
if test -f 'draw.c' -a "${1}" != "-c" ; then 
  echo shar: Will not clobber existing file \"'draw.c'\"
else
echo shar: Extracting \"'draw.c'\" \(17080 characters\)
sed "s/^X//" >'draw.c' <<'END_OF_FILE'
X/*
X
XCopyright (C) 1988, 1989 by Juha Takala, jta at sah.vtt.fi
X
X     This program is free software; you can redistribute it and/or modify
X     it under the terms of the GNU General Public License as published by
X     the Free Software Foundation; version 1.
X
X     This program is distributed in the hope that it will be useful,
X     but WITHOUT ANY WARRANTY; without even the implied warranty of
X     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
X     the file "License for more details
X
X     You should have received a copy of the GNU General Public License
X     along with this program; if not, write to the Free Software
X     Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
X
X
X */
X
X/* Description:
X *  Non-interactive user interface program to provide interface between
X *  device manipulation subroutines and user supplied command files.
X *  This program will accept all kinds of command line arguments and
X *  pseudo comments mixed in between the data lines.  See the man page.
X *
X *  The good thing about this program is that it will accept infinite
X *  amount of input data points; they are not stored, but processed one
X *  at a time.
X */
Xstatic char *RCS_id =
X  "$Id: draw.c,v 2.7 89/12/08 14:00:39 jta Exp Locker: jta $";
X
X#include <stdio.h>
X#ifdef MSDOS
X#include <stdlib.h>
X#endif /* MSDOS */
X#include <string.h>
X#include <sys/types.h>
X#include <time.h>
X#include <math.h>
X#include <signal.h>
X#ifndef MSDOS
X#include <sys/inode.h>
X#include <sys/file.h>
X#endif
X#include <string.h>
X#include "dr_lib.h"
X
Xdouble strtod();
X
Xchar *myname;
Xfloat xlowl = 0.0;		/* relative plotting area on paper */
Xfloat ylowl = 0.0;
Xfloat xhighr= 1.0;
Xfloat yhighr = 1.0;
X
Xfloat xlow = -1.0;		/* x axis range */
Xfloat xincr = 0.5;
Xfloat xhigh = 1.0;
Xint xgrid_n = 0;
Xchar *xlabel_text = NULL;
X
Xfloat ylow = -1.0;		/* y axis range */
Xfloat yincr = 0.5;
Xfloat yhigh = 1.0;
Xint ygrid_n = 0;
Xchar *ylabel_text = NULL;
X
Xstatic float x, y;		/* current point, needs to be */
X				/* remembered, eg. when relative text is */
X				/* inserted by #$texthere */
Xchar *name_text = NULL;		/* name of picture */
Xint pennum = 1;
Xint penup = 1;
Xint no_axes = 0;
Xint init_done = 0;
Xint dateflag = 0;
Xint linetype = 7;		/* default for solid linetype */
Xfloat text_angle = 0.0;		/* radians */
Xfloat text_size = 1.0;		/* text size factor */
Xfloat text_shape = 1.0;		/* text shape factor */
Xint cset = 0;			/* char set, default is US Ascii */
Xint continuous = 1;		/* continuous line / discrete points */
X				/* 0 = discontinous */
Xchar marker = '*';		/* default marker char in discontinous mode */
Xchar informat[200] = "%f %f";	/* default input data line format */
Xint reversed = 0;		/* != 0 means input is y before x  */
Xint absolute = 0;		/* != 0 when 'absolute' coordinates are */
X				/* to be processed */
Xfloat xabs_origin = 0.0;	/* origin of 'absolute' mode data points */
Xfloat yabs_origin = 0.0;
X
Xint device = HP7475 | SIZE_A4;
X
Xstatic int setmode();
Xstatic int labelpic();
Xstatic void doit();
Xstatic void setdevice();
Xextern void exit();
X
Xstatic void (*old_int_handler)();
Xstatic void new_int_handler();
Xstatic void fatal();
Xstatic void usage();
X
X#define PUT_TEXT(x,y,xofs,yofs,text) \
X      dr_put_text((x), (y), (xofs), (yofs), text_angle, 1.0, (text))
X#define PUT_X_TEXT(x,y,text) \
X      dr_put_text((x), (y), -(float)strlen(text), -2.0, 0.0, 1.4, (text))
X#define PUT_Y_TEXT(x,y,text) \
X      dr_put_text((x), (y), 0.0, -2.0, -M_PI_2, 1.4, (text))
X#define PUT_N_TEXT(x,y,text) \
X      dr_put_text((x), (y), -(float)strlen(text)/2.0, 0.4, 0.0, 2.0, (text))
X
Xmain(argc, argv)
Xint argc;
Xchar *argv[];
X{
X    int c;
X    extern char *optarg;
X    extern int optind, opterr;
X    char *opts = "a:bc:df:g:hl:m:n:p:s:S:t:x:X:y:Y:v";
X
X    FILE *fp;
X
X    myname = *argv;
X    setdevice(getenv("DRAWDEV"));
X
X
X    while ((c = getopt(argc, argv, opts)) != -1) {
X	switch (c) {
X	case 'a':
X	    if (sscanf(optarg, "%f,%f,%f,%f",
X		       &xlowl, &ylowl, &xhighr, &yhighr) != 4)
X	      fatal("invalid plotting area definition: ", optarg);
X	    break;
X	case 'b': no_axes++; break;
X	case 'c': cset = atoi(optarg); break;
X	case 'd': dateflag++; break;
X	case 'f':
X	    /* input lines are : "-- -- y -- -- x -- --"
X	     * (dashes indicate fields, that we are not interested in)
X	     * and we want to plot y=f(x), so we specify option
X	     * "-f 6,3" (fields are numbered from 1...)
X	     */
X	    {
X		int xcol, ycol;
X		if (sscanf(optarg, "%d,%d", &xcol, &ycol) != 2)
X		  fatal("invalid input format definition: %s", optarg);
X		if (xcol == ycol)
X		  fatal("invalid input format definition: %s", optarg);
X		if (xcol > ycol) {
X		    int temp = xcol;
X		    reversed++;
X		    xcol = ycol;
X		    ycol = temp;
X		}
X		ycol -= xcol;
X		informat[0] = '\0';
X		while (--xcol)
X		  (void)strcat(informat, "%*f ");
X		(void)strcat(informat, "%f ");
X		while (--ycol)
X		  (void)strcat(informat, "%*f ");
X		(void)strcat(informat, "%f");
X	    }
X	    break;
X	case 'g':
X	    if (sscanf(optarg, "%d,%d", &xgrid_n, &ygrid_n) != 2)
X	      fatal("invalid grid spec.: ", optarg);
X	    break;
X	case 'h': case '?':
X	    usage(); break;
X	case 'l':
X	    if (sscanf(optarg, "%d", &linetype) != 1)
X	      fatal("invalid line type:", optarg);
X	    break;		/* draw axes before setting it! */
X	case 'm':
X	    continuous = 0;
X	    if (sscanf(optarg, "%1s", &marker) != 1)
X	      fatal ("no marker found.", "\0");
X	    break;
X	case 'n':
X	    name_text = optarg; break;
X	case 'p':
X	    if (sscanf(optarg, "%d", &pennum) != 1)
X	      fatal("invalid pen:", optarg);
X	    if ((pennum > 6) || (pennum < 1))
X	      pennum = 1;
X	    break;
X	case 's':
X	    if (sscanf(optarg, "%f", &text_size) != 1)
X	      fatal("invalid text size: ", optarg);
X	    break;
X	case 'S':
X	    if (sscanf(optarg, "%f", &text_shape) != 1)
X	      fatal("invalid text shape: ", optarg);
X	    break;
X	case 't': setdevice(optarg); break;
X	case 'x':
X	    if (sscanf(optarg, "%f,%f,%f",
X		       &xlow, &xincr, &xhigh) != 3)
X	      fatal("invalid x range: ", optarg);
X	    break;
X	case 'X':
X	    xlabel_text = optarg; break;
X	case 'y':
X	    if (sscanf(optarg, "%f,%f,%f",
X		       &ylow, &yincr, &yhigh) != 3)
X	      fatal("invalid y range: ", optarg);
X	    break;
X	case 'Y':
X	    ylabel_text = optarg; break;
X	case 'v':
X	    fprintf(stderr, "%s\n", RCS_id);
X	    exit (0);
X	default:
X	    exit(1);
X	}
X    }
X    argc -= optind;
X    argv += optind;
X
X    old_int_handler = signal (SIGINT, new_int_handler);
X
X    /*
X     * Options have been taken, start the work.  Initializations will be
X     * done just before first data point is used...  This is jippo to
X     * make it possible to give coordinate specifications from inside
X     * the data....
X     */
X    if (argc == 0)		/* no files, take points from stdin */
X      doit(stdin);
X    else {			/* take points from all files */
X	for ( ; argc; argc--, argv++) {
X	    if (NULL == (fp = fopen(*argv, "r"))) {
X		perror(*argv);
X		continue;
X	    }
X	    doit(fp);
X	    fclose(fp);
X	}
X    }
X    if (init_done)
X      return(dr_finish(1) ? 1 : 0); /* 1 == error(s) */
X    else
X      return 2;			/* nothing was done! */
X}
X
X
X/* gl_plot will set int handler at its initialization time, and restore
X   what it saw there previously, and finally it will
X   kill(SIGINT,getpid()).  This is what we are going to do, too, so that
X   OUR caller can clean ITS data stuctures etc.
X*/
Xstatic void new_int_handler() {
X    (void)dr_finish(0);
X    signal (SIGINT, old_int_handler);
X    kill (getpid(), SIGINT);	/* zap my self */
X    sleep(1);
X    exit(3);			/* if that is not enough */
X}
X
X
Xstatic void do_initializations()
X{
X    init_done++;		/* remember it is done */
X    dr_start(xlowl, ylowl, xhighr, yhighr, device);
X    if (cset != 0) dr_set_cset(cset);
X    dr_set_size(text_size);
X    dr_set_shape(text_shape);
X    dr_set_pen(pennum);
X    dr_area_def(xlow, xincr, xhigh, ylow, yincr, yhigh, no_axes);
X    if (!no_axes) {
X	if (xgrid_n) dr_xgrid(xgrid_n);
X	if (ygrid_n) dr_ygrid(ygrid_n);
X    }
X    if (xlabel_text) PUT_X_TEXT(xhigh, ylow, xlabel_text);
X    if (ylabel_text) PUT_Y_TEXT(xlow, yhigh, ylabel_text);
X    if (name_text) PUT_N_TEXT((xhigh+xlow)/2.0, yhigh, name_text);
X    if (dateflag) putdate(dateflag);
X    dr_set_ltype(linetype);
X}
X
Xstatic void setdevice(s)
Xchar *s;
X{
X    if (s == NULL) return;
X    if      (!strcmp(s,"a4")) device = (HP7475 | SIZE_A4);
X    else if (!strcmp(s,"a3")) device = (HP7475 | SIZE_A3);
X    else if (!strcmp(s,"a"))  device = (HP7475 | SIZE_A);
X    else if (!strcmp(s,"b"))  device = (HP7475 | SIZE_B);
X#ifdef GL_LIB
X    else if (!strcmp(s,"ega")) device = SCREEN;
X#endif
X#ifdef PLOT
X    else if (!strcmp(s,"plot")) device = PLOTLIB;
X#endif
X#ifdef TEK
X    else if (!strcmp(s,"tek"))
X      fatal("Sorry,", "Tek4010 not yet.");
X#endif
X    else fatal("Unsupported plotting device", s);
X}
X  
Xstatic void doit(fp)
XFILE *fp;
X{
X#define LINESIZE 512
X    char buf[LINESIZE+1], *p;
X    char emsg[200];
X    int lineno = 0;
X    
X    while (1) {
X	if (NULL == fgets(buf, LINESIZE, fp)) /* EOF */
X	  return;
X
X	lineno++;
X	for (p=buf; *p; p++)	/* skip white space */
X	  if (strchr(" \t",*p) == NULL)
X	    break;
X	if (*p == '\n')
X	  continue;
X
X	if (*p == '#') {
X	    if (setmode(p+1)) {
X		sprintf(emsg, "Invalid pseudo comment line (%d): %s",
X			lineno, buf);
X		if (init_done)
X		  dr_errmsg(emsg);
X		else
X		  write(2, emsg, strlen(emsg));
X	    }
X	    continue;
X	}
X				/* now we expect data lines, */
X				/* initializations must be done by now */
X	if (!init_done) do_initializations();
X
X	if (sscanf(p, informat, &x, &y) != 2) {
X	    sprintf(emsg, "Invalid data line (%d): %s", lineno, buf);
X	    dr_errmsg(emsg);
X	    continue;
X	}
X	if (reversed) {
X	    float temp = x;
X	    x = y;
X	    y = temp;
X	}
X	if (absolute) {
X	    x = xlow + (xabs_origin + x) * (xhigh - xlow);
X	    y = ylow + (yabs_origin + y) * (yhigh - ylow);
X	}
X	/* this needs to be fixed: reversed directions not handled! */
X	if (x < xlow || x > xhigh || y < ylow || y > yhigh) {
X	    sprintf (emsg, "point out of range at line (%d): %s",
X		     lineno, buf);
X	    dr_warnmsg(emsg);
X	}
X	if (continuous) {
X	    dr_goto(x, y, penup);
X	    penup = 0;
X	}
X	else {
X	    dr_put_mark(x, y, marker);
X	}
X    }
X}
X
Xstatic void remember(to, from)	/* copy string to safe place, put its */
X				/* addr into `*to' */
Xchar **to, *from;
X{
X#ifndef MSDOS
X    extern char *malloc();
X#endif
X    *to = malloc(strlen(from)+1);
X    strcpy(*to, from);
X}
X
X				/* handle '#<something>' lines */
Xstatic int setmode(s)		/* return 0 on success, != when trouble */
Xchar *s;
X{
X    for (; *s; s++)		/* skip white space */
X      if (strchr(" \t",*s) == NULL)
X	break;
X    if (*s++ != '$')		/* only '$' causes action */
X      return 0;
X
X    for (; *s; s++)		/* skip white space */
X      if (strchr(" \t",*s) == NULL)
X	break;
X    
X    if (!strncmp (s, "absolute", 8)) {
X	absolute = 1;
X	if (sscanf(s+8, "%f,%f", &xabs_origin, &yabs_origin) == 2)
X	  return 0;
X	xabs_origin = 0.0;
X	yabs_origin = 0.0;
X	return 1;
X    }
X
X    if (!strncmp (s, "endabsolute", 11)) {
X	absolute = 0;
X	return 0;
X    }
X
X    if (!strncmp (s,"penup",5)) {
X	penup = 1;
X	return 0;
X    }
X
X
X    if (!strncmp(s, "include", 7)) {
X	FILE *fp;
X	char *p;
X	
X	for (s+=8; *s; s++)	/* skip whitespace */
X	  if (strchr(" \t",*s) == NULL)
X 	    break;
X	for (p=s; *p; p++)
X 	  if (strchr(" \t\n",*s) != NULL)
X 	    break;
X 	*(p-1) = '\0';		/* put terminator */
X 	
X 	if (NULL == (fp = fopen(s, "r"))) {
X 	    perror(p);
X 	    return 1;
X 	}
X 	doit(fp);		/* recursive call */
X 	fclose(fp);
X 	return 0;
X    }
X 
X    if (!strncmp (s, "cset", 4)) {
X	int csetnum = 0;
X	if (sscanf(s+4, "%d", &csetnum) != 1)
X	  return 1;
X	if (!init_done) do_initializations();
X	dr_set_cset(csetnum);
X	return 0;
X    }
X
X    if (!strncmp (s, "pen", 3)) {
X	if (sscanf(s+3, "%d", &pennum) != 1)
X	  return 1;
X	if (!init_done) do_initializations();
X	dr_set_pen(pennum);
X	return 0;
X    }
X
X    if (!strncmp (s, "ltype", 5)) {
X	int linetype = 1;
X	if (sscanf(s+5, "%d", &linetype) != 1)
X	  return 1;
X	if (!init_done) do_initializations();
X	dr_set_ltype(linetype);
X	return 0;
X    }
X
X    if (!strncmp(s, "mark", 4)) {
X	s += 4;
X	while (1) {
X	    marker = *s++;
X	    if (strchr(" \t", marker) == NULL)
X	      break;
X	}
X	continuous = (strchr (" \n\r\t", marker) != NULL);
X	return 0;
X    }
X
X    if (!strncmp(s, "date",   4)) {
X	if (dateflag)
X	  return 1;
X	dateflag = 1;
X	if (init_done)
X	  putdate(1);
X	return 0;
X    }
X    if (!strncmp(s, "time",   4)) {
X	if (dateflag)
X	  return 1;
X	dateflag = 2;
X	if (init_done)
X	  putdate(2);
X	return 0;
X    }
X
X    if (!strncmp(s, "xlabel", 6)) {
X	if (xlabel_text) return 1;
X	s += 7;			/* skip one space/separator */
X	s[strlen(s)-1] = '\0';	/* drop newline */
X	if (!init_done) {
X	    remember(&xlabel_text, s);
X	    return 0;
X	}
X	PUT_X_TEXT(xhigh, ylow, s);
X	xlabel_text = myname;	/* some char-pointer that has non NULL */
X				/* value, so that we remember this */
X	penup = 1;
X	return 0;
X    }
X
X    if (!strncmp(s, "ylabel", 6)) {
X	if (ylabel_text) return 1;
X	s += 7;
X	s[strlen(s)-1] = '\0';
X	if (!init_done) {
X	    remember(&ylabel_text, s);
X	    return 0;
X	}
X	PUT_Y_TEXT(xlow, yhigh, s);
X	ylabel_text = myname;
X	penup = 1;
X	return 0;
X    }
X
X    if (!strncmp(s, "name", 4)) {
X	if (name_text) return 1;
X	s += 5;			/* skip the space */
X	s[strlen(s)-1] = '\0';	/* drop newline */
X	if (!init_done) {
X	    remember(&name_text, s);
X	    return 0;
X	}
X	PUT_N_TEXT((xhigh+xlow)/2.0, yhigh, s);
X	name_text = myname;	
X	penup = 1;
X	return 0;
X    }
X
X    if (!strncmp(s, "xrange", 6)) {
X	if (sscanf(s+6, "%f,%f,%f", &xlow, &xincr, &xhigh) != 3)
X	  return 1;
X	return 0;
X    }
X
X    if (!strncmp(s, "yrange", 6)) {
X	if (sscanf(s+6, "%f,%f,%f", &ylow, &yincr, &yhigh) != 3)
X	  return 1;
X	return 0;
X    }
X
X    if (!strncmp(s, "area", 4)) {
X	if (sscanf(s+4, "%f,%f,%f,%f",
X		   &xlowl, &ylowl, &xhighr, &yhighr) != 4)
X	  return 1;
X	return 0;
X    }
X
X    if (!strncmp(s, "grid", 4)) {
X	if (sscanf(s+4, "%d,%d", &xgrid_n, &ygrid_n) != 2)
X	  return 1;
X	return 0;
X    }
X
X    if (!strncmp(s, "noaxes", 6)) {
X	no_axes++;
X	return 0;
X    }
X
X    if (!strncmp(s, "init", 4)) {
X	if (!init_done) do_initializations();
X	return 0;
X    }
X
X    if (!strncmp(s, "size", 4)) {
X	if (sscanf(s+4, "%f", &text_size) != 1)
X	  return 1;
X	if (init_done)
X	  dr_set_size(text_size);
X	return 0;
X    }
X
X    if (!strncmp(s, "shape", 5)) {
X	if (sscanf(s+5, "%f", &text_shape) != 1)
X	  return 1;
X	if (init_done)
X	  dr_set_shape(text_shape);
X	return 0;
X    }
X
X    if (!strncmp(s, "direction", 9)) {
X	float xx;
X
X	if (sscanf(s+9, "%f", &xx) != 1)
X	  return 1;
X	text_angle = xx;
X	return 0;
X    }
X
X    if (!strncmp(s, "text",   4)) return (labelpic(s+4, 0, 0));
X    if (!strncmp(s, "alabel", 6)) return (labelpic(s+6, 0, 1));
X    if (!strncmp(s, "rlabel", 6)) return (labelpic(s+6, 1, 1));
X
X				/* other options are errors */
X    return 1;
X}
X
Xstatic int labelpic(s, rel, conn)
Xchar *s;
Xint rel, conn;
X{
X    float xx, yy;
X    char *s1;
X    int old_linetype;
X
X    if (!init_done) do_initializations();
X
X    old_linetype = dr_set_ltype(7); /* set solid line */
X    xx = strtod(s, &s1);
X    if (s1 == s) return 1;
X    s = s1+1;
X    yy = strtod(s, &s1);
X    if (s1 == s) return 1;
X    s1++;			/* skip the separator */
X    s1[strlen(s1)-1] = '\0';
X    if (rel) {
X	xx += x;
X	yy += y;
X    }
X    PUT_TEXT(xx, yy, 0.0, 0.15, s1);
X    if (conn) {
X	dr_goto(xx, yy, 0);
X	dr_goto(x, y, 0);
X	dr_draw_circle();
X    }
X    else
X      penup = 1;
X    (void) dr_set_ltype(old_linetype);
X    return 0;
X}
X
Xint putdate(flag)
Xint flag;
X{
X    struct tm *timp;
X    long clk;
X    char buf1[40], buf2[20];
X
X    if (time(&clk) < 0) {
X	(void) sprintf (buf1, "time() failed.");
X	(void) dr_errmsg(buf1);
X	return 1;
X    }
X    timp = localtime(&clk);
X    sprintf(buf1, "%d/%d/%d",
X	    timp->tm_mday, timp->tm_mon+1, timp->tm_year + 1900);
X    if (flag > 1)
X      sprintf(buf2, " %02d:%02d:%02d",
X	      timp->tm_hour, timp->tm_min, timp->tm_sec);
X    else
X      buf2[0] = '\0';
X    strcat(buf1, buf2);
X    PUT_TEXT(xhigh, yhigh, -(float)strlen(buf1), 0.7, buf1);
X    penup = 1;
X    return 0;
X}
X
X
Xvoid fatal(s1, s2)
Xchar *s1, *s2;
X{
X    fprintf (stderr,"%s: Fatal error: %s %s\n", myname, s1, s2);
X    exit(1);
X}
X
Xstatic char *utbl[] = {
X    "-a xlo,ylo,xhi,yhi\tspecify relative plot area on paper",
X    "-b\t\t\tdon't draw axes",
X    "-c num\t\t\tchoose character set `num'",
X    "-d\t\t\tinclude date (time) stamp",
X    "-f c1,c2\t\tuse these columns of input, higher is for y",
X    "-g xgrid,ygrid\t\tspecify x&y grids",
X    "-h\t\t\tgive this stuff",
X    "-l\t\t\tlinetype",
X    "-m marker\t\tspecify marker & discontinous mode",
X    "-n text\t\tname the picture",
X    "-p num\t\t\tspecify pen number",
X#ifdef GL_LIB
X    "-t target\t\tspecify device: `ega', `a3', `a4' etc. Defaults from env.",
X#else
X    "-t target\t\tspecify device: `a3', `a4' etc. Defaults from env.",
X#endif
X    "-x|y low,incr,high\tspecify axis range",
X    "-X|Y text\t\taxis label",
X    "-v\t\t\tprint version ID",
X    "Most (all?) of these may also be given inside the data, look for manual",
X    NULL
X};
X
Xvoid usage()
X{
X    register char **p;
X    printf ("usage: \"%s [opts] [file ...]\"\n", myname);
X    for (p=utbl; *p; p++)
X      (void)printf (" %s\n", *p);
X    exit(0);
X}
END_OF_FILE
if test 17080 -ne `wc -c <'draw.c'`; then
    echo shar: \"'draw.c'\" unpacked with wrong size!
fi
# end of 'draw.c'
fi
echo shar: End of archive 2 \(of 2\).
cp /dev/null ark2isdone
MISSING=""
for I in 1 2 ; do
    if test ! -f ark${I}isdone ; then
	MISSING="${MISSING} ${I}"
    fi
done
if test "${MISSING}" = "" ; then
    echo You have unpacked both archives.
    rm -f ark[1-9]isdone
else
    echo You still need to unpack the following archives:
    echo "        " ${MISSING}
fi
##  End of shell archive.
exit 0



More information about the Comp.sources.misc mailing list