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

« back to all changes in this revision

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

Show diffs side-by-side

added added

removed removed

Lines of Context:
 
1
#ifndef lint
 
2
static const char SCCSID[]="@(#)PJ_lask.c       4.1     94/02/15        GIE     REL";
 
3
#endif
 
4
#define PROJ_PARMS__
 
5
#define PJ_LIB__
 
6
#include        <projects.h>
 
7
PROJ_HEAD(lask, "Laskowski") "\n\tMisc Sph, no inv.";
 
8
#define a10      0.975534
 
9
#define a12     -0.119161
 
10
#define a32     -0.0143059
 
11
#define a14     -0.0547009
 
12
#define b01      1.00384
 
13
#define b21      0.0802894
 
14
#define b03      0.0998909
 
15
#define b41      0.000199025
 
16
#define b23     -0.0285500
 
17
#define b05     -0.0491032
 
18
FORWARD(s_forward); /* sphere */
 
19
        double l2, p2;
 
20
 
 
21
        l2 = lp.lam * lp.lam;
 
22
        p2 = lp.phi * lp.phi;
 
23
        xy.x = lp.lam * (a10 + p2 * (a12 + l2 * a32 + p2 * a14));
 
24
        xy.y = lp.phi * (b01 + l2 * (b21 + p2 * b23 + l2 * b41) +
 
25
                p2 * (b03 + p2 * b05));
 
26
        return (xy);
 
27
}
 
28
FREEUP; if (P) pj_dalloc(P); }
 
29
ENTRY0(lask) P->fwd = s_forward; P->inv = 0; P->es = 0.; ENDENTRY(P)