| /dports/science/cdcl/dcl-5.4.8-C/src/math1/maplib/ |
| H A D | mpnwtn.c | 43 if (*ylat > 0.f) { in mpnwtn_()
44 ylat2 = *ylat - 1e-4f; in mpnwtn_()
46 ylat2 = *ylat + 1e-4f; in mpnwtn_()
48 (*func)(xlon, ylat, &x0, &y0); in mpnwtn_()
49 (*func)(&xlon2, ylat, &x1, &y1); in mpnwtn_()
54 xb = (x2 - x0) / (ylat2 - *ylat); in mpnwtn_()
56 yb = (y2 - y0) / (ylat2 - *ylat); in mpnwtn_()
87 *ylat += cy; in mpnwtn_()
90 *ylat = atan(1.f) * 4.f * .5f - .01f; in mpnwtn_()
94 *ylat = atan(1.f) * 4.f * -.5f + .01f; in mpnwtn_()
[all …]
|
| H A D | mpfplc.c | 50 mpxplc_(&xx, ylat, x, y); in mpfplc_0_()
78 *ylat = ((j - .5f) / 6 * 180.f - 90.f) * .017453292519943292f; in mpfplc_0_()
85 mpnwtn_(x, y, xlon, ylat, (S_fp)mpxplc_, &c_true); in mpfplc_0_()
89 /* Subroutine */ int mpfplc_(real *xlon, real *ylat, real *x, real *y) in mpfplc_() argument
91 return mpfplc_0_(0, xlon, ylat, x, y); in mpfplc_()
96 return mpfplc_0_(1, xlon, ylat, x, y); in mpiplc_()
107 if (abs(*ylat) < 1e-4f) { in mpxplc_()
109 *y = *ylat; in mpxplc_()
111 th = *xlon * sin(*ylat); in mpxplc_()
112 tanz = 1.f / tan(*ylat); in mpxplc_()
[all …]
|
| H A D | mpfmwl.c | 20 /* Subroutine */ int mpfmwl_0_(int n__, real *xlon, real *ylat, real *x, real in mpfmwl_0_() argument
37 *x = xmplon_(xlon) * 2.8284272000000001f * cos(*ylat) / pi; in mpfmwl_0_()
38 *y = sin(*ylat) * 1.4142136f; in mpfmwl_0_()
45 *ylat = asin(*y / 1.4142136f); in mpfmwl_0_()
46 *xlon = *x / 1.4142136f / cos(*ylat) * pi / 2; in mpfmwl_0_()
52 *ylat = *y / 1.4142136f * pi / 2; in mpfmwl_0_()
56 *ylat = rna; in mpfmwl_0_()
60 /* Subroutine */ int mpfmwl_(real *xlon, real *ylat, real *x, real *y) in mpfmwl_() argument
62 return mpfmwl_0_(0, xlon, ylat, x, y); in mpfmwl_()
65 /* Subroutine */ int mpimwl_(real *x, real *y, real *xlon, real *ylat) in mpimwl_() argument
[all …]
|
| H A D | mpfotg.c | 20 /* Subroutine */ int mpfotg_0_(int n__, real *xlon, real *ylat, real *x, real in mpfotg_0_() argument
43 r__ = cos(*ylat) / (1 - c__ * sin(*ylat)); in mpfotg_0_()
45 r__ = cos(*ylat); in mpfotg_0_()
58 *ylat = acos(r2) - atan(c__ * r1); in mpfotg_0_()
71 *ylat = acos(r__); in mpfotg_0_()
75 *ylat = acos(r__); in mpfotg_0_()
81 *ylat = rna; in mpfotg_0_()
93 /* Subroutine */ int mpfotg_(real *xlon, real *ylat, real *x, real *y) in mpfotg_() argument
95 return mpfotg_0_(0, xlon, ylat, x, y, (real *)0); in mpfotg_()
98 /* Subroutine */ int mpiotg_(real *x, real *y, real *xlon, real *ylat) in mpiotg_() argument
[all …]
|
| H A D | mpfglb.c | 20 /* Subroutine */ int mpfglb_0_(int n__, real *xlon, real *ylat, real *x, real in mpfglb_0_() argument
41 *y = sin(*ylat); in mpfglb_0_()
42 *x = cos(*ylat) * xx / (atan(1.f) * 2.f); in mpfglb_0_()
54 *ylat = asin(*y); in mpfglb_0_()
55 coslat = cos(*ylat); in mpfglb_0_()
65 *ylat = rna; in mpfglb_0_()
69 /* Subroutine */ int mpfglb_(real *xlon, real *ylat, real *x, real *y) in mpfglb_() argument
71 return mpfglb_0_(0, xlon, ylat, x, y); in mpfglb_()
74 /* Subroutine */ int mpiglb_(real *x, real *y, real *xlon, real *ylat) in mpiglb_() argument
76 return mpfglb_0_(1, xlon, ylat, x, y); in mpiglb_()
|
| H A D | mpfrbs.c | 41 r__1 = -pi * .5f + 1e-5f, r__2 = min(r__3,*ylat); in mpfrbs_0_()
89 r__1 = *ylat; in mpfrbs_0_()
91 r__2 = *ylat, r__2 *= r__2; in mpfrbs_0_()
93 r__3 = *ylat, r__3 *= r__3; in mpfrbs_0_()
95 r__4 = *ylat, r__4 *= r__4, r__4 *= r__4; in mpfrbs_0_()
97 r__5 = *ylat, r__5 *= r__5, r__6 = r__5, r__5 *= r__5; in mpfrbs_0_()
99 r__7 = *ylat, r__7 *= r__7, r__7 *= r__7; in mpfrbs_0_()
103 if (abs(*xlon) <= pi + .1f && abs(*ylat) <= pi / 2) { in mpfrbs_0_()
108 *ylat = rna; in mpfrbs_0_()
114 return mpfrbs_0_(0, xlon, ylat, x, y); in mpfrbs_()
[all …]
|
| H A D | mpfmwd.c | 20 /* Subroutine */ int mpfmwd_0_(int n__, real *xlon, real *ylat, real *x, real in mpfmwd_0_() argument
40 alpha = *ylat; in mpfmwd_0_()
41 mpznwt_((U_fp)mpzmwd_, ylat, &alpha); in mpfmwd_0_()
53 *ylat = asin((alpha * 2 + sin(alpha * 2)) / pi); in mpfmwd_0_()
58 *ylat = *y / 1.4142136f / 2 * pi; in mpfmwd_0_()
63 *ylat = rna; in mpfmwd_0_()
67 /* Subroutine */ int mpfmwd_(real *xlon, real *ylat, real *x, real *y) in mpfmwd_() argument
69 return mpfmwd_0_(0, xlon, ylat, x, y); in mpfmwd_()
72 /* Subroutine */ int mpimwd_(real *x, real *y, real *xlon, real *ylat) in mpimwd_() argument
74 return mpfmwd_0_(1, xlon, ylat, x, y); in mpimwd_()
|
| H A D | mpfcyl.c | 20 /* Subroutine */ int mpfcyl_0_(int n__, real *xlon, real *ylat, real *x, real in mpfcyl_0_() argument
34 *y = *ylat; in mpfcyl_0_()
47 *ylat = *y; in mpfcyl_0_()
48 if (abs(*xlon) <= pi && abs(*ylat) <= pi / 2) { in mpfcyl_0_()
53 *ylat = rna; in mpfcyl_0_()
57 /* Subroutine */ int mpfcyl_(real *xlon, real *ylat, real *x, real *y) in mpfcyl_() argument
59 return mpfcyl_0_(0, xlon, ylat, x, y); in mpfcyl_()
62 /* Subroutine */ int mpicyl_(real *x, real *y, real *xlon, real *ylat) in mpicyl_() argument
64 return mpfcyl_0_(1, xlon, ylat, x, y); in mpicyl_()
|
| H A D | mpfcya.c | 20 /* Subroutine */ int mpfcya_0_(int n__, real *xlon, real *ylat, real *x, real in mpfcya_0_() argument
41 *y = sin(*ylat); in mpfcya_0_()
51 *ylat = asin(*y); in mpfcya_0_()
52 if (abs(*xlon) <= pi && abs(*ylat) <= pi / 2) { in mpfcya_0_()
57 *ylat = rna; in mpfcya_0_()
61 /* Subroutine */ int mpfcya_(real *xlon, real *ylat, real *x, real *y) in mpfcya_() argument
63 return mpfcya_0_(0, xlon, ylat, x, y); in mpfcya_()
66 /* Subroutine */ int mpicya_(real *x, real *y, real *xlon, real *ylat) in mpicya_() argument
68 return mpfcya_0_(1, xlon, ylat, x, y); in mpicya_()
|
| H A D | mpfsin.c | 20 /* Subroutine */ int mpfsin_0_(int n__, real *xlon, real *ylat, real *x, real in mpfsin_0_() argument
45 r__1 = -pi * .5f, r__2 = min(r__3,*ylat); in mpfsin_0_()
52 *ylat = *y; in mpfsin_0_()
59 if (abs(*xlon) <= pi + .1f && abs(*ylat) <= pi / 2) { in mpfsin_0_()
64 *ylat = rna; in mpfsin_0_()
68 /* Subroutine */ int mpfsin_(real *xlon, real *ylat, real *x, real *y) in mpfsin_() argument
70 return mpfsin_0_(0, xlon, ylat, x, y); in mpfsin_()
73 /* Subroutine */ int mpisin_(real *x, real *y, real *xlon, real *ylat) in mpisin_() argument
75 return mpfsin_0_(1, xlon, ylat, x, y); in mpisin_()
|
| H A D | mpfcyc.c | 20 /* Subroutine */ int mpfcyc_0_(int n__, real *xlon, real *ylat, real *x, real in mpfcyc_0_() argument
44 r__1 = -pi * .5f + 1e-5f, r__2 = min(r__3,*ylat); in mpfcyc_0_()
53 *ylat = atan(*y); in mpfcyc_0_()
54 if (abs(*xlon) <= pi && abs(*ylat) <= pi / 2) { in mpfcyc_0_()
59 *ylat = rna; in mpfcyc_0_()
63 /* Subroutine */ int mpfcyc_(real *xlon, real *ylat, real *x, real *y) in mpfcyc_() argument
65 return mpfcyc_0_(0, xlon, ylat, x, y); in mpfcyc_()
68 /* Subroutine */ int mpicyc_(real *x, real *y, real *xlon, real *ylat) in mpicyc_() argument
70 return mpfcyc_0_(1, xlon, ylat, x, y); in mpicyc_()
|
| H A D | mpfktd.c | 20 /* Subroutine */ int mpfktd_0_(int n__, real *xlon, real *ylat, real *x, real in mpfktd_0_() argument
40 alpha = *ylat * .9f; in mpfktd_0_()
41 mpznwt_((U_fp)mpzktd_, ylat, &alpha); in mpfktd_0_()
53 *ylat = asin((alpha * 2 + sin(alpha * 2)) / 2.960421f); in mpfktd_0_()
59 *ylat = rna; in mpfktd_0_()
63 /* Subroutine */ int mpfktd_(real *xlon, real *ylat, real *x, real *y) in mpfktd_() argument
65 return mpfktd_0_(0, xlon, ylat, x, y); in mpfktd_()
68 /* Subroutine */ int mpiktd_(real *x, real *y, real *xlon, real *ylat) in mpiktd_() argument
70 return mpfktd_0_(1, xlon, ylat, x, y); in mpiktd_()
|
| H A D | mpfmil.c | 20 /* Subroutine */ int mpfmil_0_(int n__, real *xlon, real *ylat, real *x, real in mpfmil_0_() argument
45 r__1 = -pi * .5f, r__2 = min(r__3,*ylat); in mpfmil_0_()
54 *ylat = atan(sinh(*y * .8f)) * 1.25f; in mpfmil_0_()
55 if (abs(*xlon) <= pi && abs(*ylat) <= pi / 2) { in mpfmil_0_()
60 *ylat = rna; in mpfmil_0_()
64 /* Subroutine */ int mpfmil_(real *xlon, real *ylat, real *x, real *y) in mpfmil_() argument
66 return mpfmil_0_(0, xlon, ylat, x, y); in mpfmil_()
69 /* Subroutine */ int mpimil_(real *x, real *y, real *xlon, real *ylat) in mpimil_() argument
71 return mpfmil_0_(1, xlon, ylat, x, y); in mpimil_()
|
| H A D | mpfazm.c | 20 /* Subroutine */ int mpfazm_0_(int n__, real *xlon, real *ylat, real *x, real in mpfazm_0_() argument
37 r__ = pi / 2 - *ylat; in mpfazm_0_()
54 *ylat = pi / 2.f; in mpfazm_0_()
58 *ylat = pi / 2 - r__; in mpfazm_0_()
63 *ylat = rna; in mpfazm_0_()
67 /* Subroutine */ int mpfazm_(real *xlon, real *ylat, real *x, real *y) in mpfazm_() argument
69 return mpfazm_0_(0, xlon, ylat, x, y); in mpfazm_()
72 /* Subroutine */ int mpiazm_(real *x, real *y, real *xlon, real *ylat) in mpiazm_() argument
74 return mpfazm_0_(1, xlon, ylat, x, y); in mpiazm_()
|
| H A D | mpfek6.c | 20 /* Subroutine */ int mpfek6_0_(int n__, real *xlon, real *ylat, real *x, real in mpfek6_0_() argument
40 alpha = *ylat; in mpfek6_0_()
41 mpznwt_((U_fp)mpzek6_, ylat, &alpha); in mpfek6_0_()
55 *ylat = phi; in mpfek6_0_()
61 *ylat = rna; in mpfek6_0_()
65 /* Subroutine */ int mpfek6_(real *xlon, real *ylat, real *x, real *y) in mpfek6_() argument
67 return mpfek6_0_(0, xlon, ylat, x, y); in mpfek6_()
70 /* Subroutine */ int mpiek6_(real *x, real *y, real *xlon, real *ylat) in mpiek6_() argument
72 return mpfek6_0_(1, xlon, ylat, x, y); in mpiek6_()
|
| H A D | mpfaza.c | 20 /* Subroutine */ int mpfaza_0_(int n__, real *xlon, real *ylat, real *x, real in mpfaza_0_() argument
37 th = pi / 2 - *ylat; in mpfaza_0_()
55 *ylat = pi / 2.f; in mpfaza_0_()
59 *ylat = pi / 2 - asin(r__ / 2) * 2; in mpfaza_0_()
64 *ylat = rna; in mpfaza_0_()
68 /* Subroutine */ int mpfaza_(real *xlon, real *ylat, real *x, real *y) in mpfaza_() argument
70 return mpfaza_0_(0, xlon, ylat, x, y); in mpfaza_()
73 /* Subroutine */ int mpiaza_(real *x, real *y, real *xlon, real *ylat) in mpiaza_() argument
75 return mpfaza_0_(1, xlon, ylat, x, y); in mpiaza_()
|
| H A D | mpfbon.c | 24 /* Subroutine */ int mpfbon_0_(int n__, real *xlon, real *ylat, real *x, real in mpfbon_0_() argument
42 th = pi / 2 - s * *ylat; in mpfbon_0_()
58 *ylat = s * (pi / 2 - th); in mpfbon_0_()
62 } else if (abs(*ylat) < pi / 2) { in mpfbon_0_()
67 } else if (abs(*ylat) == pi / 2) { in mpfbon_0_()
72 *ylat = rna; in mpfbon_0_()
85 /* Subroutine */ int mpfbon_(real *xlon, real *ylat, real *x, real *y) in mpfbon_() argument
87 return mpfbon_0_(0, xlon, ylat, x, y, (real *)0); in mpfbon_()
90 /* Subroutine */ int mpibon_(real *x, real *y, real *xlon, real *ylat) in mpibon_() argument
92 return mpfbon_0_(1, xlon, ylat, x, y, (real *)0); in mpibon_()
|
| H A D | mpfhmr.c | 20 /* Subroutine */ int mpfhmr_0_(int n__, real *xlon, real *ylat, real *x, real in mpfhmr_0_() argument
40 cphi = cos(*ylat); in mpfhmr_0_()
43 *y = sin(*ylat) * 1.4142136f / cf; in mpfhmr_0_()
61 *ylat = phi; in mpfhmr_0_()
68 *ylat = rna; in mpfhmr_0_()
72 /* Subroutine */ int mpfhmr_(real *xlon, real *ylat, real *x, real *y) in mpfhmr_() argument
74 return mpfhmr_0_(0, xlon, ylat, x, y); in mpfhmr_()
77 /* Subroutine */ int mpihmr_(real *x, real *y, real *xlon, real *ylat) in mpihmr_() argument
79 return mpfhmr_0_(1, xlon, ylat, x, y); in mpihmr_()
|
| H A D | mpfmer.c | 20 /* Subroutine */ int mpfmer_0_(int n__, real *xlon, real *ylat, real *x, real in mpfmer_0_() argument
39 theta = (*ylat + pi / 2) / 2.f; in mpfmer_0_()
55 *ylat = atan(exp(*y)) * 2 - pi / 2; in mpfmer_0_()
60 *ylat = rna; in mpfmer_0_()
64 /* Subroutine */ int mpfmer_(real *xlon, real *ylat, real *x, real *y) in mpfmer_() argument
66 return mpfmer_0_(0, xlon, ylat, x, y); in mpfmer_()
69 /* Subroutine */ int mpimer_(real *x, real *y, real *xlon, real *ylat) in mpimer_() argument
71 return mpfmer_0_(1, xlon, ylat, x, y); in mpimer_()
|
| H A D | mpfgno.c | 20 /* Subroutine */ int mpfgno_0_(int n__, real *xlon, real *ylat, real *x, real in mpfgno_0_() argument
39 r__3 = atan(1.f) * 2.f - *ylat; in mpfgno_0_()
52 *ylat = atan(1.f) * 2.f; in mpfgno_0_()
56 *ylat = atan(1.f) * 2.f - atan(r__); in mpfgno_0_()
62 /* Subroutine */ int mpfgno_(real *xlon, real *ylat, real *x, real *y) in mpfgno_() argument
64 return mpfgno_0_(0, xlon, ylat, x, y); in mpfgno_()
67 /* Subroutine */ int mpigno_(real *x, real *y, real *xlon, real *ylat) in mpigno_() argument
69 return mpfgno_0_(1, xlon, ylat, x, y); in mpigno_()
|
| H A D | mpfcyb.c | 20 /* Subroutine */ int mpfcyb_0_(int n__, real *xlon, real *ylat, real *x, real in mpfcyb_0_() argument
43 *y = c__ * 2.f * tan(*ylat * .5f); in mpfcyb_0_()
50 *ylat = atan(*y * .5f / c__) * 2.f; in mpfcyb_0_()
56 *ylat = rna; in mpfcyb_0_()
65 /* Subroutine */ int mpfcyb_(real *xlon, real *ylat, real *x, real *y) in mpfcyb_() argument
67 return mpfcyb_0_(0, xlon, ylat, x, y, (real *)0); in mpfcyb_()
70 /* Subroutine */ int mpicyb_(real *x, real *y, real *xlon, real *ylat) in mpicyb_() argument
72 return mpfcyb_0_(1, xlon, ylat, x, y, (real *)0); in mpicyb_()
|
| H A D | mpfcon.c | 24 /* Subroutine */ int mpfcon_0_(int n__, real *xlon, real *ylat, real *x, real in mpfcon_0_() argument
42 th = pi / 2 - s * *ylat; in mpfcon_0_()
57 *ylat = s * (pi / 2 - r__ + cc); in mpfcon_0_()
58 if (abs(*xlon) > pi || abs(*ylat) > pi / 2) { in mpfcon_0_()
60 *ylat = rna; in mpfcon_0_()
75 /* Subroutine */ int mpfcon_(real *xlon, real *ylat, real *x, real *y) in mpfcon_() argument
77 return mpfcon_0_(0, xlon, ylat, x, y, (real *)0); in mpfcon_()
80 /* Subroutine */ int mpicon_(real *x, real *y, real *xlon, real *ylat) in mpicon_() argument
82 return mpfcon_0_(1, xlon, ylat, x, y, (real *)0); in mpicon_()
|
| H A D | mpfpst.c | 20 /* Subroutine */ int mpfpst_0_(int n__, real *xlon, real *ylat, real *x, real in mpfpst_0_() argument
36 th = pi / 2 - *ylat; in mpfpst_0_()
49 *ylat = pi / 2 - atan(r__ / 2) * 2; in mpfpst_0_()
58 /* Subroutine */ int mpfpst_(real *xlon, real *ylat, real *x, real *y) in mpfpst_() argument
60 return mpfpst_0_(0, xlon, ylat, x, y); in mpfpst_()
63 /* Subroutine */ int mpipst_(real *x, real *y, real *xlon, real *ylat) in mpipst_() argument
65 return mpfpst_0_(1, xlon, ylat, x, y); in mpipst_()
|
| H A D | mpfcct.c | 24 /* Subroutine */ int mpscct_0_(int n__, real *ylat0, real *xlon, real *ylat, in mpscct_0_() argument
57 r__ = r0 + tan(phi0 - s * *ylat); in mpscct_0_()
70 *ylat = s * (phi0 - atan(r__ - r0)); in mpscct_0_()
73 *ylat = rna; in mpscct_0_()
82 /* Subroutine */ int mpfcct_(real *xlon, real *ylat, real *x, real *y) in mpfcct_() argument
84 return mpscct_0_(1, (real *)0, xlon, ylat, x, y); in mpfcct_()
87 /* Subroutine */ int mpicct_(real *x, real *y, real *xlon, real *ylat) in mpicct_() argument
89 return mpscct_0_(2, (real *)0, xlon, ylat, x, y); in mpicct_()
|
| /dports/science/cdcl/dcl-5.4.8-C/src/grph1/stpack/ |
| H A D | stfusr.c | 87 mpfcyc_(xlon, ylat, x, y); in stfusr_0_()
89 mpfcya_(xlon, ylat, x, y); in stfusr_0_()
91 mpfcyb_(xlon, ylat, x, y); in stfusr_0_()
93 mpfglb_(xlon, ylat, x, y); in stfusr_0_()
95 mpfcct_(xlon, ylat, x, y); in stfusr_0_()
108 mpicyc_(x, y, xlon, ylat); in stfusr_0_()
110 mpicya_(x, y, xlon, ylat); in stfusr_0_()
112 mpicyb_(x, y, xlon, ylat); in stfusr_0_()
114 mpiglb_(x, y, xlon, ylat); in stfusr_0_()
116 mpicct_(x, y, xlon, ylat); in stfusr_0_()
[all …]
|