diff options
Diffstat (limited to 'libtellurian_sea_level.c')
-rw-r--r-- | libtellurian_sea_level.c | 31 |
1 files changed, 31 insertions, 0 deletions
diff --git a/libtellurian_sea_level.c b/libtellurian_sea_level.c index 76f817f..03893e2 100644 --- a/libtellurian_sea_level.c +++ b/libtellurian_sea_level.c @@ -1,5 +1,6 @@ /* See LICENSE file for copyright and license details. */ #include "common.h" +#ifndef TEST double @@ -8,3 +9,33 @@ libtellurian_sea_level(double latitude) latitude = radians(latitude); return libtellurian_sea_level_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_sea_level(0) == libtellurian_sea_level_radians(0)); + ASSERT(approx(libtellurian_sea_level(180), libtellurian_sea_level_radians(0))); + ASSERT(approx(libtellurian_sea_level(360), libtellurian_sea_level_radians(0))); + ASSERT(approx(libtellurian_sea_level(90), libtellurian_sea_level(-90))); + ASSERT(approx(libtellurian_sea_level(90), libtellurian_sea_level(270))); + ASSERT(approx(libtellurian_sea_level(45), libtellurian_sea_level(-45))); + ASSERT(approx(libtellurian_sea_level(135), libtellurian_sea_level(45))); + ASSERT(approx(libtellurian_sea_level(-135), libtellurian_sea_level(135))); + ASSERT(libtellurian_sea_level(90) == libtellurian_sea_level_radians(D90)); + ASSERT(libtellurian_sea_level(45) == libtellurian_sea_level_radians(D45)); + ASSERT(libtellurian_sea_level(30) == libtellurian_sea_level_radians(D30)); + return 0; +} + + +#endif |