#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>

#include "tap.h"
#include "isea.h"

#include "testing.h"

#include "cities.pt"

void iso(struct hex *h) {
	hex_iso(h);
	printf("iso(%d,%d,%d)", h->x, h->y, h->z);
}
void xy(struct hex *h) {
	hex_xy(h);
	printf("xy(%d,%d)", h->x, h->y);
}

struct dgg {
	int triangle;
	double x, y;
	char *data;
}  dgg[] = {
{5, 0.5485949, 0.1666408, "49900 Corvallis"},
{5, 0.5498729, 0.1533820, "42300 Aloha"},
{5, 0.5485337, 0.1538048, "41700 Tigard"},
{5, 0.5469778, 0.1653373, "41400 Albany"},
{5, 0.5477250, 0.1538168, "35700 LakeOswego"},
{5, 0.5484942, 0.1602033, "32600 Keizer"},
{5, 0.5514087, 0.1580878, "26800 MacMinnville"},
{5, 0.5463249, 0.1543639, "26100 OregonCity"},
{5, 0.5366604, 0.1942501, "23300 GrantsPass"},
{5, 0.5482398, 0.1544470, "23100 Tualatin"},
{5, 0.5468969, 0.1541197, "22500 WestLinn"},
{5, 0.5471126, 0.1531509, "20700 Milwaukie"},
{5, 0.5338688, 0.1458412, "20600 CitrusPark"},
{5, 0.5477809, 0.1577168, "20400 Woodburn"},
{5, 0.5415903, 0.1843270, "20300 Roseburg"},
{0, 0.0, 0.0, NULL}
};

int main(void) {
	int i, j;

	char isea[64], sahr[64];
	struct isea_dgg g;
	struct isea_pt xy;
	struct isea_geo ll;

	isea_grid_init(&g);
	g.output = ISEA_PROJTRI;
       	g.radius = ISEA_SCALE;

	i=0;
	while (cities[i].data) i++;

	plan_tests(i * 3);

	for (j=0; j<i; j++) {
		ll.lon = cities[j].lon * M_PI / 180.0;
		ll.lat = cities[j].lat * M_PI / 180.0;
		xy = isea_forward(&g, &ll);
		sprintf(isea, "%.7f", xy.x);
		sprintf(sahr, "%.7f", dgg[j].x);
		ok(!strcmp(isea,sahr), "projtri X (%f, %f %s) %s != %s",
			       	cities[j].lon, cities[j].lat, cities[j].data,
				isea, sahr,
				xy.x, dgg[j].x);
		ok(!strcmp(isea,sahr), "projtri Y (%f, %f %s) %f != %f",
			       	cities[j].lon, cities[j].lat, cities[j].data,
				xy.y, dgg[j].y);
		ok(g.triangle-1 == dgg[j].triangle,
			       	"projtri T (%f, %f %s) %d != %d",
			       	cities[j].lon, cities[j].lat, cities[j].data,
				g.triangle-1, dgg[j].triangle);
	}

	return exit_status();

}