git.sesse.net Git - vlc/blob - src/ac3_decoder/ac3_srfft.h

   1 static complex_t delta16[4] =
   2  { {1.00000000000000,  0.00000000000000},
   3    {0.92387953251129, -0.38268343236509},
   4    {0.70710678118655, -0.70710678118655},
   5    {0.38268343236509, -0.92387953251129}};
   6
   7 static complex_t delta16_3[4] =
   8  { {1.00000000000000,  0.00000000000000},
   9    {0.38268343236509, -0.92387953251129},
  10    {-0.70710678118655, -0.70710678118655},
  11    {-0.92387953251129, 0.38268343236509}};
  12
  13 static complex_t delta32[8] =
  14  { {1.00000000000000,  0.00000000000000},
  15    {0.98078528040323, -0.19509032201613},
  16    {0.92387953251129, -0.38268343236509},
  17    {0.83146961230255, -0.55557023301960},
  18    {0.70710678118655, -0.70710678118655},
  19    {0.55557023301960, -0.83146961230255},
  20    {0.38268343236509, -0.92387953251129},
  21    {0.19509032201613, -0.98078528040323}};
  22
  23 static complex_t delta32_3[8] =
  24  { {1.00000000000000,  0.00000000000000},
  25    {0.83146961230255, -0.55557023301960},
  26    {0.38268343236509, -0.92387953251129},
  27    {-0.19509032201613, -0.98078528040323},
  28    {-0.70710678118655, -0.70710678118655},
  29    {-0.98078528040323, -0.19509032201613},
  30    {-0.92387953251129, 0.38268343236509},
  31    {-0.55557023301960, 0.83146961230255}};
  32
  33 static complex_t delta64[16] =
  34  { {1.00000000000000,  0.00000000000000},
  35    {0.99518472667220, -0.09801714032956},
  36    {0.98078528040323, -0.19509032201613},
  37    {0.95694033573221, -0.29028467725446},
  38    {0.92387953251129, -0.38268343236509},
  39    {0.88192126434836, -0.47139673682600},
  40    {0.83146961230255, -0.55557023301960},
  41    {0.77301045336274, -0.63439328416365},
  42    {0.70710678118655, -0.70710678118655},
  43    {0.63439328416365, -0.77301045336274},
  44    {0.55557023301960, -0.83146961230255},
  45    {0.47139673682600, -0.88192126434835},
  46    {0.38268343236509, -0.92387953251129},
  47    {0.29028467725446, -0.95694033573221},
  48    {0.19509032201613, -0.98078528040323},
  49    {0.09801714032956, -0.99518472667220}};
  50
  51 static complex_t delta64_3[16] =
  52  { {1.00000000000000,  0.00000000000000},
  53    {0.95694033573221, -0.29028467725446},
  54    {0.83146961230255, -0.55557023301960},
  55    {0.63439328416365, -0.77301045336274},
  56    {0.38268343236509, -0.92387953251129},
  57    {0.09801714032956, -0.99518472667220},
  58    {-0.19509032201613, -0.98078528040323},
  59    {-0.47139673682600, -0.88192126434836},
  60    {-0.70710678118655, -0.70710678118655},
  61    {-0.88192126434835, -0.47139673682600},
  62    {-0.98078528040323, -0.19509032201613},
  63    {-0.99518472667220, 0.09801714032956},
  64    {-0.92387953251129, 0.38268343236509},
  65    {-0.77301045336274, 0.63439328416365},
  66    {-0.55557023301960, 0.83146961230255},
  67    {-0.29028467725446, 0.95694033573221}};
  68
  69 static complex_t delta128[32] =
  70  { {1.00000000000000,  0.00000000000000},
  71    {0.99879545620517, -0.04906767432742},
  72    {0.99518472667220, -0.09801714032956},
  73    {0.98917650996478, -0.14673047445536},
  74    {0.98078528040323, -0.19509032201613},
  75    {0.97003125319454, -0.24298017990326},
  76    {0.95694033573221, -0.29028467725446},
  77    {0.94154406518302, -0.33688985339222},
  78    {0.92387953251129, -0.38268343236509},
  79    {0.90398929312344, -0.42755509343028},
  80    {0.88192126434836, -0.47139673682600},
  81    {0.85772861000027, -0.51410274419322},
  82    {0.83146961230255, -0.55557023301960},
  83    {0.80320753148064, -0.59569930449243},
  84    {0.77301045336274, -0.63439328416365},
  85    {0.74095112535496, -0.67155895484702},
  86    {0.70710678118655, -0.70710678118655},
  87    {0.67155895484702, -0.74095112535496},
  88    {0.63439328416365, -0.77301045336274},
  89    {0.59569930449243, -0.80320753148064},
  90    {0.55557023301960, -0.83146961230255},
  91    {0.51410274419322, -0.85772861000027},
  92    {0.47139673682600, -0.88192126434835},
  93    {0.42755509343028, -0.90398929312344},
  94    {0.38268343236509, -0.92387953251129},
  95    {0.33688985339222, -0.94154406518302},
  96    {0.29028467725446, -0.95694033573221},
  97    {0.24298017990326, -0.97003125319454},
  98    {0.19509032201613, -0.98078528040323},
  99    {0.14673047445536, -0.98917650996478},
 100    {0.09801714032956, -0.99518472667220},
 101    {0.04906767432742, -0.99879545620517}};
 102
 103 static complex_t delta128_3[32] =
 104  { {1.00000000000000,  0.00000000000000},
 105    {0.98917650996478, -0.14673047445536},
 106    {0.95694033573221, -0.29028467725446},
 107    {0.90398929312344, -0.42755509343028},
 108    {0.83146961230255, -0.55557023301960},
 109    {0.74095112535496, -0.67155895484702},
 110    {0.63439328416365, -0.77301045336274},
 111    {0.51410274419322, -0.85772861000027},
 112    {0.38268343236509, -0.92387953251129},
 113    {0.24298017990326, -0.97003125319454},
 114    {0.09801714032956, -0.99518472667220},
 115    {-0.04906767432742, -0.99879545620517},
 116    {-0.19509032201613, -0.98078528040323},
 117    {-0.33688985339222, -0.94154406518302},
 118    {-0.47139673682600, -0.88192126434836},
 119    {-0.59569930449243, -0.80320753148065},
 120    {-0.70710678118655, -0.70710678118655},
 121    {-0.80320753148065, -0.59569930449243},
 122    {-0.88192126434835, -0.47139673682600},
 123    {-0.94154406518302, -0.33688985339222},
 124    {-0.98078528040323, -0.19509032201613},
 125    {-0.99879545620517, -0.04906767432742},
 126    {-0.99518472667220, 0.09801714032956},
 127    {-0.97003125319454, 0.24298017990326},
 128    {-0.92387953251129, 0.38268343236509},
 129    {-0.85772861000027, 0.51410274419322},
 130    {-0.77301045336274, 0.63439328416365},
 131    {-0.67155895484702, 0.74095112535496},
 132    {-0.55557023301960, 0.83146961230255},
 133    {-0.42755509343028, 0.90398929312344},
 134    {-0.29028467725446, 0.95694033573221},
 135    {-0.14673047445536, 0.98917650996478}};
 136
 137 #define HSQRT2 0.707106781188;
 138
 139 #define TRANSZERO(A0,A4,A8,A12) { \
 140   u_r = wTB[0].real; \
 141   v_i = u_r - wTB[k*2].real; \
 142   u_r += wTB[k*2].real; \
 143   u_i = wTB[0].imag; \
 144   v_r = wTB[k*2].imag - u_i; \
 145   u_i += wTB[k*2].imag; \
 146   a_r = A0.real; \
 147   a_i = A0.imag; \
 148   a1_r = a_r; \
 149   a1_r += u_r; \
 150   A0.real = a1_r; \
 151   a_r -= u_r; \
 152   A8.real = a_r; \
 153   a1_i = a_i; \
 154   a1_i += u_i; \
 155   A0.imag = a1_i; \
 156   a_i -= u_i; \
 157   A8.imag = a_i; \
 158   a1_r = A4.real; \
 159   a1_i = A4.imag; \
 160   a_r = a1_r; \
 161   a_r -= v_r; \
 162   A4.real = a_r; \
 163   a1_r += v_r; \
 164   A12.real = a1_r; \
 165   a_i = a1_i; \
 166   a_i -= v_i; \
 167   A4.imag = a_i; \
 168   a1_i += v_i; \
 169   A12.imag = a1_i; \
 170   }
 171
 172 #define TRANSHALF_16(A2,A6,A10,A14) {\
 173   u_r = wTB[2].real; \
 174   a_r = u_r; \
 175   u_i = wTB[2].imag; \
 176   u_r += u_i; \
 177   u_i -= a_r; \
 178   a_r = wTB[6].real; \
 179   a1_r = a_r; \
 180   a_i = wTB[6].imag; \
 181   a_r = a_i - a_r; \
 182   a_i += a1_r; \
 183   v_i = u_r - a_r; \
 184   u_r += a_r; \
 185   v_r = u_i + a_i; \
 186   u_i -= a_i; \
 187   v_i *= HSQRT2; \
 188   v_r *= HSQRT2; \
 189   u_r *= HSQRT2; \
 190   u_i *= HSQRT2; \
 191   a_r = A2.real; \
 192   a_i = A2.imag; \
 193   a1_r = a_r; \
 194   a1_r += u_r; \
 195   A2.real = a1_r; \
 196   a_r -= u_r; \
 197   A10.real = a_r; \
 198   a1_i = a_i; \
 199   a1_i += u_i; \
 200   A2.imag = a1_i; \
 201   a_i -= u_i; \
 202   A10.imag = a_i; \
 203   a1_r = A6.real; \
 204   a1_i = A6.imag;  \
 205   a_r = a1_r; \
 206   a1_r += v_r; \
 207   A6.real = a1_r; \
 208   a_r -= v_r; \
 209   A14.real = a_r; \
 210   a_i = a1_i; \
 211   a1_i -= v_i; \
 212   A6.imag = a1_i; \
 213   a_i += v_i; \
 214   A14.imag = a_i; \
 215   }
 216
 217 #define TRANS(A1,A5,A9,A13,WT,WB,D,D3) { \
 218   u_r = WT.real; \
 219   a_r = u_r; \
 220   a_r *= D.imag; \
 221   u_r *= D.real; \
 222   a_i = WT.imag; \
 223   a1_i = a_i; \
 224   a1_i *= D.real; \
 225   a_i *= D.imag; \
 226   u_r -= a_i; \
 227   u_i = a_r; \
 228   u_i += a1_i; \
 229   a_r = WB.real; \
 230   a1_r = a_r; \
 231   a1_r *= D3.real; \
 232   a_r *= D3.imag; \
 233   a_i = WB.imag; \
 234   a1_i = a_i; \
 235   a_i *= D3.real; \
 236   a1_i *= D3.imag; \
 237   a1_r -= a1_i; \
 238   a_r += a_i; \
 239   v_i = u_r - a1_r; \
 240   u_r += a1_r; \
 241   v_r = a_r - u_i; \
 242   u_i += a_r; \
 243   a_r = A1.real; \
 244   a_i = A1.imag; \
 245   a1_r = a_r; \
 246   a1_r += u_r; \
 247   A1.real = a1_r; \
 248   a_r -= u_r; \
 249   A9.real = a_r; \
 250   a1_i = a_i; \
 251   a1_i += u_i; \
 252   A1.imag = a1_i; \
 253   a_i -= u_i; \
 254   A9.imag = a_i; \
 255   a1_r = A5.real; \
 256   a1_i = A5.imag;  \
 257   a_r = a1_r; \
 258   a1_r -= v_r; \
 259   A5.real = a1_r; \
 260   a_r += v_r; \
 261   A13.real = a_r; \
 262   a_i = a1_i; \
 263   a1_i -= v_i; \
 264   A5.imag = a1_i; \
 265   a_i += v_i; \
 266   A13.imag = a_i; \
 267   }