aboutsummaryrefslogtreecommitdiff
path: root/src/PJ_tcc.cpp
diff options
context:
space:
mode:
authorEven Rouault <even.rouault@spatialys.com>2018-12-18 20:24:11 +0100
committerEven Rouault <even.rouault@spatialys.com>2018-12-26 10:08:53 +0100
commit610957f7035242f15743c399ffd429b92bc36206 (patch)
tree73f0d51147e2f4860c4bfc875f7a4bf9359386d4 /src/PJ_tcc.cpp
parent355d681ed88019e97742344bd642c2fd97e700a1 (diff)
downloadPROJ-610957f7035242f15743c399ffd429b92bc36206.tar.gz
PROJ-610957f7035242f15743c399ffd429b92bc36206.zip
cpp conversion: minimal steps to fix compilation errors, not warnings
Diffstat (limited to 'src/PJ_tcc.cpp')
-rw-r--r--src/PJ_tcc.cpp34
1 files changed, 34 insertions, 0 deletions
diff --git a/src/PJ_tcc.cpp b/src/PJ_tcc.cpp
new file mode 100644
index 00000000..60ded63e
--- /dev/null
+++ b/src/PJ_tcc.cpp
@@ -0,0 +1,34 @@
+#define PJ_LIB__
+
+#include <math.h>
+
+#include "proj.h"
+#include "projects.h"
+
+PROJ_HEAD(tcc, "Transverse Central Cylindrical") "\n\tCyl, Sph, no inv";
+
+#define EPS10 1.e-10
+
+
+static XY s_forward (LP lp, PJ *P) { /* Spheroidal, forward */
+ XY xy = {0.0, 0.0};
+ double b, bt;
+
+ b = cos (lp.phi) * sin (lp.lam);
+ if ((bt = 1. - b * b) < EPS10) {
+ proj_errno_set(P, PJD_ERR_TOLERANCE_CONDITION);
+ return xy;
+ }
+ xy.x = b / sqrt(bt);
+ xy.y = atan2 (tan (lp.phi) , cos (lp.lam));
+ return xy;
+}
+
+
+PJ *PROJECTION(tcc) {
+ P->es = 0.;
+ P->fwd = s_forward;
+ P->inv = 0;
+
+ return P;
+}