src/px/psum.c

#ifndef lint
static const char       RCSid[] = "$Id$";
#endif
/*
 *  psum.c - program to sum pictures.
 *
 *     10/14/86
 */

#include  <stdio.h>

#include  <math.h>

#include  "color.h"


#define  MAXFILE        8

int  xsiz, ysiz;

char  *progname;

char  *fname[MAXFILE];                  /* the file names */
FILE  *fptr[MAXFILE];                   /* the file pointers */
COLOR  scale[MAXFILE];                  /* scaling factors */
int  nfile;                             /* number of files */


int
tabputs(s)                      /* print line preceded by a tab */
char  *s;
{
        putc('\t', stdout);
        return(fputs(s, stdout));
}


main(argc, argv)
int  argc;
char  *argv[];
{
        double  d;
        int  xres, yres;
        int  an;

        progname = argv[0];

        nfile = 0;
        setcolor(scale[0], 1.0, 1.0, 1.0);

        for (an = 1; an < argc; an++) {
                if (nfile >= MAXFILE) {
                        fprintf(stderr, "%s: too many files\n", progname);
                        quit(1);
                }
                if (argv[an][0] == '-')
                        switch (argv[an][1]) {
                        case 's':
                                d = atof(argv[an+1]);
                                switch (argv[an][2]) {
                                case '\0':
                                        scalecolor(scale[nfile], d);
                                        break;
                                case 'r':
                                        colval(scale[nfile],RED) *= d;
                                        break;
                                case 'g':
                                        colval(scale[nfile],GRN) *= d;
                                        break;
                                case 'b':
                                        colval(scale[nfile],BLU) *= d;
                                        break;
                                default:
                                        goto unkopt;
                                }
                                an++;
                                continue;
                        case '\0':
                                fptr[nfile] = stdin;
                                fname[nfile] = "<stdin>";
                                break;
                        default:;
                        unkopt:
                                fprintf(stderr, "%s: unknown option: %s\n",
                                                progname, argv[an]);
                                quit(1);
                        }
                else if ((fptr[nfile] = fopen(argv[an], "r")) == NULL) {
                        fprintf(stderr, "%s: can't open file: %s\n",
                                                progname, argv[an]);
                        quit(1);
                } else
                        fname[nfile] = argv[an];
                                                /* get header */
                fputs(fname[nfile], stdout);
                fputs(":\n", stdout);
                getheader(fptr[nfile], tabputs);
                                                /* get picture size */
                if (fgetresolu(&xres, &yres, fptr[nfile]) != (YMAJOR|YDECR)) {
                        fprintf(stderr, "%s: bad picture size\n", progname);
                        quit(1);
                } else if (nfile == 0) {
                        xsiz = xres;
                        ysiz = yres;
                } else if (xres != xsiz || yres != ysiz) {
                        fprintf(stderr, "%s: pictures are different sizes\n",
                                        progname);
                        quit(1);
                }
                nfile++;
                setcolor(scale[nfile], 1.0, 1.0, 1.0);
        }
                                        /* add new header info. */
        printargs(argc, argv, stdout);
        putchar('\n');
        fputresolu(YMAJOR|YDECR, xsiz, ysiz, stdout);

        psum();
        
        quit(0);
}


psum()                          /* sum the files */
{
        COLOR  *scanin, *scanout;
        int  y, i;
        register int  x;

        scanin = (COLOR *)malloc(xsiz*sizeof(COLOR));
        scanout = (COLOR *)malloc(xsiz*sizeof(COLOR));
        if (scanin == NULL || scanout == NULL) {
                fprintf(stderr, "%s: out of memory\n", progname);
                quit(1);
        }
        for (y = ysiz-1; y >= 0; y--) {
                for (x = 0; x < xsiz; x++)
                        setcolor(scanout[x], 0.0, 0.0, 0.0);
                for (i = 0; i < nfile; i++) {
                        if (freadscan(scanin, xsiz, fptr[i]) < 0) {
                                fprintf(stderr, "%s: read error on file: %s\n",
                                                progname, fname[i]);
                                quit(1);
                        }
                        for (x = 0; x < xsiz; x++)
                                multcolor(scanin[x], scale[i]);
                        for (x = 0; x < xsiz; x++)
                                addcolor(scanout[x], scanin[x]);
                }
                if (fwritescan(scanout, xsiz, stdout) < 0) {
                        fprintf(stderr, "%s: write error\n", progname);
                        quit(1);
                }
        }
        free((void *)scanin);
        free((void *)scanout);
}


