Category

Development Platform
Home > SourceCode/Document > Multimedia program > DVD > View Code
subband_layer_2.cc
Package: dvd-munitions.tar.gz　[view]
Upload User: aoeyumen
Upload Date: 2007-01-06
Package Size: 3329k
Code Size: 19k
Category:

DVD

Development Platform:

Unix_Linux

subband_layer_2.cc：Code Content
  1. /*
  2.    File: subband_layer_2.cc
  3. */
  5. #include <stdio.h>
  6. #include <stdlib.h>
  7. #include <fstream.h>
  8. #ifdef IRIX
  9. #include <dmedia/audio.h>
  10. #endif
  11. #ifdef SOLARIS
  12. #include <sys/audioio.h>
  13. #endif
  15. #include "athread.hh"
  17. #include "error.hh"
  18. #include "debug.hh"
  19. #include "util.hh"
  21. #include "sync.hh"
  22. #include "mpeg2const.hh"
  23. #include "mpeg2buff.hh"
  25. #include "astream.hh"
  26. #include "crc.hh"
  27. #include "header.hh"
  28. #include "obuffer.hh"
  29. #include "synthesis_filter.hh"
  30. #include "subband.hh"
  31. #include "subband_layer_1.hh"
  32. #include "subband_layer_2.hh"
  33. #include "scalefactors.hh"
  36. // data taken from ISO/IEC DIS 11172, Annexes 3-B.2[abcd] and 3-B.4:
  38. // subbands 0-2 in tables 3-B.2a and 2b: (index is allocation)
  39. static const uint32 table_ab1_codelength[16] =
  40.   // bits per codeword
  41. { 0, 5, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16 };
  43. static const uint32 table_ab1_quantizationsteps[16] =
  44.   // number of steps if grouped, 0 if ungrouped
  45. { 0, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 };
  47. static const real table_ab1_factor[16] =
  48.   // factor for requantization: (real)sample * factor - 1.0 gives requantized sample
  49. { 0.0, 1.0/2.0, 1.0/4.0, 1.0/8.0, 1.0/16.0, 1.0/32.0, 1.0/64.0,
  50.   1.0/128.0, 1.0/256.0, 1.0/512.0, 1.0/1024.0, 1.0/2048.0,
  51.   1.0/4096.0, 1.0/8192.0, 1.0/16384.0, 1.0/32768.0 };
  53. static const real table_ab1_c[16] =
  54.   // factor c for requantization from table 3-B.4
  55. { 0.0,           1.33333333333, 1.14285714286, 1.06666666666, 1.03225806452,
  56.   1.01587301587, 1.00787401575, 1.00392156863, 1.00195694716, 1.00097751711,
  57.   1.00048851979, 1.00024420024, 1.00012208522, 1.00006103888, 1.00003051851,
  58.   1.00001525902 };
  59. static const real table_ab1_d[16] =
  60.   // addend d for requantization from table 3-B.4
  61. { 0.0,           0.50000000000, 0.25000000000, 0.12500000000, 0.06250000000,
  62.   0.03125000000, 0.01562500000, 0.00781250000, 0.00390625000, 0.00195312500,
  63.   0.00097656250, 0.00048828125, 0.00024414063, 0.00012207031, 0.00006103516,
  64.   0.00003051758 };
  66. // subbands 3-... tables 3-B.2a and 2b:
  67. static const uint32 table_ab234_quantizationsteps[16] =
  68. { 0, 3, 5, 0, 9, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 };
  70. // subbands 3-10 in tables 3-B.2a and 2b:
  71. static const uint32 table_ab2_codelength[16] =
  72. { 0, 5, 7, 3, 10, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 16 };
  73. static const real table_ab2_factor[16] =
  74. { 0.0, 1.0/2.0, 1.0/4.0, 1.0/4.0, 1.0/8.0, 1.0/8.0, 1.0/16.0,
  75.   1.0/32.0, 1.0/64.0, 1.0/128.0, 1.0/256.0, 1.0/512.0,
  76.   1.0/1024.0, 1.0/2048.0, 1.0/4096.0, 1.0/32768.0 };
  77. static const real table_ab2_c[16] =
  78. { 0.0,           1.33333333333, 1.60000000000, 1.14285714286, 1.77777777777,
  79.   1.06666666666, 1.03225806452, 1.01587301587, 1.00787401575, 1.00392156863,
  80.   1.00195694716, 1.00097751711, 1.00048851979, 1.00024420024, 1.00012208522,
  81.   1.00001525902 };
  82. static const real table_ab2_d[16] =
  83. { 0.0,           0.50000000000, 0.50000000000, 0.25000000000, 0.50000000000,
  84.   0.12500000000, 0.06250000000, 0.03125000000, 0.01562500000, 0.00781250000,
  85.   0.00390625000, 0.00195312500, 0.00097656250, 0.00048828125, 0.00024414063,
  86.   0.00003051758 };
  88. // subbands 11-22 in tables 3-B.2a and 2b:
  89. static const uint32 table_ab3_codelength[8] = { 0, 5, 7, 3, 10, 4, 5, 16 };
  90. static const real table_ab3_factor[8] =
  91. { 0.0, 1.0/2.0, 1.0/4.0, 1.0/4.0, 1.0/8.0, 1.0/8.0, 1.0/16.0, 1.0/32768.0 };
  92. static const real table_ab3_c[8] =
  93. { 0.0,           1.33333333333, 1.60000000000, 1.14285714286, 1.77777777777,
  94.   1.06666666666, 1.03225806452, 1.00001525902 };
  95. static const real table_ab3_d[8] =
  96. { 0.0,           0.50000000000, 0.50000000000, 0.25000000000, 0.50000000000,
  97.   0.12500000000, 0.06250000000, 0.00003051758 };
  99. // subbands 23-... in tables 3-B.2a and 2b:
  100. static const uint32 table_ab4_codelength[4] = { 0, 5, 7, 16 };
  101. static const real table_ab4_factor[8] = { 0.0, 1.0/2.0, 1.0/4.0, 1.0/32768.0 };
  102. static const real table_ab4_c[4] = { 0.0, 1.33333333333, 1.60000000000, 1.00001525902 };
  103. static const real table_ab4_d[4] = { 0.0, 0.50000000000, 0.50000000000, 0.00003051758 };
  105. // subbands in tables 3-B.2c and 2d:
  106. static const uint32 table_cd_codelength[16] =
  107. { 0, 5, 7, 10, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15 };
  108. static const uint32 table_cd_quantizationsteps[16] =
  109. { 0, 3, 5, 9, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 };
  110. static const real table_cd_factor[16] =
  111. { 0.0, 1.0/2.0, 1.0/4.0, 1.0/8.0, 1.0/8.0, 1.0/16.0, 1.0/32.0, 1.0/64.0,
  112.   1.0/128.0, 1.0/256.0, 1.0/512.0, 1.0/1024.0, 1.0/2048.0, 1.0/4096.0,
  113.   1.0/8192.0, 1.0/16384.0 };
  114. static const real table_cd_c[16] =
  115. { 0.0,           1.33333333333, 1.60000000000, 1.77777777777, 1.06666666666,
  116.   1.03225806452, 1.01587301587, 1.00787401575, 1.00392156863, 1.00195694716,
  117.   1.00097751711, 1.00048851979, 1.00024420024, 1.00012208522, 1.00006103888,
  118.   1.00003051851 };
  119. static const real table_cd_d[16] =
  120. { 0.0,           0.50000000000, 0.50000000000, 0.50000000000, 0.12500000000,
  121.   0.06250000000, 0.03125000000, 0.01562500000, 0.00781250000, 0.00390625000,
  122.   0.00195312500, 0.00097656250, 0.00048828125, 0.00024414063, 0.00012207031,
  123.   0.00006103516 };
  127. /**********************/ // used for single channel mode
  128. /*** Standard Class ***/ // and in derived class for intensity
  129. /**********************/ // stereo mode
  131. SubbandLayer2::SubbandLayer2(uint32 subbandnumber){
  132.   this->subbandnumber=subbandnumber;
  133.   groupnumber=samplenumber=0;
  134. }
  137. uint32 SubbandLayer2::get_allocationlength (Header *header){
  138.   uint32 channel_bitrate = header->bitrate_index ();
  140.   // calculate bitrate per channel:
  141.   if (header->mode () != single_channel)
  142.     if (channel_bitrate == 4)
  143.       channel_bitrate = 1;
  144.     else
  145.       channel_bitrate -= 4;
  147.   if (channel_bitrate == 1 || channel_bitrate == 2)
  148.     // table 3-B.2c or 3-B.2d
  149.     if (subbandnumber <= 1)
  150.       return 4;
  151.     else
  152.       return 3;
  153.   else
  154.     // tables 3-B.2a or 3-B.2b
  155.     if (subbandnumber <= 10)
  156.       return 4;
  157.     else if (subbandnumber <= 22)
  158.       return 3;
  159.     else
  160.       return 2;
  161. }
  164. uint32 SubbandLayer2::prepare_sample_reading(Header *header, uint32 allocation, bool *grouping,
  165.  uint32 *quantizationsteps, real *factor, uint32 *codelength, real *c, real *d){
  166.   uint32 channel_bitrate = header->bitrate_index ();
  168.   // calculate bitrate per channel:
  169.   if (header->mode () != single_channel)
  170.     if (channel_bitrate == 4)  channel_bitrate = 1;
  171.     else                       channel_bitrate -= 4;
  173.   *grouping = False;
  174.   if (channel_bitrate == 1 || channel_bitrate == 2){
  175.     // table 3-B.2c or 3-B.2d
  176.     if ((*quantizationsteps=table_cd_quantizationsteps[allocation])!=0)
  177.       *grouping=True;
  178.     *factor=table_cd_factor[allocation];
  179.     *codelength=table_cd_codelength[allocation];
  180.     *c=table_cd_c[allocation];
  181.     *d=table_cd_d[allocation];
  182.   }
  183.   else {
  184.     // tables 3-B.2a or 3-B.2b
  185.     if (subbandnumber <= 2){
  186.       if ((*quantizationsteps=table_ab1_quantizationsteps[allocation])!=0)
  187.         *grouping = True;
  188.       *factor = table_ab1_factor[allocation];
  189.       *codelength = table_ab1_codelength[allocation];
  190.       *c = table_ab1_c[allocation];
  191.       *d = table_ab1_d[allocation];
  192.     }
  193.     else {
  194.       if ((*quantizationsteps=table_ab234_quantizationsteps[allocation])!=0)
  195.         *grouping = True;
  196.       if (subbandnumber <= 10){
  197.         *factor = table_ab2_factor[allocation];
  198.         *codelength = table_ab2_codelength[allocation];
  199.         *c = table_ab2_c[allocation];
  200.         *d = table_ab2_d[allocation];
  201.       }
  202.       else if (subbandnumber <= 22){
  203.         *factor = table_ab3_factor[allocation];
  204.         *codelength = table_ab3_codelength[allocation];
  205.         *c = table_ab3_c[allocation];
  206.         *d = table_ab3_d[allocation];
  207.       }
  208.       else {
  209.         *factor = table_ab4_factor[allocation];
  210.         *codelength = table_ab4_codelength[allocation];
  211.         *c = table_ab4_c[allocation];
  212.         *d = table_ab4_d[allocation];
  213.       }
  214.     }
  215.   }
  216.   
  217.   return channel_bitrate;
  218. }
  221. void SubbandLayer2::read_allocation(AudioStream *stream, Header *header, Crc16 *crc){
  222.   register uint32 length = get_allocationlength (header);
  223.   allocation = stream->get_bits (length);
  224.   if (crc)
  225.     crc->add_bits (allocation, length);
  226. }
  229. void SubbandLayer2::read_scalefactor_selection(AudioStream *stream, Crc16 *crc){
  230.   if (allocation){
  231.     scfsi = stream->get_bits (2);
  232.     if (crc) crc->add_bits (scfsi, 2);
  233.   }
  234. }
  237. void SubbandLayer2::read_scalefactor(AudioStream *stream, Header *header){
  238.   if (allocation){
  239.     switch (scfsi){
  240.       case 0:
  241. scalefactor1 = stream->get_bits (6);
  242. scalefactor2 = stream->get_bits (6);
  243. scalefactor3 = stream->get_bits (6);
  244. break;
  245.       case 1:
  246. scalefactor1 = scalefactor2 = stream->get_bits (6);
  247. scalefactor3 = stream->get_bits (6);
  248. break;
  249.       case 2:
  250. scalefactor1 = scalefactor2 = scalefactor3 = stream->get_bits (6);
  251. break;
  252.       case 3:
  253. scalefactor1 = stream->get_bits (6);
  254. scalefactor2 = scalefactor3 = stream->get_bits (6);
  255. break;
  256.     }
  257.     if (scalefactor1 == 63 || scalefactor2 == 63 || scalefactor3 == 63)
  258.       cerr << "WARNING: stream contains an illegal scalefactor!n"; // MPEG-stream is corrupted!
  260.     prepare_sample_reading (header, allocation, &grouping, &quantizationsteps,
  261.     &factor, &codelength, &c, &d);
  262.   }
  263. }
  266. bool SubbandLayer2::read_sampledata(AudioStream *stream){
  267.   if (allocation)
  268.     if (grouping){
  269.       uint32 samplecode = stream->get_bits (codelength);
  270. #ifdef DEBUG
  271.       if (samplecode == (1 << codelength) - 1)
  272.         cerr << "WARNING: stream contains an illegal subband sample!n";  // MPEG-stream is corrupted!
  273. #endif
  274.       samples[0] = real (samplecode % quantizationsteps) * factor - 1.0;
  275.       samplecode /= quantizationsteps;
  276.       samples[1] = real (samplecode % quantizationsteps) * factor - 1.0;
  277.       samplecode /= quantizationsteps;
  278.       samples[2] = real (samplecode % quantizationsteps) * factor - 1.0;
  279.     }
  280.     else {
  281.       samples[0] = real (stream->get_bits (codelength)) * factor - 1.0;
  282. #ifdef DEBUG
  283.       if (samples[0] == (1 << codelength) - 1)
  284.         cerr << "WARNING: stream contains an illegal subband sample!n";  // MPEG-stream is corrupted!
  285. #endif
  286.       samples[1] = real (stream->get_bits (codelength)) * factor - 1.0;
  287. #ifdef DEBUG
  288.       if (samples[1] == (1 << codelength) - 1)
  289.         cerr << "WARNING: stream contains an illegal subband sample!n";  // MPEG-stream is corrupted!
  290. #endif
  291.       samples[2] = real (stream->get_bits (codelength)) * factor - 1.0;
  292. #ifdef DEBUG
  293.       if (samples[2] == (1 << codelength) - 1)
  294.         cerr << "WARNING: stream contains an illegal subband sample!n";  // MPEG-stream is corrupted!
  295. #endif
  296.     }
  298.   samplenumber = 0;
  299.   return (++groupnumber == 12) ? True : False;
  300. }
  303. bool SubbandLayer2::put_next_sample(e_channels channels,
  304.      SynthesisFilter *filter1, SynthesisFilter*){
  305. #ifdef DEBUG
  306.   if (samplenumber >= 3){
  307.     // don't call put_next_sample if the previvious call returned True;
  308.     cerr << "Illegal call to SubbandLayer2::put_next_sample()!n";
  309.     exit (1);
  310.   }
  311. #endif
  312.   if (allocation && channels != right){
  313.     register real sample = (samples[samplenumber] + d) * c;
  314.     if (groupnumber <= 4)
  315.       sample *= scalefactors[scalefactor1];
  316.     else if (groupnumber <= 8)
  317.       sample *= scalefactors[scalefactor2];
  318.     else
  319.       sample *= scalefactors[scalefactor3];
  320. #ifdef DEBUG
  321.     if (sample < -1.0 || sample > 1.0)
  322.       cerr << "WARNING: rescaled subband sample is not in [-1.0, 1.0]n";
  323.       // this should never occur
  324. #endif
  325.     if (sample < -1.0E-7 || sample > 1.0E-7)
  326.       filter1->input_sample (sample, subbandnumber);
  327.   }
  329.   return (++samplenumber == 3) ? True : False;
  330. }
  334. /******************************/
  335. /*** Intensity Stereo Class ***/
  336. /******************************/
  338. SubbandLayer2IntensityStereo::SubbandLayer2IntensityStereo (uint32 subbandnumber)
  339. : SubbandLayer2 (subbandnumber){
  340. }
  343. void SubbandLayer2IntensityStereo::read_scalefactor_selection (AudioStream *stream, Crc16 *crc){
  344.   if (allocation){
  345.     scfsi = stream->get_bits (2);
  346.     channel2_scfsi = stream->get_bits (2);
  347.     if (crc){
  348.       crc->add_bits (scfsi, 2);
  349.       crc->add_bits (channel2_scfsi, 2);
  350.     }
  351.   }
  352. }
  355. void SubbandLayer2IntensityStereo::read_scalefactor (AudioStream *stream, Header *header){
  356.   if (allocation){
  357.     SubbandLayer2::read_scalefactor (stream, header);
  358.     switch (channel2_scfsi)
  359.     {
  360.       case 0:
  361. channel2_scalefactor1 = stream->get_bits (6);
  362. channel2_scalefactor2 = stream->get_bits (6);
  363. channel2_scalefactor3 = stream->get_bits (6);
  364. break;
  365.       case 1:
  366. channel2_scalefactor1 = channel2_scalefactor2 = stream->get_bits (6);
  367. channel2_scalefactor3 = stream->get_bits (6);
  368. break;
  369.       case 2:
  370. channel2_scalefactor1 = channel2_scalefactor2 =
  371. channel2_scalefactor3 = stream->get_bits (6);
  372. break;
  373.       case 3:
  374. channel2_scalefactor1 = stream->get_bits (6);
  375. channel2_scalefactor2 = channel2_scalefactor3 = stream->get_bits (6);
  376. break;
  377.     }
  378.     if (channel2_scalefactor1 == 63 || channel2_scalefactor2 == 63 || channel2_scalefactor3 == 63)
  379.       cerr << "WARNING: stream contains an illegal scalefactor!n"; // MPEG-stream is corrupted!
  380.   }
  381. }
  384. bool SubbandLayer2IntensityStereo::put_next_sample (e_channels channels,
  385. SynthesisFilter *filter1, SynthesisFilter *filter2){
  386. #ifdef DEBUG
  387.   if (samplenumber >= 3)
  388.   {
  389.     // don't call put_next_sample if the previvious call returned True;
  390.     cerr << "Illegal call to SubbandLayer2::put_next_sample()!n";
  391.     exit (1);
  392.   }
  393. #endif
  394.   if (allocation)
  395.   {
  396.     register real sample = (samples[samplenumber] + d) * c;
  397.     if (channels == both)
  398.     {
  399.       register float sample2 = sample;
  400.       if (groupnumber <= 4)
  401.       {
  402. sample *= scalefactors[scalefactor1];
  403. sample2 *= scalefactors[channel2_scalefactor1];
  404.       }
  405.       else if (groupnumber <= 8)
  406.       {
  407. sample *= scalefactors[scalefactor2];
  408. sample2 *= scalefactors[channel2_scalefactor2];
  409.       }
  410.       else
  411.       {
  412. sample *= scalefactors[scalefactor3];
  413. sample2 *= scalefactors[channel2_scalefactor3];
  414.       }
  415. #ifdef DEBUG
  416.       if (sample < -1.0 || sample > 1.0 || sample2 < -1.0 || sample2 > 1.0)
  417. cerr << "WARNING: rescaled subband sample is not in [-1.0, 1.0]n";
  418. // this should never occur
  419. #endif
  420.       if (sample < -1.0E-7 || sample > 1.0E-7)
  421. filter1->input_sample (sample, subbandnumber);
  422.       if (sample2 < -1.0E-7 || sample2 > 1.0E-7)
  423. filter2->input_sample (sample2, subbandnumber);
  424.     }
  425.     else if (channels == left)
  426.     {
  427.       if (groupnumber <= 4)
  428. sample *= scalefactors[scalefactor1];
  429.       else if (groupnumber <= 8)
  430. sample *= scalefactors[scalefactor2];
  431.       else
  432. sample *= scalefactors[scalefactor3];
  433. #ifdef DEBUG
  434.       if (sample < -1.0 || sample > 1.0)
  435. cerr << "WARNING: rescaled subband sample is not in [-1.0, 1.0]n";
  436. // this should never occur
  437. #endif
  438.       if (sample < -1.0E-7 || sample > 1.0E-7)
  439. filter1->input_sample (sample, subbandnumber);
  440.     }
  441.     else
  442.     {
  443.       if (groupnumber <= 4)
  444. sample *= scalefactors[channel2_scalefactor1];
  445.       else if (groupnumber <= 8)
  446. sample *= scalefactors[channel2_scalefactor2];
  447.       else
  448. sample *= scalefactors[channel2_scalefactor3];
  449. #ifdef DEBUG
  450.       if (sample < -1.0 || sample > 1.0)
  451. cerr << "WARNING: rescaled subband sample is not in [-1.0, 1.0]n";
  452. // this should never occur
  453. #endif
  454.       if (sample < -1.0E-7 || sample > 1.0E-7)
  455. filter1->input_sample (sample, subbandnumber);
  456.     }
  457.   }
  459.   if (++samplenumber == 3)
  460.     return True;
  461.   else
  462.     return False;
  463. }
  467. /********************/
  468. /*** Stereo Class ***/
  469. /********************/
  471. SubbandLayer2Stereo::SubbandLayer2Stereo (uint32 subbandnumber)
  472. : SubbandLayer2 (subbandnumber)
  473. {
  474. }
  477. void SubbandLayer2Stereo::read_allocation (AudioStream *stream, Header *header, Crc16 *crc)
  478. {
  479.   uint32 length = get_allocationlength (header);
  480.   allocation = stream->get_bits (length);
  481.   channel2_allocation = stream->get_bits (length);
  482.   if (crc)
  483.   {
  484.     crc->add_bits (allocation, length);
  485.     crc->add_bits (channel2_allocation, length);
  486.   }
  487. }
  490. void SubbandLayer2Stereo::read_scalefactor_selection (AudioStream *stream, Crc16 *crc)
  491. {
  492.   if (allocation)
  493.   {
  494.     scfsi = stream->get_bits (2);
  495.     if (crc)
  496.       crc->add_bits (scfsi, 2);
  497.   }
  498.   if (channel2_allocation)
  499.   {
  500.     channel2_scfsi = stream->get_bits (2);
  501.     if (crc)
  502.       crc->add_bits (channel2_scfsi, 2);
  503.   }
  504. }
  507. void SubbandLayer2Stereo::read_scalefactor (AudioStream *stream, Header *header)
  508. {
  509.   SubbandLayer2::read_scalefactor (stream, header);
  510.   if (channel2_allocation)
  511.   {
  512.     switch (channel2_scfsi)
  513.     {
  514.       case 0:
  515. channel2_scalefactor1 = stream->get_bits (6);
  516. channel2_scalefactor2 = stream->get_bits (6);
  517. channel2_scalefactor3 = stream->get_bits (6);
  518. break;
  519.       case 1:
  520. channel2_scalefactor1 = channel2_scalefactor2 = stream->get_bits (6);
  521. channel2_scalefactor3 = stream->get_bits (6);
  522. break;
  523.       case 2:
  524. channel2_scalefactor1 = channel2_scalefactor2 =
  525. channel2_scalefactor3 = stream->get_bits (6);
  526. break;
  527.       case 3:
  528. channel2_scalefactor1 = stream->get_bits (6);
  529. channel2_scalefactor2 = channel2_scalefactor3 = stream->get_bits (6);
  530. break;
  531.     }
  532.     if (channel2_scalefactor1 == 63 || channel2_scalefactor2 == 63 || channel2_scalefactor3 == 63)
  533.       cerr << "WARNING: stream contains an illegal scalefactor!n"; // MPEG-stream is corrupted!
  535.     prepare_sample_reading (header, channel2_allocation, &channel2_grouping,
  536.     &channel2_quantizationsteps, &channel2_factor,
  537.     &channel2_codelength, &channel2_c, &channel2_d);
  538.   }
  539. }
  542. bool SubbandLayer2Stereo::read_sampledata (AudioStream *stream){
  543.   bool returnvalue = SubbandLayer2::read_sampledata (stream);
  545.   if (channel2_allocation)
  546.     if (channel2_grouping)
  547.     {
  548.       uint32 samplecode = stream->get_bits (channel2_codelength);
  549. #ifdef DEBUG
  550.       if (samplecode == (1 << channel2_codelength) - 1)
  551. cerr << "WARNING: stream contains an illegal subband sample!n";  // MPEG-stream is corrupted!
  552. #endif
  553.       channel2_samples[0] = real (samplecode % channel2_quantizationsteps) * channel2_factor - 1.0;
  554.       samplecode /= channel2_quantizationsteps;
  555.       channel2_samples[1] = real (samplecode % channel2_quantizationsteps) * channel2_factor - 1.0;
  556.       samplecode /= channel2_quantizationsteps;
  557.       channel2_samples[2] = real (samplecode % channel2_quantizationsteps) * channel2_factor - 1.0;
  558.     }
  559.     else
  560.     {
  561.       channel2_samples[0] = real (stream->get_bits (channel2_codelength)) * channel2_factor - 1.0;
  562. #ifdef DEBUG
  563.       if (channel2_samples[0] == (1 << channel2_codelength) - 1)
  564. cerr << "WARNING: stream contains an illegal subband sample!n";  // MPEG-stream is corrupted!
  565. #endif
  566.       channel2_samples[1] = real (stream->get_bits (channel2_codelength)) * channel2_factor - 1.0;
  567. #ifdef DEBUG
  568.       if (channel2_samples[1] == (1 << channel2_codelength) - 1)
  569. cerr << "WARNING: stream contains an illegal subband sample!n";  // MPEG-stream is corrupted!
  570. #endif
  571.       channel2_samples[2] = real (stream->get_bits (channel2_codelength)) * channel2_factor - 1.0;
  572. #ifdef DEBUG
  573.       if (channel2_samples[2] == (1 << channel2_codelength) - 1)
  574. cerr << "WARNING: stream contains an illegal subband sample!n";  // MPEG-stream is corrupted!
  575. #endif
  576.     }
  577.   return returnvalue;
  578. }
  581. bool SubbandLayer2Stereo::put_next_sample (e_channels channels,
  582.    SynthesisFilter *filter1, SynthesisFilter *filter2)
  583. {
  584.   bool returnvalue = SubbandLayer2::put_next_sample (channels, filter1, filter2);
  585.   if (channel2_allocation && channels != left)
  586.   {
  587.     register real sample = (channel2_samples[samplenumber - 1] + channel2_d) * channel2_c;
  588.     if (groupnumber <= 4)
  589.       sample *= scalefactors[channel2_scalefactor1];
  590.     else if (groupnumber <= 8)
  591.       sample *= scalefactors[channel2_scalefactor2];
  592.     else
  593.       sample *= scalefactors[channel2_scalefactor3];
  594. #ifdef DEBUG
  595.     if (sample < -1.0 || sample > 1.0)
  596.       cerr << "WARNING: rescaled subband sample is not in [-1.0, 1.0]n";
  597.       // this should never occur
  598. #endif
  599.     if (sample < -1.0E-7 || sample > 1.0E-7)
  600.       if (channels == both)
  601. filter2->input_sample (sample, subbandnumber);
  602.       else
  603. filter1->input_sample (sample, subbandnumber);
  604.   }
  605.   return returnvalue;
  606. }