From 09e6857ae73249bc7433f2971dcf291c70e4c766 Mon Sep 17 00:00:00 2001 From: Mattias Andrée Date: Sun, 20 Oct 2024 17:28:46 +0200 Subject: Fourth commit MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit Signed-off-by: Mattias Andrée --- libtellurian_normal_gravity.c | 31 +++++++++++++++++++++++++++++++ 1 file changed, 31 insertions(+) (limited to 'libtellurian_normal_gravity.c') diff --git a/libtellurian_normal_gravity.c b/libtellurian_normal_gravity.c index 27722b1..a0d3e35 100644 --- a/libtellurian_normal_gravity.c +++ b/libtellurian_normal_gravity.c @@ -1,5 +1,6 @@ /* See LICENSE file for copyright and license details. */ #include "common.h" +#ifndef TEST double @@ -8,3 +9,33 @@ libtellurian_normal_gravity(double latitude) latitude = radians(latitude); return libtellurian_normal_gravity_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_normal_gravity(0) == libtellurian_normal_gravity_radians(0)); + ASSERT(libtellurian_normal_gravity(180) == libtellurian_normal_gravity_radians(0)); + ASSERT(libtellurian_normal_gravity(360) == libtellurian_normal_gravity_radians(0)); + ASSERT(approx(libtellurian_normal_gravity(90), libtellurian_normal_gravity(-90))); + ASSERT(approx(libtellurian_normal_gravity(90), libtellurian_normal_gravity(270))); + ASSERT(approx(libtellurian_normal_gravity(45), libtellurian_normal_gravity(-45))); + ASSERT(approx(libtellurian_normal_gravity(135), libtellurian_normal_gravity(45))); + ASSERT(approx(libtellurian_normal_gravity(-135), libtellurian_normal_gravity(135))); + ASSERT(libtellurian_normal_gravity(90) == libtellurian_normal_gravity_radians(D90)); + ASSERT(libtellurian_normal_gravity(45) == libtellurian_normal_gravity_radians(D45)); + ASSERT(libtellurian_normal_gravity(30) == libtellurian_normal_gravity_radians(D30)); + return 0; +} + + +#endif -- cgit v1.2.3-70-g09d2