~ubuntu-branches/ubuntu/vivid/proj/vivid

« back to all changes in this revision

Viewing changes to src/PJ_poly.c

Committer: Bazaar Package Importer
Author(s): Peter S Galbraith
Date: 2002-01-11 10:27:12 UTC
Revision ID: james.westby@ubuntu.com-20020111102712-ayi18r8y2eesv0y9

Tags: upstream-4.4.5

Import upstream version 4.4.5

files added:

AUTHORS

COPYING

ChangeLog

INSTALL

Makefile.am

Makefile.in

NEWS

README

aclocal.m4

config.guess

config.sub

configure

configure.in

install-sh

ltconfig

ltmain.sh

man/Makefile.am

man/Makefile.in

man/man1

man/man1/Makefile.am

man/man1/Makefile.in

man/man1/cs2cs.1

man/man1/geod.1

man/man1/nad2nad.1

man/man1/proj.1

man/man3

man/man3/Makefile.am

man/man3/Makefile.in

man/man3/pj_init.3

missing

mkinstalldirs

nad/GL27

nad/Makefile.am

nad/Makefile.in

nad/README

nad/epsg

nad/nad.lst

nad/nad27

nad/nad83

nad/pj_out27.dist

nad/pj_out83.dist

nad/proj_def.dat

nad/test27

nad/test83

nad/world

src/Makefile.am

src/Makefile.in

src/PJ_aea.c

src/PJ_aeqd.c

src/PJ_airy.c

src/PJ_aitoff.c

src/PJ_august.c

src/PJ_bacon.c

src/PJ_bipc.c

src/PJ_boggs.c

src/PJ_bonne.c

src/PJ_cass.c

src/PJ_cc.c

src/PJ_cea.c

src/PJ_chamb.c

src/PJ_collg.c

src/PJ_crast.c

src/PJ_denoy.c

src/PJ_eck1.c

src/PJ_eck2.c

src/PJ_eck3.c

src/PJ_eck4.c

src/PJ_eck5.c

src/PJ_eqc.c

src/PJ_eqdc.c

src/PJ_fahey.c

src/PJ_fouc_s.c

src/PJ_gall.c

src/PJ_gins8.c

src/PJ_gn_sinu.c

src/PJ_gnom.c

src/PJ_goode.c

src/PJ_hammer.c

src/PJ_hatano.c

src/PJ_imw_p.c

src/PJ_labrd.c

src/PJ_laea.c

src/PJ_lagrng.c

src/PJ_larr.c

src/PJ_lask.c

src/PJ_lcc.c

src/PJ_loxim.c

src/PJ_lsat.c

src/PJ_mbt_fps.c

src/PJ_mbtfpp.c

src/PJ_mbtfpq.c

src/PJ_merc.c

src/PJ_mill.c

src/PJ_mod_ster.c

src/PJ_moll.c

src/PJ_mpoly.c

src/PJ_nell.c

src/PJ_nell_h.c

src/PJ_nocol.c

src/PJ_nsper.c

src/PJ_nzmg.c

src/PJ_ob_tran.c

src/PJ_ocea.c

src/PJ_oea.c

src/PJ_omerc.c

src/PJ_ortho.c

src/PJ_poly.c

src/PJ_putp2.c

src/PJ_putp3.c

src/PJ_putp4p.c

src/PJ_putp5.c

src/PJ_putp6.c

src/PJ_robin.c

src/PJ_rpoly.c

src/PJ_sconics.c

src/PJ_somerc.c

src/PJ_stere.c

src/PJ_sts.c

src/PJ_tcc.c

src/PJ_tcea.c

src/PJ_tmerc.c

src/PJ_tpeqd.c

src/PJ_urm5.c

src/PJ_urmfps.c

src/PJ_vandg.c

src/PJ_vandg2.c

src/PJ_vandg4.c

src/PJ_wag2.c

src/PJ_wag3.c

src/PJ_wag7.c

src/PJ_wink1.c

src/PJ_wink2.c

src/aasincos.c

src/adjlon.c

src/bch2bps.c

src/bchgen.c

src/biveval.c

src/cs2cs.c

src/dmstor.c

src/emess.c

src/emess.h

src/gen_cheb.c

src/geocent.c

src/geocent.h

src/geod.c

src/geod_for.c

src/geod_inv.c

src/geod_set.c

src/geodesic.h

src/makefile.vc

src/mk_cheby.c

src/nad2bin.c

src/nad2nad.c

src/nad_cvt.c

src/nad_init.c

src/nad_intr.c

src/nad_list.h

src/p_series.c

src/pj_apply_gridshift.c

src/pj_auth.c

src/pj_datum_set.c

src/pj_datums.c

src/pj_deriv.c

src/pj_ell_set.c

src/pj_ellps.c

src/pj_errno.c

src/pj_factors.c

src/pj_fwd.c

src/pj_init.c

src/pj_inv.c

src/pj_latlong.c

src/pj_list.c

