/dports/graphics/proj/proj-7.2.1/src/projections/ |
H A D | qsc.cpp | 98 theta -= M_HALFPI; in qsc_fwd_equat_face_theta() 99 } else if (theta > M_HALFPI + M_FORTPI || theta <= -(M_HALFPI + M_FORTPI)) { in qsc_fwd_equat_face_theta() 104 theta += M_HALFPI; in qsc_fwd_equat_face_theta() 146 phi = M_HALFPI - lat; in qsc_e_forward() 150 } else if (lon > M_HALFPI + M_FORTPI || lon <= -(M_HALFPI + M_FORTPI)) { in qsc_e_forward() 161 phi = M_HALFPI + lat; in qsc_e_forward() 223 mu += M_HALFPI; in qsc_e_forward() 254 mu -= M_HALFPI; in qsc_e_inverse() 260 mu += M_HALFPI; in qsc_e_inverse() 286 lp.phi = M_HALFPI - phi; in qsc_e_inverse() [all …]
|
H A D | healpix.cpp | 172 {-3*M_FORTPI, M_HALFPI + EPS}, in in_image() 173 {-M_HALFPI, M_FORTPI + EPS}, in in_image() 266 lamc = -3*M_FORTPI + M_HALFPI*cn; in healpix_sphere() 287 } else if (fabs(y) < M_HALFPI) { in healpix_spherhealpix_e_inverse() 293 xc = -3*M_FORTPI + M_HALFPI*cn; in healpix_spherhealpix_e_inverse() 299 lp.phi = sign(y)*M_HALFPI; in healpix_spherhealpix_e_inverse() 364 c = M_HALFPI; in get_cap() 367 c = -M_HALFPI; in get_cap() 374 if (x < -M_HALFPI) { in get_cap() 395 capmap.y = M_HALFPI; in get_cap() [all …]
|
H A D | nicol.cpp | 23 } else if (fabs(fabs(lp.lam) - M_HALFPI) < EPS) { in nicol_s_forward() 25 xy.y = M_HALFPI * sin(lp.phi); in nicol_s_forward() 26 } else if (fabs(fabs(lp.phi) - M_HALFPI) < EPS) { in nicol_s_forward() 32 tb = M_HALFPI / lp.lam - lp.lam / M_HALFPI; in nicol_s_forward() 33 c = lp.phi / M_HALFPI; in nicol_s_forward() 41 xy.x = M_HALFPI * ( m + (lp.lam < 0. ? -xy.x : xy.x)); in nicol_s_forward() 44 xy.y = M_HALFPI * ( n + (lp.phi < 0. ? xy.y : -xy.y )); in nicol_s_forward()
|
H A D | stere.cpp | 39 return (tan (.5 * (M_HALFPI + phit)) * in ssfn_() 130 if (fabs (lp.phi - M_HALFPI) < TOL) { in stere_s_forward() 161 tp = tan (.5 * (M_HALFPI + phi_l)); in stere_e_inverse() 164 halfpi = M_HALFPI; in stere_e_inverse() 172 phi_l = M_HALFPI - 2. * atan (tp); in stere_e_inverse() 173 halfpi = -M_HALFPI; in stere_e_inverse() 243 if (fabs ((t = fabs (P->phi0)) - M_HALFPI) < EPS10) in setup() 255 if (fabs (Q->phits - M_HALFPI) < EPS10) in setup() 288 Q->akm1 = fabs (Q->phits - M_HALFPI) >= EPS10 ? in setup() 321 P->phi0 = pj_param(P->ctx, P->params, "bsouth").i ? - M_HALFPI: M_HALFPI; in PROJECTION() [all …]
|
H A D | gnom.cpp | 97 lp.phi = lp.phi > 0. ? M_HALFPI : - M_HALFPI; in gnom_s_inverse() 106 lp.phi = lp.phi > 0. ? M_HALFPI : - M_HALFPI; in gnom_s_inverse() 113 lp.phi -= M_HALFPI; in gnom_s_inverse() 116 lp.phi = M_HALFPI - lp.phi; in gnom_s_inverse() 132 if (fabs(fabs(P->phi0) - M_HALFPI) < EPS10) { in PROJECTION()
|
H A D | aea.cpp | 146 lp.phi = lp.phi < 0. ? -M_HALFPI : M_HALFPI; in aea_e_inverse() 152 lp.phi = lp.phi < 0. ? -M_HALFPI : M_HALFPI; in aea_e_inverse() 157 lp.phi = Q->n > 0. ? M_HALFPI : - M_HALFPI; in aea_e_inverse() 170 if (fabs(Q->phi1) > M_HALFPI || fabs(Q->phi2) > M_HALFPI) in setup() 242 Q->phi1 = pj_param(P->ctx, P->params, "bsouth").i ? - M_HALFPI: M_HALFPI; in PROJECTION()
|
H A D | lcc.cpp | 28 if (fabs(fabs(lp.phi) - M_HALFPI) < EPS10) { in lcc_e_forward() 70 lp.phi = 2. * atan(pow(Q->c / rho, 1./Q->n)) - M_HALFPI; in lcc_e_inverse() 74 lp.phi = Q->n > 0. ? M_HALFPI : -M_HALFPI; in lcc_e_inverse() 97 if (fabs(Q->phi1) > M_HALFPI || fabs(Q->phi2) > M_HALFPI) in PROJECTION() 133 Q->rho0 *= (fabs(fabs(P->phi0) - M_HALFPI) < EPS10) ? 0. : in PROJECTION() 149 Q->rho0 = (fabs(fabs(P->phi0) - M_HALFPI) < EPS10) ? 0. : in PROJECTION()
|
H A D | mbtfpp.cpp | 38 lp.phi = (lp.phi < 0.) ? -M_HALFPI : M_HALFPI; in mbtfpp_s_inverse() 51 lp.phi = (lp.phi < 0.) ? -M_HALFPI : M_HALFPI; in mbtfpp_s_inverse()
|
H A D | vandg4.cpp | 21 } else if (fabs(lp.lam) < TOL || fabs(fabs(lp.phi) - M_HALFPI) < TOL) { in vandg4_s_forward() 33 if ((fabs(lp.lam) - M_HALFPI) < 0.) dt = -dt; in vandg4_s_forward() 43 xy.x = M_HALFPI * x1; in vandg4_s_forward() 44 xy.y = M_HALFPI * sqrt(1. + dt * fabs(x1) - x1 * x1); in vandg4_s_forward()
|
H A D | hatano.cpp | 53 th = th > 0. ? M_HALFPI : - M_HALFPI; in hatano_s_inverse() 67 lp.phi = lp.phi > 0. ? M_HALFPI : - M_HALFPI; in hatano_s_inverse()
|
/dports/math/vtk9/VTK-9.1.0/ThirdParty/libproj/vtklibproj/src/ |
H A D | PJ_qsc.c | 83 theta -= M_HALFPI; in qsc_fwd_equat_face_theta() 84 } else if (theta > M_HALFPI + M_FORTPI || theta <= -(M_HALFPI + M_FORTPI)) { in qsc_fwd_equat_face_theta() 89 theta += M_HALFPI; in qsc_fwd_equat_face_theta() 131 phi = M_HALFPI - lat; in e_forward() 135 } else if (lon > M_HALFPI + M_FORTPI || lon <= -(M_HALFPI + M_FORTPI)) { in e_forward() 146 phi = M_HALFPI + lat; in e_forward() 208 mu += M_HALFPI; in e_forward() 239 mu -= M_HALFPI; in e_inverse() 245 mu += M_HALFPI; in e_inverse() 271 lp.phi = M_HALFPI - phi; in e_inverse() [all …]
|
H A D | PJ_healpix.c | 267 lamc = -3*M_FORTPI + M_HALFPI*cn; in healpix_sphere() 288 } else if (fabsl(y) < M_HALFPI) { in healpix_sphere_inverse() 294 xc = -3*M_FORTPI + M_HALFPI*cn; in healpix_sphere_inverse() 300 lp.phi = pj_sign(y)*M_HALFPI; in healpix_sphere_inverse() 365 c = M_HALFPI; in get_cap() 368 c = -M_HALFPI; in get_cap() 375 if (x < -M_HALFPI) { in get_cap() 396 capmap.y = M_HALFPI; in get_cap() 401 capmap.y = -M_HALFPI; in get_cap() 471 a[1] = M_HALFPI; in combine_caps() [all …]
|
H A D | PJ_nocol.c | 19 } else if (fabs(fabs(lp.lam) - M_HALFPI) < EPS) { in s_forward() 21 xy.y = M_HALFPI * sin(lp.phi); in s_forward() 22 } else if (fabs(fabs(lp.phi) - M_HALFPI) < EPS) { in s_forward() 28 tb = M_HALFPI / lp.lam - lp.lam / M_HALFPI; in s_forward() 29 c = lp.phi / M_HALFPI; in s_forward() 37 xy.x = M_HALFPI * ( m + (lp.lam < 0. ? -xy.x : xy.x)); in s_forward() 40 xy.y = M_HALFPI * ( n + (lp.phi < 0. ? xy.y : -xy.y )); in s_forward()
|
H A D | PJ_stere.c | 29 return (tan (.5 * (M_HALFPI + phit)) * in ssfn_() 102 if (fabs (lp.phi - M_HALFPI) < TOL) F_ERROR; in s_forward() 129 tp = tan (.5 * (M_HALFPI + phi_l)); in e_inverse() 132 halfpi = M_HALFPI; in e_inverse() 138 phi_l = M_HALFPI - 2. * atan (tp = - rho / Q->akm1); in e_inverse() 139 halfpi = -M_HALFPI; in e_inverse() 215 if (fabs ((t = fabs (P->phi0)) - M_HALFPI) < EPS10) in setup() 227 if (fabs (Q->phits - M_HALFPI) < EPS10) in setup() 259 Q->akm1 = fabs (Q->phits - M_HALFPI) >= EPS10 ? in setup() 292 P->phi0 = pj_param(P->ctx, P->params, "bsouth").i ? - M_HALFPI: M_HALFPI; in PROJECTION() [all …]
|
H A D | PJ_gnom.c | 80 lp.phi = lp.phi > 0. ? M_HALFPI : - M_HALFPI; in s_inverse() 89 lp.phi = lp.phi > 0. ? M_HALFPI : - M_HALFPI; in s_inverse() 96 lp.phi -= M_HALFPI; in s_inverse() 99 lp.phi = M_HALFPI - lp.phi; in s_inverse() 132 if (fabs(fabs(P->phi0) - M_HALFPI) < EPS10) { in PROJECTION()
|
H A D | PJ_ortho.c | 39 if (fabs(lp.phi - P->phi0) - EPS10 > M_HALFPI) F_ERROR; in s_forward() 81 lp.phi = lp.phi < 0. ? -M_HALFPI : M_HALFPI; in s_inverse() 87 ? (xy.x == 0. ? 0. : xy.x < 0. ? -M_HALFPI : M_HALFPI) in s_inverse() 116 if (fabs(fabs(P->phi0) - M_HALFPI) <= EPS10) in PROJECTION()
|
H A D | PJ_aea.c | 110 lp.phi = lp.phi < 0. ? -M_HALFPI : M_HALFPI; in e_inverse() 114 lp.phi = lp.phi < 0. ? -M_HALFPI : M_HALFPI; in e_inverse() 118 lp.phi = Q->n > 0. ? M_HALFPI : - M_HALFPI; in e_inverse() 208 Q->phi1 = pj_param(P->ctx, P->params, "bsouth").i ? - M_HALFPI: M_HALFPI; in PROJECTION()
|
/dports/math/vtk8/VTK-8.2.0/ThirdParty/libproj/vtklibproj/src/ |
H A D | PJ_qsc.c | 83 theta -= M_HALFPI; in qsc_fwd_equat_face_theta() 84 } else if (theta > M_HALFPI + M_FORTPI || theta <= -(M_HALFPI + M_FORTPI)) { in qsc_fwd_equat_face_theta() 89 theta += M_HALFPI; in qsc_fwd_equat_face_theta() 131 phi = M_HALFPI - lat; in e_forward() 135 } else if (lon > M_HALFPI + M_FORTPI || lon <= -(M_HALFPI + M_FORTPI)) { in e_forward() 146 phi = M_HALFPI + lat; in e_forward() 208 mu += M_HALFPI; in e_forward() 239 mu -= M_HALFPI; in e_inverse() 245 mu += M_HALFPI; in e_inverse() 271 lp.phi = M_HALFPI - phi; in e_inverse() [all …]
|
H A D | PJ_healpix.c | 267 lamc = -3*M_FORTPI + M_HALFPI*cn; in healpix_sphere() 288 } else if (fabsl(y) < M_HALFPI) { in healpix_sphere_inverse() 294 xc = -3*M_FORTPI + M_HALFPI*cn; in healpix_sphere_inverse() 300 lp.phi = pj_sign(y)*M_HALFPI; in healpix_sphere_inverse() 365 c = M_HALFPI; in get_cap() 368 c = -M_HALFPI; in get_cap() 375 if (x < -M_HALFPI) { in get_cap() 396 capmap.y = M_HALFPI; in get_cap() 401 capmap.y = -M_HALFPI; in get_cap() 471 a[1] = M_HALFPI; in combine_caps() [all …]
|
H A D | PJ_nocol.c | 19 } else if (fabs(fabs(lp.lam) - M_HALFPI) < EPS) { in s_forward() 21 xy.y = M_HALFPI * sin(lp.phi); in s_forward() 22 } else if (fabs(fabs(lp.phi) - M_HALFPI) < EPS) { in s_forward() 28 tb = M_HALFPI / lp.lam - lp.lam / M_HALFPI; in s_forward() 29 c = lp.phi / M_HALFPI; in s_forward() 37 xy.x = M_HALFPI * ( m + (lp.lam < 0. ? -xy.x : xy.x)); in s_forward() 40 xy.y = M_HALFPI * ( n + (lp.phi < 0. ? xy.y : -xy.y )); in s_forward()
|
H A D | PJ_stere.c | 29 return (tan (.5 * (M_HALFPI + phit)) * in ssfn_() 102 if (fabs (lp.phi - M_HALFPI) < TOL) F_ERROR; in s_forward() 129 tp = tan (.5 * (M_HALFPI + phi_l)); in e_inverse() 132 halfpi = M_HALFPI; in e_inverse() 138 phi_l = M_HALFPI - 2. * atan (tp = - rho / Q->akm1); in e_inverse() 139 halfpi = -M_HALFPI; in e_inverse() 215 if (fabs ((t = fabs (P->phi0)) - M_HALFPI) < EPS10) in setup() 227 if (fabs (Q->phits - M_HALFPI) < EPS10) in setup() 259 Q->akm1 = fabs (Q->phits - M_HALFPI) >= EPS10 ? in setup() 292 P->phi0 = pj_param(P->ctx, P->params, "bsouth").i ? - M_HALFPI: M_HALFPI; in PROJECTION() [all …]
|
H A D | PJ_gnom.c | 80 lp.phi = lp.phi > 0. ? M_HALFPI : - M_HALFPI; in s_inverse() 89 lp.phi = lp.phi > 0. ? M_HALFPI : - M_HALFPI; in s_inverse() 96 lp.phi -= M_HALFPI; in s_inverse() 99 lp.phi = M_HALFPI - lp.phi; in s_inverse() 132 if (fabs(fabs(P->phi0) - M_HALFPI) < EPS10) { in PROJECTION()
|
H A D | PJ_ortho.c | 39 if (fabs(lp.phi - P->phi0) - EPS10 > M_HALFPI) F_ERROR; in s_forward() 81 lp.phi = lp.phi < 0. ? -M_HALFPI : M_HALFPI; in s_inverse() 87 ? (xy.x == 0. ? 0. : xy.x < 0. ? -M_HALFPI : M_HALFPI) in s_inverse() 116 if (fabs(fabs(P->phi0) - M_HALFPI) <= EPS10) in PROJECTION()
|
H A D | PJ_aea.c | 110 lp.phi = lp.phi < 0. ? -M_HALFPI : M_HALFPI; in e_inverse() 114 lp.phi = lp.phi < 0. ? -M_HALFPI : M_HALFPI; in e_inverse() 118 lp.phi = Q->n > 0. ? M_HALFPI : - M_HALFPI; in e_inverse() 208 Q->phi1 = pj_param(P->ctx, P->params, "bsouth").i ? - M_HALFPI: M_HALFPI; in PROJECTION()
|
H A D | PJ_mbtfpp.c | 33 lp.phi = (lp.phi < 0.) ? -M_HALFPI : M_HALFPI; in s_inverse() 42 lp.phi = (lp.phi < 0.) ? -M_HALFPI : M_HALFPI; in s_inverse()
|