Home | History | Annotate | Line # | Download | only in jot
jot.c revision 1.14
      1  1.14   garbled /*	$NetBSD: jot.c,v 1.14 2006/01/07 07:09:01 garbled 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   garbled __RCSID("$NetBSD: jot.c,v 1.14 2006/01/07 07:09:01 garbled 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.7  jdolecek 				strcpy(format, *argv + 2);
    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.7  jdolecek 				strcpy(format, *++argv);
    141   1.1       jtc 			break;
    142   1.1       jtc 		case 's':
    143   1.7  jdolecek 			if ((*argv)[2])
    144   1.7  jdolecek 				strcpy(sepstring, *argv + 2);
    145   1.7  jdolecek 			else if (!--argc)
    146   1.7  jdolecek 				errx(1, "Need string after -s");
    147   1.1       jtc 			else
    148   1.7  jdolecek 				strcpy(sepstring, *++argv);
    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