Home | History | Annotate | Line # | Download | only in jot
jot.c revision 1.14.2.1
      1  1.14.2.1      chap /*	$NetBSD: jot.c,v 1.14.2.1 2006/06/19 04:17:06 chap Exp $	*/
      2       1.2       jtc 
      3       1.1       jtc /*-
      4       1.1       jtc  * Copyright (c) 1993
      5       1.1       jtc  *	The Regents of the University of California.  All rights reserved.
      6       1.1       jtc  *
      7       1.1       jtc  * Redistribution and use in source and binary forms, with or without
      8       1.1       jtc  * modification, are permitted provided that the following conditions
      9       1.1       jtc  * are met:
     10       1.1       jtc  * 1. Redistributions of source code must retain the above copyright
     11       1.1       jtc  *    notice, this list of conditions and the following disclaimer.
     12       1.1       jtc  * 2. Redistributions in binary form must reproduce the above copyright
     13       1.1       jtc  *    notice, this list of conditions and the following disclaimer in the
     14       1.1       jtc  *    documentation and/or other materials provided with the distribution.
     15      1.10       agc  * 3. Neither the name of the University nor the names of its contributors
     16       1.1       jtc  *    may be used to endorse or promote products derived from this software
     17       1.1       jtc  *    without specific prior written permission.
     18       1.1       jtc  *
     19       1.1       jtc  * THIS SOFTWARE IS PROVIDED BY THE REGENTS AND CONTRIBUTORS ``AS IS'' AND
     20       1.1       jtc  * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
     21       1.1       jtc  * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
     22       1.1       jtc  * ARE DISCLAIMED.  IN NO EVENT SHALL THE REGENTS OR CONTRIBUTORS BE LIABLE
     23       1.1       jtc  * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
     24       1.1       jtc  * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
     25       1.1       jtc  * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
     26       1.1       jtc  * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
     27       1.1       jtc  * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
     28       1.1       jtc  * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
     29       1.1       jtc  * SUCH DAMAGE.
     30       1.1       jtc  */
     31       1.1       jtc 
     32       1.4     lukem #include <sys/cdefs.h>
     33       1.1       jtc #ifndef lint
     34       1.4     lukem __COPYRIGHT("@(#) Copyright (c) 1993\n\
     35       1.4     lukem 	The Regents of the University of California.  All rights reserved.\n");
     36       1.1       jtc #endif /* not lint */
     37       1.1       jtc 
     38       1.1       jtc #ifndef lint
     39       1.2       jtc #if 0
     40       1.1       jtc static char sccsid[] = "@(#)jot.c	8.1 (Berkeley) 6/6/93";
     41       1.2       jtc #endif
     42  1.14.2.1      chap __RCSID("$NetBSD: jot.c,v 1.14.2.1 2006/06/19 04:17:06 chap Exp $");
     43       1.1       jtc #endif /* not lint */
     44       1.1       jtc 
     45       1.1       jtc /*
     46       1.1       jtc  * jot - print sequential or random data
     47       1.1       jtc  *
     48       1.1       jtc  * Author:  John Kunze, Office of Comp. Affairs, UCB
     49       1.1       jtc  */
     50       1.1       jtc 
     51       1.1       jtc #include <ctype.h>
     52       1.4     lukem #include <err.h>
     53       1.1       jtc #include <limits.h>
     54      1.14   garbled #include <math.h>
     55       1.1       jtc #include <stdio.h>
     56       1.1       jtc #include <stdlib.h>
     57       1.1       jtc #include <string.h>
     58       1.1       jtc #include <time.h>
     59       1.9    atatat #include <unistd.h>
     60       1.1       jtc 
     61       1.1       jtc #define	REPS_DEF	100
     62       1.1       jtc #define	BEGIN_DEF	1
     63       1.1       jtc #define	ENDER_DEF	100
     64       1.1       jtc #define	STEP_DEF	1
     65       1.1       jtc 
     66       1.7  jdolecek #define	is_default(s)	(strcmp((s), "-") == 0)
     67       1.1       jtc 
     68       1.1       jtc double	begin;
     69       1.1       jtc double	ender;
     70       1.1       jtc double	s;
     71       1.1       jtc long	reps;
     72       1.1       jtc int	randomize;
     73       1.1       jtc int	infinity;
     74       1.1       jtc int	boring;
     75       1.1       jtc int	prec;
     76       1.1       jtc int	dox;
     77       1.1       jtc int	chardata;
     78       1.1       jtc int	nofinalnl;
     79       1.3        pk char	sepstring[BUFSIZ] = "\n";
     80       1.1       jtc char	format[BUFSIZ];
     81       1.1       jtc 
     82      1.13     perry void	getargs(int, char *[]);
     83      1.13     perry void	getformat(void);
     84      1.13     perry int	getprec(char *);
     85      1.13     perry int	main(int, char **);
     86      1.13     perry void	putdata(double, long);
     87      1.13     perry static void	usage(void);
     88       1.1       jtc 
     89       1.1       jtc int
     90      1.13     perry main(int argc, char *argv[])
     91       1.1       jtc {
     92       1.1       jtc 	double	xd, yd;
     93       1.1       jtc 	long	id;
     94       1.4     lukem 	double	*x = &xd;
     95       1.4     lukem 	double	*y = &yd;
     96       1.4     lukem 	long	*i = &id;
     97       1.1       jtc 
     98       1.1       jtc 	getargs(argc, argv);
     99       1.1       jtc 	if (randomize) {
    100       1.1       jtc 		*x = (ender - begin) * (ender > begin ? 1 : -1);
    101       1.9    atatat 		srandom((unsigned long) s);
    102       1.1       jtc 		for (*i = 1; *i <= reps || infinity; (*i)++) {
    103       1.1       jtc 			*y = (double) random() / INT_MAX;
    104       1.1       jtc 			putdata(*y * *x + begin, reps - *i);
    105       1.1       jtc 		}
    106       1.1       jtc 	}
    107       1.1       jtc 	else
    108       1.1       jtc 		for (*i = 1, *x = begin; *i <= reps || infinity; (*i)++, *x += s)
    109       1.1       jtc 			putdata(*x, reps - *i);
    110       1.1       jtc 	if (!nofinalnl)
    111       1.1       jtc 		putchar('\n');
    112       1.1       jtc 	exit(0);
    113       1.1       jtc }
    114       1.1       jtc 
    115       1.1       jtc void
    116      1.13     perry getargs(int argc, char *argv[])
    117       1.1       jtc {
    118       1.4     lukem 	unsigned int	mask = 0;
    119       1.4     lukem 	int		n = 0;
    120       1.1       jtc 
    121       1.7  jdolecek 	while (--argc && **++argv == '-' && !is_default(*argv))
    122       1.7  jdolecek 		switch ((*argv)[1]) {
    123       1.1       jtc 		case 'r':
    124       1.1       jtc 			randomize = 1;
    125       1.1       jtc 			break;
    126       1.1       jtc 		case 'c':
    127       1.1       jtc 			chardata = 1;
    128       1.1       jtc 			break;
    129       1.1       jtc 		case 'n':
    130       1.1       jtc 			nofinalnl = 1;
    131       1.1       jtc 			break;
    132       1.1       jtc 		case 'b':
    133       1.1       jtc 			boring = 1;
    134       1.1       jtc 		case 'w':
    135       1.7  jdolecek 			if ((*argv)[2])
    136  1.14.2.1      chap 				strlcpy(format, *argv + 2, sizeof(format));
    137       1.7  jdolecek 			else if (!--argc)
    138       1.7  jdolecek 				errx(1, "Need context word after -w or -b");
    139       1.1       jtc 			else
    140  1.14.2.1      chap 				strlcpy(format, *++argv, sizeof(format));
    141       1.1       jtc 			break;
    142       1.1       jtc 		case 's':
    143       1.7  jdolecek 			if ((*argv)[2])
    144  1.14.2.1      chap 				strlcpy(sepstring, *argv + 2, sizeof(sepstring));
    145       1.7  jdolecek 			else if (!--argc)
    146       1.7  jdolecek 				errx(1, "Need string after -s");
    147       1.1       jtc 			else
    148  1.14.2.1      chap 				strlcpy(sepstring, *++argv, sizeof(sepstring));
    149       1.1       jtc 			break;
    150       1.1       jtc 		case 'p':
    151       1.7  jdolecek 			if ((*argv)[2])
    152       1.7  jdolecek 				prec = atoi(*argv + 2);
    153       1.7  jdolecek 			else if (!--argc)
    154       1.7  jdolecek 				errx(1, "Need number after -p");
    155       1.1       jtc 			else
    156       1.7  jdolecek 				prec = atoi(*++argv);
    157       1.1       jtc 			if (prec <= 0)
    158       1.7  jdolecek 				errx(1, "Bad precision value");
    159       1.1       jtc 			break;
    160       1.1       jtc 		default:
    161       1.7  jdolecek 			warnx("unknown option `%s'", *argv);
    162       1.7  jdolecek 			usage();
    163       1.1       jtc 		}
    164       1.1       jtc 
    165       1.7  jdolecek 	switch (argc) {	/* examine args right to left, falling thru cases */
    166       1.1       jtc 	case 4:
    167       1.7  jdolecek 		if (!is_default(argv[3])) {
    168       1.7  jdolecek 			if (!sscanf(argv[3], "%lf", &s))
    169       1.7  jdolecek 				errx(1, "Bad s value:  %s", argv[3]);
    170       1.1       jtc 			mask |= 01;
    171       1.1       jtc 		}
    172       1.1       jtc 	case 3:
    173       1.7  jdolecek 		if (!is_default(argv[2])) {
    174       1.7  jdolecek 			if (!sscanf(argv[2], "%lf", &ender))
    175       1.7  jdolecek 				ender = argv[2][strlen(argv[2])-1];
    176       1.1       jtc 			mask |= 02;
    177       1.1       jtc 			if (!prec)
    178       1.7  jdolecek 				n = getprec(argv[2]);
    179       1.1       jtc 		}
    180       1.1       jtc 	case 2:
    181       1.7  jdolecek 		if (!is_default(argv[1])) {
    182       1.7  jdolecek 			if (!sscanf(argv[1], "%lf", &begin))
    183       1.7  jdolecek 				begin = argv[1][strlen(argv[1])-1];
    184       1.1       jtc 			mask |= 04;
    185       1.1       jtc 			if (!prec)
    186       1.7  jdolecek 				prec = getprec(argv[1]);
    187       1.1       jtc 			if (n > prec)		/* maximum precision */
    188       1.1       jtc 				prec = n;
    189       1.1       jtc 		}
    190       1.1       jtc 	case 1:
    191       1.7  jdolecek 		if (!is_default(argv[0])) {
    192       1.7  jdolecek 			if (!sscanf(argv[0], "%ld", &reps))
    193       1.7  jdolecek 				errx(1, "Bad reps value:  %s", argv[0]);
    194       1.1       jtc 			mask |= 010;
    195       1.1       jtc 		}
    196       1.1       jtc 		break;
    197       1.1       jtc 	case 0:
    198       1.7  jdolecek 		usage();
    199       1.7  jdolecek 		break;
    200       1.1       jtc 	default:
    201       1.7  jdolecek 		errx(1, "Too many arguments.  What do you mean by %s?", argv[4]);
    202       1.1       jtc 	}
    203       1.1       jtc 	getformat();
    204       1.1       jtc 	while (mask)	/* 4 bit mask has 1's where last 4 args were given */
    205       1.1       jtc 		switch (mask) {	/* fill in the 0's by default or computation */
    206       1.1       jtc 		case 001:
    207       1.1       jtc 			reps = REPS_DEF;
    208       1.1       jtc 			mask = 011;
    209       1.1       jtc 			break;
    210       1.1       jtc 		case 002:
    211       1.1       jtc 			reps = REPS_DEF;
    212       1.1       jtc 			mask = 012;
    213       1.1       jtc 			break;
    214       1.1       jtc 		case 003:
    215       1.1       jtc 			reps = REPS_DEF;
    216       1.1       jtc 			mask = 013;
    217       1.1       jtc 			break;
    218       1.1       jtc 		case 004:
    219       1.1       jtc 			reps = REPS_DEF;
    220       1.1       jtc 			mask = 014;
    221       1.1       jtc 			break;
    222       1.1       jtc 		case 005:
    223       1.1       jtc 			reps = REPS_DEF;
    224       1.1       jtc 			mask = 015;
    225       1.1       jtc 			break;
    226       1.1       jtc 		case 006:
    227       1.1       jtc 			reps = REPS_DEF;
    228       1.1       jtc 			mask = 016;
    229       1.1       jtc 			break;
    230       1.1       jtc 		case 007:
    231       1.1       jtc 			if (randomize) {
    232       1.1       jtc 				reps = REPS_DEF;
    233       1.1       jtc 				mask = 0;
    234       1.1       jtc 				break;
    235       1.1       jtc 			}
    236       1.1       jtc 			if (s == 0.0) {
    237       1.1       jtc 				reps = 0;
    238       1.1       jtc 				mask = 0;
    239       1.1       jtc 				break;
    240       1.1       jtc 			}
    241       1.1       jtc 			reps = (ender - begin + s) / s;
    242       1.1       jtc 			if (reps <= 0)
    243       1.7  jdolecek 				errx(1, "Impossible stepsize");
    244       1.1       jtc 			mask = 0;
    245       1.1       jtc 			break;
    246       1.1       jtc 		case 010:
    247       1.1       jtc 			begin = BEGIN_DEF;
    248       1.1       jtc 			mask = 014;
    249       1.1       jtc 			break;
    250       1.1       jtc 		case 011:
    251       1.1       jtc 			begin = BEGIN_DEF;
    252       1.1       jtc 			mask = 015;
    253       1.1       jtc 			break;
    254       1.1       jtc 		case 012:
    255       1.9    atatat 			s = (randomize ? time(NULL) * getpid() : STEP_DEF);
    256       1.1       jtc 			mask = 013;
    257       1.1       jtc 			break;
    258       1.1       jtc 		case 013:
    259       1.1       jtc 			if (randomize)
    260       1.1       jtc 				begin = BEGIN_DEF;
    261       1.1       jtc 			else if (reps == 0)
    262       1.7  jdolecek 				errx(1, "Must specify begin if reps == 0");
    263       1.1       jtc 			begin = ender - reps * s + s;
    264       1.1       jtc 			mask = 0;
    265       1.1       jtc 			break;
    266       1.1       jtc 		case 014:
    267       1.9    atatat 			s = (randomize ? time(NULL) * getpid() : STEP_DEF);
    268       1.1       jtc 			mask = 015;
    269       1.1       jtc 			break;
    270       1.1       jtc 		case 015:
    271       1.1       jtc 			if (randomize)
    272       1.1       jtc 				ender = ENDER_DEF;
    273       1.1       jtc 			else
    274       1.1       jtc 				ender = begin + reps * s - s;
    275       1.1       jtc 			mask = 0;
    276       1.1       jtc 			break;
    277       1.1       jtc 		case 016:
    278       1.1       jtc 			if (randomize)
    279       1.9    atatat 				s = time(NULL) * getpid();
    280       1.1       jtc 			else if (reps == 0)
    281       1.7  jdolecek 				errx(1, "Infinite sequences cannot be bounded");
    282       1.1       jtc 			else if (reps == 1)
    283       1.1       jtc 				s = 0.0;
    284       1.1       jtc 			else
    285       1.1       jtc 				s = (ender - begin) / (reps - 1);
    286       1.1       jtc 			mask = 0;
    287       1.1       jtc 			break;
    288       1.1       jtc 		case 017:		/* if reps given and implied, */
    289       1.1       jtc 			if (!randomize && s != 0.0) {
    290       1.1       jtc 				long t = (ender - begin + s) / s;
    291       1.1       jtc 				if (t <= 0)
    292       1.7  jdolecek 					errx(1, "Impossible stepsize");
    293       1.1       jtc 				if (t < reps)		/* take lesser */
    294       1.1       jtc 					reps = t;
    295       1.1       jtc 			}
    296       1.1       jtc 			mask = 0;
    297       1.1       jtc 			break;
    298       1.1       jtc 		default:
    299       1.7  jdolecek 			errx(1, "bad mask");
    300       1.1       jtc 		}
    301       1.1       jtc 	if (reps == 0)
    302       1.1       jtc 		infinity = 1;
    303       1.1       jtc }
    304       1.1       jtc 
    305       1.1       jtc void
    306      1.13     perry putdata(double x, long notlast)
    307       1.1       jtc {
    308       1.4     lukem 	long	d = x;
    309       1.4     lukem 	long	*dp = &d;
    310       1.1       jtc 
    311       1.1       jtc 	if (boring)				/* repeated word */
    312       1.3        pk 		printf("%s", format);
    313       1.1       jtc 	else if (dox)				/* scalar */
    314       1.1       jtc 		printf(format, *dp);
    315      1.14   garbled 	else if (prec == 0)			/* integer */
    316      1.14   garbled 		printf(format, round(x));
    317       1.1       jtc 	else					/* real */
    318       1.1       jtc 		printf(format, x);
    319       1.1       jtc 	if (notlast != 0)
    320       1.1       jtc 		fputs(sepstring, stdout);
    321       1.1       jtc }
    322       1.1       jtc 
    323       1.7  jdolecek static void
    324       1.7  jdolecek usage(void)
    325       1.1       jtc {
    326      1.12     peter 	(void)fprintf(stderr, "usage: %s [-cnr] [-b word] [-p precision] "
    327      1.12     peter 	    "[-s string] [-w word] [reps [begin [end [s]]]]\n", getprogname());
    328       1.1       jtc 	exit(1);
    329       1.1       jtc }
    330       1.1       jtc 
    331       1.1       jtc int
    332      1.13     perry getprec(char *s)
    333       1.1       jtc {
    334       1.4     lukem 	char	*p;
    335       1.4     lukem 	char	*q;
    336       1.1       jtc 
    337       1.1       jtc 	for (p = s; *p; p++)
    338       1.1       jtc 		if (*p == '.')
    339       1.1       jtc 			break;
    340       1.1       jtc 	if (!*p)
    341       1.1       jtc 		return (0);
    342       1.1       jtc 	for (q = ++p; *p; p++)
    343       1.5  christos 		if (!isdigit((unsigned char)*p))
    344       1.1       jtc 			break;
    345       1.1       jtc 	return (p - q);
    346       1.1       jtc }
    347       1.1       jtc 
    348       1.1       jtc void
    349      1.13     perry getformat(void)
    350       1.1       jtc {
    351       1.4     lukem 	char	*p;
    352       1.7  jdolecek 	size_t	sz;
    353       1.1       jtc 
    354       1.1       jtc 	if (boring)				/* no need to bother */
    355       1.1       jtc 		return;
    356       1.1       jtc 	for (p = format; *p; p++)		/* look for '%' */
    357       1.7  jdolecek 		if (*p == '%') {
    358       1.7  jdolecek 			if (*(p+1) != '%')
    359       1.7  jdolecek 				break;
    360       1.7  jdolecek 			p++;		/* leave %% alone */
    361       1.7  jdolecek 		}
    362       1.7  jdolecek 	sz = sizeof(format) - strlen(format) - 1;
    363       1.7  jdolecek 	if (!*p) {
    364       1.7  jdolecek 		if (chardata) {
    365       1.8    simonb 			if (strlcpy(p, "%c", sz) >= sz)
    366       1.7  jdolecek 				errx(1, "-w word too long");
    367       1.8    simonb 			dox = 1;
    368       1.7  jdolecek 		} else {
    369       1.8    simonb 			if (snprintf(p, sz, "%%.%df", prec) >= (int)sz)
    370       1.7  jdolecek 				errx(1, "-w word too long");
    371       1.7  jdolecek 		}
    372       1.7  jdolecek 	} else if (!*(p+1)) {
    373       1.7  jdolecek 		if (sz <= 0)
    374       1.7  jdolecek 			errx(1, "-w word too long");
    375       1.1       jtc 		strcat(format, "%");		/* cannot end in single '%' */
    376       1.7  jdolecek 	} else {
    377       1.7  jdolecek 		p++;				/* skip leading % */
    378       1.7  jdolecek 		for(; *p && !isalpha((unsigned char)*p); p++) {
    379       1.7  jdolecek 			/* allow all valid printf(3) flags, but deny '*' */
    380       1.7  jdolecek 			if (!strchr("0123456789#-+. ", *p))
    381       1.7  jdolecek 				break;
    382       1.7  jdolecek 		}
    383       1.7  jdolecek 		/* Allow 'l' prefix, but no other. */
    384       1.7  jdolecek 		if (*p == 'l')
    385       1.7  jdolecek 			p++;
    386       1.1       jtc 		switch (*p) {
    387       1.1       jtc 		case 'f': case 'e': case 'g': case '%':
    388       1.7  jdolecek 		case 'E': case 'G':
    389       1.1       jtc 			break;
    390       1.1       jtc 		case 's':
    391       1.7  jdolecek 			errx(1, "cannot convert numeric data to strings");
    392       1.1       jtc 			break;
    393       1.7  jdolecek 		case 'd': case 'o': case 'x': case 'u':
    394       1.7  jdolecek 		case 'D': case 'O': case 'X': case 'U':
    395       1.7  jdolecek 		case 'c': case 'i':
    396       1.1       jtc 			dox = 1;
    397       1.1       jtc 			break;
    398       1.7  jdolecek 		default:
    399       1.7  jdolecek 			errx(1, "unknown or invalid format `%s'", format);
    400       1.1       jtc 		}
    401       1.7  jdolecek 		/* Need to check for trailing stuff to print */
    402       1.7  jdolecek 		for (; *p; p++)		/* look for '%' */
    403       1.7  jdolecek 			if (*p == '%') {
    404       1.7  jdolecek 				if (*(p+1) != '%')
    405       1.7  jdolecek 					break;
    406       1.7  jdolecek 				p++;		/* leave %% alone */
    407       1.7  jdolecek 			}
    408       1.7  jdolecek 		if (*p)
    409       1.7  jdolecek 			errx(1, "unknown or invalid format `%s'", format);
    410       1.1       jtc 	}
    411       1.1       jtc }
    412