#include <math.h>

#include "proj.h"

#include "proj_internal.h"

PROJ_HEAD(denoy, "Denoyer Semi-Elliptical") "\n\tPCyl, no inv, Sph";

#define C0 0.95

#define C1 -0.08333333333333333333

#define C3 0.00166666666666666666

#define D1 0.9

#define D5 0.03

static PJ_XY denoy_s_forward(PJ_LP lp, PJ *P) { /* Spheroidal, forward */

    PJ_XY xy = {0.0, 0.0};

    (void)P;

    xy.y = lp.phi;

    xy.x = lp.lam;

    lp.lam = fabs(lp.lam);

    xy.x *= cos((C0 + lp.lam * (C1 + lp.lam * lp.lam * C3)) *

                (lp.phi * (D1 + D5 * lp.phi * lp.phi * lp.phi * lp.phi)));

    return xy;

}

PJ *PJ_PROJECTION(denoy) {

    P->es = 0.0;

    P->fwd = denoy_s_forward;

    return P;

}

#undef C0

#undef C1

#undef C3

#undef D1

#undef D5