diff options
author | Mattias Andrée <m@maandree.se> | 2024-10-20 17:28:46 +0200 |
---|---|---|
committer | Mattias Andrée <m@maandree.se> | 2024-10-20 17:28:46 +0200 |
commit | 09e6857ae73249bc7433f2971dcf291c70e4c766 (patch) | |
tree | f07ef152c6372083ab87e2e4289c7ab8b1f2c1ad /libtellurian_transverse_radius.c | |
parent | Third commit (diff) | |
download | libtellurian-09e6857ae73249bc7433f2971dcf291c70e4c766.tar.gz libtellurian-09e6857ae73249bc7433f2971dcf291c70e4c766.tar.bz2 libtellurian-09e6857ae73249bc7433f2971dcf291c70e4c766.tar.xz |
Fourth commit
Signed-off-by: Mattias Andrée <m@maandree.se>
Diffstat (limited to 'libtellurian_transverse_radius.c')
-rw-r--r-- | libtellurian_transverse_radius.c | 31 |
1 files changed, 31 insertions, 0 deletions
diff --git a/libtellurian_transverse_radius.c b/libtellurian_transverse_radius.c index d94ca86..15bc2f5 100644 --- a/libtellurian_transverse_radius.c +++ b/libtellurian_transverse_radius.c @@ -1,5 +1,6 @@ /* See LICENSE file for copyright and license details. */ #include "common.h" +#ifndef TEST double @@ -8,3 +9,33 @@ libtellurian_transverse_radius(double latitude) latitude = radians(latitude); return libtellurian_transverse_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_transverse_radius(0) == libtellurian_transverse_radius_radians(0)); + ASSERT(approx(libtellurian_transverse_radius(180), libtellurian_transverse_radius_radians(0))); + ASSERT(approx(libtellurian_transverse_radius(360), libtellurian_transverse_radius_radians(0))); + ASSERT(approx(libtellurian_transverse_radius(90), libtellurian_transverse_radius(-90))); + ASSERT(approx(libtellurian_transverse_radius(90), libtellurian_transverse_radius(270))); + ASSERT(approx(libtellurian_transverse_radius(45), libtellurian_transverse_radius(-45))); + ASSERT(approx(libtellurian_transverse_radius(135), libtellurian_transverse_radius(45))); + ASSERT(approx(libtellurian_transverse_radius(-135), libtellurian_transverse_radius(135))); + ASSERT(libtellurian_transverse_radius(90) == libtellurian_transverse_radius_radians(D90)); + ASSERT(libtellurian_transverse_radius(45) == libtellurian_transverse_radius_radians(D45)); + ASSERT(libtellurian_transverse_radius(30) == libtellurian_transverse_radius_radians(D30)); + return 0; +} + + +#endif |