--- ray/src/cal/tabfunc.c	2003/06/26 00:58:09	1.3
+++ ray/src/cal/tabfunc.c	2019/08/11 17:04:09	1.7
@@ -1,5 +1,5 @@
 #ifndef lint
-static const char	RCSid[] = "$Id: tabfunc.c,v 1.3 2003/06/26 00:58:09 schorsch Exp $";
+static const char	RCSid[] = "$Id: tabfunc.c,v 1.7 2019/08/11 17:04:09 greg Exp $";
 #endif
 /*
  * Put tabular data into functions suitable for cal programs.
@@ -12,20 +12,16 @@ static const char	RCSid[] = "$Id: tabfunc.c,v 1.3 2003
 #include <string.h>
 #include <math.h>
 #include <ctype.h>
-#ifdef _WIN32
- #include <process.h> /* getpid() */
-#else
- #include <sys/types.h>
- #include <unistd.h>
-#endif
+#include <sys/types.h>
 
-#include "standard.h"
-#include "platform.h"
+#include "rtprocess.h" /* getpid() */
+#include "rtmath.h"
+#include "rtio.h"
 
 #define  isdelim(c)	(isspace(c) || (c)==',')
 
-#define  MAXTAB		1024		/* maximum number of data rows */
-#define  MAXLINE	4096		/* maximum line width (characters) */
+#define  MAXTAB		8192		/* maximum number of data rows */
+#define  MAXLINE	16384		/* maximum line width (characters) */
 #define  OUTFMT		"%.7g"		/* output format conversion string */
 
 int	interpolate = 0;
@@ -155,7 +151,9 @@ absc_exp(void)			/* produce expression for abscissa */
 				strcpy(ourexp, "x");
 			else
 				sprintf(ourexp, "x-%g", abscissa[0]-1);
-		} else
+		} else if (fabs(abscissa[0]) < eps)
+			sprintf(ourexp, "x/%g+1", step);
+		else
 			sprintf(ourexp, "(x-%g)/%g+1", abscissa[0], step);
 	} else {
 		printf("X`%s(i):select(i,", locID);
@@ -191,7 +189,7 @@ char	*xe
 
 	xelen = strlen(xe);
 	for (i = 0; i < nfuncs; i++) {
-		if (func[i][0] == '\0' | func[i][0] == '0')
+		if ((func[i][0] == '\0') | (func[i][0] == '0'))
 			continue;
 		if (interpolate) {
 			printf("%s`%s(i):select(i,", func[i], locID);