/*
 * test routine for WaveFile data file readers
 *
 * $Log: test_read.c,v $
 * Revision 1.6  2005/09/27 05:33:53  sgt
 * sweeps working decently well,
 * although wavelist selection for them is ugly.
 *
 * Revision 1.5  2003/07/30 06:18:49  sgt
 * better handling of the last point in a wavevar,
 * in particular when wv_interp_val asks for a point beyond the end of the iv range
 * enhance test_read.c to
 *
 * Revision 1.4  2000/08/09 23:37:39  sgt
 * ss_hspice.c - wrote sf_guessrows_hsbin routine.
 * others - instrumented to count reallocs and print out the number.
 *
 * Revision 1.3  2000/01/07 05:04:48  tell
 * updating with old changes as we construct the CVS
 *
 * Revision 1.2  1998/09/17 18:25:09  tell
 * prints out variable type
 *
 * Revision 1.1  1998/08/31 21:00:28  tell
 * Initial revision
 *
 */

#include <stdio.h>
#include <string.h>
#include <errno.h>
#include <glib.h>

#include "wavefile.h"

void test_interp(WvTable *wt, double mytm);
void dump_table_info(WvTable *wt);
void dump_wavevar(gpointer p, gpointer u);

int
main(int argc, char **argv)
{
	WaveFile *wf;
	WvTable *wt;
	int i, j;
	extern int optind;
	extern char *optarg;
	int v_flag = 0;
	int x_flag = 0;
	int l_flag = 0;
	int errflg = 0;
	char *filetype = NULL;
	int c;

	while ((c = getopt (argc, argv, "lt:vx")) != EOF) {
		switch(c) {
		case 'v':
			v_flag = 1;
			break;
		case 'l':
			l_flag = 1;
			break;
		case 't':
			filetype = optarg;
			break;
		case 'x':
			x_flag = 1;
			break;
		default:
			errflg = 1;
			break;
		}
	}

	if(errflg || optind >= argc)  {
		fprintf(stderr, "usage: %s [-ltvx] file\n", argv[0]);
		exit(1);
	}
	
	spicestream_msg_level = DBG;
	wf = wf_read(argv[optind], filetype);
	if(!wf) {
		if(errno)
			perror(argv[1]);
		fprintf(stderr, "test_read: unable to read file\n");
		exit(1);
	}
	printf("filename: \"%s\" #tables=%d\n", 
	       wf->wf_filename, wf->wf_ntables);

	for(i = 0; i < wf->wf_ntables; i++) {
		printf("table %d", i);
		wt = wf_wtable(wf, i);
		if(wt->name) {
			printf(" %s=%g", wt->name, wt->swval);
		}
		putchar('\n');
		dump_table_info(wt);
	}
	wf_foreach_wavevar(wf, dump_wavevar, NULL);

	if(l_flag) {
		int t;
		for(t = 0; t < wf->wf_ntables; t++) {
			printf("table %d: ", t);
			wt = wf_wtable(wf, t);
			if(wt->name) {
				printf(" %s=%g", wt->name, wt->swval);
			}
			putchar('\n');
			printf("      %10s", wt->iv->wv_name);
			for(i = 0; i < wt->wt_ndv; i++) {
				printf(" %10s", wt->dv[i].wv_name);
			}
			putchar('\n');
			for(j = 0; j < wt->nvalues; j++) {
				printf("[%3d] %10g", j, wds_get_point(wt->iv->wds, j));
				for(i = 0; i < wt->wt_ndv; i++) {
					printf(" %10g", 
					       wds_get_point(&wt->dv[i].wds[0], j));
				}
				putchar('\n');
			}
		}
	}

	wt = wf_wtable(wf, 0);
	if(v_flag) {
		double mytm;
		double delta;

		mytm = wds_get_point(wt->iv->wds, 0);
		delta = (wds_get_point(wt->iv->wds, wt->nvalues-1) - mytm) / 40.0;
		printf("40 divisions, delta=%g\n", delta);
		for(i = 0; i <= 41; i++, mytm += delta) {
			test_interp(wt, mytm);
		}
		mytm = wds_get_point(wt->iv->wds, wt->nvalues-2);
		putchar('\n');
		test_interp(wt, mytm);
   	}
	exit(0);
}

void dump_table_info(WvTable *wt)
{
	int i, j;
	WaveFile *wf = wt->wf;

	printf("independent variable:\n");
	printf("  name: \"%s\"\n", wt->iv->wv_name);
	printf("  type: %s\n", vartype_name_str(wt->iv->wv_type));
	printf("  npts: %d\n", wt->nvalues);
	printf("  min: %g\n", wt->iv->wds->min);
	printf("  max: %g\n", wt->iv->wds->max);
	printf("  blocks: %d/%d\n", wt->iv->wds->bpused, wt->iv->wds->bpsize);
	printf("  reallocs: %d\n", wt->iv->wds->nreallocs);

	printf("columns: %d\n", wf->wf_ncols);
	printf("dependent variables: %d\n", wf->wf_ndv);
}

void
dump_wavevar(gpointer p, gpointer u)
{
	WaveVar *wv = (WaveVar *)p;
	WaveFile *wf = wv->wv_file;
	int j;

	printf(" dv \"%s\" ", wv->wv_name);
	printf(" (type=%s)", vartype_name_str(wv->wv_type));

	if(wv->wv_ncols > 1)
		printf(" (%d columns)\n", wv->wv_ncols);

	for(j = 0; j < wv->wv_ncols; j++) {
		if(wv->wv_ncols > 1)
			printf("    col[%d] ", j);
		printf("blocks=%d/%d ",
		       wv->wds[j].bpused, wv->wds[j].bpsize);
		printf("min=%g ",wv->wds[j].min);
		printf("max=%g ", wv->wds[j].max);
		printf("first=%g ", wds_get_point(&wv->wds[j], 0));
		printf("last=%g\n", wds_get_point(&wv->wds[j],
						  wv->wv_nvalues-1));
	}
}

void
test_interp(WvTable *wt, double mytm)
{
	int idx;
	idx = wf_find_point(wt->iv, mytm);
	printf("last %8s < %14.8g is %14.8g at [%4d];",
	       wt->iv->wv_name,
	       mytm,
	       wds_get_point(wt->iv->wds, idx),
	       idx);
	fflush(stdout);
	printf("%8s at %8s=%14.8g is %14.8g\n",
	       wt->dv[0].wv_name,
	       wt->iv->wv_name,
	       mytm,
	       wv_interp_value(&wt->dv[0], mytm));
}