src/meta/mplot.c

#ifndef lint
static const char       RCSid[] = "$Id$";
#endif
/*
 *   Plotting routines for meta-files to line-at-a-time printers
 */


#include  <fcntl.h>

#include  "meta.h"

#include  "plot.h"

#include  "span.h"


int  minwidth = 0;

static PLIST  inqueue = {NULL, NULL};

static PRIMITIVE  nextp;


plot(infp)              /* plot meta-file */

FILE  *infp;

{

    do {
        readp(&nextp, infp);
        initplot();
        while (inqueue.ptop != NULL || isprim(nextp.com))
            plotspan(infp);
        doglobal(&nextp);
        fargs(&nextp);
    } while (nextp.com != PEOF);

}


initplot()                      /* initialize this plot */

{

    thispage();
    outspan.xleft = 0;
    outspan.xright = dxsize - 1;
    outspan.ytop = dysize + linhite - 1;
    outspan.ybot = dysize;

}


doglobal(g)                     /* execute a global command */

PRIMITIVE  *g;

{
    char  c;
    int  tty;

    switch (g->com) {

        case PEOF:
            break;

        case PDRAW:
            fflush(stdout);
            break;

        case PEOP:
            if (g->arg0 & 0200)         /* advance page */
                nextpage();
            else if (g->arg0 == 3)      /* continue down */
                contpage();
            else
                error(USER, "illegal continue direction in doglobal");
            break;

        case PPAUSE:
            fflush(stdout);
            tty = open(TTY, O_RDWR);
            if (g->args != NULL) {
                write(tty, g->args, strlen(g->args));
                write(tty, " - (hit return to continue)", 27);
            } else
                write(tty, "\007", 1);
            do {
                c = '\n';
                read(tty, &c, 1);
            } while (c != '\n');
            close(tty);
            break;

        case PSET:
            set(g->arg0, g->args);
            break;

        case PUNSET:
            unset(g->arg0);
            break;

        case PRESET:
            reset(g->arg0);
            break;

        default:
            sprintf(errmsg, "unknown command '%c' in doglobal", g->com);
            error(WARNING, errmsg);
            break;
        }

}


plotspan(infp)                  /* plot next span */

FILE  *infp;

{
    PLIST  lastinq;
    register PRIMITIVE  *p;

                                        /* clear span */
    nextspan();
                                        /* plot from queue */
    lastinq.ptop = inqueue.ptop;
    lastinq.pbot = inqueue.pbot;
    inqueue.ptop = inqueue.pbot = NULL;
    while ((p = pop(&lastinq)) != NULL) {
        doprim(p);
        pfree(p);
    }
                                        /* plot from file */
    while (isprim(nextp.com) && CONV(nextp.xy[YMX],dysize) >= outspan.ybot) {
        doprim(&nextp);
        fargs(&nextp);
        readp(&nextp, infp);
    }
                                        /* print out span */
    outputspan();

}


nextspan()              /* prepare next span */

{
    register int  i;
    register char  *colp, *tcolp;
        
    if (spanmin <= spanmax) {                   /* clear span */

        i = nrows*dxsize;
        colp = outspan.cols;
        tcolp = outspan.tcols;
        while (i--)
            *colp++ = *tcolp++ = '\0';
    }
            
    outspan.ytop -= linhite;                    /* advance to next */
    outspan.ybot -= linhite;
    spanmin = dxsize;
    spanmax = 0;

}


outputspan()            /* output span to printer */
{
    register int  i;
    register char  *colp, *tcolp;

    if (spanmin <= spanmax) {                   /* overlay spans */

        i = nrows*dxsize;
        colp = outspan.cols;
        tcolp = outspan.tcols;
        while (i--)
            *colp++ |= *tcolp++;
    }
    printspan();                        /* print span */
}


doprim(p)               /* plot primitive */

register PRIMITIVE  *p;

