/Proj4/PJ_tcc.c

 1#ifndef lint
 2static const char SCCSID[]="@(#)PJ_tcc.c	4.1	94/02/15	GIE	REL";
 3#endif
 4#define PROJ_PARMS__ \
 5	double ap;
 6#define EPS10 1.e-10
 7#define PJ_LIB__
 8#include	"projects.h"
 9PROJ_HEAD(tcc, "Transverse Central Cylindrical") "\n\tCyl, Sph, no inv.";
10FORWARD(s_forward); /* spheroid */
11	double b, bt;
12
13	b = cos(lp.phi) * sin(lp.lam);
14	if ((bt = 1. - b * b) < EPS10) F_ERROR;
15	xy.x = b / sqrt(bt);
16	xy.y = atan2(tan(lp.phi) , cos(lp.lam));
17	return (xy);
18}
19FREEUP; if (P) pj_dalloc(P); }
20ENTRY0(tcc) P->es = 0.; P->fwd = s_forward; ENDENTRY(P)