src/cal/total.c

#ifndef lint
static const char       RCSid[] = "$Id: total.c,v 1.13 2016/08/18 00:52:47 greg Exp $";
#endif
/*
 *  total.c - program to reduce columns of data.
 *
 *      5/18/88
 */

#include  <stdlib.h>
#include  <ctype.h>
#include  <math.h>
#include  "platform.h"
#include  "rtio.h"

#define  MAXCOL         8192            /* maximum number of columns */

#define  ADD            0               /* add numbers */
#define  MULT           1               /* multiply numbers */
#define  MAX            2               /* maximum */
#define  MIN            3               /* minimum */

double  init_val[] = {0., 1., -1e12, 1e12};     /* initial values */

long  incnt = 0;                        /* limit number of input records? */
long  outcnt = 0;                       /* limit number of output records? */

int  func = ADD;                        /* default function */
double  power = 0.0;                    /* power for sum */
int  mean = 0;                          /* compute mean */
int  count = 0;                         /* number to sum */
int  nbicols = 0;                       /* number of binary input columns */
int  bocols = 0;                        /* produce binary output columns */
int  tabc = '\t';                       /* default separator */
int  subtotal = 0;                      /* produce subtotals? */
int  nrecsout = 0;                      /* number of records produced */

static int execute(char *fname);

int
main(
int  argc,
char  *argv[]
)
{
        int  status;
        int  a;

        for (a = 1; a < argc; a++)
                if (argv[a][0] == '-')
                        if (argv[a][1] >= '0' && argv[a][1] <= '9')
                                count = atoi(argv[a]+1);
                        else
                                switch (argv[a][1]) {
                                case 'm':
                                        mean = !mean;
                                        break;
                                case 'p':
                                        func = MULT;
                                        break;
                                case 's':
                                        func = ADD;
                                        power = atof(argv[a]+2);
                                        break;
                                case 'u':
                                        func = MAX;
                                        break;
                                case 'l':
                                        func = MIN;
                                        break;
                                case 't':
                                        tabc = argv[a][2];
                                        break;
                                case 'r':
                                        subtotal = !subtotal;
                                        break;
                                case 'i':
                                        switch (argv[a][2]) {
                                        case 'n':
                                                incnt = atol(argv[++a]);
                                                break;
                                        case 'a':
                                                nbicols = 0;
                                                break;
                                        case 'd':
                                                if (isdigit(argv[a][3]))
                                                        nbicols = atoi(argv[a]+3);
                                                else
                                                        nbicols = 1;
                                                if (nbicols > MAXCOL) {
                                                        fprintf(stderr,
                                                        "%s: too many input columns\n",
                                                                argv[0]);
                                                        exit(1);
                                                }
                                                break;
                                        case 'f':
                                                if (isdigit(argv[a][3]))
                                                        nbicols = -atoi(argv[a]+3);
                                                else
                                                        nbicols = -1;
                                                if (-nbicols > MAXCOL) {
                                                        fprintf(stderr,
                                                        "%s: too many input columns\n",
                                                                argv[0]);
                                                        exit(1);
                                                }
                                                break;
                                        default:
                                                goto userr;
                                        }
                                        break;
                                case 'o':
                                        switch (argv[a][2]) {
                                        case 'n':
                                                outcnt = atol(argv[++a]);
                                                break;
                                        case 'a':
                                                bocols = 0;
                                                break;
                                        case 'd':
                                                bocols = 1;
                                                break;
                                        case 'f':
                                                bocols = -1;
                                                break;
                                        default:
                                                goto userr;
                                        }
                                        break;
                                default:;
                                userr:
                                        fprintf(stderr,
"Usage: %s [-m][-sE|p|u|l][-tC][-i{f|d}[N]][-o{f|d}][-count [-r]] [file..]\n",
                                                        argv[0]);
                                        exit(1);
                                }
                else
                        break;
        if (mean) {
                if (func == MAX) {
                        fprintf(stderr, "%s: average maximum?!\n", argv[0]);
                        exit(1);
                }
                if (func == MIN) {
                        fprintf(stderr, "%s: average minimum?!\n", argv[0]);
                        exit(1);
                }
        }
        if (bocols)
                SET_FILE_BINARY(stdout);
        status = 0;
        if (a >= argc)
                status = execute(NULL) == -1 ? 1 : status;
        else
                for ( ; a < argc && (outcnt <= 0 || nrecsout < outcnt); a++)
                        status = execute(argv[a]) == -1 ? 2 : status;
        exit(status);
}


