#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 x, y;
        char *data;
} dgg[] = {
{0, 0, "vertex 1"},
{1, 0, "vertex 2 r4 1 0 0"},
{1, 1, "step 1 r4 1 0 1"},
{1, 2, "step 2 r4 1 0 2"},
{1, 3, "step 3 r4 1 0 3"},
{49, 6,  "triangle center 1 r4 1 3 6"},
{2, 0, "vertex 3"},
{2, 1, "step 3.1"},
{3, 0, "vertex 4"},
{4, 0, "vertex 5"},
{5, 0, "vertex 6"},
{6, 0, "vertex 7"},
{7, 0, "vertex 8"},
{8, 0, "vertex 9"},
{9, 0, "vertex 10"},
{10, 0, "vertex 11"},
{11, 0, "vertex 12"},
{0, 0, NULL}
};

struct dggap3 {
        int x, y;
        char *data;
} dggap3[] = {
{0, 0, "vertex 1"},
{1, 0, "vertex 2 r4 1 0 0"},
{1, 1, "step 1 r4 1 0 1"},
{1, 2, "step 2 r4 1 0 2"},
{1, 3, "step 3 r4 1 0 3"},
{49, 6,  "triangle center 1 r4 1 3 6"},
{6, 9, "vertex 3"},
{6, 9, "step 3.1"},
{9, 12, "vertex 4"},
{12, 15, "vertex 5"},
{15, 18, "vertex 6"},
{9, 3, "vertex 7"},
{12, 6, "vertex 8"},
{15, 9, "vertex 9"},
{18, 12, "vertex 10"},
{21, 15, "vertex 11"},
{15, 0, "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_HEX;
       	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(xy.x == dgg[j].x, "hex X %d == %d (%f, %f %s)",
				(int)xy.x, dgg[j].x,
			       	input[j].lon, input[j].lat, input[j].data
		  );
		ok(xy.y == dgg[j].y, "hex Y %d == %d", (int)xy.y, dgg[j].y);
	}

	if (0) {
	g.resolution = 3;
	for (j=0; j<i; j++) {
		if (j == 7 || j == 2 || j == 3) continue;
		ll.lon = input[j].lon * M_PI / 180.0;
		ll.lat = input[j].lat * M_PI / 180.0;
		xy = isea_forward(&g, &ll);
		skip_start(j == 7 || j == 2 || j == 3, 2, "skipping ap3 input %d as it falls on a hex side and is numerically unstable", j);
		ok(xy.x == dggap3[j].x, "hex ap3 X %d == %d (%f, %f %s)",
				(int)xy.x, dggap3[j].x,
			       	input[j].lon, input[j].lat, input[j].data
		  );
		ok(xy.y == dggap3[j].y, "hex ap3 Y %d == %d", (int)xy.y,
				dggap3[j].y);
		skip_end;
	}
	}

	return exit_status();

}