24 |
|
my $doforw = 0; |
25 |
|
my $doback = 1; |
26 |
|
my $pctcull = 90; |
27 |
< |
my $gunit = "Meter"; |
27 |
> |
my $gunit = "meter"; |
28 |
|
my @dim; |
29 |
|
# Get options |
30 |
|
while ($#ARGV >= 0) { |
108 |
|
print "<!-- File produced by: genBSDF @savedARGV -->\n"; |
109 |
|
print |
110 |
|
'<WindowElementType>System</WindowElementType> |
111 |
+ |
<FileType>BSDF</FileType> |
112 |
|
<Optical> |
113 |
|
<Layer> |
114 |
|
<Material> |
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"; |
121 |
> |
print "\t\t<DeviceType>Other</DeviceType>\n"; |
122 |
> |
print " </Material>\n"; |
123 |
|
# Output MGF description if requested |
124 |
|
if ( $geout ) { |
125 |
< |
print "\t\t<Geometry format=\"MGF\" unit=\"$gunit\">\n"; |
125 |
> |
print "\t<Geometry format=\"MGF\">\n"; |
126 |
> |
print "\t\t<MGFblock unit=\"$gunit\">\n"; |
127 |
|
printf "xf -t %.6f %.6f 0\n", -($dim[0]+$dim[1])/2, -($dim[2]+$dim[3])/2; |
128 |
|
open(MGFSCN, "< $mgfscn"); |
129 |
|
while (<MGFSCN>) { print $_; } |
130 |
|
close MGFSCN; |
131 |
|
print "xf\n"; |
132 |
< |
print "\t\t</Geometry>\n"; |
132 |
> |
print "</MGFblock>\n"; |
133 |
> |
print "\t</Geometry>\n"; |
134 |
|
} |
131 |
– |
print " </Material>\n"; |
135 |
|
# Set up surface sampling |
136 |
|
my $nx = int(sqrt($nsamp*($dim[1]-$dim[0])/($dim[3]-$dim[2])) + .5); |
137 |
|
my $ny = int($nsamp/$nx + .5); |
277 |
|
q{-e '$1=(0.265*$1+0.670*$2+0.065*$3)/Omega' }; |
278 |
|
if ($pctcull >= 0) { |
279 |
|
$cmd .= "-of $td/" . ($bmodnm,$fmodnm)[$forw] . ".flt " . |
280 |
< |
"| rttree_reduce -a -h -ff -t $pctcull -r $tensortree -g $ttlog2"; |
280 |
> |
"| rttree_reduce -h -ff -t $pctcull -r $tensortree -g $ttlog2"; |
281 |
|
system "$cmd" || die "Failure running rttree_reduce"; |
282 |
|
} else { |
283 |
|
$cmd .= "$td/" . ($bmodnm,$fmodnm)[$forw] . ".flt"; |
306 |
|
print "\t\t\t<WavelengthDataDirection>Reflection $side</WavelengthDataDirection>\n"; |
307 |
|
print |
308 |
|
' <AngleBasis>LBNL/Shirley-Chiu</AngleBasis> |
309 |
< |
<ScatteringDataType>BRDF</ScatteringDataType> |
309 |
> |
<ScatteringDataType>BTDF</ScatteringDataType> |
310 |
|
<ScatteringData> |
311 |
|
'; |
312 |
|
$cmd = "rcalc -if3 -e 'Omega:PI/($ns*$ns)' " . |
543 |
|
<WavelengthDataDirection>Reflection Front</WavelengthDataDirection> |
544 |
|
<ColumnAngleBasis>LBNL/Klems Full</ColumnAngleBasis> |
545 |
|
<RowAngleBasis>LBNL/Klems Full</RowAngleBasis> |
546 |
< |
<ScatteringDataType>BRDF</ScatteringDataType> |
546 |
> |
<ScatteringDataType>BTDF</ScatteringDataType> |
547 |
|
<ScatteringData> |
548 |
|
'; |
549 |
|
# Output front reflection (transposed order) |
593 |
|
<WavelengthDataDirection>Reflection Back</WavelengthDataDirection> |
594 |
|
<ColumnAngleBasis>LBNL/Klems Full</ColumnAngleBasis> |
595 |
|
<RowAngleBasis>LBNL/Klems Full</RowAngleBasis> |
596 |
< |
<ScatteringDataType>BRDF</ScatteringDataType> |
596 |
> |
<ScatteringDataType>BTDF</ScatteringDataType> |
597 |
|
<ScatteringData> |
598 |
|
'; |
599 |
|
# Output back reflection (transposed order) |