/*
# This file is part of the Astrometry.net suite.
# Licensed under a 3-clause BSD style license - see LICENSE
*/

#include <unistd.h>
#include <stdio.h>

#include "starutil.h"
#include "cutest.h"

void test_ra2hmsstring(CuTest* tc) {
	char rastr[32], decstr[32];
	double hr;
	double as;

	ra2hmsstring(0, rastr);
	CuAssertStrEquals(tc, "00:00:00.000", rastr);
	printf("ra %s\n", rastr);

	hr = 15.0;

	ra2hmsstring(-hr, rastr);
	CuAssertStrEquals(tc, "23:00:00.000", rastr);
	printf("ra %s\n", rastr);

	ra2hmsstring(hr, rastr);
	CuAssertStrEquals(tc, "01:00:00.000", rastr);
	printf("ra %s\n", rastr);

	as = 15.0/3600.0;

	ra2hmsstring(60.0*as, rastr);
	CuAssertStrEquals(tc, "00:01:00.000", rastr);
	printf("ra %s\n", rastr);

	ra2hmsstring(as, rastr);
	CuAssertStrEquals(tc, "00:00:01.000", rastr);
	printf("ra %s\n", rastr);

	ra2hmsstring(0.001*as, rastr);
	CuAssertStrEquals(tc, "00:00:00.001", rastr);
	printf("ra %s\n", rastr);

	ra2hmsstring(0.000999*as, rastr);
	CuAssertStrEquals(tc, "00:00:00.001", rastr);
	printf("ra %s\n", rastr);

	ra2hmsstring(360.0 - 0.001*as, rastr);
	CuAssertStrEquals(tc, "23:59:59.999", rastr);
	printf("ra %s\n", rastr);

	as = 1.0/3600.0;

	dec2dmsstring(0, decstr);
	CuAssertStrEquals(tc, "+00:00:00.000", decstr);
	printf("dec %s\n", decstr);

	dec2dmsstring(-0, decstr);
	CuAssertStrEquals(tc, "+00:00:00.000", decstr);
	printf("dec %s\n", decstr);

	dec2dmsstring(-1, decstr);
	CuAssertStrEquals(tc, "-01:00:00.000", decstr);
	printf("dec %s\n", decstr);

	dec2dmsstring(-1*as, decstr);
	CuAssertStrEquals(tc, "-00:00:01.000", decstr);
	printf("dec %s\n", decstr);

	dec2dmsstring(-0.001*as, decstr);
	CuAssertStrEquals(tc, "-00:00:00.001", decstr);
	printf("dec %s\n", decstr);

	dec2dmsstring(-0.000999*as, decstr);
	CuAssertStrEquals(tc, "-00:00:00.001", decstr);
	printf("dec %s\n", decstr);

	dec2dmsstring(as, decstr);
	CuAssertStrEquals(tc, "+00:00:01.000", decstr);
	printf("dec %s\n", decstr);

	dec2dmsstring(0.5, decstr);
	CuAssertStrEquals(tc, "+00:30:00.000", decstr);
	printf("dec %s\n", decstr);

	dec2dmsstring(-0.5, decstr);
	CuAssertStrEquals(tc, "-00:30:00.000", decstr);
	printf("dec %s\n", decstr);



}


void test_hammer_aitoff(CuTest* tc) {
	double xyz[3];
	double px, py;
	double ras [] = { 0,  0,  0,   0,   0, 45, 90, 135, 180, 225, 270, 315, 360};
	double decs[] = { 0, 45, 90, -45, -89,  0,  0,   0,   0,   0,   0,   0,   0};
	int i;
	for (i=0; i<sizeof(ras)/sizeof(double); i++) {
		double ra, dec;
		ra  = ras [i];
		dec = decs[i];
		radecdeg2xyzarr(ra, dec, xyz);
		project_hammer_aitoff_x(xyz[0], xyz[1], xyz[2], &px, &py);
		printf("RA,Dec (%f,%f) => (%g, %g)\n", ra, dec, px, py);
	}
}

void test_atora(CuTest* tc) {
    CuAssertDblEquals(tc, 0.0, atora("00:00:00.0"), 1e-6);
    CuAssertDblEquals(tc, 15.0, atora("01:00:00.0"), 1e-6);
    CuAssertDblEquals(tc, 0.25, atora("00:01:00.0"), 1e-6);
    CuAssertDblEquals(tc, 0.25/60.0, atora("00:00:01.0"), 1e-6);
    CuAssertDblEquals(tc, 0.25/60.0/10.0, atora("00:00:00.1"), 1e-6);

    CuAssertDblEquals(tc, 0.0, atora("0:00:00.0"), 1e-6);
    CuAssertDblEquals(tc, 15.0, atora("1:00:00.0"), 1e-6);
    CuAssertDblEquals(tc, 0.25, atora("0:01:00.0"), 1e-6);
    CuAssertDblEquals(tc, 0.25/60.0, atora("0:00:01.0"), 1e-6);
    CuAssertDblEquals(tc, 0.25/60.0/10.0, atora("0:00:00.1"), 1e-6);

    CuAssertDblEquals(tc, 0.0, atora("00:0:00.0"), 1e-6);
    CuAssertDblEquals(tc, 15.0, atora("01:0:00.0"), 1e-6);
    CuAssertDblEquals(tc, 0.25, atora("00:1:00.0"), 1e-6);
    CuAssertDblEquals(tc, 0.25/60.0, atora("00:0:01.0"), 1e-6);
    CuAssertDblEquals(tc, 0.25/60.0/10.0, atora("00:0:00.1"), 1e-6);

    CuAssertDblEquals(tc, 6 * 15 + 8 * 0.25 + 41.0 * 0.25/60.0, atora("06:08:41.0"), 1e-6);
    CuAssertDblEquals(tc, 6 * 15 + 8 * 0.25 + 41.0 * 0.25/60.0, atora("6:08:41.0"), 1e-6);
    CuAssertDblEquals(tc, 6 * 15 + 8 * 0.25 + 41.0 * 0.25/60.0, atora("6:8:41.0"), 1e-6);
}
void test_atodec(CuTest* tc) {
    CuAssertDblEquals(tc, 0.0, atodec("00:00:00.0"), 1e-6);
    CuAssertDblEquals(tc, 1.0, atodec("01:00:00.0"), 1e-6);
    CuAssertDblEquals(tc, 1.0/60.0, atodec("00:01:00.0"), 1e-6);
    CuAssertDblEquals(tc, 1.0/3600.0, atodec("00:00:01.0"), 1e-6);
    CuAssertDblEquals(tc, 1.0/36000.0, atodec("00:00:00.1"), 1e-6);
    CuAssertDblEquals(tc, -1.0, atodec("-01:00:00.0"), 1e-6);
    CuAssertDblEquals(tc, -(1.0 + 1.0/60.0), atodec("-01:01:00.0"), 1e-6);
}

void test_xtodistsq(CuTest* tc) {
	double distsq = 1e-4;
	double x;

	x = distsq2rad(distsq);
	CuAssertDblEquals(tc, distsq, rad2distsq(x), 1e-8);

	x = distsq2arcsec(distsq);
	CuAssertDblEquals(tc, distsq, arcsec2distsq(x), 1e-8);
}