ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/util/genBSDF.pl
(Generate patch)

Comparing ray/src/util/genBSDF.pl (file contents):
Revision 2.9 by greg, Mon Feb 21 22:48:51 2011 UTC vs.
Revision 2.11 by greg, Thu Feb 24 20:14:26 2011 UTC

# Line 13 | Line 13 | sub userror {
13   my $td = `mktemp -d /tmp/genBSDF.XXXXXX`;
14   chomp $td;
15   my $nsamp = 1000;
16 + my $rtargs = "-w -ab 5 -ad 700 -lw 3e-6";
17   my $mgfin = 0;
18   my $geout = 1;
19   my $nproc = 1;
# Line 23 | Line 24 | my @dim;
24   while ($#ARGV >= 0) {
25          if ("$ARGV[0]" =~ /^[-+]m/) {
26                  $mgfin = ("$ARGV[0]" =~ /^\+/);
27 +        } elsif ("$ARGV[0]" eq "-r") {
28 +                $rtargs = "$rtargs $ARGV[1]";
29 +                shift @ARGV;
30          } elsif ("$ARGV[0]" =~ /^[-+]g/) {
31                  $geout = ("$ARGV[0]" =~ /^\+/);
32          } elsif ("$ARGV[0]" =~ /^[-+]f/) {
# Line 66 | Line 70 | if ($#dim != 5) {
70   }
71   print STDERR "Warning: Device extends into room\n" if ($dim[5] > 1e-5);
72   # Add receiver surfaces (rectangular)
69 – my $bmodnm="receiver_behind";
73   my $fmodnm="receiver_face";
74 + my $bmodnm="receiver_behind";
75   open(RADSCN, ">> $radscn");
76 < print RADSCN "void glow $fmodnm\n0\n0\n4 0 0 0 0\n\n";
77 < print RADSCN "$fmodnm polygon f_receiver\n0\n0\n12\n";
78 < print RADSCN "\t",$dim[0],"\t",$dim[2],"\t",$dim[5]+2e-5,"\n";
79 < print RADSCN "\t",$dim[0],"\t",$dim[3],"\t",$dim[5]+2e-5,"\n";
76 < print RADSCN "\t",$dim[1],"\t",$dim[3],"\t",$dim[5]+2e-5,"\n";
77 < print RADSCN "\t",$dim[1],"\t",$dim[2],"\t",$dim[5]+2e-5,"\n";
78 < print RADSCN "void glow $bmodnm\n0\n0\n4 0 0 0 0\n\n";
79 < print RADSCN "$bmodnm polygon b_receiver\n0\n0\n12\n";
80 < print RADSCN "\t",$dim[1],"\t",$dim[2],"\t",$dim[4]-2e-5,"\n";
81 < print RADSCN "\t",$dim[1],"\t",$dim[3],"\t",$dim[4]-2e-5,"\n";
82 < print RADSCN "\t",$dim[0],"\t",$dim[3],"\t",$dim[4]-2e-5,"\n";
83 < print RADSCN "\t",$dim[0],"\t",$dim[2],"\t",$dim[4]-2e-5,"\n";
76 > print RADSCN "void glow $fmodnm\n0\n0\n4 1 1 1 0\n\n";
77 > print RADSCN "$fmodnm source f_receiver\n0\n0\n4 0 0 1 180\n";
78 > print RADSCN "void glow $bmodnm\n0\n0\n4 1 1 1 0\n\n";
79 > print RADSCN "$bmodnm source b_receiver\n0\n0\n4 0 0 -1 180\n";
80   close RADSCN;
81   # Generate octree
82   system "oconv -w $radscn > $octree";
# Line 110 | Line 106 | KprojOmega = PI * if(Kbin-.5,
106   # Compute Klems bin from exiting ray direction (forward or backward)
107   my $kcal = '
108   DEGREE : PI/180;
109 + abs(x) : if(x, x, -x);
110   Acos(x) : 1/DEGREE * if(x-1, 0, if(-1-x, 0, acos(x)));
111   posangle(a) : if(-a, a + 2*PI, a);
112   Atan2(y,x) : 1/DEGREE * posangle(atan2(y,x));
# Line 130 | Line 127 | kbin2(pol,azi) = select(kfindrow(1, pol),
127                  kaccum(7) + kazn(azi,360/knaz(8)),
128                  kaccum(8) + kazn(azi,360/knaz(9))
129          );
130 < kbin = if(Dz, kbin2(Acos(Dz),Atan2(Dy,Dx)), kbin2(Acos(-Dz),Atan2(-Dy,-Dx)));
130 > kbin = kbin2(Acos(abs(Dz)),Atan2(Dy,Dx));
131   ';
132   my $ndiv = 145;
133   my $nx = int(sqrt($nsamp*($dim[1]-$dim[0])/($dim[3]-$dim[2])) + .5);
# Line 142 | Line 139 | my @rfarr;
139   my @tbarr;
140   my @rbarr;
141   my $cmd;
145 – my $rtargs = "-w -ab 5 -ad 700 -lw 3e-6";
142   my $rtcmd = "rtcontrib -h -ff -fo -n $nproc -c $nsamp " .
143          "-e '$kcal' -b kbin -bn $ndiv " .
144          "-o '$td/%s.flt' -m $fmodnm -m $bmodnm $rtargs $octree";
# Line 151 | Line 147 | my $rccmd = "rcalc -e '$tcal' " .
147          q{-if3 -e '$1=(0.265*$1+0.670*$2+0.065*$3)/KprojOmega'};
148   if ( $doforw ) {
149   $cmd = "cnt $ndiv $ny $nx | rcalc -of -e '$tcal' " .
150 <        "-e 'xp=(\$3+rand(.35*recno-15))*(($dim[1]-$dim[0])/$nx)+$dim[0]' " .
151 <        "-e 'yp=(\$2+rand(.86*recno+11))*(($dim[3]-$dim[2])/$ny)+$dim[2]' " .
152 <        "-e 'zp:$dim[4]-1e-5' " .
153 <        q{-e 'Kbin=$1;x1=rand(1.21*recno+2.75);x2=rand(-3.55*recno-7.57)' } .
154 <        q{-e '$1=xp;$2=yp;$3=zp;$4=Dx;$5=Dy;$6=Dz' } .
150 >        "-e 'xp=(\$3+rand(.12*recno+288))*(($dim[1]-$dim[0])/$nx)+$dim[0]' " .
151 >        "-e 'yp=(\$2+rand(.37*recno-44))*(($dim[3]-$dim[2])/$ny)+$dim[2]' " .
152 >        "-e 'zp:$dim[4]' " .
153 >        q{-e 'Kbin=$1;x1=rand(2.75*recno+3.1);x2=rand(-2.01*recno-3.37)' } .
154 >        q{-e '$1=xp-Dx;$2=yp-Dy;$3=zp-Dz;$4=Dx;$5=Dy;$6=Dz' } .
155          "| $rtcmd";
156   system "$cmd" || die "Failure running: $cmd\n";
157   @tfarr = `$rccmd $td/$fmodnm.flt`;
# Line 167 | Line 163 | if ( $doback ) {
163   $cmd = "cnt $ndiv $ny $nx | rcalc -of -e '$tcal' " .
164          "-e 'xp=(\$3+rand(.35*recno-15))*(($dim[1]-$dim[0])/$nx)+$dim[0]' " .
165          "-e 'yp=(\$2+rand(.86*recno+11))*(($dim[3]-$dim[2])/$ny)+$dim[2]' " .
166 <        "-e 'zp:$dim[5]+1e-5' " .
166 >        "-e 'zp:$dim[5]' " .
167          q{-e 'Kbin=$1;x1=rand(1.21*recno+2.75);x2=rand(-3.55*recno-7.57)' } .
168 <        q{-e '$1=xp;$2=yp;$3=zp;$4=-Dx;$5=-Dy;$6=-Dz' } .
168 >        q{-e '$1=xp-Dx;$2=yp-Dy;$3=zp+Dz;$4=Dx;$5=Dy;$6=-Dz' } .
169          "| $rtcmd";
170   system "$cmd" || die "Failure running: $cmd\n";
171   @tbarr = `$rccmd $td/$bmodnm.flt`;

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines