6 |
|
# G. Ward |
7 |
|
# |
8 |
|
use strict; |
9 |
+ |
use File::Temp qw/ :mktemp /; |
10 |
|
if ($#ARGV < 0) { |
11 |
< |
print STDERR, "Usage: genklemsamp [-c N ][-f{a|f|d}] [view opts] [geom.rad ..]\n"; |
11 |
> |
print STDERR "Usage: genklemsamp [-c N ][-f{a|f|d}] [view opts] [geom.rad ..]\n"; |
12 |
|
exit 1; |
13 |
|
} |
14 |
< |
my $td = `mktemp -d /tmp/genklemsamp.XXXXXX`; |
14 |
> |
my $td = mkdtemp("/tmp/genklemsamp.XXXXXX"); |
15 |
|
chomp $td; |
16 |
|
my $nsamp = 1000; |
17 |
|
my $fmt = "a"; |
29 |
|
push @vopts, "@ARGV[0..1]"; |
30 |
|
shift @ARGV; |
31 |
|
} elsif ("$ARGV[0]" =~ /^-v./) { |
32 |
< |
print STDERR, "Unsupported view option: $ARGV[0]\n"; |
32 |
> |
print STDERR "Unsupported view option: $ARGV[0]\n"; |
33 |
|
exit 1; |
34 |
|
} elsif ("$ARGV[0]" =~ /^-./) { |
35 |
< |
print STDERR, "Unknown option: $ARGV[0]\n"; |
35 |
> |
print STDERR "Unknown option: $ARGV[0]\n"; |
36 |
|
exit 1; |
37 |
|
} else { |
38 |
|
last; |
64 |
|
sin_kpol = sin(Kpol); |
65 |
|
K0 = cos(Kazi)*sin_kpol; |
66 |
|
K1 = sin(Kazi)*sin_kpol; |
67 |
< |
K2 = cos(Kpol); |
67 |
> |
K2 = sqrt(1 - sin_kpol*sin_kpol); |
68 |
|
'; |
69 |
|
my $ndiv = 145; |
70 |
|
# Do we have any Radiance input files? |
77 |
|
push @vd, $_; |
78 |
|
$_ = $vwset; s/^.*Vdz://; s/;.*$//; |
79 |
|
push @vd, $_; |
80 |
< |
my @bcube = split /\s+/, `getinfo -d < $td/surf.oct`; |
80 |
> |
my @bcube = split ' ', `getinfo -d < $td/surf.oct`; |
81 |
|
$width = $bcube[3]*sqrt(3); |
82 |
|
$height = $width; |
83 |
|
push @vopts, ("-vp", $bcube[0]+$bcube[3]/2-$width/2*$vd[0], |
84 |
|
$bcube[1]+$bcube[3]/2-$width/2*$vd[1], |
85 |
|
$bcube[2]+$bcube[3]/2-$width/2*$vd[2]); |
86 |
< |
push @vopts, ("-vh", $width, "-vh", $height); |
86 |
> |
push @vopts, ("-vh", $width, "-vv", $height); |
87 |
|
$vwset = `vwright @vopts V`; |
88 |
+ |
$sca = sqrt($nsamp/($width*$height)); |
89 |
|
my $xres; |
90 |
|
my $yres; |
91 |
< |
my $ntot; |
91 |
> |
my $ntot = 0; |
92 |
|
# This generally passes through the loop twice to get density right |
93 |
< |
do { |
93 |
> |
while ($ntot < $nsamp) { |
94 |
|
$xres = int($width*$sca) + 1; |
95 |
|
$yres = int($height*$sca) + 1; |
96 |
|
system "vwrays -ff -x $xres -y $yres -pa 0 -pj .7 @vopts " . |
100 |
|
q{-e 'nOK=sq(Vdx*$5+Vdy*$6+Vdz*$7)-.999' } . |
101 |
|
"-if7 -of > $td/origins.flt"; |
102 |
|
$ntot = -s "$td/origins.flt"; |
103 |
< |
$ntot /= 3*4; |
103 |
> |
$ntot /= 3*4; # number of bytes per sample position |
104 |
|
if ($ntot == 0) { |
105 |
< |
print STDERR, "View direction does not correspond to any surfaces\n"; |
105 |
> |
if ($nsamp < 200) { |
106 |
> |
$sca = sqrt(200/($width*$height)); |
107 |
> |
redo; |
108 |
> |
} |
109 |
> |
print STDERR "View direction does not correspond to any surfaces\n"; |
110 |
|
exit 1; |
111 |
|
} |
112 |
|
$sca *= 1.05 * sqrt($nsamp/$ntot); |
113 |
< |
} while ($ntot < $nsamp); |
113 |
> |
} |
114 |
|
# All set to produce our samples |
115 |
< |
for ($k = 1; $k <= $ndiv; $k++) { |
115 |
> |
for (my $k = 1; $k <= $ndiv; $k++) { |
116 |
|
my $rn = rand(10); |
117 |
|
my $r1 = rand; my $r2 = rand; |
118 |
|
# Chance of using = (number_still_needed)/(number_left_avail) |
120 |
|
"-e 'cond=($nsamp-outno+1)/($ntot-recno+1)-rand($rn+recno)' " . |
121 |
|
"-e 'Kbin=$k;x1=rand($r1+recno);x2=rand($r2+recno)' " . |
122 |
|
q{-e '$1=$1+Vo*Vdx; $2=$2+Vo*Vdy; $3=$3+Vo*Vdz' } . |
123 |
< |
q{-e '$4=K0*Vhx+K1*Vvx+K2*Vdx' } . |
124 |
< |
q{-e '$5=K0*Vhy+K1*Vvy+K2*Vdy' } . |
125 |
< |
q{-e '$6=K0*Vhz+K1*Vvz+K2*Vdz' } . |
123 |
> |
q{-e '$4=-K0*Vhx-K1*Vvx+K2*Vdx' } . |
124 |
> |
q{-e '$5=-K0*Vhy-K1*Vvy+K2*Vdy' } . |
125 |
> |
q{-e '$6=-K0*Vhz-K1*Vvz+K2*Vdz' } . |
126 |
|
"-if3 $td/origins.flt"; |
127 |
|
} |
128 |
|
} else { |
130 |
|
my $xres = int($width*$sca) + 1; |
131 |
|
my $yres = int($height*$sca) + 1; |
132 |
|
my $ntot = $xres * $yres; |
133 |
< |
for ($k = 1; $k <= $ndiv; $k++) { |
133 |
> |
for (my $k = 1; $k <= $ndiv; $k++) { |
134 |
|
my $rn = rand(10); |
135 |
|
my $r1 = rand; my $r2 = rand; |
136 |
|
my $r3 = rand; my $r4 = rand; |
144 |
|
q{-e '$1=Vpx+Vo*Vdx+hpos*Vhx+vpos*Vvx' } . |
145 |
|
q{-e '$2=Vpy+Vo*Vdy+hpos*Vhy+vpos*Vvy' } . |
146 |
|
q{-e '$3=Vpz+Vo*Vdz+hpos*Vhz+vpos*Vvz' } . |
147 |
< |
q{-e '$4=K0*Vhx+K1*Vvx+K2*Vdx' } . |
148 |
< |
q{-e '$5=K0*Vhy+K1*Vvy+K2*Vdy' } . |
149 |
< |
q{-e '$6=K0*Vhz+K1*Vvz+K2*Vdz' } ; |
147 |
> |
q{-e '$4=-K0*Vhx-K1*Vvx+K2*Vdx' } . |
148 |
> |
q{-e '$5=-K0*Vhy-K1*Vvy+K2*Vdy' } . |
149 |
> |
q{-e '$6=-K0*Vhz-K1*Vvz+K2*Vdz' } ; |
150 |
|
} |
151 |
|
} |
152 |
|
system "rm -rf $td"; |