src/pj_list.h

src/pj_malloc.c

src/pj_mlfn.c

src/pj_msfn.c

src/pj_open_lib.c

src/pj_param.c

src/pj_phi2.c

src/pj_pr_list.c

src/pj_qsfn.c

src/pj_release.c

src/pj_strerrno.c

src/pj_transform.c

src/pj_tsfn.c

src/pj_units.c

src/pj_utils.c

src/pj_zpoly1.c

src/proj.c

src/proj.def

src/proj_api.h

src/proj_config.h.in

src/projects.h

src/rtodms.c

src/stamp-h.in

src/strtod.c

src/vector1.c

Show diffs side-by-side

added added

removed removed

src/PJ_poly.c

#ifndef lint

static const char SCCSID[]="@(#)PJ_poly.c 4.1 94/02/15 GIE REL";

#endif

#define PROJ_PARMS__ \

double ml0; \

double *en;

#define PJ_LIB__

#include <projects.h>

PROJ_HEAD(poly, "Polyconic (American)")

"\n\tConic, Sph&Ell";

#define TOL 1e-10

#define CONV 1e-10

#define N_ITER 10

#define I_ITER 20

#define ITOL 1.e-12

FORWARD(e_forward); /* ellipsoid */

double ms, sp, cp;

if (fabs(lp.phi) <= TOL) { xy.x = lp.lam; xy.y = -P->ml0; }

else {

sp = sin(lp.phi);

ms = fabs(cp = cos(lp.phi)) > TOL ? pj_msfn(sp, cp, P->es) / sp : 0.;

xy.x = ms * sin(lp.lam *= sp);

xy.y = (pj_mlfn(lp.phi, sp, cp, P->en) - P->ml0) + ms * (1. - cos(lp.lam));

}

return (xy);

}

FORWARD(s_forward); /* spheroid */

double cot, E;

if (fabs(lp.phi) <= TOL) { xy.x = lp.lam; xy.y = P->ml0; }

else {

cot = 1. / tan(lp.phi);

xy.x = sin(E = lp.lam * sin(lp.phi)) * cot;

xy.y = lp.phi - P->phi0 + cot * (1. - cos(E));

}

return (xy);

}

INVERSE(e_inverse); /* ellipsoid */

xy.y += P->ml0;

if (fabs(xy.y) <= TOL) { lp.lam = xy.x; lp.phi = 0.; }

else {

double r, c, sp, cp, s2ph, ml, mlb, mlp, dPhi;

int i;

r = xy.y * xy.y + xy.x * xy.x;

for (lp.phi = xy.y, i = I_ITER; i ; --i) {

sp = sin(lp.phi);

s2ph = sp * ( cp = cos(lp.phi));

if (fabs(cp) < ITOL)

I_ERROR;

c = sp * (mlp = sqrt(1. - P->es * sp * sp)) / cp;

ml = pj_mlfn(lp.phi, sp, cp, P->en);

mlb = ml * ml + r;

mlp = P->one_es / (mlp * mlp * mlp);

lp.phi += ( dPhi =

( ml + ml + c * mlb - 2. * xy.y * (c * ml + 1.) ) / (

P->es * s2ph * (mlb - 2. * xy.y * ml) / c +

2.* (xy.y - ml) * (c * mlp - 1. / s2ph) - mlp - mlp ));

if (fabs(dPhi) <= ITOL)

break;

}

if (!i)

I_ERROR;

c = sin(lp.phi);

lp.lam = asin(xy.x * tan(lp.phi) * sqrt(1. - P->es * c * c)) / sin(lp.phi);

}

return (lp);

}

INVERSE(s_inverse); /* spheroid */

double B, dphi, tp;

int i;

if (fabs(xy.y = P->phi0 + xy.y) <= TOL) { lp.lam = xy.x; lp.phi = 0.; }

else {

lp.phi = xy.y;

B = xy.x * xy.x + xy.y * xy.y;

i = N_ITER;

do {

tp = tan(lp.phi);

lp.phi -= (dphi = (xy.y * (lp.phi * tp + 1.) - lp.phi -

.5 * ( lp.phi * lp.phi + B) * tp) /

((lp.phi - xy.y) / tp - 1.));

} while (fabs(dphi) > CONV && --i);

if (! i) I_ERROR;

lp.lam = asin(xy.x * tan(lp.phi)) / sin(lp.phi);

}

return (lp);

}

FREEUP; if (P) { if (P->en) pj_dalloc(P->en); pj_dalloc(P); } }

ENTRY1(poly, en)

if (P->es) {

if (!(P->en = pj_enfn(P->es))) E_ERROR_0;

P->ml0 = pj_mlfn(P->phi0, sin(P->phi0), cos(P->phi0), P->en);

P->inv = e_inverse;

P->fwd = e_forward;

} else {

P->ml0 = -P->phi0;

P->inv = s_inverse;

100

P->fwd = s_forward;

101

}

102

ENDENTRY(P)

Older »