{
    register PRIMITIVE  *newp;
    
    switch (p->com) {

        case PLSEG:
            plotlseg(p);
            break;

        case PRFILL:
            fill((p->arg0&0103) | (pati[(p->arg0>>2)&03]<<2),
                        CONV(p->xy[XMN],dxsize),CONV(p->xy[YMN],dysize),
                        CONV(p->xy[XMX],dxsize)+(p->arg0&0100?-1:0),
                        CONV(p->xy[YMX],dysize)+(p->arg0&0100?-1:0));
            break;

        case PTFILL:
            tfill(p);
            break;
                
        case PMSTR:
            printstr(p);
            break;

        default:
            sprintf(errmsg, "unknown command '%c' in doprim", p->com);
            error(WARNING, errmsg);
            return;
    }

    if (CONV(p->xy[YMN],dysize) < outspan.ybot) {       /* save for next time */
        if ((newp = palloc()) == NULL)
            error(SYSTEM, "memory limit exceeded in doprim");
        mcopy((char *)newp, (char *)p, sizeof(PRIMITIVE));
        newp->args = savestr(p->args);
        add(newp, &inqueue);
    }
        
}


plotlseg(p)             /* plot a line segment */

register PRIMITIVE  *p;

{
    register int  ti;
    long  run2 = 0L, rise2 = 0L;
    int  x, y, run, rise, xstop, ystop, hrad, vrad, lpat, n;

                                                /* compute line pattern */
    lpat = (p->arg0 >> 2) & 014;
    if (p->arg0 & 0100 && lpat != 0)
        lpat += 014;
    lpat |= p->arg0 & 03;

    ti = (p->arg0 >> 2) & 03;                   /* compute line radius */
    ti = WIDTH(ti) / 2;
    hrad = CONV(ti, dxsize);
    vrad = CONV(ti, dysize);
    if (hrad < minwidth)
        hrad = minwidth;
    if (vrad < minwidth)
        vrad = minwidth;

    x = CONV(p->xy[XMX], dxsize);               /* start at top */
    y = CONV(p->xy[YMX], dysize);
    run = CONV(p->xy[XMN], dxsize) - x;
    rise = CONV(p->xy[YMN], dysize) - y;

    if (p->arg0 & 0100)                         /* slope < 0; reverse x */
        x -= (run = -run);

    xstop = x + run;                            /* compute end point */
    ystop = y + rise;

    if ((ti = outspan.ytop+vrad+1-y) < 0) {     /* adjust to top of span */
        run2 = rise2 = (long)ti*run;
        x += rise2/rise;
        y += ti;
    }

    if ((ti = outspan.ybot-vrad-1-ystop) > 0) { /* adjust to bottom of span */
        xstop += (long)ti*run/rise;
        ystop += ti;
    }

    if (abs(run) > -rise)
        n = abs(xstop - x);
    else
        n = y - ystop;
    
    paintline(x, y, run, rise, hrad, vrad, lpat, run2, rise2, n);

}


/*
 *  This routine paints a line with calls to fill().  The line can
 *    start and end at arbitrary points on a longer line segment.
 */

paintline(x, y, run, rise, hrad, vrad, lpat, run2, rise2, n)

register int  x, y;
int  run, rise;
int  hrad, vrad;
int  lpat;
long  run2, rise2;
int  n;

{
    int  xstep, ystep;

    if (run >= 0)
        xstep = 1;
    else {
        xstep = -1;
        run = -run;
    }
    if (rise >= 0)
        ystep = 1;
    else {
        ystep = -1;
        rise = -rise;
    }
    if (n < 0)
        n = max(run, rise);

    if (run > rise)
        while (n >= 0)
            if (run2 >= rise2) {
                fill(lpat, x-hrad, y-vrad, x+hrad, y+vrad);
                n--;
                x += xstep;
                rise2 += rise;
            } else {
                y += ystep;
                run2 += run;
            }
    else
        while (n >= 0)
            if (rise2 >= run2) {
                fill(lpat, x-hrad, y-vrad, x+hrad, y+vrad);
                n--;
                y += ystep;
                run2 += run;
            } else {
                x += xstep;
                rise2 += rise;
            }

}


