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) { |
71 |
|
} |
72 |
|
# Check that we're actually being asked to do something |
73 |
|
die "Must have at least one of +forward or +backward\n" if (!$doforw && !$doback); |
74 |
– |
# Name our own persist file? |
75 |
– |
my $persistfile; |
76 |
– |
if ( $tensortree && $nproc > 1 && "$rtargs" !~ /-PP /) { |
77 |
– |
$persistfile = "$td/pfile.txt"; |
78 |
– |
$rtargs = "-PP $persistfile $rtargs"; |
79 |
– |
} |
74 |
|
# Get scene description and dimensions |
75 |
|
my $radscn = "$td/device.rad"; |
76 |
|
my $mgfscn = "$td/device.mgf"; |
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 |
|
} |
137 |
– |
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); |
136 |
> |
my $nx = int(sqrt($nsamp*($dim[1]-$dim[0])/($dim[3]-$dim[2])) + 1); |
137 |
> |
my $ny = int($nsamp/$nx + 1); |
138 |
|
$nsamp = $nx * $ny; |
139 |
|
my $ns = 2**$ttlog2; |
140 |
|
my (@pdiv, $disk2sq, $sq2disk, $tcal, $kcal); |
155 |
|
|
156 |
|
#-------------- End of main program segment --------------# |
157 |
|
|
161 |
– |
#++++++++++++++ Kill persistent rtrace +++++++++++++++++++# |
162 |
– |
sub persist_end { |
163 |
– |
if ( $persistfile && open(PFI, "< $persistfile") ) { |
164 |
– |
while (<PFI>) { |
165 |
– |
s/^[^ ]* //; |
166 |
– |
kill('ALRM', $_); |
167 |
– |
last; |
168 |
– |
} |
169 |
– |
close PFI; |
170 |
– |
} |
171 |
– |
} |
172 |
– |
|
158 |
|
#++++++++++++++ Tensor tree BSDF generation ++++++++++++++# |
159 |
|
sub do_tree_bsdf { |
175 |
– |
# Get sampling rate and subdivide task |
176 |
– |
my $ns2 = $ns; |
177 |
– |
$ns2 /= 2 if ( $tensortree == 3 ); |
178 |
– |
my $nsplice = $nproc; |
179 |
– |
$nsplice *= 10 if ($nproc > 1); |
180 |
– |
$nsplice = $ns2 if ($nsplice > $ns2); |
181 |
– |
$nsplice = 999 if ($nsplice > 999); |
182 |
– |
@pdiv = (0, int($ns2/$nsplice)); |
183 |
– |
my $nrem = $ns2 % $nsplice; |
184 |
– |
for (my $i = 1; $i < $nsplice; $i++) { |
185 |
– |
my $nv = $pdiv[$i] + $pdiv[1]; |
186 |
– |
++$nv if ( $nrem-- > 0 ); |
187 |
– |
push @pdiv, $nv; |
188 |
– |
} |
189 |
– |
die "Script error 1" if ($pdiv[-1] != $ns2); |
160 |
|
# Shirley-Chiu mapping from unit square to disk |
161 |
|
$sq2disk = ' |
162 |
|
in_square_a = 2*in_square_x - 1; |
181 |
|
norm_radians(p) : if(-p - PI/4, p + 2*PI, p); |
182 |
|
in_disk_r = .999995*sqrt(Dx*Dx + Dy*Dy); |
183 |
|
in_disk_phi = norm_radians(atan2(Dy, Dx)); |
184 |
< |
in_disk_rgn = floor((in_disk_phi + PI/4)/(PI/2)) + 1; |
184 |
> |
in_disk_rgn = floor((.999995*in_disk_phi + PI/4)/(PI/2)) + 1; |
185 |
|
out_square_a = select(in_disk_rgn, |
186 |
|
in_disk_r, |
187 |
|
(PI/2 - in_disk_phi)*in_disk_r/(PI/4), |
199 |
|
print "\t<DataDefinition>\n"; |
200 |
|
print "\t\t<IncidentDataStructure>TensorTree$tensortree</IncidentDataStructure>\n"; |
201 |
|
print "\t</DataDefinition>\n"; |
202 |
< |
# Fork parallel rtcontrib processes to compute each side |
203 |
< |
my $npleft = $nproc; |
204 |
< |
if ( $doback ) { |
205 |
< |
for (my $splice = 0; $splice < $nsplice; $splice++) { |
206 |
< |
if (! $npleft ) { |
237 |
< |
wait(); |
238 |
< |
die "rtcontrib process reported error" if ( $? ); |
239 |
< |
$npleft++; |
240 |
< |
} |
241 |
< |
bg_tree_rtcontrib(0, $splice); |
242 |
< |
$npleft--; |
243 |
< |
} |
244 |
< |
while (wait() >= 0) { |
245 |
< |
die "rtcontrib process reported error" if ( $? ); |
246 |
< |
$npleft++; |
247 |
< |
} |
248 |
< |
persist_end(); |
249 |
< |
ttree_out(0); |
250 |
< |
} |
251 |
< |
if ( $doforw ) { |
252 |
< |
for (my $splice = 0; $splice < $nsplice; $splice++) { |
253 |
< |
if (! $npleft ) { |
254 |
< |
wait(); |
255 |
< |
die "rtcontrib process reported error" if ( $? ); |
256 |
< |
$npleft++; |
257 |
< |
} |
258 |
< |
bg_tree_rtcontrib(1, $splice); |
259 |
< |
$npleft--; |
260 |
< |
} |
261 |
< |
while (wait() >= 0) { |
262 |
< |
die "rtcontrib process reported error" if ( $? ); |
263 |
< |
$npleft++; |
264 |
< |
} |
265 |
< |
persist_end(); |
266 |
< |
ttree_out(1); |
267 |
< |
} |
202 |
> |
|
203 |
> |
# Start rcontrib processes for compute each side |
204 |
> |
do_tree_rtcontrib(0) if ( $doback ); |
205 |
> |
do_tree_rtcontrib(1) if ( $doforw ); |
206 |
> |
|
207 |
|
} # end of sub do_tree_bsdf() |
208 |
|
|
209 |
< |
# Run rtcontrib process in background to generate tensor tree samples |
210 |
< |
sub bg_tree_rtcontrib { |
272 |
< |
my $pid = fork(); |
273 |
< |
die "Cannot fork new process" unless defined $pid; |
274 |
< |
if ($pid > 0) { return $pid; } |
209 |
> |
# Run rcontrib process to generate tensor tree samples |
210 |
> |
sub do_tree_rtcontrib { |
211 |
|
my $forw = shift; |
276 |
– |
my $pn = shift; |
277 |
– |
my $pbeg = $pdiv[$pn]; |
278 |
– |
my $plen = $pdiv[$pn+1] - $pbeg; |
212 |
|
my $matargs = "-m $bmodnm"; |
213 |
< |
if ( !$forw || !$doback ) { $matargs .= " -m $fmodnm"; } |
214 |
< |
my $cmd = "rtcontrib $rtargs -h -ff -fo -c $nsamp " . |
213 |
> |
if ( !$forw || !$doback || $tensortree==3 ) { $matargs .= " -m $fmodnm"; } |
214 |
> |
my $cmd = "rcontrib $rtargs -h -ff -fo -n $nproc -c $nsamp " . |
215 |
|
"-e '$disk2sq' -bn '$ns*$ns' " . |
216 |
|
"-b '$ns*floor(out_square_x*$ns)+floor(out_square_y*$ns)' " . |
217 |
< |
"-o $td/%s_" . sprintf("%03d", $pn) . ".flt $matargs $octree"; |
217 |
> |
"-o $td/%s.flt $matargs $octree"; |
218 |
|
if ( $tensortree == 3 ) { |
219 |
|
# Isotropic BSDF |
220 |
< |
$cmd = "cnt $plen $ny $nx " . |
221 |
< |
"| rcalc -e 'r1=rand(($pn+.8681)*recno-.673892)' " . |
222 |
< |
"-e 'r2=rand(($pn-5.37138)*recno+67.1737811)' " . |
223 |
< |
"-e 'r3=rand(($pn+3.17603772)*recno+83.766771)' " . |
224 |
< |
"-e 'Dx=1-2*($pbeg+\$1+r1)/$ns;Dy:0;Dz=sqrt(1-Dx*Dx)' " . |
220 |
> |
my $ns2 = $ns / 2; |
221 |
> |
$cmd = "cnt $ns2 $ny $nx " . |
222 |
> |
"| rcalc -e 'r1=rand(.8681*recno-.673892)' " . |
223 |
> |
"-e 'r2=rand(-5.37138*recno+67.1737811)' " . |
224 |
> |
"-e 'r3=rand(+3.17603772*recno+83.766771)' " . |
225 |
> |
"-e 'Dx=1-2*(\$1+r1)/$ns;Dy:0;Dz=sqrt(1-Dx*Dx)' " . |
226 |
|
"-e 'xp=(\$3+r2)*(($dim[1]-$dim[0])/$nx)+$dim[0]' " . |
227 |
|
"-e 'yp=(\$2+r3)*(($dim[3]-$dim[2])/$ny)+$dim[2]' " . |
228 |
|
"-e 'zp=$dim[5-$forw]' -e 'myDz=Dz*($forw*2-1)' " . |
233 |
|
# Anisotropic BSDF |
234 |
|
# Sample area vertically to improve load balance, since |
235 |
|
# shading systems usually have bilateral symmetry (L-R) |
236 |
< |
$cmd = "cnt $plen $ns $ny $nx " . |
237 |
< |
"| rcalc -e 'r1=rand(($pn+.8681)*recno-.673892)' " . |
238 |
< |
"-e 'r2=rand(($pn-5.37138)*recno+67.1737811)' " . |
239 |
< |
"-e 'r3=rand(($pn+3.17603772)*recno+83.766771)' " . |
240 |
< |
"-e 'r4=rand(($pn-2.3857833)*recno-964.72738)' " . |
241 |
< |
"-e 'in_square_x=($pbeg+\$1+r1)/$ns' " . |
236 |
> |
$cmd = "cnt $ns $ns $ny $nx " . |
237 |
> |
"| rcalc -e 'r1=rand(.8681*recno-.673892)' " . |
238 |
> |
"-e 'r2=rand(-5.37138*recno+67.1737811)' " . |
239 |
> |
"-e 'r3=rand(3.17603772*recno+83.766771)' " . |
240 |
> |
"-e 'r4=rand(-2.3857833*recno-964.72738)' " . |
241 |
> |
"-e 'in_square_x=(\$1+r1)/$ns' " . |
242 |
|
"-e 'in_square_y=(\$2+r2)/$ns' -e '$sq2disk' " . |
243 |
|
"-e 'xp=(\$4+r3)*(($dim[1]-$dim[0])/$nx)+$dim[0]' " . |
244 |
|
"-e 'yp=(\$3+r4)*(($dim[3]-$dim[2])/$ny)+$dim[2]' " . |
248 |
|
"| $cmd"; |
249 |
|
} |
250 |
|
# print STDERR "Starting: $cmd\n"; |
251 |
< |
exec($cmd); # no return; status report to parent via wait |
252 |
< |
die "Cannot exec: $cmd\n"; |
253 |
< |
} # end of bg_tree_rtcontrib() |
251 |
> |
system "$cmd" || die "Failure running rcontrib"; |
252 |
> |
ttree_out($forw); |
253 |
> |
} # end of do_tree_rtcontrib() |
254 |
|
|
255 |
|
# Simplify and output tensor tree results |
256 |
|
sub ttree_out { |
257 |
|
my $forw = shift; |
258 |
|
my $side = ("Back","Front")[$forw]; |
259 |
|
my $cmd; |
260 |
< |
# Only output one transmitted distribution, preferring backwards |
261 |
< |
if ( !$forw || !$doback ) { |
260 |
> |
# Only output one transmitted anisotropic distribution, preferring backwards |
261 |
> |
if ( !$forw || !$doback || $tensortree==3 ) { |
262 |
|
print |
263 |
|
' <WavelengthData> |
264 |
|
<LayerNumber>System</LayerNumber> |
266 |
|
<SourceSpectrum>CIE Illuminant D65 1nm.ssp</SourceSpectrum> |
267 |
|
<DetectorSpectrum>ASTM E308 1931 Y.dsp</DetectorSpectrum> |
268 |
|
<WavelengthDataBlock> |
269 |
< |
<WavelengthDataDirection>Transmission</WavelengthDataDirection> |
270 |
< |
<AngleBasis>LBNL/Shirley-Chiu</AngleBasis> |
269 |
> |
'; |
270 |
> |
print "\t\t\t<WavelengthDataDirection>Transmission $side</WavelengthDataDirection>\n"; |
271 |
> |
print |
272 |
> |
' <AngleBasis>LBNL/Shirley-Chiu</AngleBasis> |
273 |
|
<ScatteringDataType>BTDF</ScatteringDataType> |
274 |
|
<ScatteringData> |
275 |
|
'; |
276 |
|
$cmd = "rcalc -if3 -e 'Omega:PI/($ns*$ns)' " . |
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"; |
279 |
> |
$cmd .= "-of $td/" . ($bmodnm,$fmodnm)[$forw] . ".flt " . |
280 |
> |
"| rttree_reduce -h -ff -t $pctcull -r $tensortree -g $ttlog2"; |
281 |
> |
$cmd .= " -a" if ($tensortree == 3); |
282 |
|
system "$cmd" || die "Failure running rttree_reduce"; |
283 |
|
} else { |
284 |
< |
$cmd .= "$td/" . ($bmodnm,$fmodnm)[$forw] . "_???.flt"; |
284 |
> |
$cmd .= "$td/" . ($bmodnm,$fmodnm)[$forw] . ".flt"; |
285 |
|
print "{\n"; |
286 |
|
system "$cmd" || die "Failure running rcalc"; |
287 |
|
for (my $i = ($tensortree==3)*$ns*$ns*$ns/2; $i-- > 0; ) { |
307 |
|
print "\t\t\t<WavelengthDataDirection>Reflection $side</WavelengthDataDirection>\n"; |
308 |
|
print |
309 |
|
' <AngleBasis>LBNL/Shirley-Chiu</AngleBasis> |
310 |
< |
<ScatteringDataType>BRDF</ScatteringDataType> |
310 |
> |
<ScatteringDataType>BTDF</ScatteringDataType> |
311 |
|
<ScatteringData> |
312 |
|
'; |
313 |
|
$cmd = "rcalc -if3 -e 'Omega:PI/($ns*$ns)' " . |
314 |
|
q{-e '$1=(0.265*$1+0.670*$2+0.065*$3)/Omega' }; |
315 |
|
if ($pctcull >= 0) { |
316 |
< |
$cmd .= "-of $td/" . ($fmodnm,$bmodnm)[$forw] . "_???.flt " . |
316 |
> |
$cmd .= "-of $td/" . ($fmodnm,$bmodnm)[$forw] . ".flt " . |
317 |
|
"| rttree_reduce -a -h -ff -t $pctcull -r $tensortree -g $ttlog2"; |
318 |
|
system "$cmd" || die "Failure running rttree_reduce"; |
319 |
|
} else { |
320 |
< |
$cmd .= "$td/" . ($fmodnm,$bmodnm)[$forw] . "_???.flt"; |
320 |
> |
$cmd .= "$td/" . ($fmodnm,$bmodnm)[$forw] . ".flt"; |
321 |
|
print "{\n"; |
322 |
|
system "$cmd" || die "Failure running rcalc"; |
323 |
|
for (my $i = ($tensortree==3)*$ns*$ns*$ns/2; $i-- > 0; ) { |
341 |
|
$tcal = ' |
342 |
|
DEGREE : PI/180; |
343 |
|
sq(x) : x*x; |
344 |
< |
Kpola(r) : select(r+1, -5, 5, 15, 25, 35, 45, 55, 65, 75, 90); |
344 |
> |
Kpola(r) : select(r+1, 0, 5, 15, 25, 35, 45, 55, 65, 75, 90); |
345 |
|
Knaz(r) : select(r, 1, 8, 16, 20, 24, 24, 24, 16, 12); |
346 |
|
Kaccum(r) : if(r-.5, Knaz(r) + Kaccum(r-1), 0); |
347 |
|
Kmax : Kaccum(Knaz(0)); |
385 |
|
kbin = kbin2(Acos(abs(Dz)),Atan2(Dy,Dx)); |
386 |
|
'; |
387 |
|
my $ndiv = 145; |
388 |
< |
# Compute scattering data using rtcontrib |
388 |
> |
# Compute scattering data using rcontrib |
389 |
|
my @tfarr; |
390 |
|
my @rfarr; |
391 |
|
my @tbarr; |
392 |
|
my @rbarr; |
393 |
|
my $cmd; |
394 |
< |
my $rtcmd = "rtcontrib $rtargs -h -ff -fo -n $nproc -c $nsamp " . |
394 |
> |
my $rtcmd = "rcontrib $rtargs -h -ff -fo -n $nproc -c $nsamp " . |
395 |
|
"-e '$kcal' -b kbin -bn $ndiv " . |
396 |
|
"-o '$td/%s.flt' -m $fmodnm -m $bmodnm $octree"; |
397 |
|
my $rccmd = "rcalc -e '$tcal' " . |
408 |
|
system "$cmd" || die "Failure running: $cmd\n"; |
409 |
|
@tfarr = `$rccmd $td/$fmodnm.flt`; |
410 |
|
die "Failure running: $rccmd $td/$fmodnm.flt\n" if ( $? ); |
411 |
+ |
chomp(@tfarr); |
412 |
|
@rfarr = `$rccmd $td/$bmodnm.flt`; |
413 |
|
die "Failure running: $rccmd $td/$bmodnm.flt\n" if ( $? ); |
414 |
+ |
chomp(@rfarr); |
415 |
|
} |
416 |
|
if ( $doback ) { |
417 |
|
$cmd = "cnt $ndiv $ny $nx | rcalc -of -e '$tcal' " . |
436 |
|
<AngleBasis> |
437 |
|
<AngleBasisName>LBNL/Klems Full</AngleBasisName> |
438 |
|
<AngleBasisBlock> |
439 |
< |
<Theta>0</Theta> |
440 |
< |
<nPhis>1</nPhis> |
441 |
< |
<ThetaBounds> |
442 |
< |
<LowerTheta>0</LowerTheta> |
443 |
< |
<UpperTheta>5</UpperTheta> |
444 |
< |
</ThetaBounds> |
506 |
< |
</AngleBasisBlock> |
507 |
< |
<AngleBasisBlock> |
508 |
< |
<Theta>10</Theta> |
509 |
< |
<nPhis>8</nPhis> |
510 |
< |
<ThetaBounds> |
511 |
< |
<LowerTheta>5</LowerTheta> |
512 |
< |
<UpperTheta>15</UpperTheta> |
513 |
< |
</ThetaBounds> |
514 |
< |
</AngleBasisBlock> |
515 |
< |
<AngleBasisBlock> |
516 |
< |
<Theta>20</Theta> |
517 |
< |
<nPhis>16</nPhis> |
518 |
< |
<ThetaBounds> |
519 |
< |
<LowerTheta>15</LowerTheta> |
520 |
< |
<UpperTheta>25</UpperTheta> |
521 |
< |
</ThetaBounds> |
522 |
< |
</AngleBasisBlock> |
523 |
< |
<AngleBasisBlock> |
524 |
< |
<Theta>30</Theta> |
525 |
< |
<nPhis>20</nPhis> |
526 |
< |
<ThetaBounds> |
527 |
< |
<LowerTheta>25</LowerTheta> |
528 |
< |
<UpperTheta>35</UpperTheta> |
529 |
< |
</ThetaBounds> |
530 |
< |
</AngleBasisBlock> |
531 |
< |
<AngleBasisBlock> |
532 |
< |
<Theta>40</Theta> |
533 |
< |
<nPhis>24</nPhis> |
534 |
< |
<ThetaBounds> |
535 |
< |
<LowerTheta>35</LowerTheta> |
536 |
< |
<UpperTheta>45</UpperTheta> |
537 |
< |
</ThetaBounds> |
538 |
< |
</AngleBasisBlock> |
539 |
< |
<AngleBasisBlock> |
540 |
< |
<Theta>50</Theta> |
541 |
< |
<nPhis>24</nPhis> |
542 |
< |
<ThetaBounds> |
543 |
< |
<LowerTheta>45</LowerTheta> |
544 |
< |
<UpperTheta>55</UpperTheta> |
545 |
< |
</ThetaBounds> |
546 |
< |
</AngleBasisBlock> |
547 |
< |
<AngleBasisBlock> |
548 |
< |
<Theta>60</Theta> |
549 |
< |
<nPhis>24</nPhis> |
550 |
< |
<ThetaBounds> |
551 |
< |
<LowerTheta>55</LowerTheta> |
552 |
< |
<UpperTheta>65</UpperTheta> |
553 |
< |
</ThetaBounds> |
554 |
< |
</AngleBasisBlock> |
555 |
< |
<AngleBasisBlock> |
556 |
< |
<Theta>70</Theta> |
557 |
< |
<nPhis>16</nPhis> |
558 |
< |
<ThetaBounds> |
559 |
< |
<LowerTheta>65</LowerTheta> |
560 |
< |
<UpperTheta>75</UpperTheta> |
561 |
< |
</ThetaBounds> |
562 |
< |
</AngleBasisBlock> |
563 |
< |
<AngleBasisBlock> |
564 |
< |
<Theta>82.5</Theta> |
565 |
< |
<nPhis>12</nPhis> |
566 |
< |
<ThetaBounds> |
567 |
< |
<LowerTheta>75</LowerTheta> |
568 |
< |
<UpperTheta>90</UpperTheta> |
569 |
< |
</ThetaBounds> |
439 |
> |
<Theta>0</Theta> |
440 |
> |
<nPhis>1</nPhis> |
441 |
> |
<ThetaBounds> |
442 |
> |
<LowerTheta>0</LowerTheta> |
443 |
> |
<UpperTheta>5</UpperTheta> |
444 |
> |
</ThetaBounds> |
445 |
|
</AngleBasisBlock> |
446 |
+ |
<AngleBasisBlock> |
447 |
+ |
<Theta>10</Theta> |
448 |
+ |
<nPhis>8</nPhis> |
449 |
+ |
<ThetaBounds> |
450 |
+ |
<LowerTheta>5</LowerTheta> |
451 |
+ |
<UpperTheta>15</UpperTheta> |
452 |
+ |
</ThetaBounds> |
453 |
+ |
</AngleBasisBlock> |
454 |
+ |
<AngleBasisBlock> |
455 |
+ |
<Theta>20</Theta> |
456 |
+ |
<nPhis>16</nPhis> |
457 |
+ |
<ThetaBounds> |
458 |
+ |
<LowerTheta>15</LowerTheta> |
459 |
+ |
<UpperTheta>25</UpperTheta> |
460 |
+ |
</ThetaBounds> |
461 |
+ |
</AngleBasisBlock> |
462 |
+ |
<AngleBasisBlock> |
463 |
+ |
<Theta>30</Theta> |
464 |
+ |
<nPhis>20</nPhis> |
465 |
+ |
<ThetaBounds> |
466 |
+ |
<LowerTheta>25</LowerTheta> |
467 |
+ |
<UpperTheta>35</UpperTheta> |
468 |
+ |
</ThetaBounds> |
469 |
+ |
</AngleBasisBlock> |
470 |
+ |
<AngleBasisBlock> |
471 |
+ |
<Theta>40</Theta> |
472 |
+ |
<nPhis>24</nPhis> |
473 |
+ |
<ThetaBounds> |
474 |
+ |
<LowerTheta>35</LowerTheta> |
475 |
+ |
<UpperTheta>45</UpperTheta> |
476 |
+ |
</ThetaBounds> |
477 |
+ |
</AngleBasisBlock> |
478 |
+ |
<AngleBasisBlock> |
479 |
+ |
<Theta>50</Theta> |
480 |
+ |
<nPhis>24</nPhis> |
481 |
+ |
<ThetaBounds> |
482 |
+ |
<LowerTheta>45</LowerTheta> |
483 |
+ |
<UpperTheta>55</UpperTheta> |
484 |
+ |
</ThetaBounds> |
485 |
+ |
</AngleBasisBlock> |
486 |
+ |
<AngleBasisBlock> |
487 |
+ |
<Theta>60</Theta> |
488 |
+ |
<nPhis>24</nPhis> |
489 |
+ |
<ThetaBounds> |
490 |
+ |
<LowerTheta>55</LowerTheta> |
491 |
+ |
<UpperTheta>65</UpperTheta> |
492 |
+ |
</ThetaBounds> |
493 |
+ |
</AngleBasisBlock> |
494 |
+ |
<AngleBasisBlock> |
495 |
+ |
<Theta>70</Theta> |
496 |
+ |
<nPhis>16</nPhis> |
497 |
+ |
<ThetaBounds> |
498 |
+ |
<LowerTheta>65</LowerTheta> |
499 |
+ |
<UpperTheta>75</UpperTheta> |
500 |
+ |
</ThetaBounds> |
501 |
+ |
</AngleBasisBlock> |
502 |
+ |
<AngleBasisBlock> |
503 |
+ |
<Theta>82.5</Theta> |
504 |
+ |
<nPhis>12</nPhis> |
505 |
+ |
<ThetaBounds> |
506 |
+ |
<LowerTheta>75</LowerTheta> |
507 |
+ |
<UpperTheta>90</UpperTheta> |
508 |
+ |
</ThetaBounds> |
509 |
+ |
</AngleBasisBlock> |
510 |
|
</AngleBasis> |
511 |
|
</DataDefinition> |
512 |
|
'; |
544 |
|
<WavelengthDataDirection>Reflection Front</WavelengthDataDirection> |
545 |
|
<ColumnAngleBasis>LBNL/Klems Full</ColumnAngleBasis> |
546 |
|
<RowAngleBasis>LBNL/Klems Full</RowAngleBasis> |
547 |
< |
<ScatteringDataType>BRDF</ScatteringDataType> |
547 |
> |
<ScatteringDataType>BTDF</ScatteringDataType> |
548 |
|
<ScatteringData> |
549 |
|
'; |
550 |
|
# Output front reflection (transposed order) |
594 |
|
<WavelengthDataDirection>Reflection Back</WavelengthDataDirection> |
595 |
|
<ColumnAngleBasis>LBNL/Klems Full</ColumnAngleBasis> |
596 |
|
<RowAngleBasis>LBNL/Klems Full</RowAngleBasis> |
597 |
< |
<ScatteringDataType>BRDF</ScatteringDataType> |
597 |
> |
<ScatteringDataType>BTDF</ScatteringDataType> |
598 |
|
<ScatteringData> |
599 |
|
'; |
600 |
|
# Output back reflection (transposed order) |