path: root/libtellurian_azimuth_radians.c

/* See LICENSE file for copyright and license details. */
#include "common.h"
#ifndef TEST


void
libtellurian_azimuth_radians(double latitude1, double longitude1,
                             double latitude2, double longitude2,
                             double *azimuth1_out, double *azimuth2_out)
{
	libtellurian_vincenty_inverse__(latitude1, longitude1,
	                                latitude2, longitude2,
	                                NULL, azimuth1_out, azimuth2_out);
}


#else


static int
eq(double a, double b)
{
	if (isfinite(a) && isfinite(b))
		return a == b;
	if (isnan(a) && isnan(b))
		return 1;
	printf("Unexpected floating-point class\n");
	return 0;
}

static int
check(double a, double b, double c, double d)
{
	double az1_ref = 1, az2_ref = 2, az1 = 3, az2 = 4, az1b = 5, az2b = 6;
	double s;

	s = libtellurian_distance_radians(a, b, c, d, &az1_ref, &az2_ref);
	(void) s;

	libtellurian_azimuth_radians(a, b, c, d, &az1, &az2);
	libtellurian_azimuth_radians(a, b, c, d, &az1b, NULL);
	libtellurian_azimuth_radians(a, b, c, d, NULL, &az2b);
	libtellurian_azimuth_radians(a, b, c, d, NULL, NULL);
	return eq(az1, az1_ref) && eq(az2, az2_ref) && eq(az1, az1b) && eq(az2, az2b);
}

int
main(void)
{
	ASSERT(check(1, 2, 3, 4));
	ASSERT(check(2, 3, 4, 5));
	ASSERT(check(3, 4, 5, 6));
	ASSERT(check(7, 8, 7, 8));
	ASSERT(check(0, 0, 0, 0));
	ASSERT(check(D90, D90, D90, D90));
	return 0;
}


#endif