tfill(p)                        /* fill a triangle */
register PRIMITIVE  *p;
{
    register int  x, txmin, txmax;      
    int  xmn, ymn, tpat;
    long  xsz, ysz;

    xmn = CONV(p->xy[XMN], dxsize);
    xsz = CONV(p->xy[XMX], dxsize) - xmn;
    ymn = CONV(p->xy[YMN], dysize);
    ysz = CONV(p->xy[YMX], dysize) - ymn;
    if (xsz <= 0 || ysz <= 0)
        return; 
    txmin = (outspan.ybot - ymn)*xsz/ysz;
    txmax = (outspan.ytop - ymn)*xsz/ysz;
    if (p->arg0 & 020) {                                /* up or down */
        x = txmin;
        txmin = xsz - txmax;
        txmax = xsz - x;
    }
    txmin += xmn;
    txmax += xmn;
    txmin = max(txmin, xmn);
    txmax = min(txmax, xmn + (int)xsz - 1);
    tpat = (p->arg0&0103) | (pati[(p->arg0>>2)&03]<<2);

    if (p->arg0 & 040) {                                /* left or down */
        fill(tpat, xmn, ymn, txmin - 1, ymn + (int)ysz - 1);
        for (x = txmin; x <= txmax; x++)
            if (p->arg0 & 020)                  /* down */
                fill(tpat, x, ymn, x, (int)(ysz-(x-xmn)*ysz/xsz) + ymn - 1);
            else                                /* left */
                fill(tpat, x, (int)((x-xmn)*ysz/xsz) + ymn, x, ymn + (int)ysz - 1);
    } else {                                            /* right or up */
        for (x = txmin; x <= txmax; x++)
            if (p->arg0 & 020)                  /* up */
                fill(tpat, x, (int)(ysz-(x-xmn)*ysz/xsz) + ymn, x, ymn + (int)ysz - 1);
            else                                /* right */
                fill(tpat, x, ymn, x, (int)((x-xmn)*ysz/xsz) + ymn - 1);
        fill(tpat, txmax + 1, ymn, xmn + (int)xsz - 1, ymn + (int)ysz - 1);
    }
}


fill(attrib, xmin, ymin, xmax, ymax)    /* fill rectangle with attribute */

int  attrib;
int  xmin, ymin, xmax, ymax;