static int
getrecord(                      /* read next input record */
        double field[MAXCOL],
        FILE *fp
)
{
        char  buf[16*MAXCOL];
        char  *cp, *num;
        int   nf;
                                                /* reading binary input? */
        if (nbicols > 0)
                return(getbinary(field, sizeof(double), nbicols, fp));
        if (nbicols < 0) {
                float   *fbuf = (float *)buf;
                int     i;
                nf = getbinary(fbuf, sizeof(float), -nbicols, fp);
                for (i = nf; i-- > 0; )
                        field[i] = fbuf[i];
                return(nf);
        }
                                                /* reading ASCII input */
        cp = fgets(buf, sizeof(buf), fp);
        if (cp == NULL)
                return(0);
        for (nf = 0; nf < MAXCOL; nf++) {
                while (isspace(*cp))
                        cp++;
                if (!*cp)
                        break;
                num = cp;
                while (*cp && !(isspace(tabc)?isspace(*cp):*cp==tabc))
                        cp++;
                if (*cp)
                        *cp++ = '\0';
                if (!*num || isspace(*num))
                        field[nf] = init_val[func];
                else
                        field[nf] = atof(num);
        }
        return(nf);
}


static void
putrecord(                      /* write out results record */
        const double *field,
        int n,
        FILE *fp
)
{
                                                /* binary output? */
        if (bocols > 0) {
                putbinary(field, sizeof(double), n, fp);
                return;
        }
        if (bocols < 0) {
                float   fv;
                while (n-- > 0) {
                        fv = *field++;
                        putbinary(&fv, sizeof(float), 1, fp);
                }
                return;
        }
                                                /* ASCII output */
        while (n-- > 0) {
                fprintf(fp, "%.9g", *field++);
                if (n) fputc(tabc, fp);
        }
        fputc('\n', fp);
}


