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_sea_level_radians.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_sea_level_radians.c')
| -rw-r--r-- | libtellurian_sea_level_radians.c | 24 |
1 files changed, 24 insertions, 0 deletions
diff --git a/libtellurian_sea_level_radians.c b/libtellurian_sea_level_radians.c index a23b222..e536025 100644 --- a/libtellurian_sea_level_radians.c +++ b/libtellurian_sea_level_radians.c @@ -1,5 +1,6 @@ /* See LICENSE file for copyright and license details. */ #include "common.h" +#ifndef TEST double @@ -15,3 +16,26 @@ libtellurian_sea_level_radians(double latitude) double denom = fma(x, c, y * s); return sqrt(num / denom); } + + +#else + + +int +main(void) +{ + double x = pow(LIBTELLURIAN_EQUATORIAL_RADIUS, 4.0) + pow(LIBTELLURIAN_POLAR_RADIUS, 4.0); + x = sqrt(x / (pow(LIBTELLURIAN_EQUATORIAL_RADIUS, 2.0) + pow(LIBTELLURIAN_POLAR_RADIUS, 2.0))); + ASSERT(libtellurian_sea_level_radians(0) == LIBTELLURIAN_EQUATORIAL_RADIUS); + ASSERT(libtellurian_sea_level_radians(D90) == LIBTELLURIAN_POLAR_RADIUS); + ASSERT(libtellurian_sea_level_radians(D180) == LIBTELLURIAN_EQUATORIAL_RADIUS); + ASSERT(libtellurian_sea_level_radians(-D90) == LIBTELLURIAN_POLAR_RADIUS); + ASSERT(libtellurian_sea_level_radians(D45) < LIBTELLURIAN_EQUATORIAL_RADIUS); + ASSERT(libtellurian_sea_level_radians(D45) > LIBTELLURIAN_POLAR_RADIUS); + ASSERT(libtellurian_sea_level_radians(-D45) == libtellurian_sea_level_radians(D45)); + ASSERT(fabs(libtellurian_sea_level_radians(D45) / x - 1.0) < 1.0e-12); + return 0; +} + + +#endif |