3 # no warranty for this script
9 use FindBin; use lib $FindBin::Bin;
10 use WeaponEncounterProfile;
12 my ($statsfile) = @ARGV;
17 $stats = WeaponEncounterProfile->new($statsfile);
27 my @bigmatrix = map { [ @{$m->[$_]}, $v->[$_] ] } 0..$n-1;
32 # first: bring the highest value to the top
37 my $v = $bigmatrix[$j]->[$i];
38 if($v*$v > $bestval*$bestval)
44 die "lindep" if $best == -1;
47 ($bigmatrix[$i], $bigmatrix[$best]) = ($bigmatrix[$best], $bigmatrix[$i]);
52 my $r = $bigmatrix[$j]->[$i];
55 $bigmatrix[$j]->[$k] -= $bigmatrix[$i]->[$k] * $r / $bestval;
61 for my $i(reverse 0..$n-1)
63 my $bestval = $bigmatrix[$i]->[$i];
66 my $r = $bigmatrix[$j]->[$i];
69 $bigmatrix[$j]->[$k] -= $bigmatrix[$i]->[$k] * $r / $bestval;
74 # 3. Read off solutions
75 return map { $bigmatrix[$_]->[$n] / $bigmatrix[$_]->[$_] } 0..($n-1);
78 sub SolveBestSquares($$)
93 my $weight = $w->[$i]->[$j];
94 $num += $weight * $d->[$i]->[$j];
97 push @result, $num / $denom;
102 # build linear equation system
104 my @matrix = map { [ map { 0 } 1..$n ] } 1..$n;
105 my @vector = map { 0 } 1..$n;
116 $matrix[$z][$i] += $w->[$i]->[$z];
117 $matrix[$z][$z] -= $w->[$i]->[$z];
118 $vector[$z] += $w->[$i]->[$z] * $d->[$i]->[$z];
122 return LinSolve(\@matrix, \@vector);
129 while(my ($k, $v) = each %$matrix)
131 for(my ($k2, $v2) = each %$v)
139 my @allweps = keys %allweps;
142 my @dmatrix = map { [ map { 0 } @allweps ] } @allweps;
143 my @wmatrix = map { [ map { 0 } @allweps ] } @allweps;
145 for my $i(0..@allweps - 1)
147 my $attackweapon = $allweps[$i];
150 for my $j(0..@allweps - 1)
152 my $defendweapon = $allweps[$j];
153 next if $attackweapon eq $defendweapon;
154 my $win = ($matrix->{$attackweapon}{$defendweapon} || 0);
155 my $lose = ($matrix->{$defendweapon}{$attackweapon} || 0);
156 my $c = ($win + $lose);
159 my $w = 1 - 1/($c * 0.1 + 1);
161 $dmatrix[$i][$j] = $p - (1 - $p); # antisymmetric
162 $wmatrix[$i][$j] = $w; # symmetric
169 @val = SolveBestSquares(\@dmatrix, \@wmatrix);
174 @val = map { undef } @allweps;
177 for my $i(0..@allweps - 1)
179 my $attackweapon = $allweps[$i];
180 $values{$attackweapon} = $val[$i];
187 my ($event, @data) = @_;
188 if($event eq 'start')
191 elsif($event eq 'startmatrix')
193 my ($addr, $map, @columns) = @data;
196 print "For server @{[$addr || 'any']} map @{[$map || 'any']}:\n";
198 elsif($event eq 'startrow')
200 my ($row, $val) = @data;
201 printf " %-30s %8s |", $stats->weaponid_to_name($row), defined $val ? sprintf("%8.5f", $val) : "N/A";
203 elsif($event eq 'cell')
205 my ($win, $lose, $p) = @data;
216 printf " %+6.3f", $p;
219 elsif($event eq 'endrow')
223 elsif($event eq 'endmatrix')
227 print " Relevance: $min\n";
230 elsif($event eq 'end')
237 my ($event, @data) = @_;
238 if($event eq 'start')
240 print "<html><body><h1>Weapon Profiling</h1>\n";
242 elsif($event eq 'startmatrix')
244 my ($addr, $map, @columns) = @data;
247 print "<h2>For server @{[$addr || 'any']} map @{[$map || 'any']}:</h2>\n";
248 print "<table><tr><th>Weapon</th><th>Rating</th>\n";
249 printf '<th><img width=70 height=80 src="http://svn.icculus.org/*checkout*/nexuiz/trunk/Docs/htmlfiles/weaponimg/thirdperson-%s.png" alt="%s"></th>', $stats->weaponid_to_model($_), $stats->weaponid_to_name($_) for @columns;
252 elsif($event eq 'startrow')
254 my ($row, $val) = @data;
255 printf '<tr><th><img width=108 height=53 src="http://svn.icculus.org/*checkout*/nexuiz/trunk/Docs/htmlfiles/weaponimg/firstperson-%s.png" alt="%s"></th><th align=right>%s</th>', $stats->weaponid_to_model($row), $stats->weaponid_to_name($row), defined $val ? sprintf("%8.5f", $val) : "N/A";
257 elsif($event eq 'cell')
259 my ($win, $lose, $p) = @data;
263 printf '<td align=center bgcolor="#808080">%d</td>', $win;
267 printf '<td align=center bgcolor="#%02x%02x%02x">%d</td>', $v - $v * $p, 255, 0, $win;
271 #printf '<td align=center bgcolor="#%02x%02x%02x">%d</td>', (255 - $v) - $v * $p, $v + $v * $p, 0, $win;
272 printf '<td align=center bgcolor="#%02x%02x%02x">%d</td>', 255, $v + $v * $p, 0, $win;
276 printf '<td align=center bgcolor="#ffff00">%d</td>', $win;
279 elsif($event eq 'endrow')
283 elsif($event eq 'endmatrix')
287 print "</table>Relevance: $min\n";
289 elsif($event eq 'end')
294 my $out = $ENV{html} ? \&out_html : \&out_text;
300 my ($addr, $map, $data) = @_;
301 my $values = Evaluate $data;
302 my $valid = defined [values %$values]->[0];
303 my @weapons_sorted = sort { $valid ? $values->{$b} <=> $values->{$a} : $a <=> $b } keys %$values;
305 $out->(startmatrix => ($addr, $map, @weapons_sorted));
306 for my $row(@weapons_sorted)
308 $out->(startrow => $row, ($valid ? $values->{$row} : undef));
309 for my $col(@weapons_sorted)
311 my $win = ($data->{$row}{$col} || 0);
312 my $lose = ($data->{$col}{$row} || 0);
314 if $row ne $col and (not defined $min or $min > $win + $lose);
315 $out->(cell => ($win, $lose, (($row ne $col) && ($win + $lose)) ? (2 * $win / ($win + $lose) - 1) : undef));
317 $out->(endrow => ());
319 $out->(endmatrix => ($min));