2
static const char SCCSID[]="@(#)PJ_bacon.c 4.1 94/02/15 GIE REL";
4
# define HLFPI2 2.46740110027233965467
11
PROJ_HEAD(apian, "Apian Globular I") "\n\tMisc Sph, no inv.";
12
PROJ_HEAD(ortel, "Ortelius Oval") "\n\tMisc Sph, no inv.";
13
PROJ_HEAD(bacon, "Bacon Globular") "\n\tMisc Sph, no inv.";
14
FORWARD(s_forward); /* spheroid */
17
xy.y = P->bacn ? HALFPI * sin(lp.phi) : lp.phi;
18
if ((ax = fabs(lp.lam)) >= EPS) {
19
if (P->ortl && ax >= HALFPI)
20
xy.x = sqrt(HLFPI2 - lp.phi * lp.phi + EPS) + ax - HALFPI;
22
f = 0.5 * (HLFPI2 / ax + ax);
23
xy.x = ax - f + sqrt(f * f - xy.y * xy.y);
25
if (lp.lam < 0.) xy.x = - xy.x;
30
FREEUP; if (P) pj_dalloc(P); }
34
P->es = 0.; P->fwd = s_forward;
37
P->bacn = P->ortl = 0;
38
P->es = 0.; P->fwd = s_forward;
43
P->es = 0.; P->fwd = s_forward;