ViewVC Help

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

root/Development/ray/src/cv/bsdf2rad.c

Comparing ray/src/cv/bsdf2rad.c (file contents):
Revision 2.16 by greg, Sun Apr 9 21:39:26 2017 UTC vs.
Revision 2.19 by greg, Mon Apr 10 01:31:37 2017 UTC

#	Line 13 | Line 13 | static const char RCSid[] = "$Id$";
13		#include "resolu.h"
14		#include "bsdfrep.h"
15
16	<	#define NINCIDENT       25              /* number of samples/hemisphere */
16	>	#define NINCIDENT       37              /* number of samples/hemisphere */
17
18		#define GRIDSTEP        2               /* our grid step size */
19		#define SAMPRES         (GRIDRES/GRIDSTEP)
#	Line 50 | Line 50 | const double   sph_xoffset = 15.;
50		char    *progname;
51
52		/* Get Fibonacci sphere vector (0 to NINCIDENT-1) */
53	<	static void
53	>	static RREAL *
54		get_ivector(FVECT iv, int i)
55		{
56	–	const double    zstep = 1./(2*NINCIDENT-1);
56		const double    phistep = PI*(3. - 2.236067978);
57		double          r;
58
59	<	iv[2] = 1. - (i+.5)*zstep;
59	>	iv[2] = 1. - (i+.5)*(1./NINCIDENT);
60		r = sqrt(1. - iv[2]*iv[2]);
61		iv[0] = r * cos((i+1.)*phistep);
62		iv[1] = r * sin((i+1.)*phistep);
63	+
64	+	return(iv);
65		}
66
67	+	/* Convert incident vector into sphere position */
68	+	static RREAL *
69	+	cvt_sposition(FVECT sp, const FVECT iv, int inc_side)
70	+	{
71	+	sp[0] = -iv[0]*sph_rad - inc_side*sph_xoffset;
72	+	sp[1] = -iv[1]*sph_rad;
73	+	sp[2] = iv[2]*sph_rad;
74	+
75	+	return(sp);
76	+	}
77	+
78		/* Get temporary file name */
79		static char *
80		tfile_name(const char *prefix, const char *suffix, int i)
#	Line 249 | Line 261 | build_wRBF(void)
261		RBFNODE *rbf;
262		get_ivector(ivec, i);
263		if (input_orient < 0) {
264	<	ivec[0] = -ivec[0]; ivec[1] = -ivec[1]; ivec[2] = -ivec[2];
264	>	ivec[0] = -ivec[0]; ivec[2] = -ivec[2];
265		}
266		rbf = advect_rbf(ivec, 15000);
267		if (!plotRBF(tfile_name(pref, dsuffix, i), rbf))
#	Line 268 | Line 280 | put_mirror_arrow(const FVECT ivec, int inc_side)
280		FVECT           origin, refl;
281		int             i;
282
283	<	for (i = 3; i--; ) origin[i] = ivec[i]*sph_rad;
272	<	origin[0] -= inc_side*sph_xoffset;
283	>	cvt_sposition(origin, ivec, inc_side);
284
285	<	refl[0] = 2.*ivec[2]*ivec[0];
286	<	refl[1] = 2.*ivec[2]*ivec[1];
285	>	refl[0] = -2.*ivec[2]*ivec[0];
286	>	refl[1] = -2.*ivec[2]*ivec[1];
287		refl[2] = 2.*ivec[2]*ivec[2] - 1.;
288
289		printf("\n# Mirror arrow\n");
#	Line 308 | Line 319 | put_trans_arrow(const FVECT ivec, int inc_side)
319		FVECT           origin;
320		int             i;
321
322	<	for (i = 3; i--; ) origin[i] = ivec[i]*sph_rad;
312	<	origin[0] -= inc_side*sph_xoffset;
322	>	cvt_sposition(origin, ivec, inc_side);
323
324		printf("\n# Transmission arrow\n");
325		printf("\narrow_mat cylinder trans_dir\n0\n0\n7");
#	Line 365 | Line 375 | static int
375		put_BSDFs(void)
376		{
377		const double    scalef = bsdf_rad/(log10(overall_max) - min_log10);
378	<	FVECT           ivec;
378	>	FVECT           ivec, sorg;
379		RREAL           vMtx[3][3];
380		char            *fname;
381		char            cmdbuf[256];
#	Line 379 | Line 389 | put_BSDFs(void)
389		printf("\nvoid glow arrow_glow\n0\n0\n4 1 0 1 0\n");
390		printf("\nvoid mixfunc arrow_mat\n4 arrow_glow void .5 .\n0\n0\n");
391
392	<	if (front_comp & SDsampR)
392	>	if (front_comp & SDsampR)                       /* front reflection */
393		for (i = 0; i < NINCIDENT; i++) {
394		get_ivector(ivec, i);
395		put_mirror_arrow(ivec, 1);
396	+	cvt_sposition(sorg, ivec, 1);
397	+	ivec[0] = -ivec[0]; ivec[1] = -ivec[1]; /* normal */
398		sprintf(xfargs, "-s %f -t %f %f %f", bsdf_rad,
399	<	ivec[0]*sph_rad - sph_xoffset,
388	<	ivec[1]*sph_rad, ivec[2]*sph_rad);
399	>	sorg[0], sorg[1], sorg[2]);
400		nxfa = 6;
401		printf("\nvoid colorfunc scale_pat\n");
402		printf("%d bsdf_red bsdf_grn bsdf_blu bsdf2rad.cal\n\t%s\n0\n0\n",
#	Line 394 | Line 405 | put_BSDFs(void)
405		SDcompXform(vMtx, ivec, Yaxis);
406		nxfa = addrot(xfargs, vMtx[0], vMtx[1], vMtx[2]);
407		sprintf(xfargs+strlen(xfargs), " -s %f -t %f %f %f",
408	<	scalef, ivec[0]*sph_rad - sph_xoffset,
398	<	ivec[1]*sph_rad, ivec[2]*sph_rad);
408	>	scalef, sorg[0], sorg[1], sorg[2]);
409		nxfa += 6;
410		fname = tfile_name(frpref, dsuffix, i);
411	<	sprintf(cmdbuf, "gensurf scale_mat %s%d %s %s %s %d %d | xform -mx -my %s",
411	>	sprintf(cmdbuf, "gensurf scale_mat %s%d %s %s %s %d %d | xform %s",
412		frpref, i+1, fname, fname, fname, SAMPRES-1, SAMPRES-1,
413		xfargs);
414		if (!run_cmd(cmdbuf))
415		return(0);
416		}
417	<	if (front_comp & SDsampT)
417	>	if (front_comp & SDsampT)                       /* front transmission */
418		for (i = 0; i < NINCIDENT; i++) {
419		get_ivector(ivec, i);
420		put_trans_arrow(ivec, 1);
421	+	cvt_sposition(sorg, ivec, 1);
422	+	ivec[0] = -ivec[0]; ivec[1] = -ivec[1]; /* normal */
423		sprintf(xfargs, "-s %f -t %f %f %f", bsdf_rad,
424	<	ivec[0]*sph_rad - sph_xoffset,
413	<	ivec[1]*sph_rad, ivec[2]*sph_rad);
424	>	sorg[0], sorg[1], sorg[2]);
425		nxfa = 6;
426		printf("\nvoid colorfunc scale_pat\n");
427		printf("%d bsdf_red bsdf_grn bsdf_blu bsdf2rad.cal\n\t%s\n0\n0\n",
#	Line 419 | Line 430 | put_BSDFs(void)
430		SDcompXform(vMtx, ivec, Yaxis);
431		nxfa = addrot(xfargs, vMtx[0], vMtx[1], vMtx[2]);
432		sprintf(xfargs+strlen(xfargs), " -s %f -t %f %f %f",
433	<	scalef, ivec[0]*sph_rad - sph_xoffset,
423	<	ivec[1]*sph_rad, ivec[2]*sph_rad);
433	>	scalef, sorg[0], sorg[1], sorg[2]);
434		nxfa += 6;
435		fname = tfile_name(ftpref, dsuffix, i);
436	<	sprintf(cmdbuf, "gensurf scale_mat %s%d %s %s %s %d %d | xform -I -mx -my %s",
436	>	sprintf(cmdbuf, "gensurf scale_mat %s%d %s %s %s %d %d | xform -I %s",
437		ftpref, i+1, fname, fname, fname, SAMPRES-1, SAMPRES-1,
438		xfargs);
439		if (!run_cmd(cmdbuf))
440		return(0);
441		}
442	<	if (back_comp & SDsampR)
442	>	if (back_comp & SDsampR)                        /* rear reflection */
443		for (i = 0; i < NINCIDENT; i++) {
444		get_ivector(ivec, i);
445		put_mirror_arrow(ivec, -1);
446	<	fname = tfile_name(brpref, dsuffix, i);
446	>	cvt_sposition(sorg, ivec, -1);
447	>	ivec[0] = -ivec[0]; ivec[1] = -ivec[1]; /* normal */
448		sprintf(xfargs, "-s %f -t %f %f %f", bsdf_rad,
449	<	ivec[0]*sph_rad + sph_xoffset,
439	<	ivec[1]*sph_rad, ivec[2]*sph_rad);
449	>	sorg[0], sorg[1], sorg[2]);
450		nxfa = 6;
451		printf("\nvoid colorfunc scale_pat\n");
452		printf("%d bsdf_red bsdf_grn bsdf_blu bsdf2rad.cal\n\t%s\n0\n0\n",
#	Line 445 | Line 455 | put_BSDFs(void)
455		SDcompXform(vMtx, ivec, Yaxis);
456		nxfa = addrot(xfargs, vMtx[0], vMtx[1], vMtx[2]);
457		sprintf(xfargs+strlen(xfargs), " -s %f -t %f %f %f",
458	<	scalef, ivec[0]*sph_rad + sph_xoffset,
449	<	ivec[1]*sph_rad, ivec[2]*sph_rad);
458	>	scalef, sorg[0], sorg[1], sorg[2]);
459		nxfa += 6;
460		fname = tfile_name(brpref, dsuffix, i);
461	<	sprintf(cmdbuf, "gensurf scale_mat %s%d %s %s %s %d %d | xform -I -ry 180 -mx -my %s",
461	>	sprintf(cmdbuf, "gensurf scale_mat %s%d %s %s %s %d %d | xform -I -ry 180 %s",
462		brpref, i+1, fname, fname, fname, SAMPRES-1, SAMPRES-1,
463		xfargs);
464		if (!run_cmd(cmdbuf))
465		return(0);
466		}
467	<	if (back_comp & SDsampT)
467	>	if (back_comp & SDsampT)                        /* rear transmission */
468		for (i = 0; i < NINCIDENT; i++) {
469		get_ivector(ivec, i);
470		put_trans_arrow(ivec, -1);
471	<	fname = tfile_name(btpref, dsuffix, i);
471	>	cvt_sposition(sorg, ivec, -1);
472	>	ivec[0] = -ivec[0]; ivec[1] = -ivec[1]; /* normal */
473		sprintf(xfargs, "-s %f -t %f %f %f", bsdf_rad,
474	<	ivec[0]*sph_rad + sph_xoffset,
465	<	ivec[1]*sph_rad, ivec[2]*sph_rad);
474	>	sorg[0], sorg[1], sorg[2]);
475		nxfa = 6;
476		printf("\nvoid colorfunc scale_pat\n");
477		printf("%d bsdf_red bsdf_grn bsdf_blu bsdf2rad.cal\n\t%s\n0\n0\n",
#	Line 471 | Line 480 | put_BSDFs(void)
480		SDcompXform(vMtx, ivec, Yaxis);
481		nxfa = addrot(xfargs, vMtx[0], vMtx[1], vMtx[2]);
482		sprintf(xfargs+strlen(xfargs), " -s %f -t %f %f %f",
483	<	scalef, ivec[0]*sph_rad + sph_xoffset,
475	<	ivec[1]*sph_rad, ivec[2]*sph_rad);
483	>	scalef, sorg[0], sorg[1], sorg[2]);
484		nxfa += 6;
485		fname = tfile_name(btpref, dsuffix, i);
486	<	sprintf(cmdbuf, "gensurf scale_mat %s%d %s %s %s %d %d | xform -ry 180 -mx -my %s",
486	>	sprintf(cmdbuf, "gensurf scale_mat %s%d %s %s %s %d %d | xform -ry 180 %s",
487		btpref, i+1, fname, fname, fname, SAMPRES-1, SAMPRES-1,
488		xfargs);
489		if (!run_cmd(cmdbuf))
#	Line 524 | Line 532 | static void
532		put_source(void)
533		{
534		printf("\n# Overhead parallel light source\n");
535	<	printf("\nvoid light bright\n0\n0\n3 1000 1000 1000\n");
535	>	printf("\nvoid light bright\n0\n0\n3 2000 2000 2000\n");
536		printf("\nbright source light\n0\n0\n4 0 0 1 2\n");
537		printf("\n# Material used for labels\n");
538		printf("\nvoid trans vellum\n0\n0\n7 1 1 1 0 0 .5 0\n");
#	Line 544 | Line 552 | put_hemispheres(void)
552		-1.01*sph_rad - sph_xoffset, -1.01*sph_rad, -1.01*sph_rad);
553		printf("\nvoid brighttext front_text\n3 helvet.fnt . FRONT\n0\n");
554		printf("12\n\t%f %f 0\n\t%f 0 0\n\t0 %f 0\n\t.01 1 -.1\n",
555	<	-.22*sph_rad - sph_xoffset, -1.2*sph_rad,
555	>	-.22*sph_rad - sph_xoffset, -1.4*sph_rad,
556		.35/5.*sph_rad, -1.6*.35/5.*sph_rad);
557		printf("\nfront_text alias front_label_mat vellum\n");
558		printf("\nfront_label_mat polygon front_label\n0\n0\n12");
559		printf("\n\t%f %f 0\n\t%f %f 0\n\t%f %f 0\n\t%f %f 0\n",
560	<	-.25*sph_rad - sph_xoffset, -1.1*sph_rad,
561	<	-.25*sph_rad - sph_xoffset, (-1.2-1.6*.35/5.-.1)*sph_rad,
562	<	.25*sph_rad - sph_xoffset, (-1.2-1.6*.35/5.-.1)*sph_rad,
563	<	.25*sph_rad - sph_xoffset, -1.1*sph_rad );
560	>	-.25*sph_rad - sph_xoffset, -1.3*sph_rad,
561	>	-.25*sph_rad - sph_xoffset, (-1.4-1.6*.35/5.-.1)*sph_rad,
562	>	.25*sph_rad - sph_xoffset, (-1.4-1.6*.35/5.-.1)*sph_rad,
563	>	.25*sph_rad - sph_xoffset, -1.3*sph_rad );
564		}
565		if (back_comp) {
566		printf("\n%s bubble Back\n0\n0\n4 %f 0 0 %f\n",
#	Line 562 | Line 570 | put_hemispheres(void)
570		-1.01*sph_rad + sph_xoffset, -1.01*sph_rad, -1.01*sph_rad);
571		printf("\nvoid brighttext back_text\n3 helvet.fnt . BACK\n0\n");
572		printf("12\n\t%f %f 0\n\t%f 0 0\n\t0 %f 0\n\t.01 1 -.1\n",
573	<	-.22*sph_rad + sph_xoffset, -1.2*sph_rad,
573	>	-.22*sph_rad + sph_xoffset, -1.4*sph_rad,
574		.35/4.*sph_rad, -1.6*.35/4.*sph_rad);
575		printf("\nback_text alias back_label_mat vellum\n");
576		printf("\nback_label_mat polygon back_label\n0\n0\n12");
577		printf("\n\t%f %f 0\n\t%f %f 0\n\t%f %f 0\n\t%f %f 0\n",
578	<	-.25*sph_rad + sph_xoffset, -1.1*sph_rad,
579	<	-.25*sph_rad + sph_xoffset, (-1.2-1.6*.35/4.-.1)*sph_rad,
580	<	.25*sph_rad + sph_xoffset, (-1.2-1.6*.35/4.-.1)*sph_rad,
581	<	.25*sph_rad + sph_xoffset, -1.1*sph_rad );
578	>	-.25*sph_rad + sph_xoffset, -1.3*sph_rad,
579	>	-.25*sph_rad + sph_xoffset, (-1.4-1.6*.35/4.-.1)*sph_rad,
580	>	.25*sph_rad + sph_xoffset, (-1.4-1.6*.35/4.-.1)*sph_rad,
581	>	.25*sph_rad + sph_xoffset, -1.3*sph_rad );
582		}
583		}
584
#	Line 582 | Line 590 | put_scale(void)
590		const double    leg_width = 2.*.75*(sph_xoffset - sph_rad);
591		const double    leg_height = 2.*sph_rad;
592		const int       text_lines = 6;
593	<	const int       text_digits = 7;
593	>	const int       text_digits = 8;
594		char            fmt[16];
595		int             i;
596
597		printf("\n# BSDF legend with falsecolor scale\n");
598		printf("\nvoid colorfunc lscale\n10 sca_red(Py) sca_grn(Py) sca_blu(Py)");
599		printf("\n\tbsdf2rad.cal -s %f -t 0 %f 0\n0\n0\n", leg_height, -.5*leg_height);
600	<	sprintf(fmt, "%%%dg", text_digits-2);
600	>	sprintf(fmt, "%%.%df", text_digits-3);
601		for (i = 0; i < text_lines; i++) {
602		char    vbuf[16];
603		sprintf(vbuf, fmt, pow(10., (i+.5)/text_lines*(max_log10-min_log10)+min_log10));
604		printf("\nlscale brighttext lscale\n");
605		printf("3 helvet.fnt . %s\n0\n12\n", vbuf);
606		printf("\t%f %f 0\n", -.45*leg_width, ((i+.9)/text_lines-.5)*leg_height);
607	<	printf("\t%f 0 0\n", .9*leg_width/strlen(vbuf));
607	>	printf("\t%f 0 0\n", .8*leg_width/strlen(vbuf));
608		printf("\t0 %f 0\n", -.9/text_lines*leg_height);
609		printf("\t.01 1 -.1\n");
610		}
#	Line 623 | Line 631 | put_scale(void)
631		return;
632		printf("\nvoid brighttext BSDFname\n3 helvet.fnt . \"%s\"\n0\n12\n", bsdf_name);
633		printf("\t%f %f 0\n", -.95*leg_width, -.6*leg_height);
634	<	printf("\t%f 0 0\n", .4/strlen(bsdf_name)*leg_width);
634	>	printf("\t%f 0 0\n", 1.8/strlen(bsdf_name)*leg_width);
635		printf("\t0 %f 0\n", -.1*leg_height);
636		printf("\t.01 1 -.1\n");
637		printf("\nBSDFname alias name_mat vellum\n");

Diff Legend

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