Remove Chebyshev polynomials from proj

docs/source/apps/proj.rst

@@ -148,26 +148,13 @@ The following control parameters can appear in any order
 .. option:: -v
 
     Causes a listing of cartographic control parameters tested for and used by
-    the program to be printed prior to input data. Should not be used with the
-    :option:`-T` option.
+    the program to be printed prior to input data.
 
 .. option:: -V
 
     This option causes an expanded annotated listing of the characteristics of
     the projected point. :option:`-v` is implied with this option.
 
-.. option:: -T <ulow,uhi,vlow,vhi,res[,umax,vmax]>
-
-    This option creates a set of bivariate Chebyshev polynomial coefficients
-    that approximate the selected cartographic projection on stdout. The values
-    *low* and *hi* denote the range of the input where the *u* or *v* prefixes apply to
-    respective longitude-x or latitude-y depending upon whether a forward or
-    inverse projection is selected. The integer *res* is a number specifying the
-    power of 10 precision of the approximation. For example, a *res* of -3
-    specifies an approximation with an accuracy better than 0.001. Optional *umax*, and *vmax*
-    specify maximum degree of the polynomials (default: 15).
-
-
 
 The *+args* run-line arguments are associated with cartographic parameters.
 Additional projection control parameters may be contained in two auxiliary

src/Makefile.am

@@ -18,9 +18,9 @@ EXTRA_DIST = bin_cct.cmake bin_gie.cmake bin_cs2cs.cmake \
 	lib_proj.cmake CMakeLists.txt bin_geodtest.cmake tests/geodtest.cpp \
 	wkt1_grammar.y wkt2_grammar.y apps/emess.h
 
-proj_SOURCES = apps/proj.cpp apps/gen_cheb.cpp apps/p_series.cpp apps/emess.cpp
+proj_SOURCES = apps/proj.cpp apps/emess.cpp
 projinfo_SOURCES = apps/projinfo.cpp
-cs2cs_SOURCES = apps/cs2cs.cpp apps/gen_cheb.cpp apps/p_series.cpp apps/emess.cpp
+cs2cs_SOURCES = apps/cs2cs.cpp apps/emess.cpp
 cct_SOURCES = apps/cct.cpp apps/proj_strtod.cpp apps/proj_strtod.h apps/optargpm.h
 nad2bin_SOURCES = apps/nad2bin.cpp
 geod_SOURCES = apps/geod.cpp apps/geod_set.cpp apps/geod_interface.cpp apps/geod_interface.h apps/emess.cpp

src/apps/proj.cpp

@@ -20,8 +20,6 @@
 #define MAX_PARGS 100
 #define PJ_INVERS(P) (P->inv ? 1 : 0)
 
-extern void gen_cheb(int, PJ_UV(*)(PJ_UV), const char *, PJ *, int, char **);
-
 static PJ *Proj;
 static union {
     PJ_UV (*generic)(PJ_UV, PJ *);
@@ -43,7 +41,6 @@ static int
     postscale = 0;
 
 static const char
-    *cheby_str,         /* string controlling Chebychev evaluation */
     *oform = nullptr; /* output format for x-y or decimal degrees */
 static char oform_buffer[16];   /* Buffer for oform when using -d */
 
@@ -56,22 +53,6 @@ static PJ_FACTORS facs;
 static double (*informat)(const char *, char **), /* input data deformatter function */
               fscale = 0.;                        /* cartesian scale factor */
 
-static PJ_UV int_proj(PJ_UV data) {
-    if (prescale) {
-        data.u *= fscale;
-        data.v *= fscale;
-    }
-
-    data = (*proj.generic)(data, Proj);
-
-    if (postscale && data.u != HUGE_VAL) {
-        data.u *= fscale;
-        data.v *= fscale;
-    }
-
-    return data;
-}
-
 /* file processing function */
 static void process(FILE *fid) {
     char line[MAX_LINE+3], *s = nullptr, pline[40];
@@ -309,10 +290,10 @@ static void vprocess(FILE *fid) {
 }
 
 int main(int argc, char **argv) {
-    char *arg, *pargv[MAX_PARGS], **iargv = argv;
+    char *arg, *pargv[MAX_PARGS];
     char **eargv = argv;
     FILE *fid;
-    int pargc = 0, iargc = argc, eargc = 0, mon = 0;
+    int pargc = 0, eargc = 0, mon = 0;
 
     if ( (emess_dat.Prog_name = strrchr(*argv,DIR_CHAR)) != nullptr)
         ++emess_dat.Prog_name;
@@ -425,10 +406,6 @@ int main(int argc, char **argv) {
                 emess(1,"missing argument for -%c",*arg);
                 oterr = *++argv;
                 continue;
-              case 'T': /* generate Chebyshev coefficients */
-                if (--argc <= 0) goto noargument;
-                cheby_str = *++argv;
-                continue;
               case 'm': /* cartesian multiplier */
                 if (--argc <= 0) goto noargument;
                 postscale = 1;
@@ -482,10 +459,9 @@ int main(int argc, char **argv) {
         } else /* assumed to be input file name(s) */
             eargv[eargc++] = *argv;
     }
-    if (eargc == 0 && !cheby_str) /* if no specific files force sysin */
+    if (eargc == 0) /* if no specific files force sysin */
         eargv[eargc++] = const_cast<char*>("-");
-    else if (eargc > 0 && cheby_str) /* warning */
-        emess(4, "data files when generating Chebychev prohibited");
+
     /* done with parameter and control input */
     if (inverse && postscale) {
         prescale = 1;
@@ -512,10 +488,6 @@ int main(int argc, char **argv) {
         proj.inv = pj_inv;
     } else
         proj.fwd = pj_fwd;
-    if (cheby_str) {
-        gen_cheb(inverse, int_proj, cheby_str, Proj, iargc, iargv);
-        exit(0);
-    }
 
     /* set input formatting control */
     if (mon) {

src/bin_cs2cs.cmake

@@ -1,6 +1,4 @@
 set(CS2CS_SRC apps/cs2cs.cpp
-              apps/gen_cheb.cpp 
-              apps/p_series.cpp
               apps/emess.cpp
 )
 

src/bin_proj.cmake

@@ -1,6 +1,4 @@
-set(PROJ_SRC apps/proj.cpp 
-             apps/gen_cheb.cpp 
-             apps/p_series.cpp
+set(PROJ_SRC apps/proj.cpp
              apps/emess.cpp
 )