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

Comparing ray/src/util/genklemsamp.pl (file contents):
Revision 2.4 by greg, Fri Jun 19 06:49:42 2009 UTC vs.
Revision 2.6 by greg, Thu Nov 5 06:30:57 2009 UTC

# Line 63 | Line 63 | Kpol = DEGREE * (x1*Kpola(Krow) + (1-x1)*Kpola(Krow-1)
63   sin_kpol = sin(Kpol);
64   K0 = cos(Kazi)*sin_kpol;
65   K1 = sin(Kazi)*sin_kpol;
66 < K2 = cos(Kpol);
66 > K2 = sqrt(1 - sin_kpol*sin_kpol);
67   ';
68   my $ndiv = 145;
69   # Do we have any Radiance input files?
# Line 84 | Line 84 | if ($#ARGV >= 0) {
84                          $bcube[2]+$bcube[3]/2-$width/2*$vd[2]);
85          push @vopts, ("-vh", $width, "-vh", $height);
86          $vwset = `vwright @vopts V`;
87 +        $sca = sqrt($nsamp/($width*$height));
88          my $xres;
89          my $yres;
90          my $ntot = 0;
# Line 100 | Line 101 | if ($#ARGV >= 0) {
101                  $ntot = -s "$td/origins.flt";
102                  $ntot /= 3*4;
103                  if ($ntot == 0) {
104 <                        if ($sca < sqrt(199/($width*$height))) {
104 >                        if ($nsamp < 200) {
105                                  $sca = sqrt(200/($width*$height));
106                                  redo;
107                          }
# Line 118 | Line 119 | if ($#ARGV >= 0) {
119                          "-e 'cond=($nsamp-outno+1)/($ntot-recno+1)-rand($rn+recno)' " .
120                          "-e 'Kbin=$k;x1=rand($r1+recno);x2=rand($r2+recno)' " .
121                          q{-e '$1=$1+Vo*Vdx; $2=$2+Vo*Vdy; $3=$3+Vo*Vdz' } .
122 <                        q{-e '$4=K0*Vhx+K1*Vvx+K2*Vdx' } .
123 <                        q{-e '$5=K0*Vhy+K1*Vvy+K2*Vdy' } .
124 <                        q{-e '$6=K0*Vhz+K1*Vvz+K2*Vdz' } .
122 >                        q{-e '$4=K0*Vhx-K1*Vvx+K2*Vdx' } .
123 >                        q{-e '$5=K0*Vhy-K1*Vvy+K2*Vdy' } .
124 >                        q{-e '$6=K0*Vhz-K1*Vvz+K2*Vdz' } .
125                          "-if3 $td/origins.flt";
126          }
127   } else {
# Line 142 | Line 143 | if ($#ARGV >= 0) {
143                          q{-e '$1=Vpx+Vo*Vdx+hpos*Vhx+vpos*Vvx' } .
144                          q{-e '$2=Vpy+Vo*Vdy+hpos*Vhy+vpos*Vvy' } .
145                          q{-e '$3=Vpz+Vo*Vdz+hpos*Vhz+vpos*Vvz' } .
146 <                        q{-e '$4=K0*Vhx+K1*Vvx+K2*Vdx' } .
147 <                        q{-e '$5=K0*Vhy+K1*Vvy+K2*Vdy' } .
148 <                        q{-e '$6=K0*Vhz+K1*Vvz+K2*Vdz' } ;
146 >                        q{-e '$4=K0*Vhx-K1*Vvx+K2*Vdx' } .
147 >                        q{-e '$5=K0*Vhy-K1*Vvy+K2*Vdy' } .
148 >                        q{-e '$6=K0*Vhz-K1*Vvz+K2*Vdz' } ;
149          }
150   }
151   system "rm -rf $td";

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines