开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > Windows编程 > 多媒体编程 > 查看源码
layer3.c
资源名称：Visual C++视频音频开发实用工程案例精选.rar　[点击查看]
上传用户：tuheem
上传日期：2007-05-01
资源大小：21889k
文件大小：65k
源码类别：

多媒体编程

开发平台：

Visual C++

layer3.c：源码内容
  1. /*
  2.  * libmad - MPEG audio decoder library
  3.  
  4.  */
  6. # ifdef HAVE_CONFIG_H
  7. #  include "config.h"
  8. # endif
  10. # include "global.h"
  12. # include <stdlib.h>
  13. # include <string.h>
  15. # ifdef HAVE_ASSERT_H
  16. #  include <assert.h>
  17. # endif
  19. # ifdef HAVE_LIMITS_H
  20. #  include <limits.h>
  21. # else
  22. #  define CHAR_BIT  8
  23. # endif
  25. # include "fixed.h"
  26. # include "bit.h"
  27. # include "stream.h"
  28. # include "frame.h"
  29. # include "huffman.h"
  30. # include "layer3.h"
  32. /* --- Layer III ----------------------------------------------------------- */
  34. enum {
  35.   count1table_select = 0x01,
  36.   scalefac_scale     = 0x02,
  37.   preflag      = 0x04,
  38.   mixed_block_flag   = 0x08
  39. };
  41. enum {
  42.   I_STEREO  = 0x1,
  43.   MS_STEREO = 0x2
  44. };
  46. struct sideinfo {
  47.   unsigned int main_data_begin;
  48.   unsigned int private_bits;
  50.   unsigned char scfsi[2];
  52.   struct granule {
  53.     struct channel {
  54.       /* from side info */
  55.       unsigned short part2_3_length;
  56.       unsigned short big_values;
  57.       unsigned short global_gain;
  58.       unsigned short scalefac_compress;
  60.       unsigned char flags;
  61.       unsigned char block_type;
  62.       unsigned char table_select[3];
  63.       unsigned char subblock_gain[3];
  64.       unsigned char region0_count;
  65.       unsigned char region1_count;
  67.       /* from main_data */
  68.       unsigned char scalefac[39]; /* scalefac_l and/or scalefac_s */
  69.     } ch[2];
  70.   } gr[2];
  71. };
  73. /*
  74.  * scalefactor bit lengths
  75.  * derived from section 2.4.2.7 of ISO/IEC 11172-3
  76.  */
  77. static
  78. struct {
  79.   unsigned char slen1;
  80.   unsigned char slen2;
  81. } const sflen_table[16] = {
  82.   { 0, 0 }, { 0, 1 }, { 0, 2 }, { 0, 3 },
  83.   { 3, 0 }, { 1, 1 }, { 1, 2 }, { 1, 3 },
  84.   { 2, 1 }, { 2, 2 }, { 2, 3 }, { 3, 1 },
  85.   { 3, 2 }, { 3, 3 }, { 4, 2 }, { 4, 3 }
  86. };
  88. /*
  89.  * number of LSF scalefactor band values
  90.  * derived from section 2.4.3.2 of ISO/IEC 13818-3
  91.  */
  92. static
  93. unsigned char const nsfb_table[6][3][4] = {
  94.   { {  6,  5,  5, 5 },
  95.     {  9,  9,  9, 9 },
  96.     {  6,  9,  9, 9 } },
  98.   { {  6,  5,  7, 3 },
  99.     {  9,  9, 12, 6 },
  100.     {  6,  9, 12, 6 } },
  102.   { { 11, 10,  0, 0 },
  103.     { 18, 18,  0, 0 },
  104.     { 15, 18,  0, 0 } },
  106.   { {  7,  7,  7, 0 },
  107.     { 12, 12, 12, 0 },
  108.     {  6, 15, 12, 0 } },
  110.   { {  6,  6,  6, 3 },
  111.     { 12,  9,  9, 6 },
  112.     {  6, 12,  9, 6 } },
  114.   { {  8,  8,  5, 0 },
  115.     { 15, 12,  9, 0 },
  116.     {  6, 18,  9, 0 } }
  117. };
  119. /*
  120.  * MPEG-1 scalefactor band widths
  121.  * derived from Table B.8 of ISO/IEC 11172-3
  122.  */
  123. static
  124. unsigned char const sfb_48000_long[] = {
  125.    4,  4,  4,  4,  4,  4,  6,  6,  6,   8,  10,
  126.   12, 16, 18, 22, 28, 34, 40, 46, 54,  54, 192
  127. };
  129. static
  130. unsigned char const sfb_44100_long[] = {
  131.    4,  4,  4,  4,  4,  4,  6,  6,  8,   8,  10,
  132.   12, 16, 20, 24, 28, 34, 42, 50, 54,  76, 158
  133. };
  135. static
  136. unsigned char const sfb_32000_long[] = {
  137.    4,  4,  4,  4,  4,  4,  6,  6,  8,  10,  12,
  138.   16, 20, 24, 30, 38, 46, 56, 68, 84, 102,  26
  139. };
  141. static
  142. unsigned char const sfb_48000_short[] = {
  143.    4,  4,  4,  4,  4,  4,  4,  4,  4,  4,  4,  4,  6,
  144.    6,  6,  6,  6,  6, 10, 10, 10, 12, 12, 12, 14, 14,
  145.   14, 16, 16, 16, 20, 20, 20, 26, 26, 26, 66, 66, 66
  146. };
  148. static
  149. unsigned char const sfb_44100_short[] = {
  150.    4,  4,  4,  4,  4,  4,  4,  4,  4,  4,  4,  4,  6,
  151.    6,  6,  8,  8,  8, 10, 10, 10, 12, 12, 12, 14, 14,
  152.   14, 18, 18, 18, 22, 22, 22, 30, 30, 30, 56, 56, 56
  153. };
  155. static
  156. unsigned char const sfb_32000_short[] = {
  157.    4,  4,  4,  4,  4,  4,  4,  4,  4,  4,  4,  4,  6,
  158.    6,  6,  8,  8,  8, 12, 12, 12, 16, 16, 16, 20, 20,
  159.   20, 26, 26, 26, 34, 34, 34, 42, 42, 42, 12, 12, 12
  160. };
  162. static
  163. unsigned char const sfb_48000_mixed[] = {
  164.   /* long */   4,  4,  4,  4,  4,  4,  6,  6,
  165.   /* short */  4,  4,  4,  6,  6,  6,  6,  6,  6, 10,
  166.               10, 10, 12, 12, 12, 14, 14, 14, 16, 16,
  167.               16, 20, 20, 20, 26, 26, 26, 66, 66, 66
  168. };
  170. static
  171. unsigned char const sfb_44100_mixed[] = {
  172.   /* long */   4,  4,  4,  4,  4,  4,  6,  6,
  173.   /* short */  4,  4,  4,  6,  6,  6,  8,  8,  8, 10,
  174.               10, 10, 12, 12, 12, 14, 14, 14, 18, 18,
  175.               18, 22, 22, 22, 30, 30, 30, 56, 56, 56
  176. };
  178. static
  179. unsigned char const sfb_32000_mixed[] = {
  180.   /* long */   4,  4,  4,  4,  4,  4,  6,  6,
  181.   /* short */  4,  4,  4,  6,  6,  6,  8,  8,  8, 12,
  182.               12, 12, 16, 16, 16, 20, 20, 20, 26, 26,
  183.               26, 34, 34, 34, 42, 42, 42, 12, 12, 12
  184. };
  186. /*
  187.  * MPEG-2 scalefactor band widths
  188.  * derived from Table B.2 of ISO/IEC 13818-3
  189.  */
  190. static
  191. unsigned char const sfb_24000_long[] = {
  192.    6,  6,  6,  6,  6,  6,  8, 10, 12,  14,  16,
  193.   18, 22, 26, 32, 38, 46, 54, 62, 70,  76,  36
  194. };
  196. static
  197. unsigned char const sfb_22050_long[] = {
  198.    6,  6,  6,  6,  6,  6,  8, 10, 12,  14,  16,
  199.   20, 24, 28, 32, 38, 46, 52, 60, 68,  58,  54
  200. };
  202. # define sfb_16000_long  sfb_22050_long
  204. static
  205. unsigned char const sfb_24000_short[] = {
  206.    4,  4,  4,  4,  4,  4,  4,  4,  4,  6,  6,  6,  8,
  207.    8,  8, 10, 10, 10, 12, 12, 12, 14, 14, 14, 18, 18,
  208.   18, 24, 24, 24, 32, 32, 32, 44, 44, 44, 12, 12, 12
  209. };
  211. static
  212. unsigned char const sfb_22050_short[] = {
  213.    4,  4,  4,  4,  4,  4,  4,  4,  4,  6,  6,  6,  6,
  214.    6,  6,  8,  8,  8, 10, 10, 10, 14, 14, 14, 18, 18,
  215.   18, 26, 26, 26, 32, 32, 32, 42, 42, 42, 18, 18, 18
  216. };
  218. static
  219. unsigned char const sfb_16000_short[] = {
  220.    4,  4,  4,  4,  4,  4,  4,  4,  4,  6,  6,  6,  8,
  221.    8,  8, 10, 10, 10, 12, 12, 12, 14, 14, 14, 18, 18,
  222.   18, 24, 24, 24, 30, 30, 30, 40, 40, 40, 18, 18, 18
  223. };
  225. static
  226. unsigned char const sfb_24000_mixed[] = {
  227.   /* long */   6,  6,  6,  6,  6,  6,
  228.   /* short */  6,  6,  6,  8,  8,  8, 10, 10, 10, 12,
  229.               12, 12, 14, 14, 14, 18, 18, 18, 24, 24,
  230.               24, 32, 32, 32, 44, 44, 44, 12, 12, 12
  231. };
  233. static
  234. unsigned char const sfb_22050_mixed[] = {
  235.   /* long */   6,  6,  6,  6,  6,  6,
  236.   /* short */  6,  6,  6,  6,  6,  6,  8,  8,  8, 10,
  237.               10, 10, 14, 14, 14, 18, 18, 18, 26, 26,
  238.               26, 32, 32, 32, 42, 42, 42, 18, 18, 18
  239. };
  241. static
  242. unsigned char const sfb_16000_mixed[] = {
  243.   /* long */   6,  6,  6,  6,  6,  6,
  244.   /* short */  6,  6,  6,  8,  8,  8, 10, 10, 10, 12,
  245.               12, 12, 14, 14, 14, 18, 18, 18, 24, 24,
  246.               24, 30, 30, 30, 40, 40, 40, 18, 18, 18
  247. };
  249. /*
  250.  * MPEG 2.5 scalefactor band widths
  251.  * derived from public sources
  252.  */
  253. # define sfb_12000_long  sfb_16000_long
  254. # define sfb_11025_long  sfb_12000_long
  256. static
  257. unsigned char const sfb_8000_long[] = {
  258.   12, 12, 12, 12, 12, 12, 16, 20, 24,  28,  32,
  259.   40, 48, 56, 64, 76, 90,  2,  2,  2,   2,   2
  260. };
  262. # define sfb_12000_short  sfb_16000_short
  263. # define sfb_11025_short  sfb_12000_short
  265. static
  266. unsigned char const sfb_8000_short[] = {
  267.    8,  8,  8,  8,  8,  8,  8,  8,  8, 12, 12, 12, 16,
  268.   16, 16, 20, 20, 20, 24, 24, 24, 28, 28, 28, 36, 36,
  269.   36,  2,  2,  2,  2,  2,  2,  2,  2,  2, 26, 26, 26
  270. };
  272. # define sfb_12000_mixed  sfb_16000_mixed
  273. # define sfb_11025_mixed  sfb_12000_mixed
  275. /* the 8000 Hz short block scalefactor bands do not break after the first 36
  276.    frequency lines, so this is probably wrong */
  277. static
  278. unsigned char const sfb_8000_mixed[] = {
  279.   /* long */  12, 12, 12,
  280.   /* short */  4,  4,  4,  8,  8,  8, 12, 12, 12, 16, 16, 16,
  281.               20, 20, 20, 24, 24, 24, 28, 28, 28, 36, 36, 36,
  282.                2,  2,  2,  2,  2,  2,  2,  2,  2, 26, 26, 26
  283. };
  285. static
  286. struct {
  287.   unsigned char const *l;
  288.   unsigned char const *s;
  289.   unsigned char const *m;
  290. } const sfbwidth_table[9] = {
  291.   { sfb_48000_long, sfb_48000_short, sfb_48000_mixed },
  292.   { sfb_44100_long, sfb_44100_short, sfb_44100_mixed },
  293.   { sfb_32000_long, sfb_32000_short, sfb_32000_mixed },
  294.   { sfb_24000_long, sfb_24000_short, sfb_24000_mixed },
  295.   { sfb_22050_long, sfb_22050_short, sfb_22050_mixed },
  296.   { sfb_16000_long, sfb_16000_short, sfb_16000_mixed },
  297.   { sfb_12000_long, sfb_12000_short, sfb_12000_mixed },
  298.   { sfb_11025_long, sfb_11025_short, sfb_11025_mixed },
  299.   {  sfb_8000_long,  sfb_8000_short,  sfb_8000_mixed }
  300. };
  302. /*
  303.  * scalefactor band preemphasis (used only when preflag is set)
  304.  * derived from Table B.6 of ISO/IEC 11172-3
  305.  */
  306. static
  307. unsigned char const pretab[22] = {
  308.   0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 2, 2, 3, 3, 3, 2, 0
  309. };
  311. /*
  312.  * table for requantization
  313.  *
  314.  * rq_table[x].mantissa * 2^(rq_table[x].exponent) = x^(4/3)
  315.  */
  316. static
  317. struct fixedfloat {
  318.   unsigned long mantissa  : 27;
  319.   unsigned short exponent :  5;
  320. } const rq_table[8207] = {
  321. # include "rq_table.dat"
  322. };
  324. /*
  325.  * fractional powers of two
  326.  * used for requantization and joint stereo decoding
  327.  *
  328.  * root_table[3 + x] = 2^(x/4)
  329.  */
  330. static
  331. mad_fixed_t const root_table[7] = {
  332.   MAD_F(0x09837f05) /* 2^(-3/4) == 0.59460355750136 */,
  333.   MAD_F(0x0b504f33) /* 2^(-2/4) == 0.70710678118655 */,
  334.   MAD_F(0x0d744fcd) /* 2^(-1/4) == 0.84089641525371 */,
  335.   MAD_F(0x10000000) /* 2^( 0/4) == 1.00000000000000 */,
  336.   MAD_F(0x1306fe0a) /* 2^(+1/4) == 1.18920711500272 */,
  337.   MAD_F(0x16a09e66) /* 2^(+2/4) == 1.41421356237310 */,
  338.   MAD_F(0x1ae89f99) /* 2^(+3/4) == 1.68179283050743 */
  339. };
  341. /*
  342.  * coefficients for aliasing reduction
  343.  * derived from Table B.9 of ISO/IEC 11172-3
  344.  *
  345.  *  c[]  = { -0.6, -0.535, -0.33, -0.185, -0.095, -0.041, -0.0142, -0.0037 }
  346.  * cs[i] =    1 / sqrt(1 + c[i]^2)
  347.  * ca[i] = c[i] / sqrt(1 + c[i]^2)
  348.  */
  349. static
  350. mad_fixed_t const cs[8] = {
  351.   +MAD_F(0x0db84a81) /* +0.857492926 */, +MAD_F(0x0e1b9d7f) /* +0.881741997 */,
  352.   +MAD_F(0x0f31adcf) /* +0.949628649 */, +MAD_F(0x0fbba815) /* +0.983314592 */,
  353.   +MAD_F(0x0feda417) /* +0.995517816 */, +MAD_F(0x0ffc8fc8) /* +0.999160558 */,
  354.   +MAD_F(0x0fff964c) /* +0.999899195 */, +MAD_F(0x0ffff8d3) /* +0.999993155 */
  355. };
  357. static
  358. mad_fixed_t const ca[8] = {
  359.   -MAD_F(0x083b5fe7) /* -0.514495755 */, -MAD_F(0x078c36d2) /* -0.471731969 */,
  360.   -MAD_F(0x05039814) /* -0.313377454 */, -MAD_F(0x02e91dd1) /* -0.181913200 */,
  361.   -MAD_F(0x0183603a) /* -0.094574193 */, -MAD_F(0x00a7cb87) /* -0.040965583 */,
  362.   -MAD_F(0x003a2847) /* -0.014198569 */, -MAD_F(0x000f27b4) /* -0.003699975 */
  363. };
  365. /*
  366.  * IMDCT coefficients for short blocks
  367.  * derived from section 2.4.3.4.10.2 of ISO/IEC 11172-3
  368.  *
  369.  * imdct_s[i/even][k] = cos((PI / 24) * (2 * (i / 2) + 7) * (2 * k + 1))
  370.  * imdct_s[i /odd][k] = cos((PI / 24) * (2 * (6 + (i-1)/2) + 7) * (2 * k + 1))
  371.  */
  372. static
  373. mad_fixed_t const imdct_s[6][6] = {
  374. # include "imdct_s.dat"
  375. };
  377. # if !defined(ASO_IMDCT)
  378. /*
  379.  * windowing coefficients for long blocks
  380.  * derived from section 2.4.3.4.10.3 of ISO/IEC 11172-3
  381.  *
  382.  * window_l[i] = sin((PI / 36) * (i + 1/2))
  383.  */
  384. static
  385. mad_fixed_t const window_l[36] = {
  386.   MAD_F(0x00b2aa3e) /* 0.043619387 */, MAD_F(0x0216a2a2) /* 0.130526192 */,
  387.   MAD_F(0x03768962) /* 0.216439614 */, MAD_F(0x04cfb0e2) /* 0.300705800 */,
  388.   MAD_F(0x061f78aa) /* 0.382683432 */, MAD_F(0x07635284) /* 0.461748613 */,
  389.   MAD_F(0x0898c779) /* 0.537299608 */, MAD_F(0x09bd7ca0) /* 0.608761429 */,
  390.   MAD_F(0x0acf37ad) /* 0.675590208 */, MAD_F(0x0bcbe352) /* 0.737277337 */,
  391.   MAD_F(0x0cb19346) /* 0.793353340 */, MAD_F(0x0d7e8807) /* 0.843391446 */,
  393.   MAD_F(0x0e313245) /* 0.887010833 */, MAD_F(0x0ec835e8) /* 0.923879533 */,
  394.   MAD_F(0x0f426cb5) /* 0.953716951 */, MAD_F(0x0f9ee890) /* 0.976296007 */,
  395.   MAD_F(0x0fdcf549) /* 0.991444861 */, MAD_F(0x0ffc19fd) /* 0.999048222 */,
  396.   MAD_F(0x0ffc19fd) /* 0.999048222 */, MAD_F(0x0fdcf549) /* 0.991444861 */,
  397.   MAD_F(0x0f9ee890) /* 0.976296007 */, MAD_F(0x0f426cb5) /* 0.953716951 */,
  398.   MAD_F(0x0ec835e8) /* 0.923879533 */, MAD_F(0x0e313245) /* 0.887010833 */,
  400.   MAD_F(0x0d7e8807) /* 0.843391446 */, MAD_F(0x0cb19346) /* 0.793353340 */,
  401.   MAD_F(0x0bcbe352) /* 0.737277337 */, MAD_F(0x0acf37ad) /* 0.675590208 */,
  402.   MAD_F(0x09bd7ca0) /* 0.608761429 */, MAD_F(0x0898c779) /* 0.537299608 */,
  403.   MAD_F(0x07635284) /* 0.461748613 */, MAD_F(0x061f78aa) /* 0.382683432 */,
  404.   MAD_F(0x04cfb0e2) /* 0.300705800 */, MAD_F(0x03768962) /* 0.216439614 */,
  405.   MAD_F(0x0216a2a2) /* 0.130526192 */, MAD_F(0x00b2aa3e) /* 0.043619387 */,
  406. };
  407. # endif  /* ASO_IMDCT */
  409. /*
  410.  * windowing coefficients for short blocks
  411.  * derived from section 2.4.3.4.10.3 of ISO/IEC 11172-3
  412.  *
  413.  * window_s[i] = sin((PI / 12) * (i + 1/2))
  414.  */
  415. static
  416. mad_fixed_t const window_s[12] = {
  417.   MAD_F(0x0216a2a2) /* 0.130526192 */, MAD_F(0x061f78aa) /* 0.382683432 */,
  418.   MAD_F(0x09bd7ca0) /* 0.608761429 */, MAD_F(0x0cb19346) /* 0.793353340 */,
  419.   MAD_F(0x0ec835e8) /* 0.923879533 */, MAD_F(0x0fdcf549) /* 0.991444861 */,
  420.   MAD_F(0x0fdcf549) /* 0.991444861 */, MAD_F(0x0ec835e8) /* 0.923879533 */,
  421.   MAD_F(0x0cb19346) /* 0.793353340 */, MAD_F(0x09bd7ca0) /* 0.608761429 */,
  422.   MAD_F(0x061f78aa) /* 0.382683432 */, MAD_F(0x0216a2a2) /* 0.130526192 */,
  423. };
  425. /*
  426.  * coefficients for intensity stereo processing
  427.  * derived from section 2.4.3.4.9.3 of ISO/IEC 11172-3
  428.  *
  429.  * is_ratio[i] = tan(i * (PI / 12))
  430.  * is_table[i] = is_ratio[i] / (1 + is_ratio[i])
  431.  */
  432. static
  433. mad_fixed_t const is_table[7] = {
  434.   MAD_F(0x00000000) /* 0.000000000 */,
  435.   MAD_F(0x0361962f) /* 0.211324865 */,
  436.   MAD_F(0x05db3d74) /* 0.366025404 */,
  437.   MAD_F(0x08000000) /* 0.500000000 */,
  438.   MAD_F(0x0a24c28c) /* 0.633974596 */,
  439.   MAD_F(0x0c9e69d1) /* 0.788675135 */,
  440.   MAD_F(0x10000000) /* 1.000000000 */
  441. };
  443. /*
  444.  * coefficients for LSF intensity stereo processing
  445.  * derived from section 2.4.3.2 of ISO/IEC 13818-3
  446.  *
  447.  * is_lsf_table[0][i] = (1 / sqrt(sqrt(2)))^(i + 1)
  448.  * is_lsf_table[1][i] = (1 / sqrt(2))^(i + 1)
  449.  */
  450. static
  451. mad_fixed_t const is_lsf_table[2][15] = {
  452.   {
  453.     MAD_F(0x0d744fcd) /* 0.840896415 */,
  454.     MAD_F(0x0b504f33) /* 0.707106781 */,
  455.     MAD_F(0x09837f05) /* 0.594603558 */,
  456.     MAD_F(0x08000000) /* 0.500000000 */,
  457.     MAD_F(0x06ba27e6) /* 0.420448208 */,
  458.     MAD_F(0x05a8279a) /* 0.353553391 */,
  459.     MAD_F(0x04c1bf83) /* 0.297301779 */,
  460.     MAD_F(0x04000000) /* 0.250000000 */,
  461.     MAD_F(0x035d13f3) /* 0.210224104 */,
  462.     MAD_F(0x02d413cd) /* 0.176776695 */,
  463.     MAD_F(0x0260dfc1) /* 0.148650889 */,
  464.     MAD_F(0x02000000) /* 0.125000000 */,
  465.     MAD_F(0x01ae89fa) /* 0.105112052 */,
  466.     MAD_F(0x016a09e6) /* 0.088388348 */,
  467.     MAD_F(0x01306fe1) /* 0.074325445 */
  468.   }, {
  469.     MAD_F(0x0b504f33) /* 0.707106781 */,
  470.     MAD_F(0x08000000) /* 0.500000000 */,
  471.     MAD_F(0x05a8279a) /* 0.353553391 */,
  472.     MAD_F(0x04000000) /* 0.250000000 */,
  473.     MAD_F(0x02d413cd) /* 0.176776695 */,
  474.     MAD_F(0x02000000) /* 0.125000000 */,
  475.     MAD_F(0x016a09e6) /* 0.088388348 */,
  476.     MAD_F(0x01000000) /* 0.062500000 */,
  477.     MAD_F(0x00b504f3) /* 0.044194174 */,
  478.     MAD_F(0x00800000) /* 0.031250000 */,
  479.     MAD_F(0x005a827a) /* 0.022097087 */,
  480.     MAD_F(0x00400000) /* 0.015625000 */,
  481.     MAD_F(0x002d413d) /* 0.011048543 */,
  482.     MAD_F(0x00200000) /* 0.007812500 */,
  483.     MAD_F(0x0016a09e) /* 0.005524272 */
  484.   }
  485. };
  487. /*
  488.  * NAME: III_sideinfo()
  489.  * DESCRIPTION: decode frame side information from a bitstream
  490.  */
  491. static
  492. enum mad_error III_sideinfo(struct mad_bitptr *ptr, unsigned int nch,
  493.     int lsf, struct sideinfo *si,
  494.     unsigned int *data_bitlen,
  495.     unsigned int *priv_bitlen)
  496. {
  497.   unsigned int ngr, gr, ch, i;
  498.   enum mad_error result = MAD_ERROR_NONE;
  500.   *data_bitlen = 0;
  501.   *priv_bitlen = lsf ? ((nch == 1) ? 1 : 2) : ((nch == 1) ? 5 : 3);
  503.   si->main_data_begin = mad_bit_read(ptr, lsf ? 8 : 9);
  504.   si->private_bits    = mad_bit_read(ptr, *priv_bitlen);
  506.   ngr = 1;
  507.   if (!lsf) {
  508.     ngr = 2;
  510.     for (ch = 0; ch < nch; ++ch)
  511.       si->scfsi[ch] = mad_bit_read(ptr, 4);
  512.   }
  514.   for (gr = 0; gr < ngr; ++gr) {
  515.     struct granule *granule = &si->gr[gr];
  517.     for (ch = 0; ch < nch; ++ch) {
  518.       struct channel *channel = &granule->ch[ch];
  520.       channel->part2_3_length    = mad_bit_read(ptr, 12);
  521.       channel->big_values        = mad_bit_read(ptr, 9);
  522.       channel->global_gain       = mad_bit_read(ptr, 8);
  523.       channel->scalefac_compress = mad_bit_read(ptr, lsf ? 9 : 4);
  525.       *data_bitlen += channel->part2_3_length;
  527.       if (channel->big_values > 288 && result == 0)
  528. result = MAD_ERROR_BADBIGVALUES;
  530.       channel->flags = 0;
  532.       /* window_switching_flag */
  533.       if (mad_bit_read(ptr, 1)) {
  534. channel->block_type = mad_bit_read(ptr, 2);
  536. if (channel->block_type == 0 && result == 0)
  537.   result = MAD_ERROR_BADBLOCKTYPE;
  539. if (!lsf && channel->block_type == 2 && si->scfsi[ch] && result == 0)
  540.   result = MAD_ERROR_BADSCFSI;
  542. channel->region0_count = 7;
  543. channel->region1_count = 36;
  545. if (mad_bit_read(ptr, 1))
  546.   channel->flags |= mixed_block_flag;
  547. else if (channel->block_type == 2)
  548.   channel->region0_count = 8;
  550. for (i = 0; i < 2; ++i)
  551.   channel->table_select[i] = mad_bit_read(ptr, 5);
  553. # if defined(DEBUG)
  554. channel->table_select[2] = 4;  /* not used */
  555. # endif
  557. for (i = 0; i < 3; ++i)
  558.   channel->subblock_gain[i] = mad_bit_read(ptr, 3);
  559.       }
  560.       else {
  561. channel->block_type = 0;
  563. for (i = 0; i < 3; ++i)
  564.   channel->table_select[i] = mad_bit_read(ptr, 5);
  566. channel->region0_count = mad_bit_read(ptr, 4);
  567. channel->region1_count = mad_bit_read(ptr, 3);
  568.       }
  570.       /* [preflag,] scalefac_scale, count1table_select */
  571.       channel->flags |= mad_bit_read(ptr, lsf ? 2 : 3);
  572.     }
  573.   }
  575.   return result;
  576. }
  578. /*
  579.  * NAME: III_scalefactors_lsf()
  580.  * DESCRIPTION: decode channel scalefactors for LSF from a bitstream
  581.  */
  582. static
  583. unsigned int III_scalefactors_lsf(struct mad_bitptr *ptr,
  584.   struct channel *channel,
  585.   struct channel *gr1ch, int mode_extension)
  586. {
  587.   struct mad_bitptr start;
  588.   unsigned int scalefac_compress, index, slen[4], part, n, i;
  589.   unsigned char const *nsfb;
  591.   start = *ptr;
  593.   scalefac_compress = channel->scalefac_compress;
  594.   index = (channel->block_type == 2) ?
  595.     ((channel->flags & mixed_block_flag) ? 2 : 1) : 0;
  597.   if (!((mode_extension & I_STEREO) && gr1ch)) {
  598.     if (scalefac_compress < 400) {
  599.       slen[0] = (scalefac_compress >> 4) / 5;
  600.       slen[1] = (scalefac_compress >> 4) % 5;
  601.       slen[2] = (scalefac_compress % 16) >> 2;
  602.       slen[3] =  scalefac_compress %  4;
  604.       nsfb = nsfb_table[0][index];
  605.     }
  606.     else if (scalefac_compress < 500) {
  607.       scalefac_compress -= 400;
  609.       slen[0] = (scalefac_compress >> 2) / 5;
  610.       slen[1] = (scalefac_compress >> 2) % 5;
  611.       slen[2] =  scalefac_compress %  4;
  612.       slen[3] = 0;
  614.       nsfb = nsfb_table[1][index];
  615.     }
  616.     else {
  617.       scalefac_compress -= 500;
  619.       slen[0] = scalefac_compress / 3;
  620.       slen[1] = scalefac_compress % 3;
  621.       slen[2] = 0;
  622.       slen[3] = 0;
  624.       channel->flags |= preflag;
  626.       nsfb = nsfb_table[2][index];
  627.     }
  629.     n = 0;
  630.     for (part = 0; part < 4; ++part) {
  631.       for (i = 0; i < nsfb[part]; ++i)
  632. channel->scalefac[n++] = mad_bit_read(ptr, slen[part]);
  633.     }
  635.     while (n < 39)
  636.       channel->scalefac[n++] = 0;
  637.   }
  638.   else {  /* (mode_extension & I_STEREO) && gr1ch (i.e. ch == 1) */
  639.     scalefac_compress >>= 1;
  641.     if (scalefac_compress < 180) {
  642.       slen[0] =  scalefac_compress / 36;
  643.       slen[1] = (scalefac_compress % 36) / 6;
  644.       slen[2] = (scalefac_compress % 36) % 6;
  645.       slen[3] = 0;
  647.       nsfb = nsfb_table[3][index];
  648.     }
  649.     else if (scalefac_compress < 244) {
  650.       scalefac_compress -= 180;
  652.       slen[0] = (scalefac_compress % 64) >> 4;
  653.       slen[1] = (scalefac_compress % 16) >> 2;
  654.       slen[2] =  scalefac_compress %  4;
  655.       slen[3] = 0;
  657.       nsfb = nsfb_table[4][index];
  658.     }
  659.     else {
  660.       scalefac_compress -= 244;
  662.       slen[0] = scalefac_compress / 3;
  663.       slen[1] = scalefac_compress % 3;
  664.       slen[2] = 0;
  665.       slen[3] = 0;
  667.       nsfb = nsfb_table[5][index];
  668.     }
  670.     n = 0;
  671.     for (part = 0; part < 4; ++part) {
  672.       unsigned int max, is_pos;
  674.       max = (1 << slen[part]) - 1;
  676.       for (i = 0; i < nsfb[part]; ++i) {
  677. is_pos = mad_bit_read(ptr, slen[part]);
  679. channel->scalefac[n] = is_pos;
  680. gr1ch->scalefac[n++] = (is_pos == max);
  681.       }
  682.     }
  684.     while (n < 39) {
  685.       channel->scalefac[n] = 0;
  686.       gr1ch->scalefac[n++] = 0;  /* apparently not illegal */
  687.     }
  688.   }
  690.   return mad_bit_length(&start, ptr);
  691. }
  693. /*
  694.  * NAME: III_scalefactors()
  695.  * DESCRIPTION: decode channel scalefactors of one granule from a bitstream
  696.  */
  697. static
  698. unsigned int III_scalefactors(struct mad_bitptr *ptr, struct channel *channel,
  699.       struct channel const *gr0ch, unsigned int scfsi)
  700. {
  701.   struct mad_bitptr start;
  702.   unsigned int slen1, slen2, sfbi;
  704.   start = *ptr;
  706.   slen1 = sflen_table[channel->scalefac_compress].slen1;
  707.   slen2 = sflen_table[channel->scalefac_compress].slen2;
  709.   if (channel->block_type == 2) {
  710.     unsigned int nsfb;
  712.     sfbi = 0;
  714.     nsfb = (channel->flags & mixed_block_flag) ? 8 + 3 * 3 : 6 * 3;
  715.     while (nsfb--)
  716.       channel->scalefac[sfbi++] = mad_bit_read(ptr, slen1);
  718.     nsfb = 6 * 3;
  719.     while (nsfb--)
  720.       channel->scalefac[sfbi++] = mad_bit_read(ptr, slen2);
  722.     nsfb = 1 * 3;
  723.     while (nsfb--)
  724.       channel->scalefac[sfbi++] = 0;
  725.   }
  726.   else {  /* channel->block_type != 2 */
  727.     if (scfsi & 0x8) {
  728.       for (sfbi = 0; sfbi < 6; ++sfbi)
  729. channel->scalefac[sfbi] = gr0ch->scalefac[sfbi];
  730.     }
  731.     else {
  732.       for (sfbi = 0; sfbi < 6; ++sfbi)
  733. channel->scalefac[sfbi] = mad_bit_read(ptr, slen1);
  734.     }
  736.     if (scfsi & 0x4) {
  737.       for (sfbi = 6; sfbi < 11; ++sfbi)
  738. channel->scalefac[sfbi] = gr0ch->scalefac[sfbi];
  739.     }
  740.     else {
  741.       for (sfbi = 6; sfbi < 11; ++sfbi)
  742. channel->scalefac[sfbi] = mad_bit_read(ptr, slen1);
  743.     }
  745.     if (scfsi & 0x2) {
  746.       for (sfbi = 11; sfbi < 16; ++sfbi)
  747. channel->scalefac[sfbi] = gr0ch->scalefac[sfbi];
  748.     }
  749.     else {
  750.       for (sfbi = 11; sfbi < 16; ++sfbi)
  751. channel->scalefac[sfbi] = mad_bit_read(ptr, slen2);
  752.     }
  754.     if (scfsi & 0x1) {
  755.       for (sfbi = 16; sfbi < 21; ++sfbi)
  756. channel->scalefac[sfbi] = gr0ch->scalefac[sfbi];
  757.     }
  758.     else {
  759.       for (sfbi = 16; sfbi < 21; ++sfbi)
  760. channel->scalefac[sfbi] = mad_bit_read(ptr, slen2);
  761.     }
  763.     channel->scalefac[21] = 0;
  764.   }
  766.   return mad_bit_length(&start, ptr);
  767. }
  769. /*
  770.  * The Layer III formula for requantization and scaling is defined by
  771.  * section 2.4.3.4.7.1 of ISO/IEC 11172-3, as follows:
  772.  *
  773.  *   long blocks:
  774.  *   xr[i] = sign(is[i]) * abs(is[i])^(4/3) *
  775.  *           2^((1/4) * (global_gain - 210)) *
  776.  *           2^-(scalefac_multiplier *
  777.  *               (scalefac_l[sfb] + preflag * pretab[sfb]))
  778.  *
  779.  *   short blocks:
  780.  *   xr[i] = sign(is[i]) * abs(is[i])^(4/3) *
  781.  *           2^((1/4) * (global_gain - 210 - 8 * subblock_gain[w])) *
  782.  *           2^-(scalefac_multiplier * scalefac_s[sfb][w])
  783.  *
  784.  *   where:
  785.  *   scalefac_multiplier = (scalefac_scale + 1) / 2
  786.  *
  787.  * The routines III_exponents() and III_requantize() facilitate this
  788.  * calculation.
  789.  */
  791. /*
  792.  * NAME: III_exponents()
  793.  * DESCRIPTION: calculate scalefactor exponents
  794.  */
  795. static
  796. void III_exponents(struct channel const *channel,
  797.    unsigned char const *sfbwidth, signed int exponents[39])
  798. {
  799.   signed int gain;
  800.   unsigned int scalefac_multiplier, sfbi;
  802.   gain = (signed int) channel->global_gain - 210;
  803.   scalefac_multiplier = (channel->flags & scalefac_scale) ? 2 : 1;
  805.   if (channel->block_type == 2) {
  806.     unsigned int l;
  807.     signed int gain0, gain1, gain2;
  809.     sfbi = l = 0;
  811.     if (channel->flags & mixed_block_flag) {
  812.       unsigned int premask;
  814.       premask = (channel->flags & preflag) ? ~0 : 0;
  816.       /* long block subbands 0-1 */
  818.       while (l < 36) {
  819. exponents[sfbi] = gain -
  820.   (signed int) ((channel->scalefac[sfbi] + (pretab[sfbi] & premask)) <<
  821. scalefac_multiplier);
  823. l += sfbwidth[sfbi++];
  824.       }
  825.     }
  827.     /* this is probably wrong for 8000 Hz short/mixed blocks */
  829.     gain0 = gain - 8 * (signed int) channel->subblock_gain[0];
  830.     gain1 = gain - 8 * (signed int) channel->subblock_gain[1];
  831.     gain2 = gain - 8 * (signed int) channel->subblock_gain[2];
  833.     while (l < 576) {
  834.       exponents[sfbi + 0] = gain0 -
  835. (signed int) (channel->scalefac[sfbi + 0] << scalefac_multiplier);
  836.       exponents[sfbi + 1] = gain1 -
  837. (signed int) (channel->scalefac[sfbi + 1] << scalefac_multiplier);
  838.       exponents[sfbi + 2] = gain2 -
  839. (signed int) (channel->scalefac[sfbi + 2] << scalefac_multiplier);
  841.       l    += 3 * sfbwidth[sfbi];
  842.       sfbi += 3;
  843.     }
  844.   }
  845.   else {  /* channel->block_type != 2 */
  846.     if (channel->flags & preflag) {
  847.       for (sfbi = 0; sfbi < 22; ++sfbi) {
  848. exponents[sfbi] = gain -
  849.   (signed int) ((channel->scalefac[sfbi] + pretab[sfbi]) <<
  850. scalefac_multiplier);
  851.       }
  852.     }
  853.     else {
  854.       for (sfbi = 0; sfbi < 22; ++sfbi) {
  855. exponents[sfbi] = gain -
  856.   (signed int) (channel->scalefac[sfbi] << scalefac_multiplier);
  857.       }
  858.     }
  859.   }
  860. }
  862. /*
  863.  * NAME: III_requantize()
  864.  * DESCRIPTION: requantize one (positive) value
  865.  */
  866. static
  867. mad_fixed_t III_requantize(unsigned int value, signed int exp)
  868. {
  869.   mad_fixed_t requantized;
  870.   signed int frac;
  871.   struct fixedfloat const *power;
  873.   frac = exp % 4;  /* assumes sign(frac) == sign(exp) */
  874.   exp /= 4;
  876.   power = &rq_table[value];
  877.   requantized = power->mantissa;
  878.   exp += power->exponent;
  880.   if (exp < 0) {
  881.     if (-exp >= sizeof(mad_fixed_t) * CHAR_BIT) {
  882.       /* underflow */
  883.       requantized = 0;
  884.     }
  885.     else {
  886.       requantized += 1L << (-exp - 1);
  887.       requantized >>= -exp;
  888.     }
  889.   }
  890.   else {
  891.     if (exp >= 5) {
  892.       /* overflow */
  893. # if defined(DEBUG)
  894.       fprintf(stderr, "requantize overflow (%f * 2^%d)n",
  895.       mad_f_todouble(requantized), exp);
  896. # endif
  897.       requantized = MAD_F_MAX;
  898.     }
  899.     else
  900.       requantized <<= exp;
  901.   }
  903.   return frac ? mad_f_mul(requantized, root_table[3 + frac]) : requantized;
  904. }
  906. /* we must take care that sz >= bits and sz < sizeof(long) lest bits == 0 */
  907. # define MASK(cache, sz, bits)
  908.     (((cache) >> ((sz) - (bits))) & ((1 << (bits)) - 1))
  909. # define MASK1BIT(cache, sz)  
  910.     ((cache) & (1 << ((sz) - 1)))
  912. /*
  913.  * NAME: III_huffdecode()
  914.  * DESCRIPTION: decode Huffman code words of one channel of one granule
  915.  */
  916. static
  917. enum mad_error III_huffdecode(struct mad_bitptr *ptr, mad_fixed_t xr[576],
  918.       struct channel *channel,
  919.       unsigned char const *sfbwidth,
  920.       unsigned int part2_length)
  921. {
  922.   signed int exponents[39], exp;
  923.   signed int const *expptr;
  924.   struct mad_bitptr peek;
  925.   signed int bits_left, cachesz;
  926.   register mad_fixed_t *xrptr;
  927.   mad_fixed_t const *sfbound;
  928.   register unsigned long bitcache;
  930.   bits_left = (signed) channel->part2_3_length - (signed) part2_length;
  931.   if (bits_left < 0)
  932.     return MAD_ERROR_BADPART3LEN;
  934.   III_exponents(channel, sfbwidth, exponents);
  936.   peek = *ptr;
  937.   mad_bit_skip(ptr, bits_left);
  939.   /* align bit reads to byte boundaries */
  940.   cachesz  = mad_bit_bitsleft(&peek);
  941.   cachesz += ((32 - 1 - 24) + (24 - cachesz)) & ~7;
  943.   bitcache   = mad_bit_read(&peek, cachesz);
  944.   bits_left -= cachesz;
  946.   xrptr = &xr[0];
  948.   /* big_values */
  949.   {
  950.     unsigned int region, rcount;
  951.     struct hufftable const *entry;
  952.     union huffpair const *table;
  953.     unsigned int linbits, startbits, big_values, reqhits;
  954.     mad_fixed_t reqcache[16];
  956.     sfbound = xrptr + *sfbwidth++;
  957.     rcount  = channel->region0_count + 1;
  959.     entry     = &mad_huff_pair_table[channel->table_select[region = 0]];
  960.     table     = entry->table;
  961.     linbits   = entry->linbits;
  962.     startbits = entry->startbits;
  964.     if (table == 0)
  965.       return MAD_ERROR_BADHUFFTABLE;
  967.     expptr  = &exponents[0];
  968.     exp     = *expptr++;
  969.     reqhits = 0;
  971.     big_values = channel->big_values;
  973.     while (big_values-- && cachesz + bits_left > 0) {
  974.       union huffpair const *pair;
  975.       unsigned int clumpsz, value;
  976.       register mad_fixed_t requantized;
  978.       if (xrptr == sfbound) {
  979. sfbound += *sfbwidth++;
  981. /* change table if region boundary */
  983. if (--rcount == 0) {
  984.   if (region == 0)
  985.     rcount = channel->region1_count + 1;
  986.   else
  987.     rcount = 0;  /* all remaining */
  989.   entry     = &mad_huff_pair_table[channel->table_select[++region]];
  990.   table     = entry->table;
  991.   linbits   = entry->linbits;
  992.   startbits = entry->startbits;
  994.   if (table == 0)
  995.     return MAD_ERROR_BADHUFFTABLE;
  996. }
  998. if (exp != *expptr) {
  999.   exp = *expptr;
  1000.   reqhits = 0;
  1001. }
  1003. ++expptr;
  1004.       }
  1006.       if (cachesz < 21) {
  1007. unsigned int bits;
  1009. bits       = ((32 - 1 - 21) + (21 - cachesz)) & ~7;
  1010. bitcache   = (bitcache << bits) | mad_bit_read(&peek, bits);
  1011. cachesz   += bits;
  1012. bits_left -= bits;
  1013.       }
  1015.       /* hcod (0..19) */
  1017.       clumpsz = startbits;
  1018.       pair    = &table[MASK(bitcache, cachesz, clumpsz)];
  1020.       while (!pair->final) {
  1021. cachesz -= clumpsz;
  1023. clumpsz = pair->ptr.bits;
  1024. pair    = &table[pair->ptr.offset + MASK(bitcache, cachesz, clumpsz)];
  1025.       }
  1027.       cachesz -= pair->value.hlen;
  1029.       if (linbits) {
  1030. /* x (0..14) */
  1032. value = pair->value.x;
  1034. switch (value) {
  1035. case 0:
  1036.   xrptr[0] = 0;
  1037.   break;
  1039. case 15:
  1040.   if (cachesz < linbits + 2) {
  1041.     bitcache   = (bitcache << 16) | mad_bit_read(&peek, 16);
  1042.     cachesz   += 16;
  1043.     bits_left -= 16;
  1044.   }
  1046.   value += MASK(bitcache, cachesz, linbits);
  1047.   cachesz -= linbits;
  1049.   requantized = III_requantize(value, exp);
  1050.   goto x_final;
  1052. default:
  1053.   if (reqhits & (1 << value))
  1054.     requantized = reqcache[value];
  1055.   else {
  1056.     reqhits |= (1 << value);
  1057.     requantized = reqcache[value] = III_requantize(value, exp);
  1058.   }
  1060. x_final:
  1061.   xrptr[0] = MASK1BIT(bitcache, cachesz--) ?
  1062.     -requantized : requantized;
  1063. }
  1065. /* y (0..14) */
  1067. value = pair->value.y;
  1069. switch (value) {
  1070. case 0:
  1071.   xrptr[1] = 0;
  1072.   break;
  1074. case 15:
  1075.   if (cachesz < linbits + 1) {
  1076.     bitcache   = (bitcache << 16) | mad_bit_read(&peek, 16);
  1077.     cachesz   += 16;
  1078.     bits_left -= 16;
  1079.   }
  1081.   value += MASK(bitcache, cachesz, linbits);
  1082.   cachesz -= linbits;
  1084.   requantized = III_requantize(value, exp);
  1085.   goto y_final;
  1087. default:
  1088.   if (reqhits & (1 << value))
  1089.     requantized = reqcache[value];
  1090.   else {
  1091.     reqhits |= (1 << value);
  1092.     requantized = reqcache[value] = III_requantize(value, exp);
  1093.   }
  1095. y_final:
  1096.   xrptr[1] = MASK1BIT(bitcache, cachesz--) ?
  1097.     -requantized : requantized;
  1098. }
  1099.       }
  1100.       else {
  1101. /* x (0..1) */
  1103. value = pair->value.x;
  1105. if (value == 0)
  1106.   xrptr[0] = 0;
  1107. else {
  1108.   if (reqhits & (1 << value))
  1109.     requantized = reqcache[value];
  1110.   else {
  1111.     reqhits |= (1 << value);
  1112.     requantized = reqcache[value] = III_requantize(value, exp);
  1113.   }
  1115.   xrptr[0] = MASK1BIT(bitcache, cachesz--) ?
  1116.     -requantized : requantized;
  1117. }
  1119. /* y (0..1) */
  1121. value = pair->value.y;
  1123. if (value == 0)
  1124.   xrptr[1] = 0;
  1125. else {
  1126.   if (reqhits & (1 << value))
  1127.     requantized = reqcache[value];
  1128.   else {
  1129.     reqhits |= (1 << value);
  1130.     requantized = reqcache[value] = III_requantize(value, exp);
  1131.   }
  1133.   xrptr[1] = MASK1BIT(bitcache, cachesz--) ?
  1134.     -requantized : requantized;
  1135. }
  1136.       }
  1138.       xrptr += 2;
  1139.     }
  1140.   }
  1142.   if (cachesz + bits_left < 0)
  1143.     return MAD_ERROR_BADHUFFDATA;  /* big_values overrun */
  1145.   /* count1 */
  1146.   {
  1147.     union huffquad const *table;
  1148.     register mad_fixed_t requantized;
  1150.     table = mad_huff_quad_table[channel->flags & count1table_select];
  1152.     requantized = III_requantize(1, exp);
  1154.     while (cachesz + bits_left > 0 && xrptr <= &xr[572]) {
  1155.       union huffquad const *quad;
  1157.       /* hcod (1..6) */
  1159.       if (cachesz < 10) {
  1160. bitcache   = (bitcache << 16) | mad_bit_read(&peek, 16);
  1161. cachesz   += 16;
  1162. bits_left -= 16;
  1163.       }
  1165.       quad = &table[MASK(bitcache, cachesz, 4)];
  1167.       /* quad tables guaranteed to have at most one extra lookup */
  1168.       if (!quad->final) {
  1169. cachesz -= 4;
  1171. quad = &table[quad->ptr.offset +
  1172.       MASK(bitcache, cachesz, quad->ptr.bits)];
  1173.       }
  1175.       cachesz -= quad->value.hlen;
  1177.       if (xrptr == sfbound) {
  1178. sfbound += *sfbwidth++;
  1180. if (exp != *expptr) {
  1181.   exp = *expptr;
  1182.   requantized = III_requantize(1, exp);
  1183. }
  1185. ++expptr;
  1186.       }
  1188.       /* v (0..1) */
  1190.       xrptr[0] = quad->value.v ?
  1191. (MASK1BIT(bitcache, cachesz--) ? -requantized : requantized) : 0;
  1193.       /* w (0..1) */
  1195.       xrptr[1] = quad->value.w ?
  1196. (MASK1BIT(bitcache, cachesz--) ? -requantized : requantized) : 0;
  1198.       xrptr += 2;
  1200.       if (xrptr == sfbound) {
  1201. sfbound += *sfbwidth++;
  1203. if (exp != *expptr) {
  1204.   exp = *expptr;
  1205.   requantized = III_requantize(1, exp);
  1206. }
  1208. ++expptr;
  1209.       }
  1211.       /* x (0..1) */
  1213.       xrptr[0] = quad->value.x ?
  1214. (MASK1BIT(bitcache, cachesz--) ? -requantized : requantized) : 0;
  1216.       /* y (0..1) */
  1218.       xrptr[1] = quad->value.y ?
  1219. (MASK1BIT(bitcache, cachesz--) ? -requantized : requantized) : 0;
  1221.       xrptr += 2;
  1222.     }
  1224.     if (cachesz + bits_left < 0) {
  1225. # if 0 && defined(DEBUG)
  1226.       fprintf(stderr, "huffman count1 overrun (%d bits)n",
  1227.       -(cachesz + bits_left));
  1228. # endif
  1230.       /* technically the bitstream is misformatted, but apparently
  1231.  some encoders are just a bit sloppy with stuffing bits */
  1233.       xrptr -= 4;
  1234.     }
  1235.   }
  1237.   assert(-bits_left <= MAD_BUFFER_GUARD * CHAR_BIT);
  1239. # if 0 && defined(DEBUG)
  1240.   if (bits_left < 0)
  1241.     fprintf(stderr, "read %d bits too manyn", -bits_left);
  1242.   else if (cachesz + bits_left > 0)
  1243.     fprintf(stderr, "%d stuffing bitsn", cachesz + bits_left);
  1244. # endif
  1246.   /* rzero */
  1247.   while (xrptr < &xr[576]) {
  1248.     xrptr[0] = 0;
  1249.     xrptr[1] = 0;
  1251.     xrptr += 2;
  1252.   }
  1254.   return MAD_ERROR_NONE;
  1255. }
  1257. # undef MASK
  1258. # undef MASK1BIT
  1260. /*
  1261.  * NAME: III_reorder()
  1262.  * DESCRIPTION: reorder frequency lines of a short block into subband order
  1263.  */
  1264. static
  1265. void III_reorder(mad_fixed_t xr[576], struct channel const *channel,
  1266.  unsigned char const sfbwidth[39])
  1267. {
  1268.   mad_fixed_t tmp[32][3][6];
  1269.   unsigned int sb, l, f, w, sbw[3], sw[3];
  1271.   /* this is probably wrong for 8000 Hz mixed blocks */
  1273.   sb = 0;
  1274.   if (channel->flags & mixed_block_flag) {
  1275.     sb = 2;
  1277.     l = 0;
  1278.     while (l < 36)
  1279.       l += *sfbwidth++;
  1280.   }
  1282.   for (w = 0; w < 3; ++w) {
  1283.     sbw[w] = sb;
  1284.     sw[w]  = 0;
  1285.   }
  1287.   f = *sfbwidth++;
  1288.   w = 0;
  1290.   for (l = 18 * sb; l < 576; ++l) {
  1291.     if (f-- == 0) {
  1292.       f = *sfbwidth++ - 1;
  1293.       w = (w + 1) % 3;
  1294.     }
  1296.     tmp[sbw[w]][w][sw[w]++] = xr[l];
  1298.     if (sw[w] == 6) {
  1299.       sw[w] = 0;
  1300.       ++sbw[w];
  1301.     }
  1302.   }
  1304.   memcpy(&xr[18 * sb], &tmp[sb], (576 - 18 * sb) * sizeof(mad_fixed_t));
  1305. }
  1307. /*
  1308.  * NAME: III_stereo()
  1309.  * DESCRIPTION: perform joint stereo processing on a granule
  1310.  */
  1311. static
  1312. enum mad_error III_stereo(mad_fixed_t xr[2][576],
  1313.   struct granule const *granule,
  1314.   struct mad_header *header,
  1315.   unsigned char const *sfbwidth)
  1316. {
  1317.   short modes[39];
  1318.   unsigned int sfbi, l, n, i;
  1320.   if (granule->ch[0].block_type !=
  1321.       granule->ch[1].block_type ||
  1322.       (granule->ch[0].flags & mixed_block_flag) !=
  1323.       (granule->ch[1].flags & mixed_block_flag))
  1324.     return MAD_ERROR_BADSTEREO;
  1326.   for (i = 0; i < 39; ++i)
  1327.     modes[i] = header->mode_extension;
  1329.   /* intensity stereo */
  1331.   if (header->mode_extension & I_STEREO) {
  1332.     struct channel const *right_ch = &granule->ch[1];
  1333.     mad_fixed_t const *right_xr = xr[1];
  1334.     unsigned int is_pos;
  1336.     header->flags |= MAD_FLAG_I_STEREO;
  1338.     /* first determine which scalefactor bands are to be processed */
  1340.     if (right_ch->block_type == 2) {
  1341.       unsigned int lower, start, max, bound[3], w;
  1343.       lower = start = max = bound[0] = bound[1] = bound[2] = 0;
  1345.       sfbi = l = 0;
  1347.       if (right_ch->flags & mixed_block_flag) {
  1348. while (l < 36) {
  1349.   n = sfbwidth[sfbi++];
  1351.   for (i = 0; i < n; ++i) {
  1352.     if (right_xr[i]) {
  1353.       lower = sfbi;
  1354.       break;
  1355.     }
  1356.   }
  1358.   right_xr += n;
  1359.   l += n;
  1360. }
  1362. start = sfbi;
  1363.       }
  1365.       w = 0;
  1366.       while (l < 576) {
  1367. n = sfbwidth[sfbi++];
  1369. for (i = 0; i < n; ++i) {
  1370.   if (right_xr[i]) {
  1371.     max = bound[w] = sfbi;
  1372.     break;
  1373.   }
  1374. }
  1376. right_xr += n;
  1377. l += n;
  1378. w = (w + 1) % 3;
  1379.       }
  1381.       if (max)
  1382. lower = start;
  1384.       /* long blocks */
  1386.       for (i = 0; i < lower; ++i)
  1387. modes[i] = header->mode_extension & ~I_STEREO;
  1389.       /* short blocks */
  1391.       w = 0;
  1392.       for (i = start; i < max; ++i) {
  1393. if (i < bound[w])
  1394.   modes[i] = header->mode_extension & ~I_STEREO;
  1396. w = (w + 1) % 3;
  1397.       }
  1398.     }
  1399.     else {  /* right_ch->block_type != 2 */
  1400.       unsigned int bound;
  1402.       bound = 0;
  1403.       for (sfbi = l = 0; l < 576; l += n) {
  1404. n = sfbwidth[sfbi++];
  1406. for (i = 0; i < n; ++i) {
  1407.   if (right_xr[i]) {
  1408.     bound = sfbi;
  1409.     break;
  1410.   }
  1411. }
  1413. right_xr += n;
  1414.       }
  1416.       for (i = 0; i < bound; ++i)
  1417. modes[i] = header->mode_extension & ~I_STEREO;
  1418.     }
  1420.     /* now do the actual processing */
  1422.     if (header->flags & MAD_FLAG_LSF_EXT) {
  1423.       unsigned char const *illegal_pos = granule[1].ch[1].scalefac;
  1424.       mad_fixed_t const *lsf_scale;
  1426.       /* intensity_scale */
  1427.       lsf_scale = is_lsf_table[right_ch->scalefac_compress & 0x1];
  1429.       for (sfbi = l = 0; l < 576; ++sfbi, l += n) {
  1430. n = sfbwidth[sfbi];
  1432. if (!(modes[sfbi] & I_STEREO))
  1433.   continue;
  1435. if (illegal_pos[sfbi]) {
  1436.   modes[sfbi] &= ~I_STEREO;
  1437.   continue;
  1438. }
  1440. is_pos = right_ch->scalefac[sfbi];
  1442. for (i = 0; i < n; ++i) {
  1443.   register mad_fixed_t left;
  1445.   left = xr[0][l + i];
  1447.   if (is_pos == 0)
  1448.     xr[1][l + i] = left;
  1449.   else {
  1450.     register mad_fixed_t opposite;
  1452.     opposite = mad_f_mul(left, lsf_scale[(is_pos - 1) / 2]);
  1454.     if (is_pos & 1) {
  1455.       xr[0][l + i] = opposite;
  1456.       xr[1][l + i] = left;
  1457.     }
  1458.     else
  1459.       xr[1][l + i] = opposite;
  1460.   }
  1461. }
  1462.       }
  1463.     }
  1464.     else {  /* !(header->flags & MAD_FLAG_LSF_EXT) */
  1465.       for (sfbi = l = 0; l < 576; ++sfbi, l += n) {
  1466. n = sfbwidth[sfbi];
  1468. if (!(modes[sfbi] & I_STEREO))
  1469.   continue;
  1471. is_pos = right_ch->scalefac[sfbi];
  1473. if (is_pos >= 7) {  /* illegal intensity position */
  1474.   modes[sfbi] &= ~I_STEREO;
  1475.   continue;
  1476. }
  1478. for (i = 0; i < n; ++i) {
  1479.   register mad_fixed_t left;
  1481.   left = xr[0][l + i];
  1483.   xr[0][l + i] = mad_f_mul(left, is_table[    is_pos]);
  1484.   xr[1][l + i] = mad_f_mul(left, is_table[6 - is_pos]);
  1485. }
  1486.       }
  1487.     }
  1488.   }
  1490.   /* middle/side stereo */
  1492.   if (header->mode_extension & MS_STEREO) {
  1493.     register mad_fixed_t invsqrt2;
  1495.     header->flags |= MAD_FLAG_MS_STEREO;
  1497.     invsqrt2 = root_table[3 + -2];
  1499.     for (sfbi = l = 0; l < 576; ++sfbi, l += n) {
  1500.       n = sfbwidth[sfbi];
  1502.       if (modes[sfbi] != MS_STEREO)
  1503. continue;
  1505.       for (i = 0; i < n; ++i) {
  1506. register mad_fixed_t m, s;
  1508. m = xr[0][l + i];
  1509. s = xr[1][l + i];
  1511. xr[0][l + i] = mad_f_mul(m + s, invsqrt2);  /* l = (m + s) / sqrt(2) */
  1512. xr[1][l + i] = mad_f_mul(m - s, invsqrt2);  /* r = (m - s) / sqrt(2) */
  1513.       }
  1514.     }
  1515.   }
  1517.   return MAD_ERROR_NONE;
  1518. }
  1520. /*
  1521.  * NAME: III_aliasreduce()
  1522.  * DESCRIPTION: perform frequency line alias reduction
  1523.  */
  1524. static
  1525. void III_aliasreduce(mad_fixed_t xr[576], int lines)
  1526. {
  1527.   mad_fixed_t const *bound;
  1528.   int i;
  1530.   bound = &xr[lines];
  1531.   for (xr += 18; xr < bound; xr += 18) {
  1532.     for (i = 0; i < 8; ++i) {
  1533.       register mad_fixed_t a, b;
  1534.       register mad_fixed64hi_t hi;
  1535.       register mad_fixed64lo_t lo;
  1537.       a = xr[-1 - i];
  1538.       b = xr[     i];
  1540. # if defined(ASO_ZEROCHECK)
  1541.       if (a | b) {
  1542. # endif
  1543. MAD_F_ML0(hi, lo,  a, cs[i]);
  1544. MAD_F_MLA(hi, lo, -b, ca[i]);
  1546. xr[-1 - i] = MAD_F_MLZ(hi, lo);
  1548. MAD_F_ML0(hi, lo,  b, cs[i]);
  1549. MAD_F_MLA(hi, lo,  a, ca[i]);
  1551. xr[     i] = MAD_F_MLZ(hi, lo);
  1552. # if defined(ASO_ZEROCHECK)
  1553.       }
  1554. # endif
  1555.     }
  1556.   }
  1557. }
  1559. # if defined(ASO_IMDCT)
  1560. void III_imdct_l(mad_fixed_t const [18], mad_fixed_t [36], unsigned int);
  1561. # else
  1562. /*
  1563.  * NAME: imdct36
  1564.  * DESCRIPTION: perform X[18]->x[36] IMDCT
  1565.  */
  1566. static inline
  1567. void imdct36(mad_fixed_t const X[18], mad_fixed_t x[36])
  1568. {
  1569.   mad_fixed_t t0, t1, t2,  t3,  t4,  t5,  t6,  t7;
  1570.   mad_fixed_t t8, t9, t10, t11, t12, t13, t14, t15;
  1571.   register mad_fixed64hi_t hi;
  1572.   register mad_fixed64lo_t lo;
  1574.   MAD_F_ML0(hi, lo, X[4],  MAD_F(0x0ec835e8));
  1575.   MAD_F_MLA(hi, lo, X[13], MAD_F(0x061f78aa));
  1577.   t6 = MAD_F_MLZ(hi, lo);
  1579.   MAD_F_MLA(hi, lo, (t14 = X[1] - X[10]), -MAD_F(0x061f78aa));
  1580.   MAD_F_MLA(hi, lo, (t15 = X[7] + X[16]), -MAD_F(0x0ec835e8));
  1582.   t0 = MAD_F_MLZ(hi, lo);
  1584.   MAD_F_MLA(hi, lo, (t8  = X[0] - X[11] - X[12]),  MAD_F(0x0216a2a2));
  1585.   MAD_F_MLA(hi, lo, (t9  = X[2] - X[9]  - X[14]),  MAD_F(0x09bd7ca0));
  1586.   MAD_F_MLA(hi, lo, (t10 = X[3] - X[8]  - X[15]), -MAD_F(0x0cb19346));
  1587.   MAD_F_MLA(hi, lo, (t11 = X[5] - X[6]  - X[17]), -MAD_F(0x0fdcf549));
  1589.   x[7]  = MAD_F_MLZ(hi, lo);
  1590.   x[10] = -x[7];
  1592.   MAD_F_ML0(hi, lo, t8,  -MAD_F(0x0cb19346));
  1593.   MAD_F_MLA(hi, lo, t9,   MAD_F(0x0fdcf549));
  1594.   MAD_F_MLA(hi, lo, t10,  MAD_F(0x0216a2a2));
  1595.   MAD_F_MLA(hi, lo, t11, -MAD_F(0x09bd7ca0));
  1597.   x[19] = x[34] = MAD_F_MLZ(hi, lo) - t0;
  1599.   t12 = X[0] - X[3] + X[8] - X[11] - X[12] + X[15];
  1600.   t13 = X[2] + X[5] - X[6] - X[9]  - X[14] - X[17];
  1602.   MAD_F_ML0(hi, lo, t12, -MAD_F(0x0ec835e8));
  1603.   MAD_F_MLA(hi, lo, t13,  MAD_F(0x061f78aa));
  1605.   x[22] = x[31] = MAD_F_MLZ(hi, lo) + t0;
  1607.   MAD_F_ML0(hi, lo, X[1],  -MAD_F(0x09bd7ca0));
  1608.   MAD_F_MLA(hi, lo, X[7],   MAD_F(0x0216a2a2));
  1609.   MAD_F_MLA(hi, lo, X[10], -MAD_F(0x0fdcf549));
  1610.   MAD_F_MLA(hi, lo, X[16],  MAD_F(0x0cb19346));
  1612.   t1 = MAD_F_MLZ(hi, lo) + t6;
  1614.   MAD_F_ML0(hi, lo, X[0],   MAD_F(0x03768962));
  1615.   MAD_F_MLA(hi, lo, X[2],   MAD_F(0x0e313245));
  1616.   MAD_F_MLA(hi, lo, X[3],  -MAD_F(0x0ffc19fd));
  1617.   MAD_F_MLA(hi, lo, X[5],  -MAD_F(0x0acf37ad));
  1618.   MAD_F_MLA(hi, lo, X[6],   MAD_F(0x04cfb0e2));
  1619.   MAD_F_MLA(hi, lo, X[8],  -MAD_F(0x0898c779));
  1620.   MAD_F_MLA(hi, lo, X[9],   MAD_F(0x0d7e8807));
  1621.   MAD_F_MLA(hi, lo, X[11],  MAD_F(0x0f426cb5));
  1622.   MAD_F_MLA(hi, lo, X[12], -MAD_F(0x0bcbe352));
  1623.   MAD_F_MLA(hi, lo, X[14],  MAD_F(0x00b2aa3e));
  1624.   MAD_F_MLA(hi, lo, X[15], -MAD_F(0x07635284));
  1625.   MAD_F_MLA(hi, lo, X[17], -MAD_F(0x0f9ee890));
  1627.   x[6]  = MAD_F_MLZ(hi, lo) + t1;
  1628.   x[11] = -x[6];
  1630.   MAD_F_ML0(hi, lo, X[0],  -MAD_F(0x0f426cb5));
  1631.   MAD_F_MLA(hi, lo, X[2],  -MAD_F(0x00b2aa3e));
  1632.   MAD_F_MLA(hi, lo, X[3],   MAD_F(0x0898c779));
  1633.   MAD_F_MLA(hi, lo, X[5],   MAD_F(0x0f9ee890));
  1634.   MAD_F_MLA(hi, lo, X[6],   MAD_F(0x0acf37ad));
  1635.   MAD_F_MLA(hi, lo, X[8],  -MAD_F(0x07635284));
  1636.   MAD_F_MLA(hi, lo, X[9],  -MAD_F(0x0e313245));
  1637.   MAD_F_MLA(hi, lo, X[11], -MAD_F(0x0bcbe352));
  1638.   MAD_F_MLA(hi, lo, X[12], -MAD_F(0x03768962));
  1639.   MAD_F_MLA(hi, lo, X[14],  MAD_F(0x0d7e8807));
  1640.   MAD_F_MLA(hi, lo, X[15],  MAD_F(0x0ffc19fd));
  1641.   MAD_F_MLA(hi, lo, X[17],  MAD_F(0x04cfb0e2));
  1643.   x[23] = x[30] = MAD_F_MLZ(hi, lo) + t1;
  1645.   MAD_F_ML0(hi, lo, X[0],  -MAD_F(0x0bcbe352));
  1646.   MAD_F_MLA(hi, lo, X[2],   MAD_F(0x0d7e8807));
  1647.   MAD_F_MLA(hi, lo, X[3],  -MAD_F(0x07635284));
  1648.   MAD_F_MLA(hi, lo, X[5],   MAD_F(0x04cfb0e2));
  1649.   MAD_F_MLA(hi, lo, X[6],   MAD_F(0x0f9ee890));
  1650.   MAD_F_MLA(hi, lo, X[8],  -MAD_F(0x0ffc19fd));
  1651.   MAD_F_MLA(hi, lo, X[9],  -MAD_F(0x00b2aa3e));
  1652.   MAD_F_MLA(hi, lo, X[11],  MAD_F(0x03768962));
  1653.   MAD_F_MLA(hi, lo, X[12], -MAD_F(0x0f426cb5));
  1654.   MAD_F_MLA(hi, lo, X[14],  MAD_F(0x0e313245));
  1655.   MAD_F_MLA(hi, lo, X[15],  MAD_F(0x0898c779));
  1656.   MAD_F_MLA(hi, lo, X[17], -MAD_F(0x0acf37ad));
  1658.   x[18] = x[35] = MAD_F_MLZ(hi, lo) - t1;
  1660.   MAD_F_ML0(hi, lo, X[4],   MAD_F(0x061f78aa));
  1661.   MAD_F_MLA(hi, lo, X[13], -MAD_F(0x0ec835e8));
  1663.   t7 = MAD_F_MLZ(hi, lo);
  1665.   MAD_F_MLA(hi, lo, X[1],  -MAD_F(0x0cb19346));
  1666.   MAD_F_MLA(hi, lo, X[7],   MAD_F(0x0fdcf549));
  1667.   MAD_F_MLA(hi, lo, X[10],  MAD_F(0x0216a2a2));
  1668.   MAD_F_MLA(hi, lo, X[16], -MAD_F(0x09bd7ca0));
  1670.   t2 = MAD_F_MLZ(hi, lo);
  1672.   MAD_F_MLA(hi, lo, X[0],   MAD_F(0x04cfb0e2));
  1673.   MAD_F_MLA(hi, lo, X[2],   MAD_F(0x0ffc19fd));
  1674.   MAD_F_MLA(hi, lo, X[3],  -MAD_F(0x0d7e8807));
  1675.   MAD_F_MLA(hi, lo, X[5],   MAD_F(0x03768962));
  1676.   MAD_F_MLA(hi, lo, X[6],  -MAD_F(0x0bcbe352));
  1677.   MAD_F_MLA(hi, lo, X[8],  -MAD_F(0x0e313245));
  1678.   MAD_F_MLA(hi, lo, X[9],   MAD_F(0x07635284));
  1679.   MAD_F_MLA(hi, lo, X[11], -MAD_F(0x0acf37ad));
  1680.   MAD_F_MLA(hi, lo, X[12],  MAD_F(0x0f9ee890));
  1681.   MAD_F_MLA(hi, lo, X[14],  MAD_F(0x0898c779));
  1682.   MAD_F_MLA(hi, lo, X[15],  MAD_F(0x00b2aa3e));
  1683.   MAD_F_MLA(hi, lo, X[17],  MAD_F(0x0f426cb5));
  1685.   x[5]  = MAD_F_MLZ(hi, lo);
  1686.   x[12] = -x[5];
  1688.   MAD_F_ML0(hi, lo, X[0],   MAD_F(0x0acf37ad));
  1689.   MAD_F_MLA(hi, lo, X[2],  -MAD_F(0x0898c779));
  1690.   MAD_F_MLA(hi, lo, X[3],   MAD_F(0x0e313245));
  1691.   MAD_F_MLA(hi, lo, X[5],  -MAD_F(0x0f426cb5));
  1692.   MAD_F_MLA(hi, lo, X[6],  -MAD_F(0x03768962));
  1693.   MAD_F_MLA(hi, lo, X[8],   MAD_F(0x00b2aa3e));
  1694.   MAD_F_MLA(hi, lo, X[9],  -MAD_F(0x0ffc19fd));
  1695.   MAD_F_MLA(hi, lo, X[11],  MAD_F(0x0f9ee890));
  1696.   MAD_F_MLA(hi, lo, X[12], -MAD_F(0x04cfb0e2));
  1697.   MAD_F_MLA(hi, lo, X[14],  MAD_F(0x07635284));
  1698.   MAD_F_MLA(hi, lo, X[15],  MAD_F(0x0d7e8807));
  1699.   MAD_F_MLA(hi, lo, X[17], -MAD_F(0x0bcbe352));
  1701.   x[0]  = MAD_F_MLZ(hi, lo) + t2;
  1702.   x[17] = -x[0];
  1704.   MAD_F_ML0(hi, lo, X[0],  -MAD_F(0x0f9ee890));
  1705.   MAD_F_MLA(hi, lo, X[2],  -MAD_F(0x07635284));
  1706.   MAD_F_MLA(hi, lo, X[3],  -MAD_F(0x00b2aa3e));
  1707.   MAD_F_MLA(hi, lo, X[5],   MAD_F(0x0bcbe352));
  1708.   MAD_F_MLA(hi, lo, X[6],   MAD_F(0x0f426cb5));
  1709.   MAD_F_MLA(hi, lo, X[8],   MAD_F(0x0d7e8807));
  1710.   MAD_F_MLA(hi, lo, X[9],   MAD_F(0x0898c779));
  1711.   MAD_F_MLA(hi, lo, X[11], -MAD_F(0x04cfb0e2));
  1712.   MAD_F_MLA(hi, lo, X[12], -MAD_F(0x0acf37ad));
  1713.   MAD_F_MLA(hi, lo, X[14], -MAD_F(0x0ffc19fd));
  1714.   MAD_F_MLA(hi, lo, X[15], -MAD_F(0x0e313245));
  1715.   MAD_F_MLA(hi, lo, X[17], -MAD_F(0x03768962));
  1717.   x[24] = x[29] = MAD_F_MLZ(hi, lo) + t2;
  1719.   MAD_F_ML0(hi, lo, X[1],  -MAD_F(0x0216a2a2));
  1720.   MAD_F_MLA(hi, lo, X[7],  -MAD_F(0x09bd7ca0));
  1721.   MAD_F_MLA(hi, lo, X[10],  MAD_F(0x0cb19346));
  1722.   MAD_F_MLA(hi, lo, X[16],  MAD_F(0x0fdcf549));
  1724.   t3 = MAD_F_MLZ(hi, lo) + t7;
  1726.   MAD_F_ML0(hi, lo, X[0],   MAD_F(0x00b2aa3e));
  1727.   MAD_F_MLA(hi, lo, X[2],   MAD_F(0x03768962));
  1728.   MAD_F_MLA(hi, lo, X[3],  -MAD_F(0x04cfb0e2));
  1729.   MAD_F_MLA(hi, lo, X[5],  -MAD_F(0x07635284));
  1730.   MAD_F_MLA(hi, lo, X[6],   MAD_F(0x0898c779));
  1731.   MAD_F_MLA(hi, lo, X[8],   MAD_F(0x0acf37ad));
  1732.   MAD_F_MLA(hi, lo, X[9],  -MAD_F(0x0bcbe352));
  1733.   MAD_F_MLA(hi, lo, X[11], -MAD_F(0x0d7e8807));
  1734.   MAD_F_MLA(hi, lo, X[12],  MAD_F(0x0e313245));
  1735.   MAD_F_MLA(hi, lo, X[14],  MAD_F(0x0f426cb5));
  1736.   MAD_F_MLA(hi, lo, X[15], -MAD_F(0x0f9ee890));
  1737.   MAD_F_MLA(hi, lo, X[17], -MAD_F(0x0ffc19fd));
  1739.   x[8] = MAD_F_MLZ(hi, lo) + t3;
  1740.   x[9] = -x[8];
  1742.   MAD_F_ML0(hi, lo, X[0],  -MAD_F(0x0e313245));
  1743.   MAD_F_MLA(hi, lo, X[2],   MAD_F(0x0bcbe352));
  1744.   MAD_F_MLA(hi, lo, X[3],   MAD_F(0x0f9ee890));
  1745.   MAD_F_MLA(hi, lo, X[5],  -MAD_F(0x0898c779));
  1746.   MAD_F_MLA(hi, lo, X[6],  -MAD_F(0x0ffc19fd));
  1747.   MAD_F_MLA(hi, lo, X[8],   MAD_F(0x04cfb0e2));
  1748.   MAD_F_MLA(hi, lo, X[9],   MAD_F(0x0f426cb5));
  1749.   MAD_F_MLA(hi, lo, X[11], -MAD_F(0x00b2aa3e));
  1750.   MAD_F_MLA(hi, lo, X[12], -MAD_F(0x0d7e8807));
  1751.   MAD_F_MLA(hi, lo, X[14], -MAD_F(0x03768962));
  1752.   MAD_F_MLA(hi, lo, X[15],  MAD_F(0x0acf37ad));
  1753.   MAD_F_MLA(hi, lo, X[17],  MAD_F(0x07635284));
  1755.   x[21] = x[32] = MAD_F_MLZ(hi, lo) + t3;
  1757.   MAD_F_ML0(hi, lo, X[0],  -MAD_F(0x0d7e8807));
  1758.   MAD_F_MLA(hi, lo, X[2],   MAD_F(0x0f426cb5));
  1759.   MAD_F_MLA(hi, lo, X[3],   MAD_F(0x0acf37ad));
  1760.   MAD_F_MLA(hi, lo, X[5],  -MAD_F(0x0ffc19fd));
  1761.   MAD_F_MLA(hi, lo, X[6],  -MAD_F(0x07635284));
  1762.   MAD_F_MLA(hi, lo, X[8],   MAD_F(0x0f9ee890));
  1763.   MAD_F_MLA(hi, lo, X[9],   MAD_F(0x03768962));
  1764.   MAD_F_MLA(hi, lo, X[11], -MAD_F(0x0e313245));
  1765.   MAD_F_MLA(hi, lo, X[12],  MAD_F(0x00b2aa3e));
  1766.   MAD_F_MLA(hi, lo, X[14],  MAD_F(0x0bcbe352));
  1767.   MAD_F_MLA(hi, lo, X[15], -MAD_F(0x04cfb0e2));
  1768.   MAD_F_MLA(hi, lo, X[17], -MAD_F(0x0898c779));
  1770.   x[20] = x[33] = MAD_F_MLZ(hi, lo) - t3;
  1772.   MAD_F_ML0(hi, lo, t14, -MAD_F(0x0ec835e8));
  1773.   MAD_F_MLA(hi, lo, t15,  MAD_F(0x061f78aa));
  1775.   t4 = MAD_F_MLZ(hi, lo) - t7;
  1777.   MAD_F_ML0(hi, lo, t12, MAD_F(0x061f78aa));
  1778.   MAD_F_MLA(hi, lo, t13, MAD_F(0x0ec835e8));
  1780.   x[4]  = MAD_F_MLZ(hi, lo) + t4;
  1781.   x[13] = -x[4];
  1783.   MAD_F_ML0(hi, lo, t8,   MAD_F(0x09bd7ca0));
  1784.   MAD_F_MLA(hi, lo, t9,  -MAD_F(0x0216a2a2));
  1785.   MAD_F_MLA(hi, lo, t10,  MAD_F(0x0fdcf549));
  1786.   MAD_F_MLA(hi, lo, t11, -MAD_F(0x0cb19346));
  1788.   x[1]  = MAD_F_MLZ(hi, lo) + t4;
  1789.   x[16] = -x[1];
  1791.   MAD_F_ML0(hi, lo, t8,  -MAD_F(0x0fdcf549));
  1792.   MAD_F_MLA(hi, lo, t9,  -MAD_F(0x0cb19346));
  1793.   MAD_F_MLA(hi, lo, t10, -MAD_F(0x09bd7ca0));
  1794.   MAD_F_MLA(hi, lo, t11, -MAD_F(0x0216a2a2));
  1796.   x[25] = x[28] = MAD_F_MLZ(hi, lo) + t4;
  1798.   MAD_F_ML0(hi, lo, X[1],  -MAD_F(0x0fdcf549));
  1799.   MAD_F_MLA(hi, lo, X[7],  -MAD_F(0x0cb19346));
  1800.   MAD_F_MLA(hi, lo, X[10], -MAD_F(0x09bd7ca0));
  1801.   MAD_F_MLA(hi, lo, X[16], -MAD_F(0x0216a2a2));
  1803.   t5 = MAD_F_MLZ(hi, lo) - t6;
  1805.   MAD_F_ML0(hi, lo, X[0],   MAD_F(0x0898c779));
  1806.   MAD_F_MLA(hi, lo, X[2],   MAD_F(0x04cfb0e2));
  1807.   MAD_F_MLA(hi, lo, X[3],   MAD_F(0x0bcbe352));
  1808.   MAD_F_MLA(hi, lo, X[5],   MAD_F(0x00b2aa3e));
  1809.   MAD_F_MLA(hi, lo, X[6],   MAD_F(0x0e313245));
  1810.   MAD_F_MLA(hi, lo, X[8],  -MAD_F(0x03768962));
  1811.   MAD_F_MLA(hi, lo, X[9],   MAD_F(0x0f9ee890));
  1812.   MAD_F_MLA(hi, lo, X[11], -MAD_F(0x07635284));
  1813.   MAD_F_MLA(hi, lo, X[12],  MAD_F(0x0ffc19fd));
  1814.   MAD_F_MLA(hi, lo, X[14], -MAD_F(0x0acf37ad));
  1815.   MAD_F_MLA(hi, lo, X[15],  MAD_F(0x0f426cb5));
  1816.   MAD_F_MLA(hi, lo, X[17], -MAD_F(0x0d7e8807));
  1818.   x[2]  = MAD_F_MLZ(hi, lo) + t5;
  1819.   x[15] = -x[2];
  1821.   MAD_F_ML0(hi, lo, X[0],   MAD_F(0x07635284));
  1822.   MAD_F_MLA(hi, lo, X[2],   MAD_F(0x0acf37ad));
  1823.   MAD_F_MLA(hi, lo, X[3],   MAD_F(0x03768962));
  1824.   MAD_F_MLA(hi, lo, X[5],   MAD_F(0x0d7e8807));
  1825.   MAD_F_MLA(hi, lo, X[6],  -MAD_F(0x00b2aa3e));
  1826.   MAD_F_MLA(hi, lo, X[8],   MAD_F(0x0f426cb5));
  1827.   MAD_F_MLA(hi, lo, X[9],  -MAD_F(0x04cfb0e2));
  1828.   MAD_F_MLA(hi, lo, X[11],  MAD_F(0x0ffc19fd));
  1829.   MAD_F_MLA(hi, lo, X[12], -MAD_F(0x0898c779));
  1830.   MAD_F_MLA(hi, lo, X[14],  MAD_F(0x0f9ee890));
  1831.   MAD_F_MLA(hi, lo, X[15], -MAD_F(0x0bcbe352));
  1832.   MAD_F_MLA(hi, lo, X[17],  MAD_F(0x0e313245));
  1834.   x[3]  = MAD_F_MLZ(hi, lo) + t5;
  1835.   x[14] = -x[3];
  1837.   MAD_F_ML0(hi, lo, X[0],  -MAD_F(0x0ffc19fd));
  1838.   MAD_F_MLA(hi, lo, X[2],  -MAD_F(0x0f9ee890));
  1839.   MAD_F_MLA(hi, lo, X[3],  -MAD_F(0x0f426cb5));
  1840.   MAD_F_MLA(hi, lo, X[5],  -MAD_F(0x0e313245));
  1841.   MAD_F_MLA(hi, lo, X[6],  -MAD_F(0x0d7e8807));
  1842.   MAD_F_MLA(hi, lo, X[8],  -MAD_F(0x0bcbe352));
  1843.   MAD_F_MLA(hi, lo, X[9],  -MAD_F(0x0acf37ad));
  1844.   MAD_F_MLA(hi, lo, X[11], -MAD_F(0x0898c779));
  1845.   MAD_F_MLA(hi, lo, X[12], -MAD_F(0x07635284));
  1846.   MAD_F_MLA(hi, lo, X[14], -MAD_F(0x04cfb0e2));
  1847.   MAD_F_MLA(hi, lo, X[15], -MAD_F(0x03768962));
  1848.   MAD_F_MLA(hi, lo, X[17], -MAD_F(0x00b2aa3e));
  1850.   x[26] = x[27] = MAD_F_MLZ(hi, lo) + t5;
  1851. }
  1853. /*
  1854.  * NAME: III_imdct_l()
  1855.  * DESCRIPTION: perform IMDCT and windowing for long blocks
  1856.  */
  1857. static
  1858. void III_imdct_l(mad_fixed_t const X[18], mad_fixed_t z[36],
  1859.  unsigned int block_type)
  1860. {
  1861.   unsigned int i;
  1863.   /* IMDCT */
  1865.   imdct36(X, z);
  1867.   /* windowing */
  1869.   switch (block_type) {
  1870.   case 0:  /* normal window */
  1871. # if defined(ASO_INTERLEAVE1)
  1872.     {
  1873.       register mad_fixed_t tmp1, tmp2;
  1875.       tmp1 = window_l[0];
  1876.       tmp2 = window_l[1];
  1878.       for (i = 0; i < 34; i += 2) {
  1879. z[i + 0] = mad_f_mul(z[i + 0], tmp1);
  1880. tmp1 = window_l[i + 2];
  1881. z[i + 1] = mad_f_mul(z[i + 1], tmp2);
  1882. tmp2 = window_l[i + 3];
  1883.       }
  1885.       z[34] = mad_f_mul(z[34], tmp1);
  1886.       z[35] = mad_f_mul(z[35], tmp2);
  1887.     }
  1888. # elif defined(ASO_INTERLEAVE2)
  1889.     {
  1890.       register mad_fixed_t tmp1, tmp2;
  1892.       tmp1 = z[0];
  1893.       tmp2 = window_l[0];
  1895.       for (i = 0; i < 35; ++i) {
  1896. z[i] = mad_f_mul(tmp1, tmp2);
  1897. tmp1 = z[i + 1];
  1898. tmp2 = window_l[i + 1];
  1899.       }
  1901.       z[35] = mad_f_mul(tmp1, tmp2);
  1902.     }
  1903. # elif 1
  1904.     for (i = 0; i < 36; i += 4) {
  1905.       z[i + 0] = mad_f_mul(z[i + 0], window_l[i + 0]);
  1906.       z[i + 1] = mad_f_mul(z[i + 1], window_l[i + 1]);
  1907.       z[i + 2] = mad_f_mul(z[i + 2], window_l[i + 2]);
  1908.       z[i + 3] = mad_f_mul(z[i + 3], window_l[i + 3]);
  1909.     }
  1910. # else
  1911.     for (i =  0; i < 36; ++i) z[i] = mad_f_mul(z[i], window_l[i]);
  1912. # endif
  1913.     break;
  1915.   case 1:  /* start block */
  1916.     for (i =  0; i < 18; ++i) z[i] = mad_f_mul(z[i], window_l[i]);
  1917.     /*  (i = 18; i < 24; ++i) z[i] unchanged */
  1918.     for (i = 24; i < 30; ++i) z[i] = mad_f_mul(z[i], window_s[i - 18]);
  1919.     for (i = 30; i < 36; ++i) z[i] = 0;
  1920.     break;
  1922.   case 3:  /* stop block */
  1923.     for (i =  0; i <  6; ++i) z[i] = 0;
  1924.     for (i =  6; i < 12; ++i) z[i] = mad_f_mul(z[i], window_s[i - 6]);
  1925.     /*  (i = 12; i < 18; ++i) z[i] unchanged */
  1926.     for (i = 18; i < 36; ++i) z[i] = mad_f_mul(z[i], window_l[i]);
  1927.     break;
  1928.   }
  1929. }
  1930. # endif  /* ASO_IMDCT */
  1932. /*
  1933.  * NAME: III_imdct_s()
  1934.  * DESCRIPTION: perform IMDCT and windowing for short blocks
  1935.  */
  1936. static
  1937. void III_imdct_s(mad_fixed_t const X[18], mad_fixed_t z[36])
  1938. {
  1939.   mad_fixed_t y[36], *yptr;
  1940.   mad_fixed_t const *wptr;
  1941.   int w, i;
  1942.   register mad_fixed64hi_t hi;
  1943.   register mad_fixed64lo_t lo;
  1945.   /* IMDCT */
  1947.   yptr = &y[0];
  1949.   for (w = 0; w < 3; ++w) {
  1950.     register mad_fixed_t const (*s)[6];
  1952.     s = imdct_s;
  1954.     for (i = 0; i < 3; ++i) {
  1955.       MAD_F_ML0(hi, lo, X[0], (*s)[0]);
  1956.       MAD_F_MLA(hi, lo, X[1], (*s)[1]);
  1957.       MAD_F_MLA(hi, lo, X[2], (*s)[2]);
  1958.       MAD_F_MLA(hi, lo, X[3], (*s)[3]);
  1959.       MAD_F_MLA(hi, lo, X[4], (*s)[4]);
  1960.       MAD_F_MLA(hi, lo, X[5], (*s)[5]);
  1962.       yptr[i + 0] = MAD_F_MLZ(hi, lo);
  1963.       yptr[5 - i] = -yptr[i + 0];
  1965.       ++s;
  1967.       MAD_F_ML0(hi, lo, X[0], (*s)[0]);
  1968.       MAD_F_MLA(hi, lo, X[1], (*s)[1]);
  1969.       MAD_F_MLA(hi, lo, X[2], (*s)[2]);
  1970.       MAD_F_MLA(hi, lo, X[3], (*s)[3]);
  1971.       MAD_F_MLA(hi, lo, X[4], (*s)[4]);
  1972.       MAD_F_MLA(hi, lo, X[5], (*s)[5]);
  1974.       yptr[ i + 6] = MAD_F_MLZ(hi, lo);
  1975.       yptr[11 - i] = yptr[i + 6];
  1977.       ++s;
  1978.     }
  1980.     yptr += 12;
  1981.     X    += 6;
  1982.   }
  1984.   /* windowing, overlapping and concatenation */
  1986.   yptr = &y[0];
  1987.   wptr = &window_s[0];
  1989.   for (i = 0; i < 6; ++i) {
  1990.     z[i +  0] = 0;
  1991.     z[i +  6] = mad_f_mul(yptr[ 0 + 0], wptr[0]);
  1993.     MAD_F_ML0(hi, lo, yptr[ 0 + 6], wptr[6]);
  1994.     MAD_F_MLA(hi, lo, yptr[12 + 0], wptr[0]);
  1996.     z[i + 12] = MAD_F_MLZ(hi, lo);
  1998.     MAD_F_ML0(hi, lo, yptr[12 + 6], wptr[6]);
  1999.     MAD_F_MLA(hi, lo, yptr[24 + 0], wptr[0]);
  2001.     z[i + 18] = MAD_F_MLZ(hi, lo);
  2003.     z[i + 24] = mad_f_mul(yptr[24 + 6], wptr[6]);
  2004.     z[i + 30] = 0;
  2006.     ++yptr;
  2007.     ++wptr;
  2008.   }
  2009. }
  2011. /*
  2012.  * NAME: III_overlap()
  2013.  * DESCRIPTION: perform overlap-add of windowed IMDCT outputs
  2014.  */
  2015. static
  2016. void III_overlap(mad_fixed_t const output[36], mad_fixed_t overlap[18],
  2017.  mad_fixed_t sample[18][32], unsigned int sb)
  2018. {
  2019.   unsigned int i;
  2021. # if defined(ASO_INTERLEAVE2)
  2022.   {
  2023.     register mad_fixed_t tmp1, tmp2;
  2025.     tmp1 = overlap[0];
  2026.     tmp2 = overlap[1];
  2028.     for (i = 0; i < 16; i += 2) {
  2029.       sample[i + 0][sb] = output[i + 0] + tmp1;
  2030.       overlap[i + 0]    = output[i + 0 + 18];
  2031.       tmp1 = overlap[i + 2];
  2033.       sample[i + 1][sb] = output[i + 1] + tmp2;
  2034.       overlap[i + 1]    = output[i + 1 + 18];
  2035.       tmp2 = overlap[i + 3];
  2036.     }
  2038.     sample[16][sb] = output[16] + tmp1;
  2039.     overlap[16]    = output[16 + 18];
  2040.     sample[17][sb] = output[17] + tmp2;
  2041.     overlap[17]    = output[17 + 18];
  2042.   }
  2043. # elif 0
  2044.   for (i = 0; i < 18; i += 2) {
  2045.     sample[i + 0][sb] = output[i + 0] + overlap[i + 0];
  2046.     overlap[i + 0]    = output[i + 0 + 18];
  2048.     sample[i + 1][sb] = output[i + 1] + overlap[i + 1];
  2049.     overlap[i + 1]    = output[i + 1 + 18];
  2050.   }
  2051. # else
  2052.   for (i = 0; i < 18; ++i) {
  2053.     sample[i][sb] = output[i] + overlap[i];
  2054.     overlap[i]    = output[i + 18];
  2055.   }
  2056. # endif
  2057. }
  2059. /*
  2060.  * NAME: III_overlap_z()
  2061.  * DESCRIPTION: perform "overlap-add" of zero IMDCT outputs
  2062.  */
  2063. static inline
  2064. void III_overlap_z(mad_fixed_t overlap[18],
  2065.    mad_fixed_t sample[18][32], unsigned int sb)
  2066. {
  2067.   unsigned int i;
  2069. # if defined(ASO_INTERLEAVE2)
  2070.   {
  2071.     register mad_fixed_t tmp1, tmp2;
  2073.     tmp1 = overlap[0];
  2074.     tmp2 = overlap[1];
  2076.     for (i = 0; i < 16; i += 2) {
  2077.       sample[i + 0][sb] = tmp1;
  2078.       overlap[i + 0]    = 0;
  2079.       tmp1 = overlap[i + 2];
  2081.       sample[i + 1][sb] = tmp2;
  2082.       overlap[i + 1]    = 0;
  2083.       tmp2 = overlap[i + 3];
  2084.     }
  2086.     sample[16][sb] = tmp1;
  2087.     overlap[16]    = 0;
  2088.     sample[17][sb] = tmp2;
  2089.     overlap[17]    = 0;
  2090.   }
  2091. # else
  2092.   for (i = 0; i < 18; ++i) {
  2093.     sample[i][sb] = overlap[i];
  2094.     overlap[i]    = 0;
  2095.   }
  2096. # endif
  2097. }
  2099. /*
  2100.  * NAME: III_freqinver()
  2101.  * DESCRIPTION: perform subband frequency inversion for odd sample lines
  2102.  */
  2103. static
  2104. void III_freqinver(mad_fixed_t sample[18][32], unsigned int sb)
  2105. {
  2106.   unsigned int i;
  2108. # if 1 || defined(ASO_INTERLEAVE1) || defined(ASO_INTERLEAVE2)
  2109.   {
  2110.     register mad_fixed_t tmp1, tmp2;
  2112.     tmp1 = sample[1][sb];
  2113.     tmp2 = sample[3][sb];
  2115.     for (i = 1; i < 13; i += 4) {
  2116.       sample[i + 0][sb] = -tmp1;
  2117.       tmp1 = sample[i + 4][sb];
  2118.       sample[i + 2][sb] = -tmp2;
  2119.       tmp2 = sample[i + 6][sb];
  2120.     }
  2122.     sample[13][sb] = -tmp1;
  2123.     tmp1 = sample[17][sb];
  2124.     sample[15][sb] = -tmp2;
  2125.     sample[17][sb] = -tmp1;
  2126.   }
  2127. # else
  2128.   for (i = 1; i < 18; i += 2)
  2129.     sample[i][sb] = -sample[i][sb];
  2130. # endif
  2131. }
  2133. /*
  2134.  * NAME: III_decode()
  2135.  * DESCRIPTION: decode frame main_data
  2136.  */
  2137. static
  2138. enum mad_error III_decode(struct mad_bitptr *ptr, struct mad_frame *frame,
  2139.   struct sideinfo *si, unsigned int nch)
  2140. {
  2141.   struct mad_header *header = &frame->header;
  2142.   unsigned int sfreqi, ngr, gr;
  2144.   {
  2145.     unsigned int sfreq;
  2147.     sfreq = header->samplerate;
  2148.     if (header->flags & MAD_FLAG_MPEG_2_5_EXT)
  2149.       sfreq *= 2;
  2151.     /* 48000 => 0, 44100 => 1, 32000 => 2,
  2152.        24000 => 3, 22050 => 4, 16000 => 5 */
  2153.     sfreqi = ((sfreq >>  7) & 0x000f) +
  2154.              ((sfreq >> 15) & 0x0001) - 8;
  2156.     if (header->flags & MAD_FLAG_MPEG_2_5_EXT)
  2157.       sfreqi += 3;
  2158.   }
  2160.   /* scalefactors, Huffman decoding, requantization */
  2162.   ngr = (header->flags & MAD_FLAG_LSF_EXT) ? 1 : 2;
  2164.   for (gr = 0; gr < ngr; ++gr) {
  2165.     struct granule *granule = &si->gr[gr];
  2166.     unsigned char const *sfbwidth[2];
  2167.     mad_fixed_t xr[2][576];
  2168.     unsigned int ch;
  2169.     enum mad_error error;
  2171.     for (ch = 0; ch < nch; ++ch) {
  2172.       struct channel *channel = &granule->ch[ch];
  2173.       unsigned int part2_length;
  2175.       sfbwidth[ch] = sfbwidth_table[sfreqi].l;
  2176.       if (channel->block_type == 2) {
  2177. sfbwidth[ch] = (channel->flags & mixed_block_flag) ?
  2178.   sfbwidth_table[sfreqi].m : sfbwidth_table[sfreqi].s;
  2179.       }
  2181.       if (header->flags & MAD_FLAG_LSF_EXT) {
  2182. part2_length = III_scalefactors_lsf(ptr, channel,
  2183.     ch == 0 ? 0 : &si->gr[1].ch[1],
  2184.     header->mode_extension);
  2185.       }
  2186.       else {
  2187. part2_length = III_scalefactors(ptr, channel, &si->gr[0].ch[ch],
  2188. gr == 0 ? 0 : si->scfsi[ch]);
  2189.       }
  2191.       error = III_huffdecode(ptr, xr[ch], channel, sfbwidth[ch], part2_length);
  2192.       if (error)
  2193. return error;
  2194.     }
  2196.     /* joint stereo processing */
  2198.     if (header->mode == MAD_MODE_JOINT_STEREO && header->mode_extension) {
  2199.       error = III_stereo(xr, granule, header, sfbwidth[0]);
  2200.       if (error)
  2201. return error;
  2202.     }
  2204.     /* reordering, alias reduction, IMDCT, overlap-add, frequency inversion */
  2206.     for (ch = 0; ch < nch; ++ch) {
  2207.       struct channel const *channel = &granule->ch[ch];
  2208.       mad_fixed_t (*sample)[32] = &frame->sbsample[ch][18 * gr];
  2209.       unsigned int sb, l, i, sblimit;
  2210.       mad_fixed_t output[36];
  2212.       if (channel->block_type == 2) {
  2213. III_reorder(xr[ch], channel, sfbwidth[ch]);
  2215. # if !defined(OPT_STRICT)
  2216. /*
  2217.  * According to ISO/IEC 11172-3, "Alias reduction is not applied for
  2218.  * granules with block_type == 2 (short block)." However, other
  2219.  * sources suggest alias reduction should indeed be performed on the
  2220.  * lower two subbands of mixed blocks. Most other implementations do
  2221.  * this, so by default we will too.
  2222.  */
  2223. if (channel->flags & mixed_block_flag)
  2224.   III_aliasreduce(xr[ch], 36);
  2225. # endif
  2226.       }
  2227.       else
  2228. III_aliasreduce(xr[ch], 576);
  2230.       l = 0;
  2232.       /* subbands 0-1 */
  2234.       if (channel->block_type != 2 || (channel->flags & mixed_block_flag)) {
  2235. unsigned int block_type;
  2237. block_type = channel->block_type;
  2238. if (channel->flags & mixed_block_flag)
  2239.   block_type = 0;
  2241. /* long blocks */
  2242. for (sb = 0; sb < 2; ++sb, l += 18) {
  2243.   III_imdct_l(&xr[ch][l], output, block_type);
  2244.   III_overlap(output, (*frame->overlap)[ch][sb], sample, sb);
  2245. }
  2246.       }
  2247.       else {
  2248. /* short blocks */
  2249. for (sb = 0; sb < 2; ++sb, l += 18) {
  2250.   III_imdct_s(&xr[ch][l], output);
  2251.   III_overlap(output, (*frame->overlap)[ch][sb], sample, sb);
  2252. }
  2253.       }
  2255.       III_freqinver(sample, 1);
  2257.       /* (nonzero) subbands 2-31 */
  2259.       i = 576;
  2260.       while (i > 36 && xr[ch][i - 1] == 0)
  2261. --i;
  2263.       sblimit = 32 - (576 - i) / 18;
  2265.       if (channel->block_type != 2) {
  2266. /* long blocks */
  2267. for (sb = 2; sb < sblimit; ++sb, l += 18) {
  2268.   III_imdct_l(&xr[ch][l], output, channel->block_type);
  2269.   III_overlap(output, (*frame->overlap)[ch][sb], sample, sb);
  2271.   if (sb & 1)
  2272.     III_freqinver(sample, sb);
  2273. }
  2274.       }
  2275.       else {
  2276. /* short blocks */
  2277. for (sb = 2; sb < sblimit; ++sb, l += 18) {
  2278.   III_imdct_s(&xr[ch][l], output);
  2279.   III_overlap(output, (*frame->overlap)[ch][sb], sample, sb);
  2281.   if (sb & 1)
  2282.     III_freqinver(sample, sb);
  2283. }
  2284.       }
  2286.       /* remaining (zero) subbands */
  2288.       for (sb = sblimit; sb < 32; ++sb) {
  2289. III_overlap_z((*frame->overlap)[ch][sb], sample, sb);
  2291. if (sb & 1)
  2292.   III_freqinver(sample, sb);
  2293.       }
  2294.     }
  2295.   }
  2297.   return MAD_ERROR_NONE;
  2298. }
  2300. /*
  2301.  * NAME: layer->III()
  2302.  * DESCRIPTION: decode a single Layer III frame
  2303.  */
  2304. int mad_layer_III(struct mad_stream *stream, struct mad_frame *frame)
  2305. {
  2306.   struct mad_header *header = &frame->header;
  2307.   unsigned int nch, priv_bitlen, next_md_begin = 0;
  2308.   unsigned int si_len, data_bitlen, md_len;
  2309.   unsigned int frame_space, frame_used, frame_free;
  2310.   struct mad_bitptr ptr;
  2311.   struct sideinfo si;
  2312.   enum mad_error error;
  2313.   int result = 0;
  2315.   /* allocate Layer III dynamic structures */
  2317.   if (stream->main_data == 0) {
  2318.     stream->main_data = malloc(MAD_BUFFER_MDLEN);
  2319.     if (stream->main_data == 0) {
  2320.       stream->error = MAD_ERROR_NOMEM;
  2321.       return -1;
  2322.     }
  2323.   }
  2325.   if (frame->overlap == 0) {
  2326.     frame->overlap = calloc(2 * 32 * 18, sizeof(mad_fixed_t));
  2327.     if (frame->overlap == 0) {
  2328.       stream->error = MAD_ERROR_NOMEM;
  2329.       return -1;
  2330.     }
  2331.   }
  2333.   nch = MAD_NCHANNELS(header);
  2334.   si_len = (header->flags & MAD_FLAG_LSF_EXT) ?
  2335.     (nch == 1 ? 9 : 17) : (nch == 1 ? 17 : 32);
  2337.   /* check frame sanity */
  2339.   if (stream->next_frame - mad_bit_nextbyte(&stream->ptr) <
  2340.       (signed int) si_len) {
  2341.     stream->error = MAD_ERROR_BADFRAMELEN;
  2342.     stream->md_len = 0;
  2343.     return -1;
  2344.   }
  2346.   /* check CRC word */
  2348.   if (header->flags & MAD_FLAG_PROTECTION) {
  2349.     header->crc_check =
  2350.       mad_bit_crc(stream->ptr, si_len * CHAR_BIT, header->crc_check);
  2352.     if (header->crc_check != header->crc_target &&
  2353. !(frame->options & MAD_OPTION_IGNORECRC)) {
  2354.       stream->error = MAD_ERROR_BADCRC;
  2355.       result = -1;
  2356.     }
  2357.   }
  2359.   /* decode frame side information */
  2361.   error = III_sideinfo(&stream->ptr, nch, header->flags & MAD_FLAG_LSF_EXT,
  2362.        &si, &data_bitlen, &priv_bitlen);
  2363.   if (error && result == 0) {
  2364.     stream->error = error;
  2365.     result = -1;
  2366.   }
  2368.   header->flags        |= priv_bitlen;
  2369.   header->private_bits |= si.private_bits;
  2371.   /* find main_data of next frame */
  2373.   {
  2374.     struct mad_bitptr peek;
  2375.     unsigned long header;
  2377.     mad_bit_init(&peek, stream->next_frame);
  2379.     header = mad_bit_read(&peek, 32);
  2380.     if ((header & 0xffe60000L) /* syncword | layer */ == 0xffe20000L) {
  2381.       if (!(header & 0x00010000L))  /* protection_bit */
  2382. mad_bit_skip(&peek, 16);  /* crc_check */
  2384.       next_md_begin =
  2385. mad_bit_read(&peek, (header & 0x00080000L) /* ID */ ? 9 : 8);
  2386.     }
  2388.     mad_bit_finish(&peek);
  2389.   }
  2391.   /* find main_data of this frame */
  2393.   frame_space = stream->next_frame - mad_bit_nextbyte(&stream->ptr);
  2395.   if (next_md_begin > si.main_data_begin + frame_space)
  2396.     next_md_begin = 0;
  2398.   md_len = si.main_data_begin + frame_space - next_md_begin;
  2400.   frame_used = 0;
  2402.   if (si.main_data_begin == 0) {
  2403.     ptr = stream->ptr;
  2404.     stream->md_len = 0;
  2406.     frame_used = md_len;
  2407.   }
  2408.   else {
  2409.     if (si.main_data_begin > stream->md_len) {
  2410.       if (result == 0) {
  2411. stream->error = MAD_ERROR_BADDATAPTR;
  2412. result = -1;
  2413.       }
  2414.     }
  2415.     else {
  2416.       mad_bit_init(&ptr,
  2417.    *stream->main_data + stream->md_len - si.main_data_begin);
  2419.       if (md_len > si.main_data_begin) {
  2420. assert(stream->md_len + md_len -
  2421.        si.main_data_begin <= MAD_BUFFER_MDLEN);
  2423. memcpy(*stream->main_data + stream->md_len,
  2424.        mad_bit_nextbyte(&stream->ptr),
  2425.        frame_used = md_len - si.main_data_begin);
  2426. stream->md_len += frame_used;
  2427.       }
  2428.     }
  2429.   }
  2431.   frame_free = frame_space - frame_used;
  2433.   /* decode main_data */
  2435.   if (result == 0) {
  2436.     error = III_decode(&ptr, frame, &si, nch);
  2437.     if (error) {
  2438.       stream->error = error;
  2439.       result = -1;
  2440.     }
  2442.     /* designate ancillary bits */
  2444.     stream->anc_ptr    = ptr;
  2445.     stream->anc_bitlen = md_len * CHAR_BIT - data_bitlen;
  2446.   }
  2448. # if 0 && defined(DEBUG)
  2449.   fprintf(stderr,
  2450.   "main_data_begin:%u, md_len:%u, frame_free:%u, "
  2451.   "data_bitlen:%u, anc_bitlen: %un",
  2452.   si.main_data_begin, md_len, frame_free,
  2453.   data_bitlen, stream->anc_bitlen);
  2454. # endif
  2456.   /* preload main_data buffer with up to 511 bytes for next frame(s) */
  2458.   if (frame_free >= next_md_begin) {
  2459.     memcpy(*stream->main_data,
  2460.    stream->next_frame - next_md_begin, next_md_begin);
  2461.     stream->md_len = next_md_begin;
  2462.   }
  2463.   else {
  2464.     if (md_len < si.main_data_begin) {
  2465.       unsigned int extra;
  2467.       extra = si.main_data_begin - md_len;
  2468.       if (extra + frame_free > next_md_begin)
  2469. extra = next_md_begin - frame_free;
  2471.       if (extra < stream->md_len) {
  2472. memmove(*stream->main_data,
  2473. *stream->main_data + stream->md_len - extra, extra);
  2474. stream->md_len = extra;
  2475.       }
  2476.     }
  2477.     else
  2478.       stream->md_len = 0;
  2480.     memcpy(*stream->main_data + stream->md_len,
  2481.    stream->next_frame - frame_free, frame_free);
  2482.     stream->md_len += frame_free;
  2483.   }
  2485.   return result;
  2486. }