/*
** table -- program to calculate tables for
**          the CORDIC algorithm
**
**          written by Joseph B. Evans, 11/15/91
**                     evans@shannon.tisl.ukans.edu
*/

#include <math.h>
long arctantab[32];

main(ac, av)
int ac;
char **av;
{
	double scale;
	int i,bits,msbits;

	if (ac == 4) {
		bits = atoi(av[2]);
		msbits = atoi(av[3]);
	} else {
		printf("usage: %s -d|r|b bits msbits\n",av[0]);
		exit(1);
	}

	switch (*(++av[1])) {
	case 'd':
		arctantab[0] = (1<<bits-msbits) * (180 / M_PI) * atan(2.0);
		arctantab[1] = (1<<bits-msbits) * (180 / M_PI) * atan(1.0);
		for (i = 2; i < bits; i++) {
			arctantab[i] = (1<<bits-msbits) * (180 / M_PI)
				* atan(1.0/(double)(1<<(i-1)));
		}
		break;
	case 'r':
		arctantab[0] = (1<<bits-msbits) * atan(2.0);
		arctantab[1] = (1<<bits-msbits) * atan(1.0);
		for (i = 2; i < bits; i++) {
			arctantab[i] = (1<<bits-msbits)
				* atan(1.0/(double)(1<<(i-1)));
		}
		break;
	case 'b':
		arctantab[0] = (1<<bits-msbits) * (2.0 / M_PI) * atan(2.0);
		arctantab[1] = (1<<bits-msbits) * (2.0 / M_PI) * atan(1.0);
		for (i = 2; i < bits; i++) {
			arctantab[i] = (1<<bits-msbits) * (2.0 / M_PI)
				* atan(1.0/(double)(1<<(i-1)));
		}
		break;
	default:
		printf("error: no angle mode\n");
		exit;
	}

	for (i = 0; i < bits; i++) printf("theta[%d] = %d\n",i,arctantab[i]);
}