diff options
Diffstat (limited to 'src/projections/august.cpp')
| -rw-r--r-- | src/projections/august.cpp | 34 |
1 files changed, 34 insertions, 0 deletions
diff --git a/src/projections/august.cpp b/src/projections/august.cpp new file mode 100644 index 00000000..b5a21ef7 --- /dev/null +++ b/src/projections/august.cpp @@ -0,0 +1,34 @@ +#define PJ_LIB__ + +#include <math.h> + +#include "projects.h" + +PROJ_HEAD(august, "August Epicycloidal") "\n\tMisc Sph, no inv"; +#define M 1.333333333333333 + + +static XY s_forward (LP lp, PJ *P) { /* Spheroidal, forward */ + XY xy = {0.0,0.0}; + double t, c1, c, x1, x12, y1, y12; + (void) P; + + t = tan(.5 * lp.phi); + c1 = sqrt(1. - t * t); + c = 1. + c1 * cos(lp.lam *= .5); + x1 = sin(lp.lam) * c1 / c; + y1 = t / c; + xy.x = M * x1 * (3. + (x12 = x1 * x1) - 3. * (y12 = y1 * y1)); + xy.y = M * y1 * (3. + 3. * x12 - y12); + return (xy); +} + + + + +PJ *PROJECTION(august) { + P->inv = nullptr; + P->fwd = s_forward; + P->es = 0.; + return P; +} |
