--- ray/src/rt/ambient.c	1992/07/16 14:03:27	2.8
+++ ray/src/rt/ambient.c	1993/01/20 10:48:23	2.15
@@ -1,4 +1,4 @@
-/* Copyright (c) 1991 Regents of the University of California */
+/* Copyright (c) 1993 Regents of the University of California */
 
 #ifndef lint
 static char SCCSid[] = "$SunId$ LBL";
@@ -6,11 +6,6 @@ static char SCCSid[] = "$SunId$ LBL";
 
 /*
  *  ambient.c - routines dealing with ambient (inter-reflected) component.
- *
- *  The macro AMBFLUSH (if defined) is the number of ambient values
- *	to wait before flushing to the ambient file.
- *
- *     5/9/86
  */
 
 #include  "ray.h"
@@ -23,34 +18,34 @@ static char SCCSid[] = "$SunId$ LBL";
 
 #include  "random.h"
 
-#define  OCTSCALE	0.5	/* ceil((valid rad.)/(cube size)) */
+#define	 OCTSCALE	0.5	/* ceil((valid rad.)/(cube size)) */
 
 typedef struct ambtree {
-	AMBVAL  *alist;		/* ambient value list */
-	struct ambtree  *kid;	/* 8 child nodes */
+	AMBVAL	*alist;		/* ambient value list */
+	struct ambtree	*kid;	/* 8 child nodes */
 }  AMBTREE;			/* ambient octree */
 
 extern CUBE  thescene;		/* contains space boundaries */
 
-#define  MAXASET	511	/* maximum number of elements in ambient set */
-OBJECT  ambset[MAXASET+1]={0};	/* ambient include/exclude set */
+#define	 MAXASET	511	/* maximum number of elements in ambient set */
+OBJECT	ambset[MAXASET+1]={0};	/* ambient include/exclude set */
 
-double  maxarad;		/* maximum ambient radius */
-double  minarad;		/* minimum ambient radius */
+double	maxarad;		/* maximum ambient radius */
+double	minarad;		/* minimum ambient radius */
 
-static AMBTREE  atrunk;		/* our ambient trunk node */
+static AMBTREE	atrunk;		/* our ambient trunk node */
 
+static char  *ambfname = NULL;	/* ambient file name */
 static FILE  *ambfp = NULL;	/* ambient file pointer */
-static char  *afname;		/* ambient file name */
-static long  ambheadlen;	/* length of ambient file header */
 
-#define  AMBFLUSH	(BUFSIZ/AMBVALSIZ)
+#define	 AMBFLUSH	(BUFSIZ/AMBVALSIZ)
 
-#define  newambval()	(AMBVAL *)bmalloc(sizeof(AMBVAL))
+#define	 newambval()	(AMBVAL *)bmalloc(sizeof(AMBVAL))
 
-#define  newambtree()	(AMBTREE *)calloc(8, sizeof(AMBTREE))
+#define	 newambtree()	(AMBTREE *)calloc(8, sizeof(AMBTREE))
 
 extern long  ftell(), lseek();
+static int  initambfile(), avsave(), avinsert(), ambsync();
 
 
 setambres(ar)				/* set ambient resolution */
