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

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

多媒体编程

开发平台：

Visual C++

layer12.c：源码内容
  1. /*
  2.  * libmad - MPEG audio decoder library
  4.  */
  6. # ifdef HAVE_CONFIG_H
  7. #  include "config.h"
  8. # endif
  10. # include "global.h"
  12. # ifdef HAVE_LIMITS_H
  13. #  include <limits.h>
  14. # else
  15. #  define CHAR_BIT  8
  16. # endif
  18. # include "fixed.h"
  19. # include "bit.h"
  20. # include "stream.h"
  21. # include "frame.h"
  22. # include "layer12.h"
  25. static
  26. mad_fixed_t const sf_table[63] = {
  27. # include "sf_table.dat"
  28. };
  30. /* --- Layer I ------------------------------------------------------------- */
  32. static
  33. mad_fixed_t const linear_table[14] = {
  34.   MAD_F(0x15555555),  /* 2^2  / (2^2  - 1) == 1.33333333333333 */
  35.   MAD_F(0x12492492),  /* 2^3  / (2^3  - 1) == 1.14285714285714 */
  36.   MAD_F(0x11111111),  /* 2^4  / (2^4  - 1) == 1.06666666666667 */
  37.   MAD_F(0x10842108),  /* 2^5  / (2^5  - 1) == 1.03225806451613 */
  38.   MAD_F(0x10410410),  /* 2^6  / (2^6  - 1) == 1.01587301587302 */
  39.   MAD_F(0x10204081),  /* 2^7  / (2^7  - 1) == 1.00787401574803 */
  40.   MAD_F(0x10101010),  /* 2^8  / (2^8  - 1) == 1.00392156862745 */
  41.   MAD_F(0x10080402),  /* 2^9  / (2^9  - 1) == 1.00195694716243 */
  42.   MAD_F(0x10040100),  /* 2^10 / (2^10 - 1) == 1.00097751710655 */
  43.   MAD_F(0x10020040),  /* 2^11 / (2^11 - 1) == 1.00048851978505 */
  44.   MAD_F(0x10010010),  /* 2^12 / (2^12 - 1) == 1.00024420024420 */
  45.   MAD_F(0x10008004),  /* 2^13 / (2^13 - 1) == 1.00012208521548 */
  46.   MAD_F(0x10004001),  /* 2^14 / (2^14 - 1) == 1.00006103888177 */
  47.   MAD_F(0x10002000)   /* 2^15 / (2^15 - 1) == 1.00003051850948 */
  48. };
  51. static
  52. mad_fixed_t I_sample(struct mad_bitptr *ptr, unsigned int nb)
  53. {
  54.   mad_fixed_t sample;
  56.   sample = mad_bit_read(ptr, nb);
  59.   sample ^= 1 << (nb - 1);
  60.   sample |= -(sample & (1 << (nb - 1)));
  62.   sample <<= MAD_F_FRACBITS - (nb - 1);
  65.   /* s'' = (2^nb / (2^nb - 1)) * (s''' + 2^(-nb + 1)) */
  67.   sample += MAD_F_ONE >> (nb - 1);
  69.   return mad_f_mul(sample, linear_table[nb - 2]);
  71.   /* s' = factor * s'' */
  72. }
  75. int mad_layer_I(struct mad_stream *stream, struct mad_frame *frame)
  76. {
  77.   struct mad_header *header = &frame->header;
  78.   unsigned int nch, bound, ch, s, sb, nb;
  79.   unsigned char allocation[2][32], scalefactor[2][32];
  81.   nch = MAD_NCHANNELS(header);
  83.   bound = 32;
  84.   if (header->mode == MAD_MODE_JOINT_STEREO) {
  85.     header->flags |= MAD_FLAG_I_STEREO;
  86.     bound = 4 + header->mode_extension * 4;
  87.   }
  90.   if (header->flags & MAD_FLAG_PROTECTION) {
  91.     header->crc_check =
  92.       mad_bit_crc(stream->ptr, 4 * (bound * nch + (32 - bound)),
  93.   header->crc_check);
  95.     if (header->crc_check != header->crc_target &&
  96. !(frame->options & MAD_OPTION_IGNORECRC)) {
  97.       stream->error = MAD_ERROR_BADCRC;
  98.       return -1;
  99.     }
  100.   }
  103.   for (sb = 0; sb < bound; ++sb) {
  104.     for (ch = 0; ch < nch; ++ch) {
  105.       nb = mad_bit_read(&stream->ptr, 4);
  107.       if (nb == 15) {
  108. stream->error = MAD_ERROR_BADBITALLOC;
  109. return -1;
  110.       }
  112.       allocation[ch][sb] = nb ? nb + 1 : 0;
  113.     }
  114.   }
  116.   for (sb = bound; sb < 32; ++sb) {
  117.     nb = mad_bit_read(&stream->ptr, 4);
  119.     if (nb == 15) {
  120.       stream->error = MAD_ERROR_BADBITALLOC;
  121.       return -1;
  122.     }
  124.     allocation[0][sb] =
  125.     allocation[1][sb] = nb ? nb + 1 : 0;
  126.   }
  129.   for (sb = 0; sb < 32; ++sb) {
  130.     for (ch = 0; ch < nch; ++ch) {
  131.       if (allocation[ch][sb]) {
  132. scalefactor[ch][sb] = mad_bit_read(&stream->ptr, 6);
  134. if (scalefactor[ch][sb] == 63) {
  135.   stream->error = MAD_ERROR_BADSCALEFACTOR;
  136.   return -1;
  137. }
  138.       }
  139.     }
  140.   }
  143.   for (s = 0; s < 12; ++s) {
  144.     for (sb = 0; sb < bound; ++sb) {
  145.       for (ch = 0; ch < nch; ++ch) {
  146. nb = allocation[ch][sb];
  147. frame->sbsample[ch][s][sb] = nb ?
  148.   mad_f_mul(I_sample(&stream->ptr, nb),
  149.     sf_table[scalefactor[ch][sb]]) : 0;
  150.       }
  151.     }
  153.     for (sb = bound; sb < 32; ++sb) {
  154.       if ((nb = allocation[0][sb])) {
  155. mad_fixed_t sample;
  157. sample = I_sample(&stream->ptr, nb);
  159. for (ch = 0; ch < nch; ++ch) {
  160.   frame->sbsample[ch][s][sb] =
  161.     mad_f_mul(sample, sf_table[scalefactor[ch][sb]]);
  162. }
  163.       }
  164.       else {
  165. for (ch = 0; ch < nch; ++ch)
  166.   frame->sbsample[ch][s][sb] = 0;
  167.       }
  168.     }
  169.   }
  171.   return 0;
  172. }
  174. /* --- Layer II ------------------------------------------------------------ */
  176. static
  177. struct {
  178.   unsigned int sblimit;
  179.   unsigned char const offsets[30];
  180. } const sbquant_table[5] = {
  181.   /* ISO/IEC 11172-3 Table B.2a */
  182.   { 27, { 7, 7, 7, 6, 6, 6, 6, 6, 6, 6, 6, 3, 3, 3, 3, 3, /* 0 */
  183.   3, 3, 3, 3, 3, 3, 3, 0, 0, 0, 0 } },
  184.   /* ISO/IEC 11172-3 Table B.2b */
  185.   { 30, { 7, 7, 7, 6, 6, 6, 6, 6, 6, 6, 6, 3, 3, 3, 3, 3, /* 1 */
  186.   3, 3, 3, 3, 3, 3, 3, 0, 0, 0, 0, 0, 0, 0 } },
  187.   /* ISO/IEC 11172-3 Table B.2c */
  188.   {  8, { 5, 5, 2, 2, 2, 2, 2, 2 } }, /* 2 */
  189.   /* ISO/IEC 11172-3 Table B.2d */
  190.   { 12, { 5, 5, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2 } }, /* 3 */
  191.   /* ISO/IEC 13818-3 Table B.1 */
  192.   { 30, { 4, 4, 4, 4, 2, 2, 2, 2, 2, 2, 2, 1, 1, 1, 1, 1, /* 4 */
  193.   1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1 } }
  194. };
  196. static
  197. struct {
  198.   unsigned short nbal;
  199.   unsigned short offset;
  200. } const bitalloc_table[8] = {
  201.   { 2, 0 },  /* 0 */
  202.   { 2, 3 },  /* 1 */
  203.   { 3, 3 },  /* 2 */
  204.   { 3, 1 },  /* 3 */
  205.   { 4, 2 },  /* 4 */
  206.   { 4, 3 },  /* 5 */
  207.   { 4, 4 },  /* 6 */
  208.   { 4, 5 }   /* 7 */
  209. };
  211. static
  212. unsigned char const offset_table[6][15] = {
  213.   { 0, 1, 16                                             },  /* 0 */
  214.   { 0, 1,  2, 3, 4, 5, 16                                },  /* 1 */
  215.   { 0, 1,  2, 3, 4, 5,  6, 7,  8,  9, 10, 11, 12, 13, 14 },  /* 2 */
  216.   { 0, 1,  3, 4, 5, 6,  7, 8,  9, 10, 11, 12, 13, 14, 15 },  /* 3 */
  217.   { 0, 1,  2, 3, 4, 5,  6, 7,  8,  9, 10, 11, 12, 13, 16 },  /* 4 */
  218.   { 0, 2,  4, 5, 6, 7,  8, 9, 10, 11, 12, 13, 14, 15, 16 }   /* 5 */
  219. };
  221. static
  222. struct quantclass {
  223.   unsigned short nlevels;
  224.   unsigned char group;
  225.   unsigned char bits;
  226.   mad_fixed_t C;
  227.   mad_fixed_t D;
  228. } const qc_table[17] = {
  229. # include "qc_table.dat"
  230. };
  233. static
  234. void II_samples(struct mad_bitptr *ptr,
  235. struct quantclass const *quantclass,
  236. mad_fixed_t output[3])
  237. {
  238.   unsigned int nb, s, sample[3];
  240.   if ((nb = quantclass->group)) {
  241.     unsigned int c, nlevels;
  243.     c = mad_bit_read(ptr, quantclass->bits);
  244.     nlevels = quantclass->nlevels;
  246.     for (s = 0; s < 3; ++s) {
  247.       sample[s] = c % nlevels;
  248.       c /= nlevels;
  249.     }
  250.   }
  251.   else {
  252.     nb = quantclass->bits;
  254.     for (s = 0; s < 3; ++s)
  255.       sample[s] = mad_bit_read(ptr, nb);
  256.   }
  258.   for (s = 0; s < 3; ++s) {
  259.     mad_fixed_t requantized;
  262.     requantized  = sample[s] ^ (1 << (nb - 1));
  263.     requantized |= -(requantized & (1 << (nb - 1)));
  265.     requantized <<= MAD_F_FRACBITS - (nb - 1);
  268.     /* s'' = C * (s''' + D) */
  270.     output[s] = mad_f_mul(requantized + quantclass->D, quantclass->C);
  272.     /* s' = factor * s'' */
  273.   }
  274. }
  277. int mad_layer_II(struct mad_stream *stream, struct mad_frame *frame)
  278. {
  279.   struct mad_header *header = &frame->header;
  280.   struct mad_bitptr start;
  281.   unsigned int index, sblimit, nbal, nch, bound, gr, ch, s, sb;
  282.   unsigned char const *offsets;
  283.   unsigned char allocation[2][32], scfsi[2][32], scalefactor[2][32][3];
  284.   mad_fixed_t samples[3];
  286.   nch = MAD_NCHANNELS(header);
  288.   if (header->flags & MAD_FLAG_LSF_EXT)
  289.     index = 4;
  290.   else {
  291.     switch (nch == 2 ? header->bitrate / 2 : header->bitrate) {
  292.     case 32000:
  293.     case 48000:
  294.       index = (header->samplerate == 32000) ? 3 : 2;
  295.       break;
  297.     case 56000:
  298.     case 64000:
  299.     case 80000:
  300.       index = 0;
  301.       break;
  303.     default:
  304.       index = (header->samplerate == 48000) ? 0 : 1;
  305.     }
  306.   }
  308.   sblimit = sbquant_table[index].sblimit;
  309.   offsets = sbquant_table[index].offsets;
  311.   bound = 32;
  312.   if (header->mode == MAD_MODE_JOINT_STEREO) {
  313.     header->flags |= MAD_FLAG_I_STEREO;
  314.     bound = 4 + header->mode_extension * 4;
  315.   }
  317.   if (bound > sblimit)
  318.     bound = sblimit;
  320.   start = stream->ptr;
  323.   for (sb = 0; sb < bound; ++sb) {
  324.     nbal = bitalloc_table[offsets[sb]].nbal;
  326.     for (ch = 0; ch < nch; ++ch)
  327.       allocation[ch][sb] = mad_bit_read(&stream->ptr, nbal);
  328.   }
  330.   for (sb = bound; sb < sblimit; ++sb) {
  331.     nbal = bitalloc_table[offsets[sb]].nbal;
  333.     allocation[0][sb] =
  334.     allocation[1][sb] = mad_bit_read(&stream->ptr, nbal);
  335.   }
  338.   for (sb = 0; sb < sblimit; ++sb) {
  339.     for (ch = 0; ch < nch; ++ch) {
  340.       if (allocation[ch][sb])
  341. scfsi[ch][sb] = mad_bit_read(&stream->ptr, 2);
  342.     }
  343.   }
  346.   if (header->flags & MAD_FLAG_PROTECTION) {
  347.     header->crc_check =
  348.       mad_bit_crc(start, mad_bit_length(&start, &stream->ptr),
  349.   header->crc_check);
  351.     if (header->crc_check != header->crc_target &&
  352. !(frame->options & MAD_OPTION_IGNORECRC)) {
  353.       stream->error = MAD_ERROR_BADCRC;
  354.       return -1;
  355.     }
  356.   }
  359.   for (sb = 0; sb < sblimit; ++sb) {
  360.     for (ch = 0; ch < nch; ++ch) {
  361.       if (allocation[ch][sb]) {
  362. scalefactor[ch][sb][0] = mad_bit_read(&stream->ptr, 6);
  364. switch (scfsi[ch][sb]) {
  365. case 2:
  366.   scalefactor[ch][sb][2] =
  367.   scalefactor[ch][sb][1] =
  368.   scalefactor[ch][sb][0];
  369.   break;
  371. case 0:
  372.   scalefactor[ch][sb][1] = mad_bit_read(&stream->ptr, 6);
  374. case 1:
  375. case 3:
  376.   scalefactor[ch][sb][2] = mad_bit_read(&stream->ptr, 6);
  377. }
  379. if (scfsi[ch][sb] & 1)
  380.   scalefactor[ch][sb][1] = scalefactor[ch][sb][scfsi[ch][sb] - 1];
  382. if (scalefactor[ch][sb][0] == 63 ||
  383.     scalefactor[ch][sb][1] == 63 ||
  384.     scalefactor[ch][sb][2] == 63) {
  385.   stream->error = MAD_ERROR_BADSCALEFACTOR;
  386.   return -1;
  387. }
  388.       }
  389.     }
  390.   }
  393.   for (gr = 0; gr < 12; ++gr) {
  394.     for (sb = 0; sb < bound; ++sb) {
  395.       for (ch = 0; ch < nch; ++ch) {
  396. if ((index = allocation[ch][sb])) {
  397.   index = offset_table[bitalloc_table[offsets[sb]].offset][index - 1];
  399.   II_samples(&stream->ptr, &qc_table[index], samples);
  401.   for (s = 0; s < 3; ++s) {
  402.     frame->sbsample[ch][3 * gr + s][sb] =
  403.       mad_f_mul(samples[s], sf_table[scalefactor[ch][sb][gr / 4]]);
  404.   }
  405. }
  406. else {
  407.   for (s = 0; s < 3; ++s)
  408.     frame->sbsample[ch][3 * gr + s][sb] = 0;
  409. }
  410.       }
  411.     }
  413.     for (sb = bound; sb < sblimit; ++sb) {
  414.       if ((index = allocation[0][sb])) {
  415. index = offset_table[bitalloc_table[offsets[sb]].offset][index - 1];
  417. II_samples(&stream->ptr, &qc_table[index], samples);
  419. for (ch = 0; ch < nch; ++ch) {
  420.   for (s = 0; s < 3; ++s) {
  421.     frame->sbsample[ch][3 * gr + s][sb] =
  422.       mad_f_mul(samples[s], sf_table[scalefactor[ch][sb][gr / 4]]);
  423.   }
  424. }
  425.       }
  426.       else {
  427. for (ch = 0; ch < nch; ++ch) {
  428.   for (s = 0; s < 3; ++s)
  429.     frame->sbsample[ch][3 * gr + s][sb] = 0;
  430. }
  431.       }
  432.     }
  434.     for (ch = 0; ch < nch; ++ch) {
  435.       for (s = 0; s < 3; ++s) {
  436. for (sb = sblimit; sb < 32; ++sb)
  437.   frame->sbsample[ch][3 * gr + s][sb] = 0;
  438.       }
  439.     }
  440.   }
  442.   return 0;
  443. }