blob: 0a832e568173e3e1b6381537725cf344cfb94d51 (
plain) (
blame)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
|
/* See LICENSE file for copyright and license details. */
#include "common.h"
double
libtellurian_meridan_radius_radians(double latitude)
{
double a = LIBTELLURIAN_EQUATORIAL_RADIUS;
double b = LIBTELLURIAN_POLAR_RADIUS;
double f = 1.0 - b / a;
double neg_e2 = (f - 2.0) * f;
double s = sin(latitude);
return fma(a, neg_e2, a) / pow(fma(neg_e2, s * s, 1.0), 1.5);
}
|