svn-gvsig-desktop / tags / v1_10_0_Build_1258 / libraries / libjni-proj4 / src / PJ_lask.c
History | View | Annotate | Download (783 Bytes)
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) |