static XY e_forward (LP lp, PJ *P) { /* Ellipsoidal, forward */
XY xy = {0.0,0.0};
struct pj_opaque *Q = P->opaque;
double cosc, sinc, cosl, k;
lp = pj_gauss(P->ctx, lp, Q->en);
sinc = sin(lp.phi);
cosc = cos(lp.phi);
cosl = cos(lp.lam);
k = P->k0 * Q->R2 / (1. + Q->sinc0 * sinc + Q->cosc0 * cosc * cosl);
xy.x = k * cosc * sin(lp.lam);
xy.y = k * (Q->cosc0 * sinc - Q->sinc0 * cosc * cosl);
return xy;
}
static PJ_XY sterea_e_forward (PJ_LP lp, PJ *P) { /* Ellipsoidal, forward */
PJ_XY xy = {0.0,0.0};
struct pj_opaque *Q = static_cast<struct pj_opaque*>(P->opaque);
double cosc, sinc, cosl, k;
lp = pj_gauss(P->ctx, lp, Q->en);
sinc = sin(lp.phi);
cosc = cos(lp.phi);
cosl = cos(lp.lam);
const double denom = 1. + Q->sinc0 * sinc + Q->cosc0 * cosc * cosl;
if( denom == 0.0 ) {
proj_errno_set(P, PJD_ERR_TOLERANCE_CONDITION);
return proj_coord_error().xy;
}
k = P->k0 * Q->R2 / denom;
xy.x = k * cosc * sin(lp.lam);
xy.y = k * (Q->cosc0 * sinc - Q->sinc0 * cosc * cosl);
return xy;
}
double cosc0, sinc0; \
double R2; \
void *en;
#define PJ_LIB__
#include <projects.h>
PROJ_HEAD(sterea, "Oblique Stereographic Alternative")
"\n\tAzimuthal, Sph&Ell";
# define DEL_TOL 1.e-14
# define MAX_ITER 10
FORWARD(e_forward); /* ellipsoid */
double cosc, sinc, cosl, k;
lp = pj_gauss(P->ctx, lp, P->en);
sinc = sin(lp.phi);
cosc = cos(lp.phi);
cosl = cos(lp.lam);
k = P->k0 * P->R2 / (1. + P->sinc0 * sinc + P->cosc0 * cosc * cosl);
xy.x = k * cosc * sin(lp.lam);
xy.y = k * (P->cosc0 * sinc - P->sinc0 * cosc * cosl);
return (xy);
}
INVERSE(e_inverse); /* ellipsoid */
double rho, c, sinc, cosc;
xy.x /= P->k0;
xy.y /= P->k0;
if((rho = hypot(xy.x, xy.y))) {
c = 2. * atan2(rho, P->R2);
double cosc0, sinc0; \
double R2; \
void *en;
#define PJ_LIB__
#include <lib_proj.h>
PROJ_HEAD(sterea, "Oblique Stereographic Alternative")
"\n\tAzimuthal, Sph&Ell";
# define DEL_TOL 1.e-14
# define MAX_ITER 10
FORWARD(e_forward); /* ellipsoid */
double cosc, sinc, cosl, k;
lp = pj_gauss(lp, P->en);
sinc = sin(lp.phi);
cosc = cos(lp.phi);
cosl = cos(lp.lam);
k = P->k0 * P->R2 / (1. + P->sinc0 * sinc + P->cosc0 * cosc * cosl);
xy.x = k * cosc * sin(lp.lam);
xy.y = k * (P->cosc0 * sinc - P->sinc0 * cosc * cosl);
return (xy);
}
INVERSE(e_inverse); /* ellipsoid */
double rho, c, sinc, cosc;
xy.x /= P->k0;
xy.y /= P->k0;
if((rho = hypot(xy.x, xy.y))) {
c = 2. * atan2(rho, P->R2);