@@ -74,19 +69,21 @@ int  ar;
 setambient(afile)			/* initialize calculation */
 char  *afile;
 {
-	AMBVAL  amb;
+	long  headlen;
+	AMBVAL	amb;
 						/* init ambient limits */
 	setambres(ambres);
 						/* open ambient file */
-	if ((afname = afile) != NULL)
+	if ((ambfname = afile) != NULL)
 		if ((ambfp = fopen(afile, "r+")) != NULL) {
 			initambfile(0);
+			headlen = ftell(ambfp);
 			while (readambval(&amb, ambfp))
 				avinsert(&amb, &atrunk, thescene.cuorg,
 						thescene.cusize);
 							/* align */
-			fseek(ambfp, -((ftell(ambfp)-ambheadlen)%AMBVALSIZ), 1);
-		} else if ((ambfp = fopen(afile, "w")) != NULL)
+			fseek(ambfp, -((ftell(ambfp)-headlen)%AMBVALSIZ), 1);
+		} else if ((ambfp = fopen(afile, "w+")) != NULL)
 			initambfile(1);
 		else {
 			sprintf(errmsg, "cannot open ambient file \"%s\"",
@@ -96,39 +93,11 @@ char  *afile;
 }
 
 
-initambfile(creat)		/* initialize ambient file */
-int  creat;
-{
-	extern char  *progname, *octname, VersionID[];
-
-	setbuf(ambfp, bmalloc(BUFSIZ));
-	if (creat) {			/* new file */
-		fprintf(ambfp, "%s -av %g %g %g -ab %d -aa %g ",
-				progname, colval(ambval,RED),
-				colval(ambval,GRN), colval(ambval,BLU),
-				ambounce, ambacc);
-		fprintf(ambfp, "-ad %d -as %d -ar %d %s\n",
-				ambdiv, ambssamp, ambres,
-				octname==NULL ? "" : octname);
-		fprintf(ambfp, "SOFTWARE= %s\n", VersionID);
-		fputformat(AMBFMT, ambfp);
-		putc('\n', ambfp);
-		putambmagic(ambfp);
-		fflush(ambfp);
-	} else if (checkheader(ambfp, AMBFMT, NULL) < 0
-			|| !hasambmagic(ambfp)) {
-		sprintf(errmsg, "\"%s\" is not an ambient file", afname);
-		error(USER, errmsg);
-	}
-	ambheadlen = ftell(ambfp);
-}
-
-
 ambnotify(obj)			/* record new modifier */
-OBJECT  obj;
+OBJECT	obj;
 {
 	static int  hitlimit = 0;
-	register OBJREC  *o = objptr(obj);
+	register OBJREC	 *o = objptr(obj);
 	register char  **amblp;
 
 	if (hitlimit || !ismodifier(o->otype))
@@ -151,7 +120,7 @@ COLOR  acol;
 register RAY  *r;
 {
 	static int  rdepth = 0;			/* ambient recursion */
-	double  d;
+	double	d;
 
 	if (ambdiv <= 0)			/* no ambient calculation */
 		goto dumbamb;
@@ -193,17 +162,16 @@ sumambient(acol, r, al, at, c0, s)	/* get interpolated
 COLOR  acol;
 register RAY  *r;
 int  al;
-AMBTREE  *at;
+AMBTREE	 *at;
 FVECT  c0;
-double  s;
+double	s;
 {
-	extern double  sqrt();
-	double  d, e1, e2, wt, wsum;
+	double	d, e1, e2, wt, wsum;
 	COLOR  ct;
 	FVECT  ck0;
 	int  i;
 	register int  j;
-	register AMBVAL  *av;
+	register AMBVAL	 *av;
 					/* do this node */
 	wsum = 0.0;
 	for (av = at->alist; av != NULL; av = av->next) {
@@ -282,7 +250,7 @@ COLOR  acol;
 register RAY  *r;
 int  al;
 {
-	AMBVAL  amb;
+	AMBVAL	amb;
 	FVECT	gp, gd;
 						/* compute weight */
 	amb.weight = pow(AVGREFL, (double)al);
@@ -307,12 +275,12 @@ int  al;
 
 extambient(cr, ap, pv, nv)		/* extrapolate value at pv, nv */
 COLOR  cr;
-register AMBVAL  *ap;
+register AMBVAL	 *ap;
 FVECT  pv, nv;
 {
 	FVECT  v1, v2;
 	register int  i;
-	double  d;
+	double	d;
 
 	d = 1.0;			/* zeroeth order */
 					/* gradient due to translation */
@@ -332,8 +300,36 @@ FVECT  pv, nv;
 
 
 static
+initambfile(creat)		/* initialize ambient file */
+int  creat;
+{
+	extern char  *progname, *octname, VersionID[];
+
+#ifdef MSDOS
+	setmode(fileno(ambfp), O_BINARY);
+#endif
+	setbuf(ambfp, bmalloc(BUFSIZ));
+	if (creat) {			/* new file */
+		fprintf(ambfp, "%s -av %g %g %g -ab %d -aa %g ",
+				progname, colval(ambval,RED),
+				colval(ambval,GRN), colval(ambval,BLU),
+				ambounce, ambacc);
+		fprintf(ambfp, "-ad %d -as %d -ar %d %s\n",
+				ambdiv, ambssamp, ambres,
+				octname==NULL ? "" : octname);
+		fprintf(ambfp, "SOFTWARE= %s\n", VersionID);
+		fputformat(AMBFMT, ambfp);
+		putc('\n', ambfp);
+		putambmagic(ambfp);
+	} else if (checkheader(ambfp, AMBFMT, NULL) < 0 || !hasambmagic(ambfp))
+		error(USER, "bad ambient file");
+	ambsync();
+}
+
+
+static
 avsave(av)				/* insert and save an ambient value */
-AMBVAL  *av;
+AMBVAL	*av;
 {
 	static int  nunflshed = 0;
 
@@ -355,14 +351,14 @@ writerr:
 
 static
 avinsert(aval, at, c0, s)		/* insert ambient value in a tree */
-AMBVAL  *aval;
+AMBVAL	*aval;
 register AMBTREE  *at;
 FVECT  c0;
-double  s;
+double	s;
 {
 	FVECT  ck0;
 	int  branch;
-	register AMBVAL  *av;
+	register AMBVAL	 *av;
 	register int  i;
 
 	if ((av = newambval()) == NULL)
@@ -390,16 +386,24 @@ memerr:
 }
 
 
-#include  <fcntl.h>
+#ifdef	NIX
 
+static
+ambsync()			/* flush ambient file */
+{
+	return(fflush(ambfp));
+}
 
+#else
+
 static
 ambsync()			/* synchronize ambient file */
 {
 	static FILE  *ambinp = NULL;
+	static long  lastpos = -1;
 	struct flock  fls;
-	AMBVAL  avs;
-	long  lastpos, flen;
+	long  flen;
+	AMBVAL	avs;
 	register int  n;
 				/* gain exclusive access */
 	fls.l_type = F_WRLCK;
@@ -408,22 +412,33 @@ ambsync()			/* synchronize ambient file */
 	fls.l_len = 0L;
 	if (fcntl(fileno(ambfp), F_SETLKW, &fls) < 0)
 		error(SYSTEM, "cannot lock ambient file");
+	if (lastpos < 0)	/* initializing */
+		goto syncend;
 				/* see if file has grown */
-	lastpos = lseek(fileno(ambfp), 0L, 2);	/* may move pointer */
-	flen = lseek(fileno(ambfp), 0L, 1);	/* new(?) file length */
-	if (n = (flen - lastpos)%AMBVALSIZ)	/* assure alignment */
-		lseek(fileno(ambfp), flen -= n, 0);
-	if (n = (flen - lastpos)/AMBVALSIZ) {	/* file has grown */
-		if (ambinp == NULL && (ambinp = fopen(afname, "r")) == NULL)
-			error(SYSTEM, "cannot reopen ambient file");
-		fseek(ambinp, lastpos, 0);	/* go to previous position */
-		while (n--) {			/* load contributed values */
+	if ((flen = lseek(fileno(ambfp), 0L, 2)) < 0)
+		error(SYSTEM, "cannot seek on ambient file");
+	if (n = flen - lastpos) {		/* file has grown */
+		if (ambinp == NULL) {
+			ambinp = fopen(ambfname, "r");
+			if (ambinp == NULL)
+				error(SYSTEM, "fopen failed in ambsync");
+		}
+		if (fseek(ambinp, lastpos, 0) < 0)
+			error(SYSTEM, "fseek failed in ambsync");
+		while (n >= AMBVALSIZ) {	/* load contributed values */
 			readambval(&avs, ambinp);
 			avinsert(&avs,&atrunk,thescene.cuorg,thescene.cusize);
+			n -= AMBVALSIZ;
 		}
+		if (n)				/* alignment */
+			lseek(fileno(ambfp), flen-n, 0);
 	}
+syncend:
 	n = fflush(ambfp);			/* calls write() at last */
+	lastpos = lseek(fileno(ambfp), 0L, 1);
 	fls.l_type = F_UNLCK;			/* release file */
 	fcntl(fileno(ambfp), F_SETLKW, &fls);
 	return(n);
 }
+
+#endif