static int
execute(                        /* compute result */
char  *fname
)
{
        double  inpval[MAXCOL];
        double  tally[MAXCOL];
        short   rsign[MAXCOL];
        double  result[MAXCOL];
        int  n;
        int  nread, ncol;
        long  nlin, ltotal;
        FILE  *fp;
                                                        /* open file */
        if (fname == NULL)
                fp = stdin;
        else if ((fp = fopen(fname, "r")) == NULL) {
                fprintf(stderr, "%s: cannot open\n", fname);
                return(-1);
        }
        if (nbicols)
                SET_FILE_BINARY(fp);
#ifdef getc_unlocked                            /* avoid lock/unlock overhead */
        flockfile(fp);
#endif

        ltotal = 0;
        while (!feof(fp)) {
                if (ltotal == 0) {                      /* initialize */
                        if (func == MULT)       /* special case */
                                for (n = 0; n < MAXCOL; n++) {
                                        tally[n] = 0.0;
                                        rsign[n] = 1;
                                }
                        else
                                for (n = 0; n < MAXCOL; n++)
                                        tally[n] = init_val[func];
                }
                ncol = 0;
                for (nlin = 0; (count <= 0 || nlin < count) &&
                                (incnt <= 0 || nlin < incnt) &&
                                (nread = getrecord(inpval, fp)) > 0;
                                nlin++) {
                                                        /* compute */
                        for (n = 0; n < nread; n++)
                                switch (func) {
                                case ADD:
                                        if (inpval[n] == 0.0)
                                                break;
                                        if (power != 0.0)
                                                tally[n] += pow(fabs(inpval[n]),power);
                                        else
                                                tally[n] += inpval[n];
                                        break;
                                case MULT:
                                        if (inpval[n] == 0.0)
                                                rsign[n] = 0;
                                        else if (inpval[n] < 0.0) {
                                                rsign[n] = -rsign[n];
                                                inpval[n] = -inpval[n];
                                        }
                                        if (rsign[n])
                                                tally[n] += log(inpval[n]);
                                        break;
                                case MAX:
                                        if (inpval[n] > tally[n])
                                                tally[n] = inpval[n];
                                        break;
                                case MIN:
                                        if (inpval[n] < tally[n])
                                                tally[n] = inpval[n];
                                        break;
                                }
                        if (nread > ncol)
                                ncol = nread;
                }
                if (nlin == 0)
                        break;
                                                /* compute and print */
                ltotal += nlin;
                for (n = 0; n < ncol; n++) {
                        result[n] = tally[n];
                        if (mean) {
                                result[n] /= (double)ltotal;
                                if (func == ADD && power != 0.0 && result[n] != 0.0)
                                        result[n] = pow(result[n], 1.0/power);
                        }
                        if (func == MULT)
                                result[n] = rsign[n] * exp(result[n]);
                }
                putrecord(result, ncol, stdout);
                ++nrecsout;
                if (outcnt > 0 && nrecsout >= outcnt)
                        break;
                if (!subtotal)
                        ltotal = 0;
                if (incnt > 0 && nlin >= incnt)
                        break;
        }
                                                        /* close input */
        return(fclose(fp));
}
Revision:	1.14
Committed:	Sat Dec 28 18:05:13 2019 UTC (4 years, 4 months ago) by greg
Content type:	text/plain
Branch:	MAIN
CVS Tags:	rad5R3
Changes since 1.13:	+1 -2 lines
Log Message:	Removed redundant include files
#	User	Rev	Content
1	greg	1.1	#ifndef lint
2	greg	1.14	static const char RCSid[] = "$Id: total.c,v 1.13 2016/08/18 00:52:47 greg Exp $";
3	greg	1.1	#endif
4			/*
5			* total.c - program to reduce columns of data.
6			*
7			* 5/18/88
8			*/
9
10			#include <stdlib.h>
11			#include <ctype.h>
12			#include <math.h>
13	greg	1.4	#include "platform.h"
14	greg	1.13	#include "rtio.h"
15	greg	1.1
16	greg	1.8	#define MAXCOL 8192 /* maximum number of columns */
17	greg	1.1
18			#define ADD 0 /* add numbers */
19			#define MULT 1 /* multiply numbers */
20			#define MAX 2 /* maximum */
21			#define MIN 3 /* minimum */
22
23			double init_val[] = {0., 1., -1e12, 1e12}; /* initial values */
24
25	greg	1.11	long incnt = 0; /* limit number of input records? */
26			long outcnt = 0; /* limit number of output records? */
27
28	greg	1.1	int func = ADD; /* default function */
29			double power = 0.0; /* power for sum */
30			int mean = 0; /* compute mean */
31			int count = 0; /* number to sum */
32	greg	1.4	int nbicols = 0; /* number of binary input columns */
33			int bocols = 0; /* produce binary output columns */
34	greg	1.1	int tabc = '\t'; /* default separator */
35			int subtotal = 0; /* produce subtotals? */
36	greg	1.11	int nrecsout = 0; /* number of records produced */
37	greg	1.1
38	schorsch	1.2	static int execute(char *fname);
39	greg	1.1
40	schorsch	1.2	int
41			main(
42			int argc,
43			char *argv[]
44			)
45	greg	1.1	{
46			int status;
47			int a;
48
49			for (a = 1; a < argc; a++)
50			if (argv[a][0] == '-')
51			if (argv[a][1] >= '0' && argv[a][1] <= '9')
52			count = atoi(argv[a]+1);
53			else
54			switch (argv[a][1]) {
55			case 'm':
56			mean = !mean;
57			break;
58			case 'p':
59			func = MULT;
60			break;
61			case 's':
62			func = ADD;
63			power = atof(argv[a]+2);
64			break;
65			case 'u':
66			func = MAX;
67			break;
68			case 'l':
69			func = MIN;
70			break;
71			case 't':
72			tabc = argv[a][2];
73			break;
74			case 'r':
75			subtotal = !subtotal;
76			break;
77	greg	1.4	case 'i':
78			switch (argv[a][2]) {
79	greg	1.11	case 'n':
80			incnt = atol(argv[++a]);
81			break;
82	greg	1.4	case 'a':
83			nbicols = 0;
84			break;
85			case 'd':
86			if (isdigit(argv[a][3]))
87			nbicols = atoi(argv[a]+3);
88			else
89			nbicols = 1;
90			if (nbicols > MAXCOL) {
91			fprintf(stderr,
92			"%s: too many input columns\n",
93			argv[0]);
94			exit(1);
95			}
96			break;
97			case 'f':
98			if (isdigit(argv[a][3]))
99			nbicols = -atoi(argv[a]+3);
100			else
101			nbicols = -1;
102			if (-nbicols > MAXCOL) {
103			fprintf(stderr,
104			"%s: too many input columns\n",
105			argv[0]);
106			exit(1);
107			}
108			break;
109			default:
110			goto userr;
111			}
112			break;
113			case 'o':
114			switch (argv[a][2]) {
115	greg	1.11	case 'n':
116			outcnt = atol(argv[++a]);
117			break;
118	greg	1.4	case 'a':
119			bocols = 0;
120			break;
121			case 'd':
122			bocols = 1;
123			break;
124			case 'f':
125			bocols = -1;
126			break;
127			default:
128			goto userr;
129			}
130			break;
131			default:;
132			userr:
133	greg	1.1	fprintf(stderr,
134	greg	1.4	"Usage: %s [-m][-sE\|p\|u\|l][-tC][-i{f\|d}[N]][-o{f\|d}][-count [-r]] [file..]\n",
135	greg	1.1	argv[0]);
136			exit(1);
137			}
138			else
139			break;
140			if (mean) {
141			if (func == MAX) {
142			fprintf(stderr, "%s: average maximum?!\n", argv[0]);
143			exit(1);
144			}
145			if (func == MIN) {
146			fprintf(stderr, "%s: average minimum?!\n", argv[0]);
147			exit(1);
148			}
149			}
150	greg	1.4	if (bocols)
151			SET_FILE_BINARY(stdout);
152	greg	1.1	status = 0;
153			if (a >= argc)
154			status = execute(NULL) == -1 ? 1 : status;
155			else
156	greg	1.11	for ( ; a < argc && (outcnt <= 0 \|\| nrecsout < outcnt); a++)
157	greg	1.1	status = execute(argv[a]) == -1 ? 2 : status;
158			exit(status);
159			}
160
161
162	schorsch	1.2	static int
163	greg	1.4	getrecord( /* read next input record */
164			double field[MAXCOL],
165			FILE *fp
166			)
167			{
168			char buf[16*MAXCOL];
169			char cp, num;
170			int nf;
171			/* reading binary input? */
172			if (nbicols > 0)
173	greg	1.13	return(getbinary(field, sizeof(double), nbicols, fp));
174	greg	1.4	if (nbicols < 0) {
175			float fbuf = (float )buf;
176			int i;
177	greg	1.13	nf = getbinary(fbuf, sizeof(float), -nbicols, fp);
178	greg	1.4	for (i = nf; i-- > 0; )
179			field[i] = fbuf[i];
180			return(nf);
181			}
182			/* reading ASCII input */
183			cp = fgets(buf, sizeof(buf), fp);
184			if (cp == NULL)
185			return(0);
186			for (nf = 0; nf < MAXCOL; nf++) {
187			while (isspace(*cp))
188			cp++;
189			if (!*cp)
190			break;
191			num = cp;
192			while (cp && !(isspace(tabc)?isspace(cp):*cp==tabc))
193			cp++;
194			if (*cp)
195			*cp++ = '\0';
196			if (!num \|\| isspace(num))
197			field[nf] = init_val[func];
198			else
199			field[nf] = atof(num);
200			}
201			return(nf);
202			}
203
204
205			static void
206			putrecord( /* write out results record */
207			const double *field,
208			int n,
209			FILE *fp
210			)
211			{
212			/* binary output? */
213			if (bocols > 0) {
214	greg	1.13	putbinary(field, sizeof(double), n, fp);
215	greg	1.4	return;
216			}
217			if (bocols < 0) {
218			float fv;
219			while (n-- > 0) {
220			fv = *field++;
221	greg	1.13	putbinary(&fv, sizeof(float), 1, fp);
222	greg	1.4	}
223			return;
224			}
225			/* ASCII output */
226	greg	1.10	while (n-- > 0) {
227			fprintf(fp, "%.9g", *field++);
228			if (n) fputc(tabc, fp);
229			}
230	greg	1.4	fputc('\n', fp);
231			}
232
233
234			static int
235	schorsch	1.2	execute( /* compute result */
236			char *fname
237			)
238	greg	1.1	{
239	greg	1.4	double inpval[MAXCOL];
240			double tally[MAXCOL];
241	greg	1.7	short rsign[MAXCOL];
242	greg	1.1	double result[MAXCOL];
243	greg	1.7	int n;
244	greg	1.4	int nread, ncol;
245	greg	1.1	long nlin, ltotal;
246			FILE *fp;
247			/* open file */
248			if (fname == NULL)
249			fp = stdin;
250			else if ((fp = fopen(fname, "r")) == NULL) {
251			fprintf(stderr, "%s: cannot open\n", fname);
252			return(-1);
253			}
254	greg	1.4	if (nbicols)
255			SET_FILE_BINARY(fp);
256	greg	1.9	#ifdef getc_unlocked /* avoid lock/unlock overhead */
257			flockfile(fp);
258			#endif
259
260	greg	1.1	ltotal = 0;
261			while (!feof(fp)) {
262	schorsch	1.3	if (ltotal == 0) { /* initialize */
263	greg	1.1	if (func == MULT) /* special case */
264	greg	1.7	for (n = 0; n < MAXCOL; n++) {
265	greg	1.4	tally[n] = 0.0;
266	greg	1.7	rsign[n] = 1;
267			}
268	greg	1.1	else
269			for (n = 0; n < MAXCOL; n++)
270	greg	1.4	tally[n] = init_val[func];
271	schorsch	1.3	}
272	greg	1.1	ncol = 0;
273			for (nlin = 0; (count <= 0 \|\| nlin < count) &&
274	greg	1.11	(incnt <= 0 \|\| nlin < incnt) &&
275	greg	1.4	(nread = getrecord(inpval, fp)) > 0;
276	greg	1.1	nlin++) {
277			/* compute */
278	greg	1.4	for (n = 0; n < nread; n++)
279	greg	1.1	switch (func) {
280			case ADD:
281	greg	1.4	if (inpval[n] == 0.0)
282	greg	1.1	break;
283			if (power != 0.0)
284	greg	1.4	tally[n] += pow(fabs(inpval[n]),power);
285	greg	1.1	else
286	greg	1.4	tally[n] += inpval[n];
287	greg	1.1	break;
288			case MULT:
289	greg	1.4	if (inpval[n] == 0.0)
290	greg	1.7	rsign[n] = 0;
291			else if (inpval[n] < 0.0) {
292			rsign[n] = -rsign[n];
293			inpval[n] = -inpval[n];
294			}
295			if (rsign[n])
296			tally[n] += log(inpval[n]);
297	greg	1.1	break;
298			case MAX:
299	greg	1.4	if (inpval[n] > tally[n])
300			tally[n] = inpval[n];
301	greg	1.1	break;
302			case MIN:
303	greg	1.4	if (inpval[n] < tally[n])
304			tally[n] = inpval[n];
305	greg	1.1	break;
306			}
307	greg	1.4	if (nread > ncol)
308			ncol = nread;
309	greg	1.1	}
310			if (nlin == 0)
311			break;
312			/* compute and print */
313			ltotal += nlin;
314			for (n = 0; n < ncol; n++) {
315	greg	1.4	result[n] = tally[n];
316	greg	1.1	if (mean) {
317	greg	1.4	result[n] /= (double)ltotal;
318			if (func == ADD && power != 0.0 && result[n] != 0.0)
319			result[n] = pow(result[n], 1.0/power);
320	greg	1.1	}
321			if (func == MULT)
322	greg	1.7	result[n] = rsign[n] * exp(result[n]);
323	greg	1.1	}
324	greg	1.4	putrecord(result, ncol, stdout);
325	greg	1.11	++nrecsout;
326			if (outcnt > 0 && nrecsout >= outcnt)
327			break;
328	greg	1.1	if (!subtotal)
329			ltotal = 0;
330	greg	1.12	if (incnt > 0 && nlin >= incnt)
331			break;
332	greg	1.1	}
333	greg	1.4	/* close input */
334	greg	1.1	return(fclose(fp));
335			}