/*
 *	General Parabolic Blending
 *
 * Looks better at 4 spaces/tab
 *
 */

#include <stdlib.h>
#include <stdio.h>
#include <math.h>
#ifdef THINK_C
#include <console.h>
#endif
#ifdef MYGETOPT
char getopt( int	argc, char	*argv[], char *format );
#endif
#ifndef true
#define	true	1
#define false	0
#endif
#include "parab.h"
#define frand(a) ((-0.5 + ((double)rand()/(double)RAND_MAX))*(a))

#define MAX_DATA	10



void	usage(char *name);
void	set_opts( int	argc, char	*argv[] );
void	DoError(char *str);
int		get_format(void);

static char __ident[] = "@(#)(c) SPDsoft, 28 diciembre 1993. v2";
#define VERS_STR        ((char*)&__ident[4])


int		NumLines;
int		Verbose;
int		SuperSampling;
unsigned short int	SFPI;
double	FPS;
int		Stochastic;
int		PutTime;
int		PutPolar;
int		AbsTime;
double	Desv;



int	get_format(void)
{
	int done=0;
	int i;
	
	sz_data=0;
	tid = -1;
	did=0;
		
	while(!done)
	{
		switch(getc(stdin))
		{
			case '%':
				if(did==MAX_DATA)
					DoError("Too many data!");
					
				switch(getc(stdin))
				{
					case 'c':
					case 'C':
						tdata[did]='c';
						sz_data+=3;
						did++;	
						break;
					case 'p':
					case 'P':
						tdata[did]='p';
						sz_data+=3;
						did++;	
						break;
					case 'x':
					case 'X':
						tdata[did]='x';
						sz_data+=1;
						did++;	
						break;
					case 't':
						AbsTime=false;
					case 'T':
						tdata[did]='t';
						tid=sz_data;
						sz_data+=1;
						did++;	
						break;
					case '\n':
						done=-2;
						DoError("Unexpected end of line");
						break;
					case EOF:
						done=-1;
						DoError("Unexpected end of file");
						break;
					default:
						done=-3;
						DoError("Unexpected format");
						break;
				}
				break;
			case '\n':
				done=1;
				break;
			case EOF:
				done=-1;
				DoError("Unexpected end of file");
				break;
			default:
				break;
		}
	}
	if(tid==-1) DoError("No time");
	if(tid+1!=sz_data) DoError("Time must be the last token");
	
	for(i=0;i<EVAL_WINDOW;i++)
		if((data[i]=malloc(sz_data))==NULL)
			DoError("No memory");

	return(done);
}

void	DoError(char *str)
{
	fprintf(stderr,"Error: %s\n",str);
	exit(2);
}

void	usage(char *name)
{
	fprintf(stderr,"Uso: %s [-tPcvnrph]\n",name);
	fprintf(stderr,"%s: -c: Num lines\n",name);
	fprintf(stderr,"%s: -t: don't print time\n",name);
	fprintf(stderr,"%s: -P: print polar as polar\n",name);
	fprintf(stderr,"%s: -v: ver\n",name);
	fprintf(stderr,"%s: -n #: Subframes/frame\n",name);
	fprintf(stderr,"%s: -r #: frames/second (default 25)\n",name);
	fprintf(stderr,"%s: -p #: Max Desviation (per unit)\n",name);
	fprintf(stderr,"%s: -h: This text\n",name);
	fprintf(stderr,
	"\n%s: Heading format: %%c point %%p polar %%x float "
	"%%T/t time abs/dif\n",name);
	exit(1);
}

void	set_opts( int	argc, char	*argv[] )
{
	int	opcion;
	extern char     *optarg;
	extern int      optind,opterr;


	NumLines		=	false;
	Verbose			=	false;
	SuperSampling	=	false;
	SFPI			=	1;
	FPS				=	25.0;
	Stochastic		=	true;
	PutTime			=	true;
	PutPolar		=	false;
	Desv			=	0.0;
	AbsTime			=	true;

	while ( (opcion=getopt(argc,argv,"tPcvn:r:p:h")) != EOF )
	{
		switch(opcion)
		{
			case 't':
					PutTime=false;
					break;
			case 'P':
					PutPolar=true;
					break;
			case 'c':
					NumLines=true;
					break;
			case 'v':
					printf("%s: %s\n",argv[0],VERS_STR);
					Verbose=true;
					break;
			case 'n':
					SuperSampling=true;
					SFPI = atoi(optarg);
					break;
			case 'r':
					FPS = atoi(optarg);
					break;
			case 'p':
					Stochastic=true;
					Desv = atof(optarg);
					break;
			default:
					usage(argv[0]);
					break;
		}
	}
}

main( int	argc, char	*argv[] )
{
	int	done=0,i;
	double time, delta;
	
#ifdef THINK_C
	argc = ccommand(&argv);
#endif
	set_opts( argc, argv );
	
	get_format();
	init_data();
	delta=1.0/((double)FPS*(double)SFPI);
	
	for(time=data[0][tid];!done;time += delta )
	{
		done = get_data( Stochastic ? time+frand(Desv*delta) : time );
		if (!done) out_data( NumLines , PutTime, PutPolar);	
	}
  	
	flush_data();
	
	for(i=0;i<EVAL_WINDOW;i++)
		if(data[i]!=NULL) free (data[i]);
		
	exit(0);
}