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 while(my ($k2, $v2) = each %$v)
139 delete $allweps{0}; # ignore the tuba
140 my @allweps = keys %allweps;
143 my @dmatrix = map { [ map { 0 } @allweps ] } @allweps;
144 my @wmatrix = map { [ map { 0 } @allweps ] } @allweps;
146 for my $i(0..@allweps - 1)
148 my $attackweapon = $allweps[$i];
151 for my $j(0..@allweps - 1)
153 my $defendweapon = $allweps[$j];
154 next if $attackweapon eq $defendweapon;
155 my $win = ($matrix->{$attackweapon}{$defendweapon} || 0);
156 my $lose = ($matrix->{$defendweapon}{$attackweapon} || 0);
157 my $c = ($win + $lose);
160 my $w = 1 - 1/($c * 0.1 + 1);
162 $dmatrix[$i][$j] = $p - (1 - $p); # antisymmetric
163 $wmatrix[$i][$j] = $w; # symmetric
170 @val = SolveBestSquares(\@dmatrix, \@wmatrix);
175 @val = map { undef } @allweps;
178 for my $i(0..@allweps - 1)
180 my $attackweapon = $allweps[$i];
181 $values{$attackweapon} = $val[$i];
188 my ($event, @data) = @_;
189 if($event eq 'start')
192 elsif($event eq 'startmatrix')
194 my ($addr, $type, $map, @columns) = @data;
198 print "For server $addr type $type map $map:\n";
200 elsif($event eq 'startrow')
202 my ($row, $val) = @data;
203 printf " %-30s %8s |", $stats->weaponid_to_name($row), defined $val ? sprintf("%8.5f", $val) : "N/A";
205 elsif($event eq 'cell')
207 my ($win, $lose, $p) = @data;
218 printf " %+6.3f", $p;
221 elsif($event eq 'endrow')
225 elsif($event eq 'endmatrix')
229 print " Relevance: $min\n";
232 elsif($event eq 'end')
239 my ($event, @data) = @_;
240 if($event eq 'start')
242 print "<html><body><h1>Weapon Profiling</h1>\n";
244 elsif($event eq 'startmatrix')
246 my ($addr, $type, $map, @columns) = @data;
250 print "<h2>For server $addr type $type map $map</h2>\n";
251 print "<table><tr><th>Weapon</th><th>Rating</th>\n";
252 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;
255 elsif($event eq 'startrow')
257 my ($row, $val) = @data;
258 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";
260 elsif($event eq 'cell')
262 my ($win, $lose, $p) = @data;
266 printf '<td align=center bgcolor="#808080">%d</td>', $win;
270 printf '<td align=center bgcolor="#%02x%02x%02x">%d</td>', $v - $v * $p, 255, 0, $win;
274 #printf '<td align=center bgcolor="#%02x%02x%02x">%d</td>', (255 - $v) - $v * $p, $v + $v * $p, 0, $win;
275 printf '<td align=center bgcolor="#%02x%02x%02x">%d</td>', 255, $v + $v * $p, 0, $win;
279 printf '<td align=center bgcolor="#ffff00">%d</td>', $win;
282 elsif($event eq 'endrow')
286 elsif($event eq 'endmatrix')
290 print "</table>Relevance: $min\n";
292 elsif($event eq 'end')
297 my $out_html_cache_fh;
298 sub out_html_cache($@)
300 my ($event, @data) = @_;
301 if($event eq 'startmatrix')
304 my ($addr, $type, $map, @columns) = @data;
310 open $out_html_cache_fh, ">", "$type/$map"
311 or warn "open $type/$map: $!";
312 select $out_html_cache_fh;
315 out_html($event, @data)
316 if defined $out_html_cache_fh;
317 if($event eq 'endmatrix')
321 close $out_html_cache_fh
322 if defined $out_html_cache_fh;
323 undef $out_html_cache_fh;
328 $ENV{html_cache} ? \&out_html_cache :
329 $ENV{html} ? \&out_html :
336 my ($addr, $type, $map, $data) = @_;
337 my $values = Evaluate $data;
338 my $valid = defined [values %$values]->[0];
339 my @weapons_sorted = sort { $valid ? $values->{$b} <=> $values->{$a} : $a <=> $b } keys %$values;
341 $out->(startmatrix => ($addr, $type, $map, @weapons_sorted));
342 for my $row(@weapons_sorted)
344 $out->(startrow => $row, ($valid ? $values->{$row} : undef));
345 for my $col(@weapons_sorted)
347 my $win = ($data->{$row}{$col} || 0);
348 my $lose = ($data->{$col}{$row} || 0);
350 if $row ne $col and (not defined $min or $min > $win + $lose);
351 $out->(cell => ($win, $lose, (($row ne $col) && ($win + $lose)) ? (2 * $win / ($win + $lose) - 1) : undef));
353 $out->(endrow => ());
355 $out->(endmatrix => ($min));