void
quit(code)              /* exit gracefully */
int  code;
{
        exit(code);
}
Revision:	2.5
Committed:	Sat Feb 22 02:07:27 2003 UTC (21 years, 2 months ago) by greg
Content type:	text/plain
Branch:	MAIN
CVS Tags:	rad3R5
Changes since 2.4:	+4 -6 lines
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	greg	2.5	static const char RCSid[] = "$Id$";
3	greg	1.1	#endif
4			/*
5			* psum.c - program to sum pictures.
6			*
7			* 10/14/86
8			*/
9
10			#include <stdio.h>
11
12	greg	2.3	#include <math.h>
13
14	greg	1.1	#include "color.h"
15
16
17			#define MAXFILE 8
18
19			int xsiz, ysiz;
20
21			char *progname;
22
23			char fname[MAXFILE]; / the file names */
24			FILE fptr[MAXFILE]; / the file pointers */
25			COLOR scale[MAXFILE]; /* scaling factors */
26			int nfile; /* number of files */
27
28
29	gwlarson	2.4	int
30	greg	1.1	tabputs(s) /* print line preceded by a tab */
31			char *s;
32			{
33			putc('\t', stdout);
34	gwlarson	2.4	return(fputs(s, stdout));
35	greg	1.1	}
36
37
38			main(argc, argv)
39			int argc;
40			char *argv[];
41			{
42	greg	2.2	double d;
43	greg	1.1	int xres, yres;
44			int an;
45
46			progname = argv[0];
47
48			nfile = 0;
49			setcolor(scale[0], 1.0, 1.0, 1.0);
50
51			for (an = 1; an < argc; an++) {
52			if (nfile >= MAXFILE) {
53			fprintf(stderr, "%s: too many files\n", progname);
54			quit(1);
55			}
56			if (argv[an][0] == '-')
57			switch (argv[an][1]) {
58			case 's':
59			d = atof(argv[an+1]);
60			switch (argv[an][2]) {
61			case '\0':
62			scalecolor(scale[nfile], d);
63			break;
64			case 'r':
65			colval(scale[nfile],RED) *= d;
66			break;
67			case 'g':
68			colval(scale[nfile],GRN) *= d;
69			break;
70			case 'b':
71			colval(scale[nfile],BLU) *= d;
72			break;
73			default:
74			goto unkopt;
75			}
76			an++;
77			continue;
78			case '\0':
79			fptr[nfile] = stdin;
80			fname[nfile] = "<stdin>";
81			break;
82			default:;
83			unkopt:
84			fprintf(stderr, "%s: unknown option: %s\n",
85			progname, argv[an]);
86			quit(1);
87			}
88			else if ((fptr[nfile] = fopen(argv[an], "r")) == NULL) {
89			fprintf(stderr, "%s: can't open file: %s\n",
90			progname, argv[an]);
91			quit(1);
92			} else
93			fname[nfile] = argv[an];
94			/* get header */
95			fputs(fname[nfile], stdout);
96			fputs(":\n", stdout);
97			getheader(fptr[nfile], tabputs);
98			/* get picture size */
99	greg	1.2	if (fgetresolu(&xres, &yres, fptr[nfile]) != (YMAJOR\|YDECR)) {
100	greg	1.1	fprintf(stderr, "%s: bad picture size\n", progname);
101			quit(1);
102			} else if (nfile == 0) {
103			xsiz = xres;
104			ysiz = yres;
105			} else if (xres != xsiz \|\| yres != ysiz) {
106			fprintf(stderr, "%s: pictures are different sizes\n",
107			progname);
108			quit(1);
109			}
110			nfile++;
111			setcolor(scale[nfile], 1.0, 1.0, 1.0);
112			}
113			/* add new header info. */
114			printargs(argc, argv, stdout);
115	greg	1.2	putchar('\n');
116			fputresolu(YMAJOR\|YDECR, xsiz, ysiz, stdout);
117	greg	1.1
118			psum();
119
120			quit(0);
121			}
122
123
124			psum() /* sum the files */
125			{
126			COLOR scanin, scanout;
127			int y, i;
128			register int x;
129
130			scanin = (COLOR )malloc(xsizsizeof(COLOR));
131			scanout = (COLOR )malloc(xsizsizeof(COLOR));
132			if (scanin == NULL \|\| scanout == NULL) {
133			fprintf(stderr, "%s: out of memory\n", progname);
134			quit(1);
135			}
136			for (y = ysiz-1; y >= 0; y--) {
137			for (x = 0; x < xsiz; x++)
138			setcolor(scanout[x], 0.0, 0.0, 0.0);
139			for (i = 0; i < nfile; i++) {
140			if (freadscan(scanin, xsiz, fptr[i]) < 0) {
141			fprintf(stderr, "%s: read error on file: %s\n",
142			progname, fname[i]);
143			quit(1);
144			}
145			for (x = 0; x < xsiz; x++)
146			multcolor(scanin[x], scale[i]);
147			for (x = 0; x < xsiz; x++)
148			addcolor(scanout[x], scanin[x]);
149			}
150			if (fwritescan(scanout, xsiz, stdout) < 0) {
151			fprintf(stderr, "%s: write error\n", progname);
152			quit(1);
153			}
154			}
155	greg	2.5	free((void *)scanin);
156			free((void *)scanout);
157	greg	1.1	}
158
159
160	greg	2.5	void
161	greg	1.1	quit(code) /* exit gracefully */
162			int code;
163			{
164			exit(code);
165			}