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

« back to all changes in this revision

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

#ifndef lint

static const char SCCSID[]="@(#)geod_for.c 4.6 95/09/23 GIE REL";

#endif

# include "projects.h"

# include "geodesic.h"

# define MERI_TOL 1e-9

static double

th1,costh1,sinth1,sina12,cosa12,M,N,c1,c2,D,P,s1;

static int

merid, signS;

void

geod_pre(void) {

al12 = adjlon(al12); /* reduce to +- 0-PI */

signS = fabs(al12) > HALFPI ? 1 : 0;

th1 = ellipse ? atan(onef * tan(phi1)) : phi1;

costh1 = cos(th1);

sinth1 = sin(th1);

if ((merid = fabs(sina12 = sin(al12)) < MERI_TOL)) {

sina12 = 0.;

cosa12 = fabs(al12) < HALFPI ? 1. : -1.;

M = 0.;

} else {

cosa12 = cos(al12);

M = costh1 * sina12;

}

N = costh1 * cosa12;

if (ellipse) {

if (merid) {

c1 = 0.;

c2 = f4;

D = 1. - c2;

D *= D;

P = c2 / D;

} else {

c1 = f * M;

c2 = f4 * (1. - M * M);

D = (1. - c2)*(1. - c2 - c1 * M);

P = (1. + .5 * c1 * M) * c2 / D;

}

if (merid) s1 = HALFPI - th1;

else {

s1 = (fabs(M) >= 1.) ? 0. : acos(M);

s1 = sinth1 / sin(s1);

s1 = (fabs(s1) >= 1.) ? 0. : acos(s1);

}

void

geod_for(void) {

double d,sind,u,V,X,ds,cosds,sinds,ss,de;

if (ellipse) {

d = S / (D * a);

if (signS) d = -d;

u = 2. * (s1 - d);

V = cos(u + d);

X = c2 * c2 * (sind = sin(d)) * cos(d) * (2. * V * V - 1.);

ds = d + X - 2. * P * V * (1. - 2. * P * cos(u)) * sind;

ss = s1 + s1 - ds;

} else {

ds = S / a;

if (signS) ds = - ds;

}

cosds = cos(ds);

sinds = sin(ds);

if (signS) sinds = - sinds;

al21 = N * cosds - sinth1 * sinds;

if (merid) {

phi2 = atan( tan(HALFPI + s1 - ds) / onef);

if (al21 > 0.) {

al21 = PI;

if (signS)

de = PI;

else {

phi2 = - phi2;

de = 0.;

}

} else {

al21 = 0.;

if (signS) {

phi2 = - phi2;

de = 0;

} else

de = PI;

}

} else {

al21 = atan(M / al21);

if (al21 > 0)

al21 += PI;

if (al12 < 0.)

al21 -= PI;

al21 = adjlon(al21);

phi2 = atan(-(sinth1 * cosds + N * sinds) * sin(al21) /

(ellipse ? onef * M : M));

de = atan2(sinds * sina12 ,

(costh1 * cosds - sinth1 * sinds * cosa12));

if (ellipse)

if (signS)

de += c1 * ((1. - c2) * ds +

100

c2 * sinds * cos(ss));

101

else

102

de -= c1 * ((1. - c2) * ds -

103

c2 * sinds * cos(ss));

104

}

105

lam2 = adjlon( lam1 + de );

106

}

Older »