ViewVC Help

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

root/Development/ray/src/util/genBSDF.pl

Comparing ray/src/util/genBSDF.pl (file contents):
Revision 2.18 by greg, Wed Jun 1 22:50:24 2011 UTC vs.
Revision 2.22 by greg, Fri Jun 24 00:41:51 2011 UTC

#	Line 13 | Line 13 | sub userror {
13		}
14		my $td = mkdtemp("/tmp/genBSDF.XXXXXX");
15		chomp $td;
16	+	my @savedARGV = @ARGV;
17		my $tensortree = 0;
18		my $ttlog2 = 4;
19		my $nsamp = 1000;
#	Line 22 | Line 23 | my $geout = 1;
23		my $nproc = 1;
24		my $doforw = 0;
25		my $doback = 1;
26	+	my $gunit = "Meter";
27		my @dim;
28		# Get options
29		while ($#ARGV >= 0) {
#	Line 32 | Line 34 | while ($#ARGV >= 0) {
34		shift @ARGV;
35		} elsif ("$ARGV[0]" =~ /^[-+]g/) {
36		$geout = ("$ARGV[0]" =~ /^\+/);
37	+	$gunit = $ARGV[1];
38	+	if ($gunit !~ /^(?i)(meter|foot|inch|centimeter|millimeter)$/) {
39	+	die "Illegal geometry unit '$gunit': must be meter, foot, inch, centimeter, or millimeter\n";
40	+	}
41	+	shift @ARGV;
42		} elsif ("$ARGV[0]" =~ /^[-+]f/) {
43		$doforw = ("$ARGV[0]" =~ /^\+/);
44		} elsif ("$ARGV[0]" =~ /^[-+]b/) {
#	Line 98 | Line 105 | die "Could not compile scene\n" if ( $? );
105		print
106		'<?xml version="1.0" encoding="UTF-8"?>
107		<WindowElement xmlns="http://windows.lbl.gov" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://windows.lbl.gov/BSDF-v1.4.xsd">
108	<	<WindowElementType>System</WindowElementType>
108	>	';
109	>	print "<!-- File produced by: genBSDF @savedARGV -->\n";
110	>	print
111	>	'<WindowElementType>System</WindowElementType>
112		<Optical>
113		<Layer>
114		<Material>
115		<Name>Name</Name>
116		<Manufacturer>Manufacturer</Manufacturer>
117		';
118	<	printf "\t\t<Thickness unit=\"Meter\">%.3f</Thickness>\n", $dim[5] - $dim[4];
119	<	printf "\t\t<Width unit=\"Meter\">%.3f</Width>\n", $dim[1] - $dim[0];
120	<	printf "\t\t<Height unit=\"Meter\">%.3f</Height>\n", $dim[3] - $dim[2];
118	>	printf "\t\t<Thickness unit=\"$gunit\">%.3f</Thickness>\n", $dim[5] - $dim[4];
119	>	printf "\t\t<Width unit=\"$gunit\">%.3f</Width>\n", $dim[1] - $dim[0];
120	>	printf "\t\t<Height unit=\"$gunit\">%.3f</Height>\n", $dim[3] - $dim[2];
121		print "\t\t<DeviceType>Integral</DeviceType>\n";
122		# Output MGF description if requested
123		if ( $geout ) {
124	<	print "\t\t<Geometry format=\"MGF\" unit=\"Meter\">\n";
124	>	print "\t\t<Geometry format=\"MGF\" unit=\"$gunit\">\n";
125		printf "xf -t %.6f %.6f 0\n", -($dim[0]+$dim[1])/2, -($dim[2]+$dim[3])/2;
126		open(MGFSCN, "< $mgfscn");
127		while (<MGFSCN>) { print $_; }
#	Line 317 | Line 327 | print
327		<WavelengthDataBlock>
328		';
329		print "\t\t\t<WavelengthDataDirection>Reflection $side</WavelengthDataDirection>\n";
330	<	print '
331	<	<AngleBasis>LBNL/Shirley-Chiu</AngleBasis>
330	>	print
331	>	'                       <AngleBasis>LBNL/Shirley-Chiu</AngleBasis>
332		<ScatteringDataType>BRDF</ScatteringDataType>
333		<ScatteringData>
334		';
#	Line 398 | Line 408 | my $rtcmd = "rtcontrib $rtargs -h -ff -fo -n $nproc -c
408		"-o '$td/%s.flt' -m $fmodnm -m $bmodnm $octree";
409		my $rccmd = "rcalc -e '$tcal' " .
410		"-e 'mod(n,d):n-floor(n/d)*d' -e 'Kbin=mod(recno-.999,$ndiv)' " .
411	<	q{-if3 -e '$1=(0.265*$1+0.670*$2+0.065*$3)/KprojOmega'};
411	>	q{-if3 -e 'oval=(0.265*$1+0.670*$2+0.065*$3)/KprojOmega' } .
412	>	q[-o '${  oval  },'];
413		if ( $doforw ) {
414		$cmd = "cnt $ndiv $ny $nx | rcalc -of -e '$tcal' " .
415		"-e 'xp=(\$3+rand(.12*recno+288))*(($dim[1]-$dim[0])/$nx)+$dim[0]' " .

Diff Legend

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