← Back to branch summary

~drizzle-trunk/drizzle/development

« back to all changes in this revision

Viewing changes to storage/myisam/ftbench/Ecompare.pl

  • Committer: brian
  • Date: 2008-06-25 05:29:13 UTC
  • Revision ID: brian@localhost.localdomain-20080625052913-6upwo0jsrl4lnapl
clean slate

Show diffs side-by-side

added added

removed removed

Lines of Context:
storage/myisam/ftbench/Ecompare.pl
 
1
#!/usr/bin/perl
 
2
 
 
3
# compares out-files (as created by Ereport.pl) from dir1/*.out and dir2/*.out
 
4
# for each effectiveness column computes the probability of the hypothesis
 
5
# "Both files have the same effectiveness"
 
6
 
 
7
# sign test is used to verify that test results are statistically
 
8
# significant to support the hypothesis. Function is computed on the fly.
 
9
 
 
10
# basic formula is \sum_{r=0}^R C_N^r 2^{-N}
 
11
# As N can be big, we'll work with logarithms
 
12
$log2=log(2);
 
13
sub probab {
 
14
  my $N=shift, $R=shift;
 
15
 
 
16
  my $r, $sum=0;
 
17
 
 
18
  for $r (0..$R) {
 
19
    $sum+=exp(logfac($N)-logfac($r)-logfac($N-$r)-$N*$log2);
 
20
  }
 
21
  return $sum;
 
22
}
 
23
 
 
24
# log(N!)
 
25
# for N<20 exact value from the table (below) is taken
 
26
# otherwise, Stirling approximation for N! is used
 
27
sub logfac {
 
28
  my $n=shift; die "n=$n<0" if $n<0;
 
29
  return $logfactab[$n] if $n<=$#logfactab;
 
30
  return $n*log($n)-$n+log(2*3.14159265358*$n)/2;
 
31
}
 
32
@logfactab=(
 
33
0,    0,         0.693147180559945, 1.79175946922805, 3.17805383034795,
 
34
4.78749174278205, 6.57925121201010, 8.52516136106541, 10.6046029027453,
 
35
12.8018274800815, 15.1044125730755, 17.5023078458739, 19.9872144956619,
 
36
22.5521638531234, 25.1912211827387, 27.8992713838409, 30.6718601060807,
 
37
33.5050734501369, 36.3954452080331, 39.3398841871995, 42.3356164607535,
 
38
);
 
39
 
 
40
############################# main () ###############################
 
41
#$p=shift; $m=shift; $p-=$m;
 
42
#if($p>$m) {
 
43
#  print "1 > 2 [+$p-$m]: ", probab($p+$m, $m), "\n";
 
44
#} elsif($p<$m) {
 
45
#  print "1 < 2 [+$p-$m]: ", probab($p+$m, $p), "\n";
 
46
#} else {
 
47
#  print "1 = 2 [+$p-$m]: ", probab($p+$m, $m), "\n";
 
48
#}
 
49
#exit;
 
50
 
 
51
die "Use: $0 dir1 dir2\n" unless @ARGV==2 &&
 
52
                             -d ($dir1=shift) && -d ($dir2=shift);
 
53
$_=`cd $dir1; echo *.out`;
 
54
s/\.out\b//g;
 
55
$total="";
 
56
 
 
57
for $file (split) {
 
58
  open(OUT1,$out1="$dir1/$file.out") || die "Cannot open $out1: $!";
 
59
  open(OUT2,$out2="$dir2/$file.out") || die "Cannot open $out2: $!";
 
60
 
 
61
  @p=@m=();
 
62
  while(!eof(OUT1) || !eof(OUT2)) {
 
63
    $_=<OUT1>; @l1=split; shift @l1;
 
64
    $_=<OUT2>; @l2=split; shift @l2;
 
65
 
 
66
    die "Number of columns differ in line $.\n" unless $#l1 == $#l2;
 
67
 
 
68
    for (0..$#l1) {
 
69
      $p[$_]+= $l1[$_] > $l2[$_];
 
70
      $m[$_]+= $l1[$_] < $l2[$_];
 
71
    }
 
72
  }
 
73
 
 
74
  for (0..$#l1) {
 
75
    $pp[$_]+=$p[$_]; $mm[$_]+=$m[$_];
 
76
    $total[$_].=rep($file, ($#l1 ? $_ : undef), $p[$_], $m[$_]);
 
77
  }
 
78
  close OUT1;
 
79
  close OUT2;
 
80
}
 
81
 
 
82
for (0..$#l1) {
 
83
  rep($total[$_], ($#l1 ? $_ : undef), $pp[$_], $mm[$_]);
 
84
}
 
85
 
 
86
sub rep {
 
87
  my ($test, $n, $p, $m, $c, $r)=@_;
 
88
  
 
89
  if   ($p>$m) { $c=">"; $r="+"; }
 
90
  elsif($p<$m) { $c="<"; $r="-"; }
 
91
  else         { $c="="; $r="="; }
 
92
  $n=" $n: " if defined $n;
 
93
  printf "%-8s $n $dir1 $c $dir2 [+%03d-%03d]: %16.15f\n",
 
94
          $test, $p, $m, probab($p+$m, ($p>=$m ? $m : $p));
 
95
  $r;
 
96
}