~ubuntu-branches/ubuntu/wily/proj/wily

« back to all changes in this revision

Viewing changes to src/PJ_somerc.c

  • Committer: Package Import Robot
  • Author(s): Francesco Paolo Lovergine, Jerome Villeneuve Larouche, Francesco Paolo Lovergine
  • Date: 2013-11-25 15:11:25 UTC
  • mfrom: (1.2.7)
  • Revision ID: package-import@ubuntu.com-20131125151125-mvcw144wvgep1hev
Tags: 4.8.0-1
[ Jerome Villeneuve Larouche ]
* New upstream release
* Modified libproj-dev.install to remove nad_list.h and projects.h
* Modified proj-bin.install to remove nad2nad
* Modified proj-bin.manpages to remove nad2nad man
* Added symbols for libproj0

[ Francesco Paolo Lovergine ]
* Properly merged with current git master and sid version 4.7.0-2.
* Removed proj transitional package and obsolete conflicts/replaces against
  series 4.6 and older.

Show diffs side-by-side

added added

removed removed

Lines of Context:
14
14
                + P->K)) - HALFPI;
15
15
        lamp = P->c * lp.lam;
16
16
        cp = cos(phip);
17
 
        phipp = aasin(P->cosp0 * sin(phip) - P->sinp0 * cp * cos(lamp));
18
 
        lampp = aasin(cp * sin(lamp) / cos(phipp));
 
17
        phipp = aasin(P->ctx,P->cosp0 * sin(phip) - P->sinp0 * cp * cos(lamp));
 
18
        lampp = aasin(P->ctx,cp * sin(lamp) / cos(phipp));
19
19
        xy.x = P->kR * lampp;
20
20
        xy.y = P->kR * log(tan(FORTPI + 0.5 * phipp));
21
21
        return (xy);
27
27
        phipp = 2. * (atan(exp(xy.y / P->kR)) - FORTPI);
28
28
        lampp = xy.x / P->kR;
29
29
        cp = cos(phipp);
30
 
        phip = aasin(P->cosp0 * sin(phipp) + P->sinp0 * cp * cos(lampp));
31
 
        lamp = aasin(cp * sin(lampp) / cos(phip));
 
30
        phip = aasin(P->ctx,P->cosp0 * sin(phipp) + P->sinp0 * cp * cos(lampp));
 
31
        lamp = aasin(P->ctx,cp * sin(lampp) / cos(phip));
32
32
        con = (P->K - log(tan(FORTPI + 0.5 * phip)))/P->c;
33
33
        for (i = NITER; i ; --i) {
34
34
                esp = P->e * sin(phip);
55
55
        cp *= cp;
56
56
        P->c = sqrt(1 + P->es * cp * cp * P->rone_es);
57
57
        sp = sin(P->phi0);
58
 
        P->cosp0 = cos( phip0 = aasin(P->sinp0 = sp / P->c) );
 
58
        P->cosp0 = cos( phip0 = aasin(P->ctx, P->sinp0 = sp / P->c) );
59
59
        sp *= P->e;
60
60
        P->K = log(tan(FORTPI + 0.5 * phip0)) - P->c * (
61
61
                log(tan(FORTPI + 0.5 * P->phi0)) - P->hlf_e *