diff options
Diffstat (limited to '')
-rw-r--r-- | libtellurian_meridian_radius.c | 41 |
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 |