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

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

#include "testing.h"

#include "hexes.pt"

struct dgg {
	int seq_a3r4, seq_a3r3;
	char *data;
} dgg[] = {
{1, 1, "vertex 1"},
{2, 2, "vertex 2 r4 1 0 0"},
{3, 5, "step 1 r4 1 0 1"},
{4, 3, "step 2 r4 1 0 2"},
{5, 3, "step 3 r4 1 0 3"},
{35, 13,  "triangle center 1 r4 1 3 6"},
{83, 29, "vertex 3"},
{84, 32, "step 3.1"},
{164, 56, "vertex 4"},
{245, 83, "vertex 5"},
{326, 110, "vertex 6"},
{407, 137, "vertex 7"},
{488, 164, "vertex 8"},
{569, 191, "vertex 9"},
{650, 218, "vertex 10"},
{731, 245, "vertex 11"},
{812, 272, "vertex 12"},
{0, 0, NULL}
};

int main(void) {
	int i, j;

	struct isea_dgg g;
	struct isea_pt xy;
	struct isea_geo ll;

	isea_grid_init(&g);
	isea_orient_pole(&g);
	g.output = ISEA_SEQNUM;
       	g.radius = ISEA_SCALE;
	g.aperture = 3;
	g.resolution = 4;

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

	plan_tests(i * 2 - 1);

	for (j=0; j<i; j++) {
		ll.lon = input[j].lon * M_PI / 180.0;
		ll.lat = input[j].lat * M_PI / 180.0;
		xy = isea_forward(&g, &ll);
		ok(g.serial == dgg[j].seq_a3r4,
				"seqnum a3r4 %d == %d (%f, %f %s)",
				g.serial, dgg[j].seq_a3r4,
			       	input[j].lon, input[j].lat, input[j].data
		  );
	}

	g.resolution = 3;
	for (j=0; j<i; j++) {
		if (j==7) continue; /* falls on hex edge, so
				       translated hex depends on rounding
				       or precision */
		ll.lon = input[j].lon * M_PI / 180.0;
		ll.lat = input[j].lat * M_PI / 180.0;
		xy = isea_forward(&g, &ll);
		ok(g.serial == dgg[j].seq_a3r3,
				"seqnum a3r3 %d == %d (%f, %f %s)",
				g.serial, dgg[j].seq_a3r3,
			       	input[j].lon, input[j].lat, input[j].data
		  );
	}

	return exit_status();

}