| 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"; |
| 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"; |
| 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], |
| 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 { |
| 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"; |