diff options
| author | Even Rouault <even.rouault@spatialys.com> | 2017-02-26 13:50:07 +0100 |
|---|---|---|
| committer | Even Rouault <even.rouault@spatialys.com> | 2017-02-26 13:50:07 +0100 |
| commit | ac66e57366f11cd76d187dbc1e482d3f131dcd72 (patch) | |
| tree | d326ab624bef2cc7ebc87b55a4af44bcb8acc35d | |
| parent | abb21adb39de752b427e1e6fc66d94cfc9c3ee04 (diff) | |
| download | PROJ-ac66e57366f11cd76d187dbc1e482d3f131dcd72.tar.gz PROJ-ac66e57366f11cd76d187dbc1e482d3f131dcd72.zip | |
PJ_lsat.c: restore d__1 computation that was lost in conversion done in c9f12e0033474518fa460444b9948f36ce47d51f. Issue raised by VS12
| -rw-r--r-- | src/PJ_lsat.c | 1 |
1 files changed, 1 insertions, 0 deletions
diff --git a/src/PJ_lsat.c b/src/PJ_lsat.c index 292a5827..a58bb3cb 100644 --- a/src/PJ_lsat.c +++ b/src/PJ_lsat.c @@ -22,6 +22,7 @@ static void seraz0(double lam, double mult, PJ *P) { s = Q->p22 * Q->sa * cos(lam) * sqrt((1. + Q->t * sdsq) / ((1. + Q->w * sdsq) * (1. + Q->q * sdsq))); + d__1 = 1. + Q->q * sdsq; h = sqrt((1. + Q->q * sdsq) / (1. + Q->w * sdsq)) * ((1. + Q->w * sdsq) / (d__1 * d__1) - Q->p22 * Q->ca); |