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

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

#include "testing.h"

#include "cities.pt"

struct dgg {
	int quad;
	double x, y;
	char *data;
}  dgg[] = {
{1, 0.3700177, 0.3076079, "49900 Corvallis"},
{1, 0.3578962, 0.3131305, "42300 Aloha"},
{1, 0.3589320, 0.3140789, "41700 Tigard"},
{1, 0.3696974, 0.3096601, "41400 Albany"},
{1, 0.3593468, 0.3147732, "35700 LakeOswego"},
{1, 0.3644930, 0.3109138, "32600 Keizer"},
{1, 0.3612037, 0.3094476, "26800 MacMinnville"},
{1, 0.3605206, 0.3157122, "26100 OregonCity"},
{1, 0.3998953, 0.3041388, "23300 GrantsPass"},
{1, 0.3596352, 0.3140123, "23100 Tualatin"},
{1, 0.3600231, 0.3153389, "22500 WestLinn"},
{1, 0.3590763, 0.3156365, "20700 Milwaukie"},
{1, 0.3593678, 0.3307608, "20600 CitrusPark"},
{1, 0.3626963, 0.3127749, "20400 Woodburn"},
{1, 0.3888368, 0.3048309, "20300 Roseburg"},
{1, 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_Q2DD;
       	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), "q2dd 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), "q2dd Y (%f, %f %s) %f == %f",
			       	cities[j].lon, cities[j].lat, cities[j].data,
				xy.y, dgg[j].y);
		ok(g.quad == dgg[j].quad,
			       	"q2dd Q (%f, %f %s) %d == %d",
			       	cities[j].lon, cities[j].lat, cities[j].data,
				g.quad+1, dgg[j].quad);
	}

	return exit_status();

}