Fourth commit

Signed-off-by: Mattias Andrée <m@maandree.se>

1 files changed, 41 insertions, 0 deletions

diff --git a/libtellurian_meridian_radius.c b/libtellurian_meridian_radius.c
new file mode 100644
index 0000000..f75862a
--- /dev/null
+++ b/libtellurian_meridian_radius.c
@@ -0,0 +1,41 @@
+/* See LICENSE file for copyright and license details. */
+#include "common.h"
+#ifndef TEST
+
+
+double
+libtellurian_meridian_radius(double latitude)
+{
+	latitude = radians(latitude);
+	return libtellurian_meridian_radius_radians(latitude);
+}
+
+
+#else
+
+
+static int
+approx(double a, double b)
+{
+        return fabs(a - b) <= 1e-8 * (0.5 * (a + b));
+}
+
+int
+main(void)
+{
+	ASSERT(libtellurian_meridian_radius(0) == libtellurian_meridian_radius_radians(0));
+	ASSERT(approx(libtellurian_meridian_radius(180), libtellurian_meridian_radius_radians(0)));
+	ASSERT(approx(libtellurian_meridian_radius(360), libtellurian_meridian_radius_radians(0)));
+	ASSERT(approx(libtellurian_meridian_radius(90), libtellurian_meridian_radius(-90)));
+	ASSERT(approx(libtellurian_meridian_radius(90), libtellurian_meridian_radius(270)));
+	ASSERT(approx(libtellurian_meridian_radius(45), libtellurian_meridian_radius(-45)));
+	ASSERT(approx(libtellurian_meridian_radius(135), libtellurian_meridian_radius(45)));
+	ASSERT(approx(libtellurian_meridian_radius(-135), libtellurian_meridian_radius(135)));
+	ASSERT(libtellurian_meridian_radius(90) == libtellurian_meridian_radius_radians(D90));
+	ASSERT(libtellurian_meridian_radius(45) == libtellurian_meridian_radius_radians(D45));
+	ASSERT(libtellurian_meridian_radius(30) == libtellurian_meridian_radius_radians(D30));
+	return 0;
+}
+
+
+#endif