23 |
|
my $nproc = 1; |
24 |
|
my $doforw = 0; |
25 |
|
my $doback = 1; |
26 |
< |
my $pctcull = 95; |
26 |
> |
my $pctcull = 90; |
27 |
|
my $gunit = "Meter"; |
28 |
|
my @dim; |
29 |
|
# Get options |
45 |
|
} elsif ("$ARGV[0]" =~ /^[-+]b/) { |
46 |
|
$doback = ("$ARGV[0]" =~ /^\+/); |
47 |
|
} elsif ("$ARGV[0]" eq "-t") { |
48 |
+ |
# Use value < 0 for rttree_reduce bypass |
49 |
|
$pctcull = $ARGV[1]; |
50 |
|
shift @ARGV; |
51 |
|
} elsif ("$ARGV[0]" =~ /^-t[34]$/) { |
322 |
|
sub ttree_out { |
323 |
|
my $forw = shift; |
324 |
|
my $side = ("Back","Front")[$forw]; |
325 |
+ |
my $cmd; |
326 |
|
# Only output one transmitted distribution, preferring backwards |
327 |
|
if ( !$forw || !$doback ) { |
328 |
|
print |
337 |
|
<ScatteringDataType>BTDF</ScatteringDataType> |
338 |
|
<ScatteringData> |
339 |
|
'; |
340 |
< |
system "rcalc -if3 -e 'Omega:PI/($ns*$ns)' " . |
341 |
< |
q{-e '$1=(0.265*$1+0.670*$2+0.065*$3)/Omega' -of } . |
342 |
< |
"$td/" . ($bmodnm,$fmodnm)[$forw] . "_???.flt " . |
340 |
> |
$cmd = "rcalc -if3 -e 'Omega:PI/($ns*$ns)' " . |
341 |
> |
q{-e '$1=(0.265*$1+0.670*$2+0.065*$3)/Omega' }; |
342 |
> |
if ($pctcull >= 0) { |
343 |
> |
$cmd .= "-of $td/" . ($bmodnm,$fmodnm)[$forw] . "_???.flt " . |
344 |
|
"| rttree_reduce -a -h -ff -t $pctcull -r $tensortree -g $ttlog2"; |
345 |
< |
die "Failure running rttree_reduce" if ( $? ); |
345 |
> |
system "$cmd" || die "Failure running rttree_reduce"; |
346 |
> |
} else { |
347 |
> |
$cmd .= "$td/" . ($bmodnm,$fmodnm)[$forw] . "_???.flt"; |
348 |
> |
print "{\n"; |
349 |
> |
system "$cmd" || die "Failure running rcalc"; |
350 |
> |
for (my $i = ($tensortree==3)*$ns*$ns*$ns/2; $i-- > 0; ) { |
351 |
> |
print "0\n"; |
352 |
> |
} |
353 |
> |
print "}\n"; |
354 |
> |
} |
355 |
|
print |
356 |
|
' </ScatteringData> |
357 |
|
</WavelengthDataBlock> |
373 |
|
<ScatteringDataType>BRDF</ScatteringDataType> |
374 |
|
<ScatteringData> |
375 |
|
'; |
376 |
< |
system "rcalc -if3 -e 'Omega:PI/($ns*$ns)' " . |
377 |
< |
q{-e '$1=(0.265*$1+0.670*$2+0.065*$3)/Omega' -of } . |
378 |
< |
"$td/" . ($fmodnm,$bmodnm)[$forw] . "_???.flt " . |
376 |
> |
$cmd = "rcalc -if3 -e 'Omega:PI/($ns*$ns)' " . |
377 |
> |
q{-e '$1=(0.265*$1+0.670*$2+0.065*$3)/Omega' }; |
378 |
> |
if ($pctcull >= 0) { |
379 |
> |
$cmd .= "-of $td/" . ($fmodnm,$bmodnm)[$forw] . "_???.flt " . |
380 |
|
"| rttree_reduce -a -h -ff -t $pctcull -r $tensortree -g $ttlog2"; |
381 |
< |
die "Failure running rttree_reduce" if ( $? ); |
381 |
> |
system "$cmd" || die "Failure running rttree_reduce"; |
382 |
> |
} else { |
383 |
> |
$cmd .= "$td/" . ($fmodnm,$bmodnm)[$forw] . "_???.flt"; |
384 |
> |
print "{\n"; |
385 |
> |
system "$cmd" || die "Failure running rcalc"; |
386 |
> |
for (my $i = ($tensortree==3)*$ns*$ns*$ns/2; $i-- > 0; ) { |
387 |
> |
print "0\n"; |
388 |
> |
} |
389 |
> |
print "}\n"; |
390 |
> |
} |
391 |
|
print |
392 |
|
' </ScatteringData> |
393 |
|
</WavelengthDataBlock> |
404 |
|
$tcal = ' |
405 |
|
DEGREE : PI/180; |
406 |
|
sq(x) : x*x; |
407 |
< |
Kpola(r) : select(r+1, -5, 5, 15, 25, 35, 45, 55, 65, 75, 90); |
407 |
> |
Kpola(r) : select(r+1, 0, 5, 15, 25, 35, 45, 55, 65, 75, 90); |
408 |
|
Knaz(r) : select(r, 1, 8, 16, 20, 24, 24, 24, 16, 12); |
409 |
|
Kaccum(r) : if(r-.5, Knaz(r) + Kaccum(r-1), 0); |
410 |
|
Kmax : Kaccum(Knaz(0)); |
485 |
|
system "$cmd" || die "Failure running: $cmd\n"; |
486 |
|
@tbarr = `$rccmd $td/$bmodnm.flt`; |
487 |
|
die "Failure running: $rccmd $td/$bmodnm.flt\n" if ( $? ); |
488 |
+ |
chomp(@tbarr); |
489 |
|
@rbarr = `$rccmd $td/$fmodnm.flt`; |
490 |
|
die "Failure running: $rccmd $td/$fmodnm.flt\n" if ( $? ); |
491 |
+ |
chomp(@rbarr); |
492 |
|
} |
493 |
|
# Output angle basis |
494 |
|
print |
497 |
|
<AngleBasis> |
498 |
|
<AngleBasisName>LBNL/Klems Full</AngleBasisName> |
499 |
|
<AngleBasisBlock> |
500 |
< |
<Theta>0</Theta> |
501 |
< |
<nPhis>1</nPhis> |
502 |
< |
<ThetaBounds> |
503 |
< |
<LowerTheta>0</LowerTheta> |
504 |
< |
<UpperTheta>5</UpperTheta> |
505 |
< |
</ThetaBounds> |
482 |
< |
</AngleBasisBlock> |
483 |
< |
<AngleBasisBlock> |
484 |
< |
<Theta>10</Theta> |
485 |
< |
<nPhis>8</nPhis> |
486 |
< |
<ThetaBounds> |
487 |
< |
<LowerTheta>5</LowerTheta> |
488 |
< |
<UpperTheta>15</UpperTheta> |
489 |
< |
</ThetaBounds> |
490 |
< |
</AngleBasisBlock> |
491 |
< |
<AngleBasisBlock> |
492 |
< |
<Theta>20</Theta> |
493 |
< |
<nPhis>16</nPhis> |
494 |
< |
<ThetaBounds> |
495 |
< |
<LowerTheta>15</LowerTheta> |
496 |
< |
<UpperTheta>25</UpperTheta> |
497 |
< |
</ThetaBounds> |
498 |
< |
</AngleBasisBlock> |
499 |
< |
<AngleBasisBlock> |
500 |
< |
<Theta>30</Theta> |
501 |
< |
<nPhis>20</nPhis> |
502 |
< |
<ThetaBounds> |
503 |
< |
<LowerTheta>25</LowerTheta> |
504 |
< |
<UpperTheta>35</UpperTheta> |
505 |
< |
</ThetaBounds> |
506 |
< |
</AngleBasisBlock> |
507 |
< |
<AngleBasisBlock> |
508 |
< |
<Theta>40</Theta> |
509 |
< |
<nPhis>24</nPhis> |
510 |
< |
<ThetaBounds> |
511 |
< |
<LowerTheta>35</LowerTheta> |
512 |
< |
<UpperTheta>45</UpperTheta> |
513 |
< |
</ThetaBounds> |
514 |
< |
</AngleBasisBlock> |
515 |
< |
<AngleBasisBlock> |
516 |
< |
<Theta>50</Theta> |
517 |
< |
<nPhis>24</nPhis> |
518 |
< |
<ThetaBounds> |
519 |
< |
<LowerTheta>45</LowerTheta> |
520 |
< |
<UpperTheta>55</UpperTheta> |
521 |
< |
</ThetaBounds> |
522 |
< |
</AngleBasisBlock> |
523 |
< |
<AngleBasisBlock> |
524 |
< |
<Theta>60</Theta> |
525 |
< |
<nPhis>24</nPhis> |
526 |
< |
<ThetaBounds> |
527 |
< |
<LowerTheta>55</LowerTheta> |
528 |
< |
<UpperTheta>65</UpperTheta> |
529 |
< |
</ThetaBounds> |
530 |
< |
</AngleBasisBlock> |
531 |
< |
<AngleBasisBlock> |
532 |
< |
<Theta>70</Theta> |
533 |
< |
<nPhis>16</nPhis> |
534 |
< |
<ThetaBounds> |
535 |
< |
<LowerTheta>65</LowerTheta> |
536 |
< |
<UpperTheta>75</UpperTheta> |
537 |
< |
</ThetaBounds> |
538 |
< |
</AngleBasisBlock> |
539 |
< |
<AngleBasisBlock> |
540 |
< |
<Theta>82.5</Theta> |
541 |
< |
<nPhis>12</nPhis> |
542 |
< |
<ThetaBounds> |
543 |
< |
<LowerTheta>75</LowerTheta> |
544 |
< |
<UpperTheta>90</UpperTheta> |
545 |
< |
</ThetaBounds> |
500 |
> |
<Theta>0</Theta> |
501 |
> |
<nPhis>1</nPhis> |
502 |
> |
<ThetaBounds> |
503 |
> |
<LowerTheta>0</LowerTheta> |
504 |
> |
<UpperTheta>5</UpperTheta> |
505 |
> |
</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> |
570 |
+ |
</AngleBasisBlock> |
571 |
|
</AngleBasis> |
572 |
|
</DataDefinition> |
573 |
|
'; |
608 |
|
<ScatteringDataType>BRDF</ScatteringDataType> |
609 |
|
<ScatteringData> |
610 |
|
'; |
611 |
< |
# Output front reflection (reciprocity averaging) |
611 |
> |
# Output front reflection (transposed order) |
612 |
|
for (my $od = 0; $od < $ndiv; $od++) { |
613 |
|
for (my $id = 0; $id < $ndiv; $id++) { |
614 |
< |
print .5*($rfarr[$ndiv*$id + $od] + $rfarr[$ndiv*$od + $id]), ",\n"; |
614 |
> |
print $rfarr[$ndiv*$id + $od], ",\n"; |
615 |
|
} |
616 |
|
print "\n"; |
617 |
|
} |
658 |
|
<ScatteringDataType>BRDF</ScatteringDataType> |
659 |
|
<ScatteringData> |
660 |
|
'; |
661 |
< |
# Output back reflection (reciprocity averaging) |
661 |
> |
# Output back reflection (transposed order) |
662 |
|
for (my $od = 0; $od < $ndiv; $od++) { |
663 |
|
for (my $id = 0; $id < $ndiv; $id++) { |
664 |
< |
print .5*($rbarr[$ndiv*$id + $od] + $rbarr[$ndiv*$od + $id]), ",\n"; |
664 |
> |
print $rbarr[$ndiv*$id + $od], ",\n"; |
665 |
|
} |
666 |
|
print "\n"; |
667 |
|
} |