--- ray/src/util/rcrop.c	2022/03/15 20:21:18	1.7
+++ ray/src/util/rcrop.c	2022/04/07 16:58:30	1.13
@@ -1,5 +1,5 @@
 #ifndef lint
-static const char RCSid[] = "$Id: rcrop.c,v 1.7 2022/03/15 20:21:18 greg Exp $";
+static const char RCSid[] = "$Id: rcrop.c,v 1.13 2022/04/07 16:58:30 greg Exp $";
 #endif
 /*
  * rcrop.c - crop a Radiance picture or matrix data
@@ -12,13 +12,11 @@ static const char RCSid[] = "$Id: rcrop.c,v 1.7 2022/0
 #include "fvect.h"
 #include "view.h"
 
-#define MAXWORD		64	/* maximum word (number) length */
-
 char	*progname;		/* global argv[0] */
 
 VIEW	vw = STDVIEW;
 int	gotvw = 0;
-char	fmt[MAXFMTLEN] = "Unknown";
+char	fmt[MAXFMTLEN] = "ascii";	/* assumed when unspecified */
 int	ncomp = 0;
 RESOLU	res;
 int	rmin, cmin, nrows, ncols;
@@ -100,9 +98,11 @@ binary_copyf(FILE *fp, int asize)
 					/* check if fseek() useful */
 	if (skip_len > skip_thresh &&
 			fseek(fp, (rmin*width + cmin)*elsiz, SEEK_CUR) == 0) {
+		off_t	curpos;
 		buf = (char *)malloc(ncols*elsiz);
 		if (!buf)
 			goto memerr;
+#ifdef NON_POSIX
 		for (y = nrows; y-- > 0; ) {
 			if (getbinary(buf, elsiz, ncols, fp) != ncols)
 				goto readerr;
@@ -114,8 +114,20 @@ binary_copyf(FILE *fp, int asize)
 				return(0);
 			}
 		}
-		free(buf);		/* success! */
-		return(1);
+#else
+		curpos = ftello(fp);
+		for (y = nrows; y-- > 0; curpos += width*elsiz) {
+			if (pread(fileno(fp), buf, ncols*elsiz,
+						curpos) != ncols*elsiz)
+				goto readerr;
+			if (putbinary(buf, elsiz, ncols, stdout) != ncols)
+				goto writerr;
+		}
+#endif
+		free(buf);
+		if (fflush(stdout) == EOF)
+			goto writerr;
+		return(1);		/* success! */
 	}				/* else need to read it all... */
 	buf = (char *)malloc(width*elsiz);
 	if (!buf)
@@ -199,6 +211,45 @@ io_err:
 	return(0);
 }
 
+/* Adjust (crop) our view */
+static int
+adjust_view(void)
+{
+	double		p0[2], p1[2];
+	const char	*err;
+
+	if (res.rt & YMAJOR) {
+		p0[0] = cmin/(double)res.xr;
+		p0[1] = rmin/(double)res.yr;
+		p1[0] = (cmin+ncols)/(double)res.xr;
+		p1[1] = (rmin+nrows)/(double)res.yr;
+	} else {
+		p0[0] = rmin/(double)res.xr;
+		p0[1] = cmin/(double)res.yr;
+		p1[0] = (rmin+nrows)/(double)res.xr;
+		p1[1] = (cmin+ncols)/(double)res.yr;
+	}
+	if (res.rt & XDECR) {
+		p0[0] = 1. - p0[0];
+		p1[0] = 1. - p1[0];
+	}
+	if (res.rt & YDECR) {
+		p0[1] = 1. - p0[1];
+		p1[1] = 1. - p1[1];
+	}
+	err = cropview(&vw, p0[0], p0[1], p1[0], p1[1]);
+
+	if (!err)
+		return(1);	/* success! */
+
+	fputs(progname, stderr);
+	fputs(": view error - ", stderr);
+	fputs(err, stderr);
+	fputc('\n', stderr);
+	return(0);		/* something went wrong */
+}
+
+
 /* Main routine -- load header and call processor */
 int
 main(int argc, char *argv[])
@@ -261,43 +312,13 @@ main(int argc, char *argv[])
 		fputs(": illegal crop\n", stderr);
 		return(1);
 	}
-	printargs(argc, argv, stdout);
-	if (gotvw) {		/* adjust view? */
-		double		p0[2], p1[2];
-		const char	*err;
-		if (res.rt & YMAJOR) {
-			p0[0] = cmin/(double)res.xr;
-			p0[1] = rmin/(double)res.yr;
-			p1[0] = (cmin+ncols)/(double)res.xr;
-			p1[1] = (rmin+nrows)/(double)res.yr;
-		} else {
-			p0[1] = cmin/(double)res.xr;
-			p0[0] = rmin/(double)res.yr;
-			p1[1] = (cmin+ncols)/(double)res.xr;
-			p1[0] = (rmin+nrows)/(double)res.yr;
-		}
-		if (res.rt & XDECR) {
-			p0[0] = 1. - p0[0];
-			p1[0] = 1. - p1[0];
-		}
-		if (res.rt & YDECR) {
-			p0[1] = 1. - p0[1];
-			p1[1] = 1. - p1[1];
-		}
-		err = cropview(&vw, p0[0], p0[1], p1[0], p1[1]);
-		if (err) {
-			fputs(progname, stderr);
-			fputs(": view error - ", stderr);
-			fputs(err, stderr);
-			fputc('\n', stderr);
-			return(1);
-		} else {
-			fputs(VIEWSTR, stdout);
-			fprintview(&vw, stdout);
-			fputc('\n', stdout);
-		}
+	printargs(5, argv, stdout);	/* add to header */
+	if (gotvw && adjust_view()) {
+		fputs(VIEWSTR, stdout);	/* write adjusted view */
+		fprintview(&vw, stdout);
+		fputc('\n', stdout);
 	}
-	if (gotdims)
+	if (gotdims)			/* dimensions + format */
 		printf("NROWS=%d\nNCOLS=%d\n", nrows, ncols);
 	if (ncomp)
 		printf("NCOMP=%d\n", ncomp);
@@ -343,7 +364,7 @@ main(int argc, char *argv[])
 		return(1);
 	}
 	if (!(asiz < 0 ? colr_copyf(fp) :
-			!asiz ? ascii_copyf(fp) : binary_copyf(fp, asiz)))
+			asiz ? binary_copyf(fp, asiz) : ascii_copyf(fp)))
 		return(1);
 					/* need to consume the rest? */
 	if (fp == stdin && rmin+nrows < numscans(&res) &&