197 |
|
print "\t\t<IncidentDataStructure>TensorTree$tensortree</IncidentDataStructure>\n"; |
198 |
|
print "\t</DataDefinition>\n"; |
199 |
|
|
200 |
< |
# Start rtcontrib processes for compute each side |
200 |
> |
# Start rcontrib processes for compute each side |
201 |
|
do_tree_rtcontrib(0) if ( $doback ); |
202 |
|
do_tree_rtcontrib(1) if ( $doforw ); |
203 |
|
|
204 |
|
} # end of sub do_tree_bsdf() |
205 |
|
|
206 |
< |
# Run rtcontrib process in background to generate tensor tree samples |
206 |
> |
# Run rcontrib process to generate tensor tree samples |
207 |
|
sub do_tree_rtcontrib { |
208 |
|
my $forw = shift; |
209 |
|
my $matargs = "-m $bmodnm"; |
210 |
< |
if ( !$forw || !$doback ) { $matargs .= " -m $fmodnm"; } |
210 |
> |
if ( !$forw || !$doback || $tensortree==3 ) { $matargs .= " -m $fmodnm"; } |
211 |
|
my $cmd = "rcontrib $rtargs -h -ff -fo -n $nproc -c $nsamp " . |
212 |
|
"-e '$disk2sq' -bn '$ns*$ns' " . |
213 |
|
"-b '$ns*floor(out_square_x*$ns)+floor(out_square_y*$ns)' " . |
245 |
|
"| $cmd"; |
246 |
|
} |
247 |
|
# print STDERR "Starting: $cmd\n"; |
248 |
< |
system "$cmd" || die "Failure running rtcontrib"; |
248 |
> |
system "$cmd" || die "Failure running rcontrib"; |
249 |
|
ttree_out($forw); |
250 |
|
} # end of do_tree_rtcontrib() |
251 |
|
|
254 |
|
my $forw = shift; |
255 |
|
my $side = ("Back","Front")[$forw]; |
256 |
|
my $cmd; |
257 |
< |
# Only output one transmitted distribution, preferring backwards |
258 |
< |
if ( !$forw || !$doback ) { |
257 |
> |
# Only output one transmitted anisotropic distribution, preferring backwards |
258 |
> |
if ( !$forw || !$doback || $tensortree==3 ) { |
259 |
|
print |
260 |
|
' <WavelengthData> |
261 |
|
<LayerNumber>System</LayerNumber> |
263 |
|
<SourceSpectrum>CIE Illuminant D65 1nm.ssp</SourceSpectrum> |
264 |
|
<DetectorSpectrum>ASTM E308 1931 Y.dsp</DetectorSpectrum> |
265 |
|
<WavelengthDataBlock> |
266 |
< |
<WavelengthDataDirection>Transmission</WavelengthDataDirection> |
267 |
< |
<AngleBasis>LBNL/Shirley-Chiu</AngleBasis> |
266 |
> |
'; |
267 |
> |
print "\t\t\t<WavelengthDataDirection>Transmission $side</WavelengthDataDirection>\n"; |
268 |
> |
print |
269 |
> |
' <AngleBasis>LBNL/Shirley-Chiu</AngleBasis> |
270 |
|
<ScatteringDataType>BTDF</ScatteringDataType> |
271 |
|
<ScatteringData> |
272 |
|
'; |
381 |
|
kbin = kbin2(Acos(abs(Dz)),Atan2(Dy,Dx)); |
382 |
|
'; |
383 |
|
my $ndiv = 145; |
384 |
< |
# Compute scattering data using rtcontrib |
384 |
> |
# Compute scattering data using rcontrib |
385 |
|
my @tfarr; |
386 |
|
my @rfarr; |
387 |
|
my @tbarr; |
404 |
|
system "$cmd" || die "Failure running: $cmd\n"; |
405 |
|
@tfarr = `$rccmd $td/$fmodnm.flt`; |
406 |
|
die "Failure running: $rccmd $td/$fmodnm.flt\n" if ( $? ); |
407 |
+ |
chomp(@tfarr); |
408 |
|
@rfarr = `$rccmd $td/$bmodnm.flt`; |
409 |
|
die "Failure running: $rccmd $td/$bmodnm.flt\n" if ( $? ); |
410 |
+ |
chomp(@rfarr); |
411 |
|
} |
412 |
|
if ( $doback ) { |
413 |
|
$cmd = "cnt $ndiv $ny $nx | rcalc -of -e '$tcal' " . |