aboutsummaryrefslogtreecommitdiff
path: root/src/projections/tmerc.cpp
diff options
context:
space:
mode:
authorEven Rouault <even.rouault@spatialys.com>2020-03-06 20:38:38 +0100
committerEven Rouault <even.rouault@spatialys.com>2020-03-08 19:58:35 +0100
commit60d3df673ca224107eb63e459073fc11ab5f4f16 (patch)
tree52c871c17f0c98c90ce6415cdb79dc1e73c14f4d /src/projections/tmerc.cpp
parent38ec5e662a74d40e02e38dc5dca553c3ecb04356 (diff)
downloadPROJ-60d3df673ca224107eb63e459073fc11ab5f4f16.tar.gz
PROJ-60d3df673ca224107eb63e459073fc11ab5f4f16.zip
src/projections/: remove assignments in expression and multiple statements per line
Should hopefully result in no change in results, and hopefully more readable code...
Diffstat (limited to 'src/projections/tmerc.cpp')
-rw-r--r--src/projections/tmerc.cpp14
1 files changed, 7 insertions, 7 deletions
diff --git a/src/projections/tmerc.cpp b/src/projections/tmerc.cpp
index 5c76b141..4b2a96f0 100644
--- a/src/projections/tmerc.cpp
+++ b/src/projections/tmerc.cpp
@@ -152,21 +152,21 @@ static PJ_XY approx_s_fwd (PJ_LP lp, PJ *P) {
static PJ_LP approx_e_inv (PJ_XY xy, PJ *P) {
PJ_LP lp = {0.0,0.0};
struct pj_opaque_approx *Q = static_cast<struct pj_opaque_approx*>(P->opaque);
- double n, con, cosphi, d, ds, sinphi, t;
lp.phi = pj_inv_mlfn(P->ctx, Q->ml0 + xy.y / P->k0, P->es, Q->en);
if (fabs(lp.phi) >= M_HALFPI) {
lp.phi = xy.y < 0. ? -M_HALFPI : M_HALFPI;
lp.lam = 0.;
} else {
- sinphi = sin(lp.phi);
- cosphi = cos(lp.phi);
- t = fabs (cosphi) > 1e-10 ? sinphi/cosphi : 0.;
- n = Q->esp * cosphi * cosphi;
- d = xy.x * sqrt (con = 1. - P->es * sinphi * sinphi) / P->k0;
+ const double sinphi = sin(lp.phi);
+ const double cosphi = cos(lp.phi);
+ double t = fabs (cosphi) > 1e-10 ? sinphi/cosphi : 0.;
+ const double n = Q->esp * cosphi * cosphi;
+ double con = 1. - P->es * sinphi * sinphi;
+ const double d = xy.x * sqrt (con) / P->k0;
con *= t;
t *= t;
- ds = d * d;
+ const double ds = d * d;
lp.phi -= (con * ds / (1.-P->es)) * FC2 * (1. -
ds * FC4 * (5. + t * (3. - 9. * n) + n * (1. - 4 * n) -
ds * FC6 * (61. + t * (90. - 252. * n +