| /third_party/boost/boost/geometry/srs/projections/proj/ |
| D | imw_p.hpp | 84 par_imw_p<T> & proj_parm, T *del, T *sig) in phi12() argument
88 … if (!pj_param_r<srs::spar::lat_1>(params, "lat_1", srs::dpar::lat_1, proj_parm.phi_1) || in phi12()
89 … !pj_param_r<srs::spar::lat_2>(params, "lat_2", srs::dpar::lat_2, proj_parm.phi_2)) { in phi12()
94 *del = 0.5 * (proj_parm.phi_2 - proj_parm.phi_1); in phi12()
95 *sig = 0.5 * (proj_parm.phi_2 + proj_parm.phi_1); in phi12()
103 par_imw_p<T> const& proj_parm, in loc_for() argument
115 m = pj_mlfn(lp_phi, sp, cos(lp_phi), proj_parm.en); in loc_for()
116 xa = proj_parm.Pp + proj_parm.Qp * m; in loc_for()
117 ya = proj_parm.P + proj_parm.Q * m; in loc_for()
122 if (proj_parm.mode == phi_2_is_zero) { in loc_for()
[all …]
|
| D | sconics.hpp | 87 inline int phi12(Params const& params, par_sconics<T>& proj_parm, T *del) in phi12() argument
99 proj_parm.sig = 0.5 * (p2 + p1); in phi12()
100 err = (fabs(*del) < epsilon || fabs(proj_parm.sig) < epsilon) ? -42 : 0; in phi12()
166 …inline void setup(Params const& params, Parameters& par, par_sconics<T>& proj_parm, proj_type type… in setup() argument
173 proj_parm.type = type; in setup()
175 err = phi12(params, proj_parm, &del); in setup()
179 switch (proj_parm.type) { in setup()
181 proj_parm.n = sin(proj_parm.sig); in setup()
183 proj_parm.rho_c = proj_parm.n / cs + cs / proj_parm.n; in setup()
184 proj_parm.rho_0 = sqrt((proj_parm.rho_c - 2 * sin(par.phi0))/proj_parm.n); in setup()
[all …]
|
| D | tpeqd.hpp | 123 inline void setup_tpeqd(Params const& params, Parameters& par, par_tpeqd<T>& proj_parm) in setup_tpeqd() argument
137 proj_parm.dlam2 = adjlon(lam_2 - lam_1); in setup_tpeqd()
139 proj_parm.cp1 = cos(phi_1); in setup_tpeqd()
140 proj_parm.cp2 = cos(phi_2); in setup_tpeqd()
141 proj_parm.sp1 = sin(phi_1); in setup_tpeqd()
142 proj_parm.sp2 = sin(phi_2); in setup_tpeqd()
143 proj_parm.cs = proj_parm.cp1 * proj_parm.sp2; in setup_tpeqd()
144 proj_parm.sc = proj_parm.sp1 * proj_parm.cp2; in setup_tpeqd()
145 proj_parm.ccs = proj_parm.cp1 * proj_parm.cp2 * sin(proj_parm.dlam2); in setup_tpeqd()
146 …proj_parm.z02 = aacos(proj_parm.sp1 * proj_parm.sp2 + proj_parm.cp1 * proj_parm.cp2 * cos(proj_par… in setup_tpeqd()
[all …]
|
| D | aea.hpp | 177 inline void setup(Parameters const& par, par_aea<T>& proj_parm) in setup() argument
182 if (fabs(proj_parm.phi1 + proj_parm.phi2) < epsilon10) in setup()
184 proj_parm.n = sinphi = sin(proj_parm.phi1); in setup()
185 cosphi = cos(proj_parm.phi1); in setup()
186 secant = fabs(proj_parm.phi1 - proj_parm.phi2) >= epsilon10; in setup()
187 if( (proj_parm.ellips = (par.es > 0.))) { in setup()
190 proj_parm.en = pj_enfn<T>(par.es); in setup()
196 sinphi = sin(proj_parm.phi2); in setup()
197 cosphi = cos(proj_parm.phi2); in setup()
203 proj_parm.n = (m1 * m1 - m2 * m2) / (ml2 - ml1); in setup()
[all …]
|
| D | eqdc.hpp | 128 inline void setup_eqdc(Params const& params, Parameters& par, par_eqdc<T>& proj_parm) in setup_eqdc() argument
133 … proj_parm.phi1 = pj_get_param_r<T, srs::spar::lat_1>(params, "lat_1", srs::dpar::lat_1); in setup_eqdc()
134 … proj_parm.phi2 = pj_get_param_r<T, srs::spar::lat_2>(params, "lat_2", srs::dpar::lat_2); in setup_eqdc()
136 if (fabs(proj_parm.phi1 + proj_parm.phi2) < epsilon10) in setup_eqdc()
139 proj_parm.en = pj_enfn<T>(par.es); in setup_eqdc()
141 proj_parm.n = sinphi = sin(proj_parm.phi1); in setup_eqdc()
142 cosphi = cos(proj_parm.phi1); in setup_eqdc()
143 secant = fabs(proj_parm.phi1 - proj_parm.phi2) >= epsilon10; in setup_eqdc()
144 if( (proj_parm.ellips = (par.es > 0.)) ) { in setup_eqdc()
148 ml1 = pj_mlfn(proj_parm.phi1, sinphi, cosphi, proj_parm.en); in setup_eqdc()
[all …]
|
| D | lcc.hpp | 148 inline void setup_lcc(Params const& params, Parameters& par, par_lcc<T>& proj_parm) in setup_lcc() argument
156 proj_parm.phi1 = 0.0; in setup_lcc()
157 proj_parm.phi2 = 0.0; in setup_lcc()
158 …bool is_phi1_set = pj_param_r<srs::spar::lat_1>(params, "lat_1", srs::dpar::lat_1, proj_parm.phi1); in setup_lcc()
159 …bool is_phi2_set = pj_param_r<srs::spar::lat_2>(params, "lat_2", srs::dpar::lat_2, proj_parm.phi2); in setup_lcc()
166 proj_parm.phi1 = 33; in setup_lcc()
170 proj_parm.phi2 = 45; in setup_lcc()
177 proj_parm.phi2 = proj_parm.phi1; in setup_lcc()
179 par.phi0 = proj_parm.phi1; in setup_lcc()
181 if (fabs(proj_parm.phi1 + proj_parm.phi2) < epsilon10) in setup_lcc()
[all …]
|
| D | lsat.hpp | 72 seraz0(T lam, T const& mult, par_lsat<T>& proj_parm) in seraz0() argument
79 s = proj_parm.p22 * proj_parm.sa * cos(lam) * sqrt((1. + proj_parm.t * sdsq) in seraz0()
80 / ((1. + proj_parm.w * sdsq) * (1. + proj_parm.q * sdsq))); in seraz0()
82 d__1 = 1. + proj_parm.q * sdsq; in seraz0()
83 … h = sqrt((1. + proj_parm.q * sdsq) / (1. + proj_parm.w * sdsq)) * ((1. + proj_parm.w * sdsq) in seraz0()
84 / (d__1 * d__1) - proj_parm.p22 * proj_parm.ca); in seraz0()
86 sq = sqrt(proj_parm.xj * proj_parm.xj + s * s); in seraz0()
87 fc = mult * (h * proj_parm.xj - s * s) / sq; in seraz0()
88 proj_parm.b += fc; in seraz0()
89 proj_parm.a2 += fc * cos(lam + lam); in seraz0()
[all …]
|
| D | chamb.hpp | 201 inline void setup_chamb(Params const& params, Parameters& par, par_chamb<T>& proj_parm) in setup_chamb() argument
208 proj_parm.c[i].phi = chamb_init_lat<T>(params, i); in setup_chamb()
209 proj_parm.c[i].lam = chamb_init_lon<T>(params, i); in setup_chamb()
210 proj_parm.c[i].lam = adjlon(proj_parm.c[i].lam - par.lam0); in setup_chamb()
211 proj_parm.c[i].cosphi = cos(proj_parm.c[i].phi); in setup_chamb()
212 proj_parm.c[i].sinphi = sin(proj_parm.c[i].phi); in setup_chamb()
216 …proj_parm.c[i].v = vect(proj_parm.c[j].phi - proj_parm.c[i].phi, proj_parm.c[i].cosphi, proj_parm.… in setup_chamb()
217 … proj_parm.c[j].cosphi, proj_parm.c[j].sinphi, proj_parm.c[j].lam - proj_parm.c[i].lam); in setup_chamb()
218 if (proj_parm.c[i].v.r == 0.0) in setup_chamb()
222 proj_parm.beta_0 = lc(proj_parm.c[0].v.r, proj_parm.c[2].v.r, proj_parm.c[1].v.r); in setup_chamb()
[all …]
|
| D | rouss.hpp | 121 inline void setup_rouss(Parameters const& par, par_rouss<T>& proj_parm) in setup_rouss() argument
125 if (!proj_mdist_ini(par.es, proj_parm.en)) in setup_rouss()
129 proj_parm.s0 = proj_mdist(par.phi0, es2, cos(par.phi0), proj_parm.en); in setup_rouss()
136 proj_parm.C1 = proj_parm.A1 = R_R0_2 / 4.; in setup_rouss()
137 proj_parm.C2 = proj_parm.A2 = R_R0_2 * (2 * t2 - 1. - 2. * es2) / 12.; in setup_rouss()
138 proj_parm.A3 = R_R0_2 * t * (1. + 4. * t2)/ ( 12. * N0); in setup_rouss()
139 proj_parm.A4 = R_R0_4 / 24.; in setup_rouss()
140 proj_parm.A5 = R_R0_4 * ( -1. + t2 * (11. + 12. * t2))/24.; in setup_rouss()
141 proj_parm.A6 = R_R0_4 * ( -2. + t2 * (11. - 2. * t2))/240.; in setup_rouss()
142 proj_parm.B1 = t / (2. * N0); in setup_rouss()
[all …]
|
| D | aeqd.hpp | 100 …ward(T const& lp_lon, T const& lp_lat, T& xy_x, T& xy_y, Par const& par, ProjParm const& proj_parm) in e_forward() argument
109 switch (proj_parm.mode) { in e_forward()
114 … xy_x = (rho = fabs(proj_parm.Mp - pj_mlfn(lp_lat, sinphi, cosphi, proj_parm.en))) * in e_forward()
132 … >::apply(par.lam0, par.phi0, lp_lon + par.lam0, lp_lat, srs::spheroid<T>(par.a, proj_parm.b)); in e_forward()
141 …erse(T const& xy_x, T const& xy_y, T& lp_lon, T& lp_lat, Par const& par, ProjParm const& proj_parm) in e_inverse() argument
150 if (proj_parm.mode == obliq || proj_parm.mode == equit) { in e_inverse()
161 … >::apply(par.lam0, par.phi0, s12, azi1, srs::spheroid<T>(par.a, proj_parm.b)); in e_inverse()
166 … lp_lat = pj_inv_mlfn(proj_parm.mode == n_pole ? proj_parm.Mp - c : proj_parm.Mp + c, in e_inverse()
167 par.es, proj_parm.en); in e_inverse()
168 lp_lon = atan2(xy_x, proj_parm.mode == n_pole ? -xy_y : xy_y); in e_inverse()
[all …]
|
| D | omerc.hpp | 160 inline void setup_omerc(Params const& params, Parameters & par, par_omerc<T>& proj_parm) in setup_omerc() argument
171 … proj_parm.no_rot = pj_get_param_b<srs::spar::no_rot>(params, "no_rot", srs::dpar::no_rot); in setup_omerc()
205 proj_parm.B = cosph0 * cosph0; in setup_omerc()
206 proj_parm.B = sqrt(1. + par.es * proj_parm.B * proj_parm.B / par.one_es); in setup_omerc()
207 proj_parm.A = proj_parm.B * par.k0 * com / con; in setup_omerc()
208 D = proj_parm.B * com / (cosph0 * sqrt(con)); in setup_omerc()
216 proj_parm.E = F += D; in setup_omerc()
217 proj_parm.E *= math::pow(pj_tsfn(par.phi0, sinph0, par.e), proj_parm.B); in setup_omerc()
219 proj_parm.B = 1. / com; in setup_omerc()
220 proj_parm.A = par.k0; in setup_omerc()
[all …]
|
| D | nsper.hpp | 205 inline void setup(Params const& params, Parameters& par, par_nsper<T>& proj_parm) in setup() argument
207 proj_parm.height = pj_get_param_f<T, srs::spar::h>(params, "h", srs::dpar::h); in setup()
208 if (proj_parm.height <= 0.) in setup()
212 proj_parm.mode = par.phi0 < 0. ? s_pole : n_pole; in setup()
214 proj_parm.mode = equit; in setup()
216 proj_parm.mode = obliq; in setup()
217 proj_parm.sinph0 = sin(par.phi0); in setup()
218 proj_parm.cosph0 = cos(par.phi0); in setup()
220 proj_parm.pn1 = proj_parm.height / par.a; /* normalize by radius */ in setup()
221 proj_parm.p = 1. + proj_parm.pn1; in setup()
[all …]
|
| D | eck3.hpp | 109 inline void setup_eck3(Parameters& par, par_eck3<T>& proj_parm) in setup_eck3() argument
111 proj_parm.C_x = 0.42223820031577120149; in setup_eck3()
112 proj_parm.C_y = 0.84447640063154240298; in setup_eck3()
113 proj_parm.A = 1.0; in setup_eck3()
114 proj_parm.B = 0.4052847345693510857755; in setup_eck3()
121 inline void setup_putp1(Parameters& par, par_eck3<T>& proj_parm) in setup_putp1() argument
123 proj_parm.C_x = 1.89490; in setup_putp1()
124 proj_parm.C_y = 0.94745; in setup_putp1()
125 proj_parm.A = -0.5; in setup_putp1()
126 proj_parm.B = 0.30396355092701331433; in setup_putp1()
[all …]
|
| D | laea.hpp | 304 inline void setup_laea(Parameters& par, par_laea<T>& proj_parm) in setup_laea() argument
312 proj_parm.mode = par.phi0 < 0. ? s_pole : n_pole; in setup_laea()
314 proj_parm.mode = equit; in setup_laea()
316 proj_parm.mode = obliq; in setup_laea()
321 proj_parm.qp = pj_qsfn(1., par.e, par.one_es); in setup_laea()
322 proj_parm.mmf = .5 / (1. - par.es); in setup_laea()
323 proj_parm.apa = pj_authset<T>(par.es); in setup_laea()
324 switch (proj_parm.mode) { in setup_laea()
327 proj_parm.dd = 1.; in setup_laea()
330 proj_parm.dd = 1. / (proj_parm.rq = sqrt(.5 * proj_parm.qp)); in setup_laea()
[all …]
|
| D | gn_sinu.hpp | 168 inline void setup(Parameters& par, par_gn_sinu_s<T>& proj_parm) in setup() argument
172 … proj_parm.C_x = (proj_parm.C_y = sqrt((proj_parm.m + 1.) / proj_parm.n))/(proj_parm.m + 1.); in setup()
178 … inline void setup_gn_sinu(Params const& params, Parameters& par, par_gn_sinu_s<T>& proj_parm) in setup_gn_sinu() argument
180 if (pj_param_f<srs::spar::n>(params, "n", srs::dpar::n, proj_parm.n) in setup_gn_sinu()
181 && pj_param_f<srs::spar::m>(params, "m", srs::dpar::m, proj_parm.m)) { in setup_gn_sinu()
182 if (proj_parm.n <= 0 || proj_parm.m < 0) in setup_gn_sinu()
187 setup(par, proj_parm); in setup_gn_sinu()
192 inline void setup_sinu(Parameters const& par, par_gn_sinu_e<T>& proj_parm) in setup_sinu() argument
194 proj_parm.en = pj_enfn<T>(par.es); in setup_sinu()
199 inline void setup_sinu(Parameters& par, par_gn_sinu_s<T>& proj_parm) in setup_sinu() argument
[all …]
|
| D | etmerc.hpp | 239 inline void setup(Parameters& par, par_etmerc<T>& proj_parm) in setup() argument
257 proj_parm.cgb[0] = n*( 2 + n*(-2/3.0 + n*(-2 + n*(116/45.0 + n*(26/45.0 + in setup()
259 proj_parm.cbg[0] = n*(-2 + n*( 2/3.0 + n*( 4/3.0 + n*(-82/45.0 + n*(32/45.0 + in setup()
262 proj_parm.cgb[1] = np*(7/3.0 + n*( -8/5.0 + n*(-227/45.0 + n*(2704/315.0 + in setup()
264 proj_parm.cbg[1] = np*(5/3.0 + n*(-16/15.0 + n*( -13/9.0 + n*( 904/315.0 + in setup()
268 proj_parm.cgb[2] = np*( 56/15.0 + n*(-136/35.0 + n*(-1262/105.0 + in setup()
270 proj_parm.cbg[2] = np*(-26/15.0 + n*( 34/21.0 + n*( 8/5.0 + in setup()
274 proj_parm.cgb[3] = np*(4279/630.0 + n*(-332/35.0 + n*(-399572/14175.0))); in setup()
275 proj_parm.cbg[3] = np*(1237/630.0 + n*( -12/5.0 + n*( -24832/14175.0))); in setup()
277 proj_parm.cgb[4] = np*(4174/315.0 + n*(-144838/6237.0 )); in setup()
[all …]
|
| D | ocea.hpp | 115 inline void setup_ocea(Params const& params, Parameters& par, par_ocea<T>& proj_parm) in setup_ocea() argument
121 proj_parm.rok = 1. / par.k0; in setup_ocea()
122 proj_parm.rtk = par.k0; in setup_ocea()
129 proj_parm.singam = atan(-cos(alpha)/(-sin(phi_0) * sin(alpha))) + lonz; in setup_ocea()
131 proj_parm.sinphi = asin(cos(phi_0) * sin(alpha)); in setup_ocea()
140 proj_parm.singam = atan2(cos(phi_1) * sin(phi_2) * cos(lam_1) - in setup_ocea()
147 proj_parm.singam = -proj_parm.singam; in setup_ocea()
150 proj_parm.sinphi = atan(-cos(proj_parm.singam - lam_1) / tan(phi_1)); in setup_ocea()
152 par.lam0 = proj_parm.singam + half_pi; in setup_ocea()
153 proj_parm.cosphi = cos(proj_parm.sinphi); in setup_ocea()
[all …]
|
| D | oea.hpp | 119 inline void setup_oea(Params const& params, Parameters& par, par_oea<T>& proj_parm) in setup_oea() argument
121 … if (((proj_parm.n = pj_get_param_f<T, srs::spar::n>(params, "n", srs::dpar::n)) <= 0.) || in setup_oea()
122 … ((proj_parm.m = pj_get_param_f<T, srs::spar::m>(params, "m", srs::dpar::m)) <= 0.)) { in setup_oea()
125 … proj_parm.theta = pj_get_param_r<T, srs::spar::theta>(params, "theta", srs::dpar::theta); in setup_oea()
126 proj_parm.sp0 = sin(par.phi0); in setup_oea()
127 proj_parm.cp0 = cos(par.phi0); in setup_oea()
128 proj_parm.rn = 1./ proj_parm.n; in setup_oea()
129 proj_parm.rm = 1./ proj_parm.m; in setup_oea()
130 proj_parm.two_r_n = 2. * proj_parm.rn; in setup_oea()
131 proj_parm.two_r_m = 2. * proj_parm.rm; in setup_oea()
[all …]
|
| D | ob_tran.hpp | 241 inline void o_forward(T lp_lon, T lp_lat, T& xy_x, T& xy_y, Par const& proj_parm) in o_forward() argument
248 lp_lon = adjlon(aatan2(cosphi * sin(lp_lon), proj_parm.sphip * cosphi * coslam + in o_forward()
249 proj_parm.cphip * sinphi) + proj_parm.lamp); in o_forward()
250 lp_lat = aasin(proj_parm.sphip * sinphi - proj_parm.cphip * cosphi * coslam); in o_forward()
252 proj_parm.fwd(lp_lon, lp_lat, xy_x, xy_y); in o_forward()
256 … inline void o_inverse(T const& xy_x, T const& xy_y, T& lp_lon, T& lp_lat, Par const& proj_parm) in o_inverse() argument
260 proj_parm.inv(xy_x, xy_y, lp_lon, lp_lat); in o_inverse()
262 coslam = cos(lp_lon -= proj_parm.lamp); in o_inverse()
265 lp_lat = aasin(proj_parm.sphip * sinphi + proj_parm.cphip * cosphi * coslam); in o_inverse()
266 lp_lon = aatan2(cosphi * sin(lp_lon), proj_parm.sphip * cosphi * coslam - in o_inverse()
[all …]
|
| D | bonne.hpp | 166 … inline void setup_bonne(Params const& params, Parameters const& par, par_bonne<T>& proj_parm) in setup_bonne() argument
172 … proj_parm.phi1 = pj_get_param_r<T, srs::spar::lat_1>(params, "lat_1", srs::dpar::lat_1); in setup_bonne()
173 if (fabs(proj_parm.phi1) < epsilon10) in setup_bonne()
177 proj_parm.en = pj_enfn<T>(par.es); in setup_bonne()
178 proj_parm.m1 = pj_mlfn(proj_parm.phi1, proj_parm.am1 = sin(proj_parm.phi1), in setup_bonne()
179 c = cos(proj_parm.phi1), proj_parm.en); in setup_bonne()
180 … proj_parm.am1 = c / (sqrt(1. - par.es * proj_parm.am1 * proj_parm.am1) * proj_parm.am1); in setup_bonne()
182 if (fabs(proj_parm.phi1) + epsilon10 >= half_pi) in setup_bonne()
183 proj_parm.cphi1 = 0.; in setup_bonne()
185 proj_parm.cphi1 = 1. / tan(proj_parm.phi1); in setup_bonne()
|
| D | labrd.hpp | 166 … inline void setup_labrd(Params const& params, Parameters const& par, par_labrd<T>& proj_parm) in setup_labrd() argument
177 proj_parm.kRg = par.k0 * sqrt( N * R ); in setup_labrd()
178 proj_parm.p0s = atan( sqrt(R / N) * tan(par.phi0) ); in setup_labrd()
179 proj_parm.A = sinp / sin(proj_parm.p0s); in setup_labrd()
181 proj_parm.C = .5 * par.e * proj_parm.A * log((1. + t)/(1. - t)) + in setup_labrd()
182 - proj_parm.A * log( tan(fourth_pi + .5 * par.phi0)) in setup_labrd()
183 + log( tan(fourth_pi + .5 * proj_parm.p0s)); in setup_labrd()
185 … proj_parm.Ca = (1. - cos(t)) * ( proj_parm.Cb = 1. / (12. * proj_parm.kRg * proj_parm.kRg) ); in setup_labrd()
186 proj_parm.Cb *= sin(t); in setup_labrd()
187 proj_parm.Cc = 3. * (proj_parm.Ca * proj_parm.Ca - proj_parm.Cb * proj_parm.Cb); in setup_labrd()
[all …]
|
| D | stere.hpp | 297 … inline void setup(Parameters const& par, par_stere<T>& proj_parm) /* general initialization */ in setup() argument
305 proj_parm.mode = par.phi0 < 0. ? s_pole : n_pole; in setup()
307 proj_parm.mode = t > epsilon10 ? obliq : equit; in setup()
308 proj_parm.phits = fabs(proj_parm.phits); in setup()
313 switch (proj_parm.mode) { in setup()
316 if (fabs(proj_parm.phits - half_pi) < epsilon10) in setup()
317 proj_parm.akm1 = 2. * par.k0 / in setup()
320 proj_parm.akm1 = cos(proj_parm.phits) / in setup()
321 pj_tsfn(proj_parm.phits, t = sin(proj_parm.phits), par.e); in setup()
323 proj_parm.akm1 /= sqrt(1. - t * t); in setup()
[all …]
|
| D | mod_ster.hpp | 168 … inline void setup(Parameters& par, par_mod_ster<T>& proj_parm) /* general initialization */ in setup() argument
180 proj_parm.schio = sin(chio); in setup()
181 proj_parm.cchio = cos(chio); in setup()
187 inline void setup_mil_os(Parameters& par, par_mod_ster<T>& proj_parm) in setup_mil_os() argument
197 proj_parm.n = 2; in setup_mil_os()
200 proj_parm.zcoeff = AB; in setup_mil_os()
203 setup(par, proj_parm); in setup_mil_os()
208 inline void setup_lee_os(Parameters& par, par_mod_ster<T>& proj_parm) in setup_lee_os() argument
218 proj_parm.n = 2; in setup_lee_os()
221 proj_parm.zcoeff = AB; in setup_lee_os()
[all …]
|
| D | airy.hpp | 148 inline void setup_airy(Params const& params, Parameters& par, par_airy<T>& proj_parm) in setup_airy() argument
154 … proj_parm.no_cut = pj_get_param_b<srs::spar::no_cut>(params, "no_cut", srs::dpar::no_cut); in setup_airy()
157 proj_parm.Cb = -0.5; in setup_airy()
159 proj_parm.Cb = 1./tan(beta); in setup_airy()
160 proj_parm.Cb *= proj_parm.Cb * log(cos(beta)); in setup_airy()
165 proj_parm.p_halfpi = -half_pi; in setup_airy()
166 proj_parm.mode = s_pole; in setup_airy()
168 proj_parm.p_halfpi = half_pi; in setup_airy()
169 proj_parm.mode = n_pole; in setup_airy()
173 proj_parm.mode = equit; in setup_airy()
[all …]
|
| D | putp6.hpp | 117 inline void setup_putp6(Parameters& par, par_putp6<T>& proj_parm) in setup_putp6() argument
119 proj_parm.C_x = 1.01346; in setup_putp6()
120 proj_parm.C_y = 0.91910; in setup_putp6()
121 proj_parm.A = 4.; in setup_putp6()
122 proj_parm.B = 2.1471437182129378784; in setup_putp6()
123 proj_parm.D = 2.; in setup_putp6()
130 inline void setup_putp6p(Parameters& par, par_putp6<T>& proj_parm) in setup_putp6p() argument
132 proj_parm.C_x = 0.44329; in setup_putp6p()
133 proj_parm.C_y = 0.80404; in setup_putp6p()
134 proj_parm.A = 6.; in setup_putp6p()
[all …]
|