1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
|
#!/usr/bin/perl
use warnings;
use strict;
package Utils;
sub getLabel
{
my ($file,$label)=@_;
open(FILE,"$file") or die "$0: Cannot open $file: $!\n";
my $value;
while(<FILE>) {
chomp;
if (/$label(.*$)/) {
$value=$1;
last;
}
}
close(FILE);
defined($value) or die "$0: undefined $label in $file\n";
return $value;
}
sub myswap
{
my ($x,$y) = @_;
my $tmp = $$x;
$$x = $$y;
$$y = $tmp;
}
sub setVector
{
my ($a,$value) = @_;
my $total = scalar(@$a);
for (my $i = 0; $i < $total; $i++) {
$a->[$i] = $value;
}
}
sub checkRange
{
my ($val,@list) = @_;
foreach my $item (@list) {
return if ($val eq $item);
}
die "$0: checkRange: $val does not belong to @list\n";
}
sub loadParams
{
my ($file) = @_;
my %params;
open(FILE,"$file") or die "$0: Cannot open $file\n";
while(<FILE>) {
next if (/^#/);
if (/(^[a-zA-Z]+)=(.*$)/) {
$params{"$1"}=$2;
}
}
close(FILE);
return %params;
}
sub scale
{
my ($file,$f,$c) = @_;
open(FILE,"$file") or die "$0: Cannot open $file: $!\n";
my $counter = 0;
my @saved;
while(<FILE>) {
my @temp = split;
my $n = scalar(@temp);
if ($n == $c) {
for (my $i = 1; $i < scalar(@temp); ++$i) {
$temp[$i] *= $f;
}
}
$saved[$counter++] = \@temp;
}
close(FILE);
open(FOUT,">$file") or die "$0: Cannot open $file: $!\n";
for (my $i = 0; $i < $counter; ++$i) {
my $temp = $saved[$i];
print FOUT "@$temp\n";
}
close(FOUT);
}
sub findReflection
{
my ($s,$L) = @_;
return $L - 1 - $s;
}
sub reflected
{
my ($site1,$site2,$L) = @_;
my $rsite1 = findReflection($site1,$L);
my $rsite2 = findReflection($site2,$L);
die "$0: reflected $site1 $site2\n" unless ($rsite1 >= $rsite2);
return ($rsite2 >= $site1) ? 0 : 1;
}
1;
|
