Home | History | Annotate | Download | only in lib
      1 #!/usr/bin/perl
      2 print <<'EOD';
      3 /********************************************************************
      4  *                                                                  *
      5  * THIS FILE IS PART OF THE OggVorbis SOFTWARE CODEC SOURCE CODE.   *
      6  * USE, DISTRIBUTION AND REPRODUCTION OF THIS LIBRARY SOURCE IS     *
      7  * GOVERNED BY A BSD-STYLE SOURCE LICENSE INCLUDED WITH THIS SOURCE *
      8  * IN 'COPYING'. PLEASE READ THESE TERMS BEFORE DISTRIBUTING.       *
      9  *                                                                  *
     10  * THE OggVorbis SOURCE CODE IS (C) COPYRIGHT 1994-2007             *
     11  * by the Xiph.Org Foundation http://www.xiph.org/                  *
     12  *                                                                  *
     13  ********************************************************************
     14 
     15   function: lookup data; generated by lookups.pl; edit there
     16   last mod: $Id: lookups.pl 13293 2007-07-24 00:09:47Z xiphmont $
     17 
     18  ********************************************************************/
     19 
     20 #ifndef _V_LOOKUP_DATA_H_
     21 
     22 #ifdef FLOAT_LOOKUP
     23 EOD
     24 
     25 $cos_sz=128;
     26 $invsq_sz=32;
     27 $invsq2exp_min=-32;
     28 $invsq2exp_max=32;
     29 
     30 $fromdB_sz=35;
     31 $fromdB_shift=5;
     32 $fromdB2_shift=3;
     33 
     34 $invsq_i_shift=10;
     35 $cos_i_shift=9;
     36 $delta_shift=6;
     37 
     38 print "#define COS_LOOKUP_SZ $cos_sz\n";
     39 print "static float COS_LOOKUP[COS_LOOKUP_SZ+1]={\n";
     40 
     41 for($i=0;$i<=$cos_sz;){
     42     print "\t";
     43     for($j=0;$j<4 && $i<=$cos_sz;$j++){
     44 	printf "%+.13f,", cos(3.14159265358979323846*($i++)/$cos_sz) ;
     45     }
     46     print "\n";
     47 }
     48 print "};\n\n";
     49 
     50 print "#define INVSQ_LOOKUP_SZ $invsq_sz\n";
     51 print "static float INVSQ_LOOKUP[INVSQ_LOOKUP_SZ+1]={\n";
     52 
     53 for($i=0;$i<=$invsq_sz;){
     54     print "\t";
     55     for($j=0;$j<4 && $i<=$invsq_sz;$j++){
     56 	my$indexmap=$i++/$invsq_sz*.5+.5;
     57 	printf "%.12f,", 1./sqrt($indexmap);
     58     }
     59     print "\n";
     60 }
     61 print "};\n\n";
     62 
     63 print "#define INVSQ2EXP_LOOKUP_MIN $invsq2exp_min\n";
     64 print "#define INVSQ2EXP_LOOKUP_MAX $invsq2exp_max\n";
     65 print "static float INVSQ2EXP_LOOKUP[INVSQ2EXP_LOOKUP_MAX-\\\n".
     66       "                              INVSQ2EXP_LOOKUP_MIN+1]={\n";
     67 
     68 for($i=$invsq2exp_min;$i<=$invsq2exp_max;){
     69     print "\t";
     70     for($j=0;$j<4 && $i<=$invsq2exp_max;$j++){
     71 	printf "%15.10g,", 2**($i++*-.5);
     72     }
     73     print "\n";
     74 }
     75 print "};\n\n#endif\n\n";
     76 
     77 
     78 # 0 to -140 dB
     79 $fromdB2_sz=1<<$fromdB_shift;
     80 $fromdB_gran=1<<($fromdB_shift-$fromdB2_shift);
     81 print "#define FROMdB_LOOKUP_SZ $fromdB_sz\n";
     82 print "#define FROMdB2_LOOKUP_SZ $fromdB2_sz\n";
     83 print "#define FROMdB_SHIFT $fromdB_shift\n";
     84 print "#define FROMdB2_SHIFT $fromdB2_shift\n";
     85 print "#define FROMdB2_MASK ".((1<<$fromdB_shift)-1)."\n";
     86 
     87 print "static float FROMdB_LOOKUP[FROMdB_LOOKUP_SZ]={\n";
     88 
     89 for($i=0;$i<$fromdB_sz;){
     90     print "\t";
     91     for($j=0;$j<4 && $i<$fromdB_sz;$j++){
     92 	printf "%15.10g,", 10**(.05*(-$fromdB_gran*$i++));
     93     }
     94     print "\n";
     95 }
     96 print "};\n\n";
     97 
     98 print "static float FROMdB2_LOOKUP[FROMdB2_LOOKUP_SZ]={\n";
     99 
    100 for($i=0;$i<$fromdB2_sz;){
    101     print "\t";
    102     for($j=0;$j<4 && $i<$fromdB_sz;$j++){
    103 	printf "%15.10g,", 10**(.05*(-$fromdB_gran/$fromdB2_sz*(.5+$i++)));
    104     }
    105     print "\n";
    106 }
    107 print "};\n\n#ifdef INT_LOOKUP\n\n";
    108 
    109 
    110 $iisz=0x10000>>$invsq_i_shift;
    111 print "#define INVSQ_LOOKUP_I_SHIFT $invsq_i_shift\n";
    112 print "#define INVSQ_LOOKUP_I_MASK ".(0x0ffff>>(16-$invsq_i_shift))."\n";
    113 print "static long INVSQ_LOOKUP_I[$iisz+1]={\n";
    114 for($i=0;$i<=$iisz;){
    115     print "\t";
    116     for($j=0;$j<4 && $i<=$iisz;$j++){
    117 	my$indexmap=$i++/$iisz*.5+.5;
    118 	printf "%8d,", int(1./sqrt($indexmap)*65536.+.5);
    119     }
    120     print "\n";
    121 }
    122 print "};\n\n";
    123 
    124 $cisz=0x10000>>$cos_i_shift;
    125 print "#define COS_LOOKUP_I_SHIFT $cos_i_shift\n";
    126 print "#define COS_LOOKUP_I_MASK ".(0x0ffff>>(16-$cos_i_shift))."\n";
    127 print "#define COS_LOOKUP_I_SZ $cisz\n";
    128 print "static long COS_LOOKUP_I[COS_LOOKUP_I_SZ+1]={\n";
    129 
    130 for($i=0;$i<=$cisz;){
    131     print "\t";
    132     for($j=0;$j<4 && $i<=$cisz;$j++){
    133 	printf "%8d,", int(cos(3.14159265358979323846*($i++)/$cos_sz)*16384.+.5) ;
    134     }
    135     print "\n";
    136 }
    137 print "};\n\n";
    138 
    139 
    140 print "#endif\n\n#endif\n";
    141 
    142 
    143