src/meta/tcurve.c

#ifndef lint
static const char       RCSid[] = "$Id: tcurve.c,v 1.4 2003/08/01 14:14:24 schorsch Exp $";
#endif
/*
 *       PROGRAM TO PLOT TEL-A-GRAF CURVES TO METAFILE
 *
 *              Greg Ward
 *              12/12/84
 *
 *      cc -o ../tcurve tcurve.c tgraph.o primout.o mfio.o syscalls.o misc.o -lm
 */

#include "rtprocess.h" /* getpid() */
#include  "tgraph.h"
#include  "plot.h"

#define  XLEGEND  (XBEG+XSIZ+4*TSIZ)    /* x start of legend */

#define  YLEGEND  (YBEG+2*YSIZ/3)       /* y start of legend */

short  usecurve[NCUR];          /* booleans for curve usage */

double  xmin, ymin, xmax, ymax;         /* domain */

double  xsize, ysize;                   /* axis dimensions */

double  xmnset = -FHUGE, xmxset = FHUGE,        /* domain settings */
        ymnset = -FHUGE, ymxset = FHUGE;

short  logx = FALSE, logy = FALSE;      /* flags for log plots */

short  polar = FALSE;                   /* flag for polar plot */
       
short  grid = FALSE;                    /* flag for grid */

int    curtype[NCUR] = {0, 020, 040, 060, 01, 021, 041, 061,
                         02, 022, 042, 062, 03, 023, 043, 063};

char   *sym[NCUR] = {"ex", "triangle", "square", "triangle2", "diamond",
                        "cross", "octagon", "crosssquare", "exsquare",
                        "trianglesquare", "triangle2square", "crossdiamond",
                        "crossoctagon", "exoctagon", "block", "bullet"};
                       
int    ncurves;

int    xlegend,
       ylegend;                         /* current legend position */

int    symrad = SYMRAD;                 /* symbol radius */ 

char  *progname;


static void limline(int a0, double x,double y,double xout,double yout, char *s);


int
main(
        int  argc,
        char  **argv
)
/*
 *     Take Tel-A-Graf runnable files and convert them to
 *  metafile primitives to send to standard output
 */

{
 char  tfname[MAXFNAME];
 FILE  *fp;
 int  axflag;

 progname = *argv++;
 argc--;

 initialize();

 for (; argc && (**argv == '-' || **argv == '+'); argc--, argv++)
    option(*argv);

 if (polar)                     /* avoid dumb choices */
    logx = FALSE;

 axflag = XTICS|XNUMS|YTICS|YNUMS;
 if (grid)
    axflag |= XGRID|YGRID;
 if (polar)
    axflag |= ORIGIN;
 else
    axflag |= BOX;

 pglob(PINCL, 2, "symbols.mta");

 if (argc)

    for ( ; argc--; argv++)  {

       fp = efopen(*argv, "r");
       normalize(fp, NULL);
       makeaxis(axflag);
       fseek(fp, 0L, 0);
       plot(fp);
       fclose(fp);
       }
 else  {

    sprintf(tfname, "%stc%d", TDIR, getpid());
    fp = efopen(tfname, "w+");
    normalize(stdin, fp);
    makeaxis(axflag);
    fseek(fp, 0L, 0);
    plot(fp);
    fclose(fp);
    unlink(tfname);
    }

 pglob(PEOF, 0200, NULL);

 return(0);
 }


void
plot(                   /* read file and generate plot */
        FILE  *fp
)
{
 int  ncur = 0;                 /* curves seen so far */
 int  cur = 0;                  /* current curve pattern */
 char  line[255], *s;
 double  x, y;
 double lastx = 0, lasty = 0;
 short  oobounds = FALSE, firstpoint = TRUE;

 xlegend = XLEGEND;
 ylegend = YLEGEND;

 if (ncurves > 0) {
    pprim(PMSTR, 0100, xlegend, ylegend+800, xlegend, ylegend+800, "Legend:");
    pprim(PMSTR, 0100, xlegend, ylegend+800, xlegend, ylegend+800, "______");
    }

 while (fgets(line, sizeof line, fp) != NULL)

    if (istitle(line)) {
       s = snagquo(line);
       boxstring(0, 0, YBEG+YSIZ+1000, XYSIZE-1, YBEG+YSIZ+1500, s);
       }

    else if (isxlabel(line)) {
       s = snagquo(line);
       boxstring(0, XBEG, YBEG-1250, XBEG+XSIZ, YBEG-900, s);
       }

    else if (isylabel(line)) {
       s = snagquo(line);
       boxstring(020, XBEG-1900, YBEG, XBEG-1550, YBEG+YSIZ, s);
       }

    else if (islabel(line)) {
       if (++ncur < NCUR && usecurve[ncur]) {
          oobounds = FALSE;
          firstpoint = TRUE;
          cur++;
          plseg(curtype[cur], xlegend, ylegend, xlegend+2000, ylegend);
          symout(curtype[cur] & 03, xlegend, ylegend, sym[ncur]);
          pprim(PMSTR, 020, xlegend+400, ylegend+200,
                      xlegend+400, ylegend+200, snagquo(line));
          ylegend -= 500;
          }
       }
       
    else if (isdata(line) && usecurve[ncur])  {

       if (getdata(line, &x, &y) >= 0) {

          if (oobounds) {
             oobounds = FALSE;
             limline(curtype[cur], x, y, lastx, lasty, sym[ncur]);
             }
          else if (firstpoint) {
             firstpoint = FALSE;
             limline(curtype[cur], x, y, x, y, sym[ncur]);
             }
          else
             plseg(curtype[cur], XCONV(lastx), YCONV(lasty),
                                 XCONV(x), YCONV(y));

          lastx = x;
          lasty = y;

          }
       else {
          if (!oobounds) {
             oobounds = TRUE;
             if (firstpoint)
                firstpoint = FALSE;
             else
                limline(curtype[cur], lastx, lasty, x, y, NULL);
             }
          lastx = x;
          lasty = y;
          }

       }

 pglob(PEOP, 0200, NULL);

 }


