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

« back to all changes in this revision

Viewing changes to src/PJ_lsat.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:
66
66
        }
67
67
        if (l) {
68
68
                sp = sin(lp.phi);
69
 
                phidp = aasin((P->one_es * P->ca * sp - P->sa * cos(lp.phi) * 
 
69
                phidp = aasin(P->ctx,(P->one_es * P->ca * sp - P->sa * cos(lp.phi) * 
70
70
                        sin(lamt)) / sqrt(1. - P->es * sp * sp));
71
71
                tanph = log(tan(FORTPI + .5 * phidp));
72
72
                sd = sin(lamdp);
116
116
        lamt -= HALFPI * (1. - scl) * sl;
117
117
        lp.lam = lamt - P->p22 * lamdp;
118
118
        if (fabs(P->sa) < TOL)
119
 
            lp.phi = aasin(spp / sqrt(P->one_es * P->one_es + P->es * sppsq));
 
119
            lp.phi = aasin(P->ctx,spp / sqrt(P->one_es * P->one_es + P->es * sppsq));
120
120
        else
121
121
                lp.phi = atan((tan(lamdp) * cos(lamt) - P->ca * sin(lamt)) /
122
122
                        (P->one_es * P->sa));
127
127
    int land, path;
128
128
    double lam, alf, esc, ess;
129
129
 
130
 
        land = pj_param(P->params, "ilsat").i;
 
130
        land = pj_param(P->ctx, P->params, "ilsat").i;
131
131
        if (land <= 0 || land > 5) E_ERROR(-28);
132
 
        path = pj_param(P->params, "ipath").i;
 
132
        path = pj_param(P->ctx, P->params, "ipath").i;
133
133
        if (path <= 0 || path > (land <= 3 ? 251 : 233)) E_ERROR(-29);
134
134
        if (land <= 3) {
135
135
                P->lam0 = DEG_TO_RAD * 128.87 - TWOPI / 251. * path;