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

« back to all changes in this revision

Viewing changes to src/PJ_eck4.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_eck4.c

#ifndef lint

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

#endif

#define PJ_LIB__

#include <projects.h>

PROJ_HEAD(eck4, "Eckert IV") "\n\tPCyl, Sph.";

#define C_x .42223820031577120149

#define C_y 1.32650042817700232218

#define RC_y .75386330736002178205

#define C_p 3.57079632679489661922

#define RC_p .28004957675577868795

#define EPS 1e-7

#define NITER 6

FORWARD(s_forward); /* spheroid */

double p, V, s, c;

int i;

p = C_p * sin(lp.phi);

V = lp.phi * lp.phi;

lp.phi *= 0.895168 + V * ( 0.0218849 + V * 0.00826809 );

for (i = NITER; i ; --i) {

c = cos(lp.phi);

s = sin(lp.phi);

lp.phi -= V = (lp.phi + s * (c + 2.) - p) /

(1. + c * (c + 2.) - s * s);

if (fabs(V) < EPS)

break;

}

if (!i) {

xy.x = C_x * lp.lam;

xy.y = lp.phi < 0. ? -C_y : C_y;

} else {

xy.x = C_x * lp.lam * (1. + cos(lp.phi));

xy.y = C_y * sin(lp.phi);

}

return (xy);

}

INVERSE(s_inverse); /* spheroid */

double c;

lp.phi = aasin(xy.y / C_y);

lp.lam = xy.x / (C_x * (1. + (c = cos(lp.phi))));

lp.phi = aasin((lp.phi + sin(lp.phi) * (c + 2.)) / C_p);

return (lp);

}

FREEUP; if (P) pj_dalloc(P); }

ENTRY0(eck4); P->es = 0.; P->inv = s_inverse; P->fwd = s_forward; ENDENTRY(P)

Older »