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
|
#!/usr/bin/perl
use strict;
use warnings;
use Math::Complex;
my ($beta,$label,$start)=@ARGV;
defined($label) or die "USAGE: $0 beta label [start] < observerOutput.txt\n";
$start = 0 if (!defined($start));
my $flag=0;
my @matrix;
my @sum;
my $counter=0;
my $effectiveCounter=0;
while (<STDIN>) {
if (/^\#Sites=0 1 2/) {
$flag=1;
next;
}
if (/^\#Time=$beta/ and $flag==1) {
$flag++;
next;
}
if (/^\Q$label/ and $flag==2) {
matrixRead(\@matrix);
if ($counter >= $start) {
$effectiveCounter++;
matrixAdd(\@sum,\@matrix);
}
$counter++;
}
$flag=0;
}
print STDERR "#Counter=$counter EffectiveCounter=$effectiveCounter\n";
($effectiveCounter>0) or die "$0: effectiveCounter==0\n";
matrixDivide(\@sum,$effectiveCounter);
print "$label\n";
matrixPrint(\@sum);
sub matrixRead
{
my ($matrix)=@_;
$_=<STDIN>;
my @temp=split;
my $rows=$temp[0];
my $cols=$temp[1];
my $i=0;
while (<STDIN>) {
my @temp=split;
for (my $j=0;$j<scalar(@temp);$j++) {
my $index = $i+$j*$rows;
$matrix->[$index]=$temp[$j];
}
$i++;
last if ($i==$rows);
}
}
sub matrixAdd
{
my ($matrix1,$matrix2)=@_;
my $cols = sqrt(scalar(@$matrix2));
my $rows = $cols;
for (my $i=0;$i<$rows;$i++) {
for (my $j=0;$j<$cols;$j++) {
my $index = $i+$j*$rows;
if (!defined($matrix1->[$index])) {
$matrix1->[$index] = Math::Complex->make($matrix2->[$index]);
} else {
$matrix1->[$index] += Math::Complex->make($matrix2->[$index]);
}
}
}
}
sub matrixDivide
{
my ($matrix,$divisor)=@_;
my $cols = sqrt(scalar(@$matrix));
my $rows = $cols;
defined($rows) or die "$0: rows undefined\n";
for (my $i=0;$i<$rows;$i++) {
for (my $j=0;$j<$cols;$j++) {
my $index = $i+$j*$rows;
defined($matrix->[$index]) or die "$0: matrix $i $j undefined\n";
$matrix->[$index] /= $divisor;
}
}
}
sub matrixPrint
{
my ($matrix)=@_;
my $cols = sqrt(scalar(@$matrix));
my $rows = $cols;
print "$rows $cols\n";
for (my $i=0;$i<$rows;$i++) {
for (my $j=0;$j<$cols;$j++) {
my $index = $i+$j*$rows;
print "$matrix->[$index] ";
}
print "\n";
}
}
sub toReal
{
my ($t) = @_;
$t=~s/,.*$//;
$t=~s/\(//;
return $t;
}
|
