michael@0: #!/usr/bin/perl
michael@0: #
michael@0: # This Source Code Form is subject to the terms of the Mozilla Public
michael@0: # License, v. 2.0. If a copy of the MPL was not distributed with this
michael@0: # file, You can obtain one at http://mozilla.org/MPL/2.0/.
michael@0: sub GenNPL {
michael@0:   my($ret) = << "END_NPL";
michael@0: /* -*- Mode: C; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*- */
michael@0: /* This Source Code Form is subject to the terms of the Mozilla Public
michael@0:  * License, v. 2.0. If a copy of the MPL was not distributed with this
michael@0:  * file, You can obtain one at http://mozilla.org/MPL/2.0/. */
michael@0: END_NPL
michael@0: 
michael@0:   return $ret;
michael@0: }
michael@0: 
michael@0: print GenNPL();
michael@0: $total=0;
michael@0: @h;
michael@0: @l;
michael@0: 
michael@0: while(<STDIN>)
michael@0: {
michael@0:    @k = split(/\s+/, $_);
michael@0:   @i = unpack("CCCC", $k[0]);
michael@0: #  printf("%x %x %s",$i[0] ,  $i[1] , "[" . $k[0] . "]   " . $i . "  " . $j . "  " . $k[1]  ."\n");
michael@0:   if((0xA1 <= $i[0]) && (0xA1 <= $i[1])){
michael@0:   $total += $k[1];
michael@0:      $v = $i[0] - 0x00A1;
michael@0:      $h[$v] += $k[1];
michael@0:      $u = $i[1] - 0x00A1;
michael@0:      $l[$u] += $k[1];
michael@0: #     print "hello $v $h[$v] $u $l[$u]\n";
michael@0:   }
michael@0: }
michael@0: 
michael@0: 
michael@0: $ffh = 0.0;
michael@0: $ffl = 0.0;
michael@0: for($i=0x00A1;$i< 0x00FF ; $i++)
michael@0: {
michael@0:      $fh[$i - 0x00a1] = $h[$i- 0x00a1] / $total;
michael@0:      $ffh += $fh[$i - 0x00a1];
michael@0: 
michael@0:      $fl[$i - 0x00a1] = $l[$i- 0x00a1] / $total;
michael@0:      $ffl += $fl[$i - 0x00a1];
michael@0: }
michael@0: $mh = $ffh / 94.0;
michael@0: $ml = $ffl / 94.0;
michael@0: 
michael@0: $sumh=0.0;
michael@0: $suml=0.0;
michael@0: for($i=0x00A1;$i< 0x00FF ; $i++)
michael@0: {
michael@0:      $sh = $fh[$i - 0x00a1] - $mh;
michael@0:      $sh *= $sh;
michael@0:      $sumh += $sh;
michael@0: 
michael@0:      $sl = $fl[$i - 0x00a1] - $ml;
michael@0:      $sl *= $sl;
michael@0:      $suml += $sl;
michael@0: }
michael@0: $sumh /= 94.0;
michael@0: $suml /= 94.0;
michael@0: $stdh = sqrt($sumh);
michael@0: $stdl = sqrt($suml);
michael@0: 
michael@0: print "{\n";
michael@0: print "  {\n";
michael@0: for($i=0x00A1;$i< 0x00FF ; $i++)
michael@0: {
michael@0:    if($i eq 0xfe) {
michael@0:      printf("   %.6ff  \/\/ FreqH[%2x]\n", $fh[$i - 0x00a1] , $i);
michael@0:    } else {
michael@0:      printf("   %.6ff, \/\/ FreqH[%2x]\n", $fh[$i - 0x00a1] , $i);
michael@0:    }
michael@0: }
michael@0: print "  },\n";
michael@0: printf ("%.6ff, \/\/ Lead Byte StdDev\n", $stdh);
michael@0: printf ("%.6ff, \/\/ Lead Byte Mean\n", $mh);
michael@0: printf ("%.6ff, \/\/ Lead Byte Weight\n", $stdh / ($stdh + $stdl));
michael@0: print "  {\n";
michael@0: for($i=0x00A1;$i< 0x00FF ; $i++)
michael@0: {
michael@0:    if($i eq 0xfe) {
michael@0:      printf("  %.6ff  \/\/ FreqL[%2x]\n", $fl[$i - 0x00a1] ,  $i);
michael@0:    } else {
michael@0:      printf("  %.6ff, \/\/ FreqL[%2x]\n", $fl[$i - 0x00a1] ,  $i);
michael@0:    }
michael@0: }
michael@0: print "  },\n";
michael@0: printf ("%.6ff, \/\/ Trail Byte StdDev\n", $stdl);
michael@0: printf ("%.6ff, \/\/ Trail Byte Mean\n", $ml);
michael@0: printf ("%.6ff  \/\/ Trial Byte Weight\n", $stdl / ($stdh + $stdl));
michael@0: print "};\n";