aboutsummaryrefslogtreecommitdiffstats
path: root/libtellurian_gaussian_radius_radians.c
diff options
context:
space:
mode:
authorMattias Andrée <m@maandree.se>2024-10-22 22:22:41 +0200
committerMattias Andrée <m@maandree.se>2024-10-22 22:22:41 +0200
commitbb5de3aa2ee118df78f0347cffd4e58f846dc1fb (patch)
treeb313a20b79d67f35a5e4e05bc9106b22dc17f451 /libtellurian_gaussian_radius_radians.c
parent... (diff)
downloadlibtellurian-bb5de3aa2ee118df78f0347cffd4e58f846dc1fb.tar.gz
libtellurian-bb5de3aa2ee118df78f0347cffd4e58f846dc1fb.tar.bz2
libtellurian-bb5de3aa2ee118df78f0347cffd4e58f846dc1fb.tar.xz
Signed-off-by: Mattias Andrée <m@maandree.se>
Diffstat (limited to 'libtellurian_gaussian_radius_radians.c')
-rw-r--r--libtellurian_gaussian_radius_radians.c2
1 files changed, 1 insertions, 1 deletions
diff --git a/libtellurian_gaussian_radius_radians.c b/libtellurian_gaussian_radius_radians.c
index 02d0586..7b14ed5 100644
--- a/libtellurian_gaussian_radius_radians.c
+++ b/libtellurian_gaussian_radius_radians.c
@@ -11,7 +11,7 @@ libtellurian_gaussian_radius_radians(double latitude)
double c = b / a;
double neg_e2 = fma(c, c, -1.0);
double neg_e2_plus_1 = c * c;
- double s = sin(latitude);
+ double s = sin(geodetic(latitude));
double denom = fma(neg_e2, s * s, 1.0);
return a * sqrt(neg_e2_plus_1) / denom;
}