{
    int  filpat;
    int  rpos;
    unsigned char  *pattr;
    register char  *colp;
    register int  i;

    xmin = max(xmin, outspan.xleft) - outspan.xleft;
    ymin = max(ymin, outspan.ybot) - outspan.ybot;
    xmax = min(xmax, outspan.xright) - outspan.xleft;
    ymax = min(ymax, outspan.ytop) - outspan.ybot;

    for (rpos = 0; rpos < nrows; rpos++)

        if (rpos >= ymin >> 3 && rpos <= ymax >> 3) {

            filpat = 0377;
            if (rpos == ymin >> 3) {
                i = ymin & 07;
                filpat = (filpat >> i) << i;
            }
            if (rpos == ymax >> 3) {
                i = ~ymax & 07;
                filpat = ((filpat << i) & 0377) >> i;
            }

            pattr = pattern[(attrib&074)>>2]
                                [((outspan.ybot>>3)+rpos)%(PATSIZE>>3)];

            if (attrib & 0100) {
                colp = &outspan.tcols[rpos*dxsize + xmin];
                for (i = xmin; i <= xmax; i++)
                    *colp++ ^= filpat & pattr[i%PATSIZE];
            } else {
                colp = &outspan.cols[rpos*dxsize + xmin];
                for (i = xmin; i <= xmax; i++)
                    *colp++ |= filpat & pattr[i%PATSIZE];
            }

            spanmin = min(xmin, spanmin);
            spanmax = max(xmax, spanmax);
        }

}
Revision:	1.1
Committed:	Sat Feb 22 02:07:26 2003 UTC (21 years, 7 months ago) by greg
Content type:	text/plain
Branch:	MAIN
CVS Tags:	rad3R5
Log Message:	Changes and check-in for 3.5 release Includes new source files and modifications not recorded for many years See ray/doc/notes/ReleaseNotes for notes between 3.1 and 3.5 release
#	User	Rev	Content
1	greg	1.1	#ifndef lint
2			static const char RCSid[] = "$Id$";
3			#endif
4			/*
5			* Plotting routines for meta-files to line-at-a-time printers
6			*/
7
8
9			#include <fcntl.h>
10
11			#include "meta.h"
12
13			#include "plot.h"
14
15			#include "span.h"
16
17
18
19			int minwidth = 0;
20
21			static PLIST inqueue = {NULL, NULL};
22
23			static PRIMITIVE nextp;
24
25
26
27			plot(infp) /* plot meta-file */
28
29			FILE *infp;
30
31			{
32
33			do {
34			readp(&nextp, infp);
35			initplot();
36			while (inqueue.ptop != NULL \|\| isprim(nextp.com))
37			plotspan(infp);
38			doglobal(&nextp);
39			fargs(&nextp);
40			} while (nextp.com != PEOF);
41
42			}
43
44
45
46
47
48			initplot() /* initialize this plot */
49
50			{
51
52			thispage();
53			outspan.xleft = 0;
54			outspan.xright = dxsize - 1;
55			outspan.ytop = dysize + linhite - 1;
56			outspan.ybot = dysize;
57
58			}
59
60
61
62
63
64			doglobal(g) /* execute a global command */
65
66			PRIMITIVE *g;
67
68			{
69			char c;
70			int tty;
71
72			switch (g->com) {
73
74			case PEOF:
75			break;
76
77			case PDRAW:
78			fflush(stdout);
79			break;
80
81			case PEOP:
82			if (g->arg0 & 0200) /* advance page */
83			nextpage();
84			else if (g->arg0 == 3) /* continue down */
85			contpage();
86			else
87			error(USER, "illegal continue direction in doglobal");
88			break;
89
90			case PPAUSE:
91			fflush(stdout);
92			tty = open(TTY, O_RDWR);
93			if (g->args != NULL) {
94			write(tty, g->args, strlen(g->args));
95			write(tty, " - (hit return to continue)", 27);
96			} else
97			write(tty, "\007", 1);
98			do {
99			c = '\n';
100			read(tty, &c, 1);
101			} while (c != '\n');
102			close(tty);
103			break;
104
105			case PSET:
106			set(g->arg0, g->args);
107			break;
108
109			case PUNSET:
110			unset(g->arg0);
111			break;
112
113			case PRESET:
114			reset(g->arg0);
115			break;
116
117			default:
118			sprintf(errmsg, "unknown command '%c' in doglobal", g->com);
119			error(WARNING, errmsg);
120			break;
121			}
122
123			}
124
125
126
127
128
129			plotspan(infp) /* plot next span */
130
131			FILE *infp;
132
133			{
134			PLIST lastinq;
135			register PRIMITIVE *p;
136
137			/* clear span */
138			nextspan();
139			/* plot from queue */
140			lastinq.ptop = inqueue.ptop;
141			lastinq.pbot = inqueue.pbot;
142			inqueue.ptop = inqueue.pbot = NULL;
143			while ((p = pop(&lastinq)) != NULL) {
144			doprim(p);
145			pfree(p);
146			}
147			/* plot from file */
148			while (isprim(nextp.com) && CONV(nextp.xy[YMX],dysize) >= outspan.ybot) {
149			doprim(&nextp);
150			fargs(&nextp);
151			readp(&nextp, infp);
152			}
153			/* print out span */
154			outputspan();
155
156			}
157
158
159
160
161
162			nextspan() /* prepare next span */
163
164			{
165			register int i;
166			register char colp, tcolp;
167
168			if (spanmin <= spanmax) { /* clear span */
169
170			i = nrows*dxsize;
171			colp = outspan.cols;
172			tcolp = outspan.tcols;
173			while (i--)
174			colp++ = tcolp++ = '\0';
175			}
176
177			outspan.ytop -= linhite; /* advance to next */
178			outspan.ybot -= linhite;
179			spanmin = dxsize;
180			spanmax = 0;
181
182			}
183
184
185
186			outputspan() /* output span to printer */
187			{
188			register int i;
189			register char colp, tcolp;
190
191			if (spanmin <= spanmax) { /* overlay spans */
192
193			i = nrows*dxsize;
194			colp = outspan.cols;
195			tcolp = outspan.tcols;
196			while (i--)
197			colp++ \|= tcolp++;
198			}
199			printspan(); /* print span */
200			}
201
202
203
204			doprim(p) /* plot primitive */
205
206			register PRIMITIVE *p;
207
208			{
209			register PRIMITIVE *newp;
210
211			switch (p->com) {
212
213			case PLSEG:
214			plotlseg(p);
215			break;
216
217			case PRFILL:
218			fill((p->arg0&0103) \| (pati[(p->arg0>>2)&03]<<2),
219			CONV(p->xy[XMN],dxsize),CONV(p->xy[YMN],dysize),
220			CONV(p->xy[XMX],dxsize)+(p->arg0&0100?-1:0),
221			CONV(p->xy[YMX],dysize)+(p->arg0&0100?-1:0));
222			break;
223
224			case PTFILL:
225			tfill(p);
226			break;
227
228			case PMSTR:
229			printstr(p);
230			break;
231
232			default:
233			sprintf(errmsg, "unknown command '%c' in doprim", p->com);
234			error(WARNING, errmsg);
235			return;
236			}
237
238			if (CONV(p->xy[YMN],dysize) < outspan.ybot) { /* save for next time */
239			if ((newp = palloc()) == NULL)
240			error(SYSTEM, "memory limit exceeded in doprim");
241			mcopy((char )newp, (char )p, sizeof(PRIMITIVE));
242			newp->args = savestr(p->args);
243			add(newp, &inqueue);
244			}
245
246			}
247
248
249
250
251
252			plotlseg(p) /* plot a line segment */
253
254			register PRIMITIVE *p;
255
256			{
257			register int ti;
258			long run2 = 0L, rise2 = 0L;
259			int x, y, run, rise, xstop, ystop, hrad, vrad, lpat, n;
260
261			/* compute line pattern */
262			lpat = (p->arg0 >> 2) & 014;
263			if (p->arg0 & 0100 && lpat != 0)
264			lpat += 014;
265			lpat \|= p->arg0 & 03;
266
267			ti = (p->arg0 >> 2) & 03; /* compute line radius */
268			ti = WIDTH(ti) / 2;
269			hrad = CONV(ti, dxsize);
270			vrad = CONV(ti, dysize);
271			if (hrad < minwidth)
272			hrad = minwidth;
273			if (vrad < minwidth)
274			vrad = minwidth;
275
276			x = CONV(p->xy[XMX], dxsize); /* start at top */
277			y = CONV(p->xy[YMX], dysize);
278			run = CONV(p->xy[XMN], dxsize) - x;
279			rise = CONV(p->xy[YMN], dysize) - y;
280
281			if (p->arg0 & 0100) /* slope < 0; reverse x */
282			x -= (run = -run);
283
284			xstop = x + run; /* compute end point */
285			ystop = y + rise;
286
287			if ((ti = outspan.ytop+vrad+1-y) < 0) { /* adjust to top of span */
288			run2 = rise2 = (long)ti*run;
289			x += rise2/rise;
290			y += ti;
291			}
292
293			if ((ti = outspan.ybot-vrad-1-ystop) > 0) { /* adjust to bottom of span */
294			xstop += (long)ti*run/rise;
295			ystop += ti;
296			}
297
298			if (abs(run) > -rise)
299			n = abs(xstop - x);
300			else
301			n = y - ystop;
302
303			paintline(x, y, run, rise, hrad, vrad, lpat, run2, rise2, n);
304
305			}
306
307
308
309			/*
310			* This routine paints a line with calls to fill(). The line can
311			* start and end at arbitrary points on a longer line segment.
312			*/
313
314			paintline(x, y, run, rise, hrad, vrad, lpat, run2, rise2, n)
315
316			register int x, y;
317			int run, rise;
318			int hrad, vrad;
319			int lpat;
320			long run2, rise2;
321			int n;
322
323			{
324			int xstep, ystep;
325
326			if (run >= 0)
327			xstep = 1;
328			else {
329			xstep = -1;
330			run = -run;
331			}
332			if (rise >= 0)
333			ystep = 1;
334			else {
335			ystep = -1;
336			rise = -rise;
337			}
338			if (n < 0)
339			n = max(run, rise);
340
341			if (run > rise)
342			while (n >= 0)
343			if (run2 >= rise2) {
344			fill(lpat, x-hrad, y-vrad, x+hrad, y+vrad);
345			n--;
346			x += xstep;
347			rise2 += rise;
348			} else {
349			y += ystep;
350			run2 += run;
351			}
352			else
353			while (n >= 0)
354			if (rise2 >= run2) {
355			fill(lpat, x-hrad, y-vrad, x+hrad, y+vrad);
356			n--;
357			y += ystep;
358			run2 += run;
359			} else {
360			x += xstep;
361			rise2 += rise;
362			}
363
364			}
365
366
367
368			tfill(p) /* fill a triangle */
369			register PRIMITIVE *p;
370			{
371			register int x, txmin, txmax;
372			int xmn, ymn, tpat;
373			long xsz, ysz;
374
375			xmn = CONV(p->xy[XMN], dxsize);
376			xsz = CONV(p->xy[XMX], dxsize) - xmn;
377			ymn = CONV(p->xy[YMN], dysize);
378			ysz = CONV(p->xy[YMX], dysize) - ymn;
379			if (xsz <= 0 \|\| ysz <= 0)
380			return;
381			txmin = (outspan.ybot - ymn)*xsz/ysz;
382			txmax = (outspan.ytop - ymn)*xsz/ysz;
383			if (p->arg0 & 020) { /* up or down */
384			x = txmin;
385			txmin = xsz - txmax;
386			txmax = xsz - x;
387			}
388			txmin += xmn;
389			txmax += xmn;
390			txmin = max(txmin, xmn);
391			txmax = min(txmax, xmn + (int)xsz - 1);
392			tpat = (p->arg0&0103) \| (pati[(p->arg0>>2)&03]<<2);
393
394			if (p->arg0 & 040) { /* left or down */
395			fill(tpat, xmn, ymn, txmin - 1, ymn + (int)ysz - 1);
396			for (x = txmin; x <= txmax; x++)
397			if (p->arg0 & 020) /* down */
398			fill(tpat, x, ymn, x, (int)(ysz-(x-xmn)*ysz/xsz) + ymn - 1);
399			else /* left */
400			fill(tpat, x, (int)((x-xmn)*ysz/xsz) + ymn, x, ymn + (int)ysz - 1);
401			} else { /* right or up */
402			for (x = txmin; x <= txmax; x++)
403			if (p->arg0 & 020) /* up */
404			fill(tpat, x, (int)(ysz-(x-xmn)*ysz/xsz) + ymn, x, ymn + (int)ysz - 1);
405			else /* right */
406			fill(tpat, x, ymn, x, (int)((x-xmn)*ysz/xsz) + ymn - 1);
407			fill(tpat, txmax + 1, ymn, xmn + (int)xsz - 1, ymn + (int)ysz - 1);
408			}
409			}
410
411
412
413
414			fill(attrib, xmin, ymin, xmax, ymax) /* fill rectangle with attribute */
415
416			int attrib;
417			int xmin, ymin, xmax, ymax;
418
419			{
420			int filpat;
421			int rpos;
422			unsigned char *pattr;
423			register char *colp;
424			register int i;
425
426			xmin = max(xmin, outspan.xleft) - outspan.xleft;
427			ymin = max(ymin, outspan.ybot) - outspan.ybot;
428			xmax = min(xmax, outspan.xright) - outspan.xleft;
429			ymax = min(ymax, outspan.ytop) - outspan.ybot;
430
431			for (rpos = 0; rpos < nrows; rpos++)
432
433			if (rpos >= ymin >> 3 && rpos <= ymax >> 3) {
434
435			filpat = 0377;
436			if (rpos == ymin >> 3) {
437			i = ymin & 07;
438			filpat = (filpat >> i) << i;
439			}
440			if (rpos == ymax >> 3) {
441			i = ~ymax & 07;
442			filpat = ((filpat << i) & 0377) >> i;
443			}
444
445			pattr = pattern[(attrib&074)>>2]
446			[((outspan.ybot>>3)+rpos)%(PATSIZE>>3)];
447
448			if (attrib & 0100) {
449			colp = &outspan.tcols[rpos*dxsize + xmin];
450			for (i = xmin; i <= xmax; i++)
451			*colp++ ^= filpat & pattr[i%PATSIZE];
452			} else {
453			colp = &outspan.cols[rpos*dxsize + xmin];
454			for (i = xmin; i <= xmax; i++)
455			*colp++ \|= filpat & pattr[i%PATSIZE];
456			}
457
458			spanmin = min(xmin, spanmin);
459			spanmax = max(xmax, spanmax);
460			}
461
462			}