|
2
|
1 |
#!/usr/bin/perl -w
|
|
|
2 |
|
|
|
3 |
use strict;
|
|
|
4 |
use warnings;
|
|
|
5 |
use File::Temp qw/ tempdir /;
|
|
|
6 |
use File::Basename;
|
|
|
7 |
|
|
|
8 |
#my $td = tempdir( CLEANUP => 0 );
|
|
|
9 |
my $td = "/tmp/objline";
|
|
|
10 |
mkdir $td;
|
|
|
11 |
#print "Temdir: $td\n";
|
|
|
12 |
my @rad = ();
|
|
|
13 |
my @prllmgf = ();
|
|
|
14 |
my @oblqmgf = ();
|
|
|
15 |
|
|
|
16 |
#my $name;
|
|
|
17 |
#my $path;
|
|
|
18 |
#my $extn;
|
|
|
19 |
my $fn = "";
|
|
|
20 |
|
|
|
21 |
my $oblqxf = "-rz 45 -ry 45";
|
|
|
22 |
|
|
|
23 |
if ($#ARGV < 0) {
|
|
|
24 |
# Input is from STDIN: Capture to file
|
|
|
25 |
open(FH, ">$td/stdin.rad");
|
|
|
26 |
while (<>) {
|
|
|
27 |
print FH;
|
|
|
28 |
}
|
|
|
29 |
close FH;
|
|
|
30 |
# Pretend stdin.rad was passed as a filename
|
|
|
31 |
@ARGV = ("$td/stdin.rad");
|
|
|
32 |
}
|
|
|
33 |
|
|
|
34 |
foreach (@ARGV) {
|
|
|
35 |
my ($name, $path, $extn) = fileparse($_);
|
|
|
36 |
$fn = $name.$extn;
|
|
|
37 |
push(@rad, "$td/$fn.orig");
|
|
|
38 |
push(@prllmgf, "$td/$fn.orig.mgf");
|
|
|
39 |
push(@oblqmgf, "$td/$fn.oblq.mgf");
|
|
|
40 |
|
|
32
|
41 |
# Make sure in-line commands such as !genbox are expanded
|
|
|
42 |
system "xform -c $_ > $td/$fn.orig";
|
|
2
|
43 |
system "rad2mgf $td/$fn.orig > $td/$fn.orig.mgf";
|
|
|
44 |
|
|
|
45 |
open (FH, ">$td/$fn.oblq.mgf");
|
|
|
46 |
print FH "i $fn.orig.mgf $oblqxf > $td/$fn.oblq.mgf\n";
|
|
|
47 |
close (FH);
|
|
|
48 |
}
|
|
|
49 |
|
|
|
50 |
my $rce = '\'xm=($1+$2)/2;ym=($3+$4)/2;zm=($5+$6)/2';
|
|
|
51 |
$rce .= ';max(a,b):if(a-b,a,b);r=max(max($2-$1,$4-$3),$6-$5)*.52';
|
|
|
52 |
$rce .= ';$1=xm-r;$2=xm+r;$3=ym-r;$4=ym+r;$5=zm-r;$6=zm+r\'';
|
|
|
53 |
|
|
|
54 |
my $rads = join(' ', @rad);
|
|
|
55 |
my $prlldim = `getbbox -h $rads |rcalc -t' ' -e $rce`;
|
|
|
56 |
chomp($prlldim);
|
|
|
57 |
print "prlldim: $prlldim\n";
|
|
|
58 |
my $oblqdim = `xform $oblqxf $rads |getbbox -h|rcalc -t' ' -e $rce`;
|
|
|
59 |
chomp($oblqdim);
|
|
|
60 |
print "oblqdim: $oblqdim\n";
|
|
|
61 |
|
|
|
62 |
my $prlls = join(' ', @prllmgf); #!!!
|
|
|
63 |
print "prlls: $prlls\n";
|
|
|
64 |
my $oblqs = join(' ', @oblqmgf);
|
|
|
65 |
my $cmd;
|
|
|
66 |
foreach ('x', 'y', 'z') {
|
|
|
67 |
$cmd = "mgf2meta -t .005 $_ " . $prlldim . ' ' . $prlls . " >$td/$_.mta";
|
|
|
68 |
#system "mgf2meta -t .005 $_ $prlldim $prlls >$td/$_.mta";
|
|
|
69 |
system "$cmd";
|
|
|
70 |
}
|
|
|
71 |
$cmd = "mgf2meta -t .005 x " . $oblqdim . ' ' . $oblqs . " >$td/o.mta";
|
|
|
72 |
system "$cmd";
|
|
|
73 |
|
|
|
74 |
system("plot4", "$td/x.mta", "$td/y.mta", "$td/z.mta", "$td/o.mta");
|
|
|
75 |
|
|
|
76 |
#EOF
|