ViewVC Help

View File | Revision Log | Show Annotations | Download File | Root Listing

root/Development/ray/src/rt/rpict.c

Comparing ray/src/rt/rpict.c (file contents):
Revision 2.7 by greg, Wed Mar 11 09:32:40 1992 UTC vs.
Revision 2.8 by greg, Fri Jul 10 14:58:19 1992 UTC

#	Line 1 | Line 1
1	<	/* Copyright (c) 1991 Regents of the University of California */
1	>	/* Copyright (c) 1992 Regents of the University of California */
2
3		#ifndef lint
4		static char SCCSid[] = "$SunId$ LBL";
#	Line 77 | Line 77 | extern long  nrays;                    /* number of rays traced */
77
78		#define  pixjitter()    (.5+dstrpix*(.5-frandom()))
79
80	+	#define  HFTEMPLATE     "/tmp/hfXXXXXX"
81	+
82	+	static char  *hfname = NULL;            /* header file name */
83	+	static FILE  *hfp = NULL;               /* header file pointer */
84	+
85	+	static int  hres, vres;                 /* resolution for this frame */
86	+
87	+	extern char  *mktemp();
88	+
89		double  pixvalue();
90
91
92		quit(code)                      /* quit program */
93		int  code;
94		{
95	<	if (code || ralrm > 0)          /* report status */
95	>	if (code)                       /* report status */
96		report();
97	<
97	>	if (hfname != NULL) {           /* delete header file */
98	>	if (hfp != NULL)
99	>	fclose(hfp);
100	>	unlink(hfname);
101	>	}
102		exit(code);
103		}
104
#	Line 120 | Line 133 | report()               /* report progress */
133		#endif
134
135
136	+	openheader()                    /* save standard output to header file */
137	+	{
138	+	hfname = mktemp(HFTEMPLATE);
139	+	if (freopen(hfname, "w", stdout) == NULL) {
140	+	sprintf(errmsg, "cannot open header file \"%s\"", hfname);
141	+	error(SYSTEM, errmsg);
142	+	}
143	+	}
144	+
145	+
146	+	closeheader()                   /* done with header output */
147	+	{
148	+	if (hfname == NULL)
149	+	return;
150	+	if (fflush(stdout) == EOF || (hfp = fopen(hfname, "r")) == NULL)
151	+	error(SYSTEM, "error reopening header file");
152	+	}
153	+
154	+
155	+	dupheader()                     /* repeat header on standard output */
156	+	{
157	+	register int  c;
158	+
159	+	if (fseek(hfp, 0L, 0) < 0)
160	+	error(SYSTEM, "seek error on header file");
161	+	while ((c = getc(hfp)) != EOF)
162	+	putchar(c);
163	+	}
164	+
165	+
166	+	rpict(seq, pout, zout, prvr)                    /* generate image(s) */
167	+	int  seq;
168	+	char  *pout, *zout, *prvr;
169	+	/*
170	+	* If seq is greater than zero, then we will render a sequence of
171	+	* images based on view parameter strings read from the standard input.
172	+	* If pout is NULL, then all images will be sent to the standard ouput.
173	+	* If seq is greater than zero and prvr is an integer, then it is the
174	+	* frame number at which rendering should begin.  Preceeding view parameter
175	+	* strings will be skipped in the input.
176	+	* If pout and prvr are the same, prvr is renamed to avoid overwriting.
177	+	* Note that pout and zout should contain %d format specifications for
178	+	* sequenced file naming.
179	+	*/
180	+	{
181	+	extern char  *rindex(), *strncpy(), *strcat();
182	+	char  fbuf[128], fbuf2[128], *zf;
183	+	RESOLU  rs;
184	+	double  pa;
185	+	/* finished writing header */
186	+	closeheader();
187	+	/* check sampling */
188	+	if (psample < 1)
189	+	psample = 1;
190	+	else if (psample > MAXDIV) {
191	+	sprintf(errmsg, "pixel sampling reduced from %d to %d",
192	+	psample, MAXDIV);
193	+	error(WARNING, errmsg);
194	+	psample = MAXDIV;
195	+	}
196	+	/* get starting frame */
197	+	if (seq <= 0)
198	+	seq = 0;
199	+	else if (prvr != NULL && isint(prvr)) {
200	+	int  rn;                        /* skip to specified view */
201	+	if ((rn = atoi(prvr)) < seq)
202	+	error(USER, "recover frame less than start frame");
203	+	if (pout == NULL)
204	+	error(USER, "missing output file specification");
205	+	for ( ; seq < rn; seq++)
206	+	if (nextview(stdin) == EOF)
207	+	error(USER, "unexpected EOF on view input");
208	+	prvr = fbuf;                    /* mark for renaming */
209	+	}
210	+	if (pout != NULL) {
211	+	sprintf(fbuf, pout, seq);
212	+	if (!strcmp(prvr, fbuf)) {      /* rename recover file */
213	+	fbuf2[0] = '\0';
214	+	if ((prvr = rindex(fbuf, '/')) != NULL)
215	+	strncpy(fbuf2, fbuf, prvr-fbuf+1);
216	+	strcat(fbuf2, "rfXXXXXX");
217	+	prvr = mktemp(fbuf2);
218	+	if (rename(fbuf, prvr) < 0 && errno != ENOENT) {
219	+	sprintf(errmsg,
220	+	"cannot rename \"%s\" to \"%s\"",
221	+	fbuf, prvr);
222	+	error(SYSTEM, errmsg);
223	+	}
224	+	}
225	+	}
226	+	/* render sequence */
227	+	do {
228	+	if (seq && nextview(stdin) == EOF)
229	+	break;
230	+	if (pout != NULL) {
231	+	sprintf(fbuf, pout, seq);
232	+	if (freopen(fbuf, "w", stdout) == NULL) {
233	+	sprintf(errmsg,
234	+	"cannot open output file \"%s\"", fbuf);
235	+	error(SYSTEM, errmsg);
236	+	}
237	+	dupheader();
238	+	}
239	+	hres = hresolu; vres = vresolu; pa = pixaspect;
240	+	if (prvr != NULL)
241	+	if (viewfile(prvr, &ourview, &rs) <= 0
242	+	|| rs.or != PIXSTANDARD) {
243	+	sprintf(errmsg,
244	+	"cannot recover view parameters from \"%s\"", prvr);
245	+	error(WARNING, errmsg);
246	+	} else {
247	+	char  *err;
248	+	if ((err = setview(&ourview)) != NULL)
249	+	error(USER, err);
250	+	pa = 0.0;
251	+	hres = scanlen(&rs);
252	+	vres = numscans(&rs);
253	+	}
254	+	normaspect(viewaspect(&ourview), &pa, &hres, &vres);
255	+	if (seq) {
256	+	if (ralrm > 0) {
257	+	sprintf(errmsg, "starting frame %d\n", seq);
258	+	eputs(errmsg);
259	+	}
260	+	printf("FRAME=%d\n", seq);
261	+	}
262	+	fputs(VIEWSTR, stdout);
263	+	fprintview(&ourview, stdout);
264	+	putchar('\n');
265	+	if (pa < .99 || pa > 1.01)
266	+	fputaspect(pa, stdout);
267	+	fputformat(COLRFMT, stdout);
268	+	putchar('\n');
269	+	if (zout != NULL)
270	+	sprintf(zf=fbuf, zout, seq);
271	+	else
272	+	zf = NULL;
273	+	render(zf, prvr);
274	+	prvr = NULL;
275	+	} while (seq++);
276	+	}
277	+
278	+
279	+	nextview(fp)                            /* get next view from fp */
280	+	FILE  *fp;
281	+	{
282	+	char  linebuf[256], *err;
283	+
284	+	while (fgets(linebuf, sizeof(linebuf), fp) != NULL)
285	+	if (isview(linebuf) && sscanview(&ourview, linebuf) > 0) {
286	+	if ((err = setview(&ourview)) != NULL)
287	+	error(USER, err);
288	+	return(0);
289	+	}
290	+	return(EOF);
291	+	}
292	+
293	+
294		render(zfile, oldfile)                          /* render the scene */
295		char  *zfile, *oldfile;
296		{
#	Line 134 | Line 305 | char  *zfile, *oldfile;
305		COLOR  *colptr;
306		float  *zptr;
307		register int  i;
137	–	/* check sampling */
138	–	if (psample < 1)
139	–	psample = 1;
140	–	else if (psample > MAXDIV) {
141	–	sprintf(errmsg, "pixel sampling reduced from %d to %d",
142	–	psample, MAXDIV);
143	–	error(WARNING, errmsg);
144	–	psample = MAXDIV;
145	–	}
308		/* allocate scanlines */
309		for (i = 0; i <= psample; i++) {
310	<	scanbar[i] = (COLOR *)malloc(hresolu*sizeof(COLOR));
310	>	scanbar[i] = (COLOR *)malloc(hres*sizeof(COLOR));
311		if (scanbar[i] == NULL)
312		goto memerr;
313		}
314		hstep = (psample*140+49)/99;            /* quincunx sampling */
315		ystep = (psample*99+70)/140;
316		if (hstep > 2) {
317	<	i = hresolu/hstep + 2;
317	>	i = hres/hstep + 2;
318		if ((sampdens = malloc(i)) == NULL)
319		goto memerr;
320		while (i--)
#	Line 166 | Line 328 | char  *zfile, *oldfile;
328		error(SYSTEM, errmsg);
329		}
330		for (i = 0; i <= psample; i++) {
331	<	zbar[i] = (float *)malloc(hresolu*sizeof(float));
331	>	zbar[i] = (float *)malloc(hres*sizeof(float));
332		if (zbar[i] == NULL)
333		goto memerr;
334		}
#	Line 176 | Line 338 | char  *zfile, *oldfile;
338		zbar[i] = NULL;
339		}
340		/* write out boundaries */
341	<	fprtresolu(hresolu, vresolu, stdout);
341	>	fprtresolu(hres, vres, stdout);
342		/* recover file and compute first */
343		i = salvage(oldfile);
344		if (zfd != -1 && i > 0 &&
345	<	lseek(zfd, (long)i*hresolu*sizeof(float), 0) == -1)
345	>	lseek(zfd, (long)i*hres*sizeof(float), 0) == -1)
346		error(SYSTEM, "z file seek error in render");
347	<	pctdone = 100.0*i/vresolu;
347	>	pctdone = 100.0*i/vres;
348		if (ralrm > 0)                  /* report init stats */
349		report();
350		#ifndef  BSD
351		else
352		#endif
353		signal(SIGALRM, report);
354	<	ypos = vresolu-1 - i;
355	<	fillscanline(scanbar[0], zbar[0], sampdens, hresolu, ypos, hstep);
354	>	ypos = vres-1 - i;
355	>	fillscanline(scanbar[0], zbar[0], sampdens, hres, ypos, hstep);
356		/* compute scanlines */
357		for (ypos -= ystep; ypos > -ystep; ypos -= ystep) {
358		/* bottom adjust? */
#	Line 206 | Line 368 | char  *zfile, *oldfile;
368		zbar[0] = zptr;
369		/* fill base line */
370		fillscanline(scanbar[0], zbar[0], sampdens,
371	<	hresolu, ypos, hstep);
371	>	hres, ypos, hstep);
372		/* fill bar */
373	<	fillscanbar(scanbar, zbar, hresolu, ypos, ystep);
373	>	fillscanbar(scanbar, zbar, hres, ypos, ystep);
374		/* write it out */
375		#ifndef  BSD
376		signal(SIGALRM, SIG_IGN);       /* don't interrupt writes */
377		#endif
378		for (i = ystep; i > 0; i--) {
379		if (zfd != -1 && write(zfd, (char *)zbar[i],
380	<	hresolu*sizeof(float))
381	<	< hresolu*sizeof(float))
380	>	hres*sizeof(float))
381	>	< hres*sizeof(float))
382		goto writerr;
383	<	if (fwritescan(scanbar[i], hresolu, stdout) < 0)
383	>	if (fwritescan(scanbar[i], hres, stdout) < 0)
384		goto writerr;
385		}
386		if (fflush(stdout) == EOF)
387		goto writerr;
388		/* record progress */
389	<	pctdone = 100.0*(vresolu-1-ypos)/vresolu;
389	>	pctdone = 100.0*(vres-1-ypos)/vres;
390		if (ralrm > 0 && time((long *)0) >= tlastrept+ralrm)
391		report();
392		#ifndef  BSD
#	Line 235 | Line 397 | char  *zfile, *oldfile;
397		/* clean up */
398		signal(SIGALRM, SIG_IGN);
399		if (zfd != -1) {
400	<	if (write(zfd, (char *)zbar[0], hresolu*sizeof(float))
401	<	< hresolu*sizeof(float))
400	>	if (write(zfd, (char *)zbar[0], hres*sizeof(float))
401	>	< hres*sizeof(float))
402		goto writerr;
403		if (close(zfd) == -1)
404		goto writerr;
405		for (i = 0; i <= psample; i++)
406		free((char *)zbar[i]);
407		}
408	<	fwritescan(scanbar[0], hresolu, stdout);
408	>	fwritescan(scanbar[0], hres, stdout);
409		if (fflush(stdout) == EOF)
410		goto writerr;
411		for (i = 0; i <= psample; i++)
#	Line 388 | Line 550 | int  x, y;                     /* pixel position */
550		static RAY  thisray;
551
552		if (viewray(thisray.rorg, thisray.rdir, &ourview,
553	<	(x+pixjitter())/hresolu, (y+pixjitter())/vresolu) < 0) {
553	>	(x+pixjitter())/hres, (y+pixjitter())/vres) < 0) {
554		setcolor(col, 0.0, 0.0, 0.0);
555		return(0.0);
556		}
557
558		rayorigin(&thisray, NULL, PRIMARY, 1.0);
559
560	<	samplendx = pixnumber(x,y,hresolu,vresolu);     /* set pixel index */
560	>	samplendx = pixnumber(x,y,hres,vres);   /* set pixel index */
561
562		rayvalue(&thisray);                     /* trace ray */
563
#	Line 431 | Line 593 | char  *oldfile;
593		return(0);
594		}
595
596	<	if (x != hresolu || y != vresolu) {
596	>	if (x != hres || y != vres) {
597		sprintf(errmsg, "resolution mismatch in recover file \"%s\"",
598		oldfile);
599		error(USER, errmsg);
600		}
601
602	<	scanline = (COLR *)malloc(hresolu*sizeof(COLR));
602	>	scanline = (COLR *)malloc(hres*sizeof(COLR));
603		if (scanline == NULL)
604		error(SYSTEM, "out of memory in salvage");
605	<	for (y = 0; y < vresolu; y++) {
606	<	if (freadcolrs(scanline, hresolu, fp) < 0)
605	>	for (y = 0; y < vres; y++) {
606	>	if (freadcolrs(scanline, hres, fp) < 0)
607		break;
608	<	if (fwritecolrs(scanline, hresolu, stdout) < 0)
608	>	if (fwritecolrs(scanline, hres, stdout) < 0)
609		goto writerr;
610		}
611		if (fflush(stdout) == EOF)

Diff Legend

–	Removed lines
+	Added lines
<	Changed lines (old)
>	Changed lines (new)