diff options
| author | Mattias Andrée <m@maandree.se> | 2024-10-19 18:21:54 +0200 |
|---|---|---|
| committer | Mattias Andrée <m@maandree.se> | 2024-10-19 18:21:54 +0200 |
| commit | 3ce34980d7ba2bddbf3e9a1bd6f98cbc855bddc2 (patch) | |
| tree | 7724cc9cb9c3b2ce5ee60d2e012d1fbeb80aacfb /libtellurian_coarse_distance_radians.c | |
| download | libtellurian-3ce34980d7ba2bddbf3e9a1bd6f98cbc855bddc2.tar.gz libtellurian-3ce34980d7ba2bddbf3e9a1bd6f98cbc855bddc2.tar.bz2 libtellurian-3ce34980d7ba2bddbf3e9a1bd6f98cbc855bddc2.tar.xz | |
First commit
Signed-off-by: Mattias Andrée <m@maandree.se>
Diffstat (limited to 'libtellurian_coarse_distance_radians.c')
| -rw-r--r-- | libtellurian_coarse_distance_radians.c | 13 |
1 files changed, 13 insertions, 0 deletions
diff --git a/libtellurian_coarse_distance_radians.c b/libtellurian_coarse_distance_radians.c new file mode 100644 index 0000000..eb9cac0 --- /dev/null +++ b/libtellurian_coarse_distance_radians.c @@ -0,0 +1,13 @@ +/* See LICENSE file for copyright and license details. */ +#include "common.h" + + +double +libtellurian_coarse_distance_radians(double latitude1, double longitude1, + double latitude2, double longitude2) +{ + double h = fma(haversin(longitude2 - longitude1), + cos(latitude1) * cos(latitude2), + haversin(latitude2 - latitude1)); + return 2.0 * LIBTELLURIAN_AUTHALIC_RADIUS * asin(sqrt(h)); +} |