void
limline(        /* print line from/to out of bounds */
        int  a0,
        double  x,
        double  y,
        double  xout,
        double  yout,
        char  *s
)
{

    for ( ; ; )
        if (xout < xmin) {
            yout = (yout - y)*(xmin - x)/(xout - x) + y;
            xout = xmin;
        } else if (yout < ymin) {
            xout = (xout - x)*(ymin - y)/(yout - y) + x;
            yout = ymin;
        } else if (xout > xmax) {
            yout = (yout - y)*(xmax - x)/(xout - x) + y;
            xout = xmax;
        } else if (yout > ymax) {
            xout = (xout - x)*(ymax - y)/(yout - y) + x;
            yout = ymax;
        } else {
            if (s != NULL)
                symout(a0 & 03, XCONV(xout), YCONV(yout), s);
            plseg(a0, XCONV(x), YCONV(y), XCONV(xout), YCONV(yout));
            return;
        }        

}
Revision:	1.5
Committed:	Sat Nov 15 02:13:37 2003 UTC (20 years, 5 months ago) by schorsch
Content type:	text/plain
Branch:	MAIN
CVS Tags:	rad3R6P1, rad3R6
Changes since 1.4:	+22 -22 lines
Log Message:	Continued ANSIfication, and reduced other compile warnings.
#	User	Rev	Content
1	greg	1.1	#ifndef lint
2	schorsch	1.5	static const char RCSid[] = "$Id: tcurve.c,v 1.4 2003/08/01 14:14:24 schorsch Exp $";
3	greg	1.1	#endif
4			/*
5			* PROGRAM TO PLOT TEL-A-GRAF CURVES TO METAFILE
6			*
7			* Greg Ward
8			* 12/12/84
9			*
10			* cc -o ../tcurve tcurve.c tgraph.o primout.o mfio.o syscalls.o misc.o -lm
11			*/
12	schorsch	1.2
13	schorsch	1.3	#include "rtprocess.h" /* getpid() */
14	greg	1.1	#include "tgraph.h"
15	schorsch	1.5	#include "plot.h"
16	greg	1.1
17			#define XLEGEND (XBEG+XSIZ+4TSIZ) / x start of legend */
18
19			#define YLEGEND (YBEG+2YSIZ/3) / y start of legend */
20
21			short usecurve[NCUR]; /* booleans for curve usage */
22
23			double xmin, ymin, xmax, ymax; /* domain */
24
25			double xsize, ysize; /* axis dimensions */
26
27			double xmnset = -FHUGE, xmxset = FHUGE, /* domain settings */
28			ymnset = -FHUGE, ymxset = FHUGE;
29
30			short logx = FALSE, logy = FALSE; /* flags for log plots */
31
32			short polar = FALSE; /* flag for polar plot */
33
34			short grid = FALSE; /* flag for grid */
35
36			int curtype[NCUR] = {0, 020, 040, 060, 01, 021, 041, 061,
37			02, 022, 042, 062, 03, 023, 043, 063};
38
39			char *sym[NCUR] = {"ex", "triangle", "square", "triangle2", "diamond",
40			"cross", "octagon", "crosssquare", "exsquare",
41			"trianglesquare", "triangle2square", "crossdiamond",
42			"crossoctagon", "exoctagon", "block", "bullet"};
43
44			int ncurves;
45
46			int xlegend,
47			ylegend; /* current legend position */
48
49			int symrad = SYMRAD; /* symbol radius */
50
51			char *progname;
52
53
54	schorsch	1.5	static void limline(int a0, double x,double y,double xout,double yout, char *s);
55	greg	1.1
56
57	schorsch	1.5	int
58			main(
59			int argc,
60			char **argv
61			)
62	greg	1.1	/*
63			* Take Tel-A-Graf runnable files and convert them to
64			* metafile primitives to send to standard output
65			*/
66
67			{
68			char tfname[MAXFNAME];
69			FILE *fp;
70			int axflag;
71
72			progname = *argv++;
73			argc--;
74
75			initialize();
76
77			for (; argc && (argv == '-' \|\| argv == '+'); argc--, argv++)
78			option(*argv);
79
80			if (polar) /* avoid dumb choices */
81			logx = FALSE;
82
83			axflag = XTICS\|XNUMS\|YTICS\|YNUMS;
84			if (grid)
85			axflag \|= XGRID\|YGRID;
86			if (polar)
87			axflag \|= ORIGIN;
88			else
89			axflag \|= BOX;
90
91			pglob(PINCL, 2, "symbols.mta");
92
93			if (argc)
94
95			for ( ; argc--; argv++) {
96
97			fp = efopen(*argv, "r");
98			normalize(fp, NULL);
99			makeaxis(axflag);
100			fseek(fp, 0L, 0);
101			plot(fp);
102			fclose(fp);
103			}
104			else {
105
106			sprintf(tfname, "%stc%d", TDIR, getpid());
107			fp = efopen(tfname, "w+");
108			normalize(stdin, fp);
109			makeaxis(axflag);
110			fseek(fp, 0L, 0);
111			plot(fp);
112			fclose(fp);
113			unlink(tfname);
114			}
115
116			pglob(PEOF, 0200, NULL);
117
118			return(0);
119			}
120
121
122
123
124	schorsch	1.5	void
125			plot( /* read file and generate plot */
126			FILE *fp
127			)
128	greg	1.1	{
129			int ncur = 0; /* curves seen so far */
130			int cur = 0; /* current curve pattern */
131			char line[255], *s;
132			double x, y;
133	schorsch	1.5	double lastx = 0, lasty = 0;
134	greg	1.1	short oobounds = FALSE, firstpoint = TRUE;
135
136			xlegend = XLEGEND;
137			ylegend = YLEGEND;
138
139			if (ncurves > 0) {
140			pprim(PMSTR, 0100, xlegend, ylegend+800, xlegend, ylegend+800, "Legend:");
141			pprim(PMSTR, 0100, xlegend, ylegend+800, xlegend, ylegend+800, "______");
142			}
143
144			while (fgets(line, sizeof line, fp) != NULL)
145
146			if (istitle(line)) {
147			s = snagquo(line);
148			boxstring(0, 0, YBEG+YSIZ+1000, XYSIZE-1, YBEG+YSIZ+1500, s);
149			}
150
151			else if (isxlabel(line)) {
152			s = snagquo(line);
153			boxstring(0, XBEG, YBEG-1250, XBEG+XSIZ, YBEG-900, s);
154			}
155
156			else if (isylabel(line)) {
157			s = snagquo(line);
158			boxstring(020, XBEG-1900, YBEG, XBEG-1550, YBEG+YSIZ, s);
159			}
160
161			else if (islabel(line)) {
162			if (++ncur < NCUR && usecurve[ncur]) {
163			oobounds = FALSE;
164			firstpoint = TRUE;
165			cur++;
166			plseg(curtype[cur], xlegend, ylegend, xlegend+2000, ylegend);
167			symout(curtype[cur] & 03, xlegend, ylegend, sym[ncur]);
168			pprim(PMSTR, 020, xlegend+400, ylegend+200,
169			xlegend+400, ylegend+200, snagquo(line));
170			ylegend -= 500;
171			}
172			}
173
174			else if (isdata(line) && usecurve[ncur]) {
175
176			if (getdata(line, &x, &y) >= 0) {
177
178			if (oobounds) {
179			oobounds = FALSE;
180			limline(curtype[cur], x, y, lastx, lasty, sym[ncur]);
181			}
182			else if (firstpoint) {
183			firstpoint = FALSE;
184			limline(curtype[cur], x, y, x, y, sym[ncur]);
185			}
186			else
187			plseg(curtype[cur], XCONV(lastx), YCONV(lasty),
188			XCONV(x), YCONV(y));
189
190			lastx = x;
191			lasty = y;
192
193			}
194			else {
195			if (!oobounds) {
196			oobounds = TRUE;
197			if (firstpoint)
198			firstpoint = FALSE;
199			else
200			limline(curtype[cur], lastx, lasty, x, y, NULL);
201			}
202			lastx = x;
203			lasty = y;
204			}
205
206			}
207
208			pglob(PEOP, 0200, NULL);
209
210			}
211
212
213
214	schorsch	1.5	void
215			limline( /* print line from/to out of bounds */
216			int a0,
217			double x,
218			double y,
219			double xout,
220			double yout,
221			char *s
222			)
223	greg	1.1	{
224
225			for ( ; ; )
226			if (xout < xmin) {
227			yout = (yout - y)*(xmin - x)/(xout - x) + y;
228			xout = xmin;
229			} else if (yout < ymin) {
230			xout = (xout - x)*(ymin - y)/(yout - y) + x;
231			yout = ymin;
232			} else if (xout > xmax) {
233			yout = (yout - y)*(xmax - x)/(xout - x) + y;
234			xout = xmax;
235			} else if (yout > ymax) {
236			xout = (xout - x)*(ymax - y)/(yout - y) + x;
237			yout = ymax;
238			} else {
239			if (s != NULL)
240			symout(a0 & 03, XCONV(xout), YCONV(yout), s);
241			plseg(a0, XCONV(x), YCONV(y), XCONV(xout), YCONV(yout));
242			return;
243			}
244
245			}