The JUCE cross-platform C++ framework, with DISTRHO/KXStudio specific changes
You can not select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.

3146 lines
131KB

  1. /*
  2. ==============================================================================
  3. This file is part of the JUCE library.
  4. Copyright (c) 2013 - Raw Material Software Ltd.
  5. Permission is granted to use this software under the terms of either:
  6. a) the GPL v2 (or any later version)
  7. b) the Affero GPL v3
  8. Details of these licenses can be found at: www.gnu.org/licenses
  9. JUCE is distributed in the hope that it will be useful, but WITHOUT ANY
  10. WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
  11. A PARTICULAR PURPOSE. See the GNU General Public License for more details.
  12. ------------------------------------------------------------------------------
  13. To release a closed-source product which uses JUCE, commercial licenses are
  14. available: visit www.juce.com for more information.
  15. ==============================================================================
  16. */
  17. /*
  18. IMPORTANT DISCLAIMER: By choosing to enable the JUCE_USE_MP3AUDIOFORMAT flag and
  19. to compile this MP3 code into your software, you do so AT YOUR OWN RISK! By doing so,
  20. you are agreeing that Raw Material Software is in no way responsible for any patent,
  21. copyright, or other legal issues that you may suffer as a result.
  22. The code in juce_MP3AudioFormat.cpp is NOT guaranteed to be free from infringements of 3rd-party
  23. intellectual property. If you wish to use it, please seek your own independent advice about the
  24. legality of doing so. If you are not willing to accept full responsibility for the consequences
  25. of using this code, then do not enable the JUCE_USE_MP3AUDIOFORMAT setting.
  26. */
  27. #if JUCE_USE_MP3AUDIOFORMAT
  28. namespace MP3Decoder
  29. {
  30. struct AllocationTable { int16 bits, d; };
  31. const struct AllocationTable allocTable0[] =
  32. {
  33. {4, 0}, {5, 3}, {3, -3}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {8, -127}, {9, -255}, {10, -511}, {11, -1023}, {12, -2047}, {13, -4095}, {14, -8191}, {15, -16383}, {16, -32767},
  34. {4, 0}, {5, 3}, {3, -3}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {8, -127}, {9, -255}, {10, -511}, {11, -1023}, {12, -2047}, {13, -4095}, {14, -8191}, {15, -16383}, {16, -32767},
  35. {4, 0}, {5, 3}, {3, -3}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {8, -127}, {9, -255}, {10, -511}, {11, -1023}, {12, -2047}, {13, -4095}, {14, -8191}, {15, -16383}, {16, -32767},
  36. {4, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {8, -127}, {9, -255}, {10, -511}, {11, -1023}, {12, -2047}, {13, -4095}, {16, -32767},
  37. {4, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {8, -127}, {9, -255}, {10, -511}, {11, -1023}, {12, -2047}, {13, -4095}, {16, -32767},
  38. {4, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {8, -127}, {9, -255}, {10, -511}, {11, -1023}, {12, -2047}, {13, -4095}, {16, -32767},
  39. {4, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {8, -127}, {9, -255}, {10, -511}, {11, -1023}, {12, -2047}, {13, -4095}, {16, -32767},
  40. {4, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {8, -127}, {9, -255}, {10, -511}, {11, -1023}, {12, -2047}, {13, -4095}, {16, -32767},
  41. {4, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {8, -127}, {9, -255}, {10, -511}, {11, -1023}, {12, -2047}, {13, -4095}, {16, -32767},
  42. {4, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {8, -127}, {9, -255}, {10, -511}, {11, -1023}, {12, -2047}, {13, -4095}, {16, -32767},
  43. {4, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {8, -127}, {9, -255}, {10, -511}, {11, -1023}, {12, -2047}, {13, -4095}, {16, -32767},
  44. {3, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {16, -32767}, {3, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {16, -32767},
  45. {3, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {16, -32767}, {3, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {16, -32767},
  46. {3, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {16, -32767}, {3, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {16, -32767},
  47. {3, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {16, -32767}, {3, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {16, -32767},
  48. {3, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {16, -32767}, {3, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {16, -32767},
  49. {3, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {16, -32767}, {3, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {16, -32767},
  50. {2, 0}, {5, 3}, {7, 5}, {16, -32767}, {2, 0}, {5, 3}, {7, 5}, {16, -32767}, {2, 0}, {5, 3}, {7, 5}, {16, -32767}, {2, 0}, {5, 3}, {7, 5}, {16, -32767}
  51. };
  52. const struct AllocationTable allocTable1[] =
  53. {
  54. {4, 0}, {5, 3}, {3, -3}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {8, -127}, {9, -255}, {10, -511}, {11, -1023}, {12, -2047}, {13, -4095}, {14, -8191}, {15, -16383}, {16, -32767},
  55. {4, 0}, {5, 3}, {3, -3}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {8, -127}, {9, -255}, {10, -511}, {11, -1023}, {12, -2047}, {13, -4095}, {14, -8191}, {15, -16383}, {16, -32767},
  56. {4, 0}, {5, 3}, {3, -3}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {8, -127}, {9, -255}, {10, -511}, {11, -1023}, {12, -2047}, {13, -4095}, {14, -8191}, {15, -16383}, {16, -32767},
  57. {4, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {8, -127}, {9, -255}, {10, -511}, {11, -1023}, {12, -2047}, {13, -4095}, {16, -32767},
  58. {4, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {8, -127}, {9, -255}, {10, -511}, {11, -1023}, {12, -2047}, {13, -4095}, {16, -32767},
  59. {4, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {8, -127}, {9, -255}, {10, -511}, {11, -1023}, {12, -2047}, {13, -4095}, {16, -32767},
  60. {4, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {8, -127}, {9, -255}, {10, -511}, {11, -1023}, {12, -2047}, {13, -4095}, {16, -32767},
  61. {4, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {8, -127}, {9, -255}, {10, -511}, {11, -1023}, {12, -2047}, {13, -4095}, {16, -32767},
  62. {4, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {8, -127}, {9, -255}, {10, -511}, {11, -1023}, {12, -2047}, {13, -4095}, {16, -32767},
  63. {4, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {8, -127}, {9, -255}, {10, -511}, {11, -1023}, {12, -2047}, {13, -4095}, {16, -32767},
  64. {4, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {8, -127}, {9, -255}, {10, -511}, {11, -1023}, {12, -2047}, {13, -4095}, {16, -32767},
  65. {3, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {16, -32767}, {3, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {16, -32767},
  66. {3, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {16, -32767}, {3, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {16, -32767},
  67. {3, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {16, -32767}, {3, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {16, -32767},
  68. {3, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {16, -32767}, {3, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {16, -32767},
  69. {3, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {16, -32767}, {3, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {16, -32767},
  70. {3, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {16, -32767}, {3, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {16, -32767},
  71. {2, 0}, {5, 3}, {7, 5}, {16, -32767}, {2, 0}, {5, 3}, {7, 5}, {16, -32767}, {2, 0}, {5, 3}, {7, 5}, {16, -32767}, {2, 0}, {5, 3}, {7, 5}, {16, -32767},
  72. {2, 0}, {5, 3}, {7, 5}, {16, -32767}, {2, 0}, {5, 3}, {7, 5}, {16, -32767}, {2, 0}, {5, 3}, {7, 5}, {16, -32767}
  73. };
  74. const struct AllocationTable allocTable2[] =
  75. {
  76. {4, 0}, {5, 3}, {7, 5}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {8, -127}, {9, -255}, {10, -511}, {11, -1023}, {12, -2047}, {13, -4095}, {14, -8191}, {15, -16383},
  77. {4, 0}, {5, 3}, {7, 5}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {8, -127}, {9, -255}, {10, -511}, {11, -1023}, {12, -2047}, {13, -4095}, {14, -8191}, {15, -16383},
  78. {3, 0}, {5, 3}, {7, 5}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {3, 0}, {5, 3}, {7, 5}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63},
  79. {3, 0}, {5, 3}, {7, 5}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {3, 0}, {5, 3}, {7, 5}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63},
  80. {3, 0}, {5, 3}, {7, 5}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {3, 0}, {5, 3}, {7, 5}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}
  81. };
  82. const struct AllocationTable allocTable3[] =
  83. {
  84. {4, 0}, {5, 3}, {7, 5}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {8, -127}, {9, -255}, {10, -511}, {11, -1023}, {12, -2047}, {13, -4095}, {14, -8191}, {15, -16383},
  85. {4, 0}, {5, 3}, {7, 5}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {8, -127}, {9, -255}, {10, -511}, {11, -1023}, {12, -2047}, {13, -4095}, {14, -8191}, {15, -16383},
  86. {3, 0}, {5, 3}, {7, 5}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {3, 0}, {5, 3}, {7, 5}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63},
  87. {3, 0}, {5, 3}, {7, 5}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {3, 0}, {5, 3}, {7, 5}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63},
  88. {3, 0}, {5, 3}, {7, 5}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {3, 0}, {5, 3}, {7, 5}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63},
  89. {3, 0}, {5, 3}, {7, 5}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {3, 0}, {5, 3}, {7, 5}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63},
  90. {3, 0}, {5, 3}, {7, 5}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {3, 0}, {5, 3}, {7, 5}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}
  91. };
  92. const struct AllocationTable allocTable4[] =
  93. {
  94. {4, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {8, -127}, {9, -255}, {10, -511}, {11, -1023}, {12, -2047}, {13, -4095}, {14, -8191},
  95. {4, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {8, -127}, {9, -255}, {10, -511}, {11, -1023}, {12, -2047}, {13, -4095}, {14, -8191},
  96. {4, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {8, -127}, {9, -255}, {10, -511}, {11, -1023}, {12, -2047}, {13, -4095}, {14, -8191},
  97. {4, 0}, {5, 3}, {7, 5}, {3, -3}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {8, -127}, {9, -255}, {10, -511}, {11, -1023}, {12, -2047}, {13, -4095}, {14, -8191},
  98. {3, 0}, {5, 3}, {7, 5}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {3, 0}, {5, 3}, {7, 5}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63},
  99. {3, 0}, {5, 3}, {7, 5}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {3, 0}, {5, 3}, {7, 5}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63},
  100. {3, 0}, {5, 3}, {7, 5}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {3, 0}, {5, 3}, {7, 5}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63},
  101. {3, 0}, {5, 3}, {7, 5}, {10, 9}, {4, -7}, {5, -15}, {6, -31}, {7, -63}, {2, 0}, {5, 3}, {7, 5}, {10, 9}, {2, 0}, {5, 3}, {7, 5}, {10, 9},
  102. {2, 0}, {5, 3}, {7, 5}, {10, 9}, {2, 0}, {5, 3}, {7, 5}, {10, 9}, {2, 0}, {5, 3}, {7, 5}, {10, 9}, {2, 0}, {5, 3}, {7, 5}, {10, 9},
  103. {2, 0}, {5, 3}, {7, 5}, {10, 9}, {2, 0}, {5, 3}, {7, 5}, {10, 9}, {2, 0}, {5, 3}, {7, 5}, {10, 9}, {2, 0}, {5, 3}, {7, 5}, {10, 9},
  104. {2, 0}, {5, 3}, {7, 5}, {10, 9}, {2, 0}, {5, 3}, {7, 5}, {10, 9}, {2, 0}, {5, 3}, {7, 5}, {10, 9}, {2, 0}, {5, 3}, {7, 5}, {10, 9},
  105. {2, 0}, {5, 3}, {7, 5}, {10, 9}, {2, 0}, {5, 3}, {7, 5}, {10, 9}, {2, 0}, {5, 3}, {7, 5}, {10, 9}, {2, 0}, {5, 3}, {7, 5}, {10, 9},
  106. {2, 0}, {5, 3}, {7, 5}, {10, 9}
  107. };
  108. struct BandInfoStruct
  109. {
  110. int16 longIndex[23];
  111. int16 longDiff[22];
  112. int16 shortIndex[14];
  113. int16 shortDiff[13];
  114. };
  115. const BandInfoStruct bandInfo[9] =
  116. {
  117. { {0, 4, 8, 12, 16, 20, 24, 30, 36, 44, 52, 62, 74, 90, 110, 134, 162, 196, 238, 288, 342, 418, 576},
  118. {4, 4, 4, 4, 4, 4, 6, 6, 8, 8, 10, 12, 16, 20, 24, 28, 34, 42, 50, 54, 76, 158},
  119. {0, 4 * 3, 8 * 3, 12 * 3, 16 * 3, 22 * 3, 30 * 3, 40 * 3, 52 * 3, 66 * 3, 84 * 3, 106 * 3, 136 * 3, 192 * 3},
  120. {4, 4, 4, 4, 6, 8, 10, 12, 14, 18, 22, 30, 56} },
  121. { {0, 4, 8, 12, 16, 20, 24, 30, 36, 42, 50, 60, 72, 88, 106, 128, 156, 190, 230, 276, 330, 384, 576},
  122. {4, 4, 4, 4, 4, 4, 6, 6, 6, 8, 10, 12, 16, 18, 22, 28, 34, 40, 46, 54, 54, 192},
  123. {0, 4 * 3, 8 * 3, 12 * 3, 16 * 3, 22 * 3, 28 * 3, 38 * 3, 50 * 3, 64 * 3, 80 * 3, 100 * 3, 126 * 3, 192 * 3},
  124. {4, 4, 4, 4, 6, 6, 10, 12, 14, 16, 20, 26, 66} },
  125. { {0, 4, 8, 12, 16, 20, 24, 30, 36, 44, 54, 66, 82, 102, 126, 156, 194, 240, 296, 364, 448, 550, 576},
  126. {4, 4, 4, 4, 4, 4, 6, 6, 8, 10, 12, 16, 20, 24, 30, 38, 46, 56, 68, 84, 102, 26},
  127. {0, 4 * 3, 8 * 3, 12 * 3, 16 * 3, 22 * 3, 30 * 3, 42 * 3, 58 * 3, 78 * 3, 104 * 3, 138 * 3, 180 * 3, 192 * 3},
  128. {4, 4, 4, 4, 6, 8, 12, 16, 20, 26, 34, 42, 12} },
  129. { {0, 6, 12, 18, 24, 30, 36, 44, 54, 66, 80, 96, 116, 140, 168, 200, 238, 284, 336, 396, 464, 522, 576},
  130. {6, 6, 6, 6, 6, 6, 8, 10, 12, 14, 16, 20, 24, 28, 32, 38, 46, 52, 60, 68, 58, 54 },
  131. {0, 4 * 3, 8 * 3, 12 * 3, 18 * 3, 24 * 3, 32 * 3, 42 * 3, 56 * 3, 74 * 3, 100 * 3, 132 * 3, 174 * 3, 192 * 3},
  132. {4, 4, 4, 6, 6, 8, 10, 14, 18, 26, 32, 42, 18 } },
  133. { {0, 6, 12, 18, 24, 30, 36, 44, 54, 66, 80, 96, 114, 136, 162, 194, 232, 278, 332, 394, 464, 540, 576},
  134. {6, 6, 6, 6, 6, 6, 8, 10, 12, 14, 16, 18, 22, 26, 32, 38, 46, 54, 62, 70, 76, 36 },
  135. {0, 4 * 3, 8 * 3, 12 * 3, 18 * 3, 26 * 3, 36 * 3, 48 * 3, 62 * 3, 80 * 3, 104 * 3, 136 * 3, 180 * 3, 192 * 3},
  136. {4, 4, 4, 6, 8, 10, 12, 14, 18, 24, 32, 44, 12 } },
  137. { {0, 6, 12, 18, 24, 30, 36, 44, 54, 66, 80, 96, 116, 140, 168, 200, 238, 284, 336, 396, 464, 522, 576},
  138. {6, 6, 6, 6, 6, 6, 8, 10, 12, 14, 16, 20, 24, 28, 32, 38, 46, 52, 60, 68, 58, 54 },
  139. {0, 4 * 3, 8 * 3, 12 * 3, 18 * 3, 26 * 3, 36 * 3, 48 * 3, 62 * 3, 80 * 3, 104 * 3, 134 * 3, 174 * 3, 192 * 3},
  140. {4, 4, 4, 6, 8, 10, 12, 14, 18, 24, 30, 40, 18 } },
  141. { {0, 6, 12, 18, 24, 30, 36, 44, 54, 66, 80, 96, 116, 140, 168, 200, 238, 284, 336, 396, 464, 522, 576},
  142. {6, 6, 6, 6, 6, 6, 8, 10, 12, 14, 16, 20, 24, 28, 32, 38, 46, 52, 60, 68, 58, 54},
  143. {0, 12, 24, 36, 54, 78, 108, 144, 186, 240, 312, 402, 522, 576},
  144. {4, 4, 4, 6, 8, 10, 12, 14, 18, 24, 30, 40, 18} },
  145. { {0, 6, 12, 18, 24, 30, 36, 44, 54, 66, 80, 96, 116, 140, 168, 200, 238, 284, 336, 396, 464, 522, 576},
  146. {6, 6, 6, 6, 6, 6, 8, 10, 12, 14, 16, 20, 24, 28, 32, 38, 46, 52, 60, 68, 58, 54},
  147. {0, 12, 24, 36, 54, 78, 108, 144, 186, 240, 312, 402, 522, 576},
  148. {4, 4, 4, 6, 8, 10, 12, 14, 18, 24, 30, 40, 18} },
  149. { {0, 12, 24, 36, 48, 60, 72, 88, 108, 132, 160, 192, 232, 280, 336, 400, 476, 566, 568, 570, 572, 574, 576},
  150. {12, 12, 12, 12, 12, 12, 16, 20, 24, 28, 32, 40, 48, 56, 64, 76, 90, 2, 2, 2, 2, 2},
  151. {0, 24, 48, 72, 108, 156, 216, 288, 372, 480, 486, 492, 498, 576},
  152. {8, 8, 8, 12, 16, 20, 24, 28, 36, 2, 2, 2, 26} }
  153. };
  154. const double decodeWindow[] =
  155. {
  156. 0.000000000, -0.000015259, -0.000015259, -0.000015259, -0.000015259, -0.000015259, -0.000015259, -0.000030518,
  157. -0.000030518, -0.000030518, -0.000030518, -0.000045776, -0.000045776, -0.000061035, -0.000061035, -0.000076294,
  158. -0.000076294, -0.000091553, -0.000106812, -0.000106812, -0.000122070, -0.000137329, -0.000152588, -0.000167847,
  159. -0.000198364, -0.000213623, -0.000244141, -0.000259399, -0.000289917, -0.000320435, -0.000366211, -0.000396729,
  160. -0.000442505, -0.000473022, -0.000534058, -0.000579834, -0.000625610, -0.000686646, -0.000747681, -0.000808716,
  161. -0.000885010, -0.000961304, -0.001037598, -0.001113892, -0.001205444, -0.001296997, -0.001388550, -0.001480103,
  162. -0.001586914, -0.001693726, -0.001785278, -0.001907349, -0.002014160, -0.002120972, -0.002243042, -0.002349854,
  163. -0.002456665, -0.002578735, -0.002685547, -0.002792358, -0.002899170, -0.002990723, -0.003082275, -0.003173828,
  164. -0.003250122, -0.003326416, -0.003387451, -0.003433228, -0.003463745, -0.003479004, -0.003479004, -0.003463745,
  165. -0.003417969, -0.003372192, -0.003280640, -0.003173828, -0.003051758, -0.002883911, -0.002700806, -0.002487183,
  166. -0.002227783, -0.001937866, -0.001617432, -0.001266479, -0.000869751, -0.000442505, 0.000030518, 0.000549316,
  167. 0.001098633, 0.001693726, 0.002334595, 0.003005981, 0.003723145, 0.004486084, 0.005294800, 0.006118774,
  168. 0.007003784, 0.007919312, 0.008865356, 0.009841919, 0.010848999, 0.011886597, 0.012939453, 0.014022827,
  169. 0.015121460, 0.016235352, 0.017349243, 0.018463135, 0.019577026, 0.020690918, 0.021789551, 0.022857666,
  170. 0.023910522, 0.024932861, 0.025909424, 0.026840210, 0.027725220, 0.028533936, 0.029281616, 0.029937744,
  171. 0.030532837, 0.031005859, 0.031387329, 0.031661987, 0.031814575, 0.031845093, 0.031738281, 0.031478882,
  172. 0.031082153, 0.030517578, 0.029785156, 0.028884888, 0.027801514, 0.026535034, 0.025085449, 0.023422241,
  173. 0.021575928, 0.019531250, 0.017257690, 0.014801025, 0.012115479, 0.009231567, 0.006134033, 0.002822876,
  174. -0.000686646, -0.004394531, -0.008316040, -0.012420654, -0.016708374, -0.021179199, -0.025817871, -0.030609131,
  175. -0.035552979, -0.040634155, -0.045837402, -0.051132202, -0.056533813, -0.061996460, -0.067520142, -0.073059082,
  176. -0.078628540, -0.084182739, -0.089706421, -0.095169067, -0.100540161, -0.105819702, -0.110946655, -0.115921021,
  177. -0.120697021, -0.125259399, -0.129562378, -0.133590698, -0.137298584, -0.140670776, -0.143676758, -0.146255493,
  178. -0.148422241, -0.150115967, -0.151306152, -0.151962280, -0.152069092, -0.151596069, -0.150497437, -0.148773193,
  179. -0.146362305, -0.143264771, -0.139450073, -0.134887695, -0.129577637, -0.123474121, -0.116577148, -0.108856201,
  180. -0.100311279, -0.090927124, -0.080688477, -0.069595337, -0.057617187, -0.044784546, -0.031082153, -0.016510010,
  181. -0.001068115, 0.015228271, 0.032379150, 0.050354004, 0.069168091, 0.088775635, 0.109161377, 0.130310059,
  182. 0.152206421, 0.174789429, 0.198059082, 0.221984863, 0.246505737, 0.271591187, 0.297210693, 0.323318481,
  183. 0.349868774, 0.376800537, 0.404083252, 0.431655884, 0.459472656, 0.487472534, 0.515609741, 0.543823242,
  184. 0.572036743, 0.600219727, 0.628295898, 0.656219482, 0.683914185, 0.711318970, 0.738372803, 0.765029907,
  185. 0.791213989, 0.816864014, 0.841949463, 0.866363525, 0.890090942, 0.913055420, 0.935195923, 0.956481934,
  186. 0.976852417, 0.996246338, 1.014617920, 1.031936646, 1.048156738, 1.063217163, 1.077117920, 1.089782715,
  187. 1.101211548, 1.111373901, 1.120223999, 1.127746582, 1.133926392, 1.138763428, 1.142211914, 1.144287109,
  188. 1.144989014
  189. };
  190. const int16 huffmanTab0[] = { 0 };
  191. const int16 huffmanTab1[] = { -5,-3,-1,17,1,16,0 };
  192. const int16 huffmanTab2[] = { -15,-11,-9,-5,-3,-1,34,2,18,-1,33,32,17,-1,1,16,0 };
  193. const int16 huffmanTab3[] = { -13,-11,-9,-5,-3,-1,34,2,18,-1,33,32,16,17,-1,1,0 };
  194. const int16 huffmanTab5[] = { -29,-25,-23,-15,-7,-5,-3,-1,51,35,50,49,-3,-1,19,3,-1,48,34,-3,-1,18,33,-1,2,32,17,-1,1,16,0 };
  195. const int16 huffmanTab6[] = { -25,-19,-13,-9,-5,-3,-1,51,3,35,-1,50,48,-1,19,49,-3,-1,34,2,18,-3,-1,33,32,1,-1,17,-1,16,0 };
  196. const int16 huffmanTab7[] =
  197. {
  198. -69,-65,-57,-39,-29,-17,-11,-7,-3,-1,85,69,-1,84,83,-1,53,68,-3,-1,37,82,21,-5,-1,81,-1,5,52,-1,80,-1,67,51,
  199. -5,-3,-1,36,66,20,-1,65,64,-11,-7,-3,-1,4,35,-1,50,3,-1,19,49,-3,-1,48,34,18,-5,-1,33,-1,2,32,17,-1,1,16,0
  200. };
  201. const int16 huffmanTab8[] =
  202. {
  203. -65,-63,-59,-45,-31,-19,-13,-7,-5,-3,-1,85,84,69,83,-3,-1,53,68,37,-3,-1,82,5,21,-5,-1,81,-1,52,67,-3,-1,80,
  204. 51,36,-5,-3,-1,66,20,65,-3,-1,4,64,-1,35,50,-9,-7,-3,-1,19,49,-1,3,48,34,-1,2,32,-1,18,33,17,-3,-1,1,16,0
  205. };
  206. const int16 huffmanTab9[] =
  207. {
  208. -63,-53,-41,-29,-19,-11,-5,-3,-1,85,69,53,-1,83,-1,84,5,-3,-1,68,37,-1,82,21,-3,-1,81,52,-1,67,-1,80,4,-7,-3,
  209. -1,36,66,-1,51,64,-1,20,65,-5,-3,-1,35,50,19,-1,49,-1,3,48,-5,-3,-1,34,2,18,-1,33,32,-3,-1,17,1,-1,16,0
  210. };
  211. const int16 huffmanTab10[] =
  212. {
  213. -125,-121,-111,-83,-55,-35,-21,-13,-7,-3,-1,119,103,-1,118,87,-3,-1,117,102,71,-3,-1,116,86,-1,101,55,-9,-3,
  214. -1,115,70,-3,-1,85,84,99,-1,39,114,-11,-5,-3,-1,100,7,112,-1,98,-1,69,53,-5,-1,6,-1,83,68,23,-17,-5,-1,113,
  215. -1,54,38,-5,-3,-1,37,82,21,-1,81,-1,52,67,-3,-1,22,97,-1,96,-1,5,80,-19,-11,-7,-3,-1,36,66,-1,51,4,-1,20,
  216. 65,-3,-1,64,35,-1,50,3,-3,-1,19,49,-1,48,34,-7,-3,-1,18,33,-1,2,32,17,-1,1,16,0
  217. };
  218. const int16 huffmanTab11[] =
  219. {
  220. -121,-113,-89,-59,-43,-27,-17,-7,-3,-1,119,103,-1,118,117,-3,-1,102,71,-1,116,-1,87,85,-5,-3,-1,86,101,55,
  221. -1,115,70,-9,-7,-3,-1,69,84,-1,53,83,39,-1,114,-1,100,7,-5,-1,113,-1,23,112,-3,-1,54,99,-1,96,-1,68,37,-13,
  222. -7,-5,-3,-1,82,5,21,98,-3,-1,38,6,22,-5,-1,97,-1,81,52,-5,-1,80,-1,67,51,-1,36,66,-15,-11,-7,-3,-1,20,65,
  223. -1,4,64,-1,35,50,-1,19,49,-5,-3,-1,3,48,34,33,-5,-1,18,-1,2,32,17,-3,-1,1,16,0
  224. };
  225. const int16 huffmanTab12[] =
  226. {
  227. -115,-99,-73,-45,-27,-17,-9,-5,-3,-1,119,103,118,-1,87,117,-3,-1,102,71,-1,116,101,-3,-1,86,55,-3,-1,115,
  228. 85,39,-7,-3,-1,114,70,-1,100,23,-5,-1,113,-1,7,112,-1,54,99,-13,-9,-3,-1,69,84,-1,68,-1,6,5,-1,38,98,-5,
  229. -1,97,-1,22,96,-3,-1,53,83,-1,37,82,-17,-7,-3,-1,21,81,-1,52,67,-5,-3,-1,80,4,36,-1,66,20,-3,-1,51,65,-1,
  230. 35,50,-11,-7,-5,-3,-1,64,3,48,19,-1,49,34,-1,18,33,-7,-5,-3,-1,2,32,0,17,-1,1,16
  231. };
  232. const int16 huffmanTab13[] =
  233. {
  234. -509,-503,-475,-405,-333,-265,-205,-153,-115,-83,-53,-35,-21,-13,-9,-7,-5,-3,-1,254,252,253,237,255,-1,239,223,
  235. -3,-1,238,207,-1,222,191,-9,-3,-1,251,206,-1,220,-1,175,233,-1,236,221,-9,-5,-3,-1,250,205,190,-1,235,159,-3,
  236. -1,249,234,-1,189,219,-17,-9,-3,-1,143,248,-1,204,-1,174,158,-5,-1,142,-1,127,126,247,-5,-1,218,-1,173,188,-3,
  237. -1,203,246,111,-15,-7,-3,-1,232,95,-1,157,217,-3,-1,245,231,-1,172,187,-9,-3,-1,79,244,-3,-1,202,230,243,-1,
  238. 63,-1,141,216,-21,-9,-3,-1,47,242,-3,-1,110,156,15,-5,-3,-1,201,94,171,-3,-1,125,215,78,-11,-5,-3,-1,200,214,
  239. 62,-1,185,-1,155,170,-1,31,241,-23,-13,-5,-1,240,-1,186,229,-3,-1,228,140,-1,109,227,-5,-1,226,-1,46,14,-1,30,
  240. 225,-15,-7,-3,-1,224,93,-1,213,124,-3,-1,199,77,-1,139,184,-7,-3,-1,212,154,-1,169,108,-1,198,61,-37,-21,-9,-5,
  241. -3,-1,211,123,45,-1,210,29,-5,-1,183,-1,92,197,-3,-1,153,122,195,-7,-5,-3,-1,167,151,75,209,-3,-1,13,208,-1,
  242. 138,168,-11,-7,-3,-1,76,196,-1,107,182,-1,60,44,-3,-1,194,91,-3,-1,181,137,28,-43,-23,-11,-5,-1,193,-1,152,12,
  243. -1,192,-1,180,106,-5,-3,-1,166,121,59,-1,179,-1,136,90,-11,-5,-1,43,-1,165,105,-1,164,-1,120,135,-5,-1,148,-1,
  244. 119,118,178,-11,-3,-1,27,177,-3,-1,11,176,-1,150,74,-7,-3,-1,58,163,-1,89,149,-1,42,162,-47,-23,-9,-3,-1,26,
  245. 161,-3,-1,10,104,160,-5,-3,-1,134,73,147,-3,-1,57,88,-1,133,103,-9,-3,-1,41,146,-3,-1,87,117,56,-5,-1,131,-1,
  246. 102,71,-3,-1,116,86,-1,101,115,-11,-3,-1,25,145,-3,-1,9,144,-1,72,132,-7,-5,-1,114,-1,70,100,40,-1,130,24,-41,
  247. -27,-11,-5,-3,-1,55,39,23,-1,113,-1,85,7,-7,-3,-1,112,54,-1,99,69,-3,-1,84,38,-1,98,53,-5,-1,129,-1,8,128,-3,
  248. -1,22,97,-1,6,96,-13,-9,-5,-3,-1,83,68,37,-1,82,5,-1,21,81,-7,-3,-1,52,67,-1,80,36,-3,-1,66,51,20,-19,-11,
  249. -5,-1,65,-1,4,64,-3,-1,35,50,19,-3,-1,49,3,-1,48,34,-3,-1,18,33,-1,2,32,-3,-1,17,1,16,0
  250. };
  251. const int16 huffmanTab15[] =
  252. {
  253. -495,-445,-355,-263,-183,-115,-77,-43,-27,-13,-7,-3,-1,255,239,-1,254,223,-1,238,-1,253,207,-7,-3,-1,252,222,-1,
  254. 237,191,-1,251,-1,206,236,-7,-3,-1,221,175,-1,250,190,-3,-1,235,205,-1,220,159,-15,-7,-3,-1,249,234,-1,189,219,
  255. -3,-1,143,248,-1,204,158,-7,-3,-1,233,127,-1,247,173,-3,-1,218,188,-1,111,-1,174,15,-19,-11,-3,-1,203,246,
  256. -3,-1,142,232,-1,95,157,-3,-1,245,126,-1,231,172,-9,-3,-1,202,187,-3,-1,217,141,79,-3,-1,244,63,-1,243,216,
  257. -33,-17,-9,-3,-1,230,47,-1,242,-1,110,240,-3,-1,31,241,-1,156,201,-7,-3,-1,94,171,-1,186,229,-3,-1,125,215,
  258. -1,78,228,-15,-7,-3,-1,140,200,-1,62,109,-3,-1,214,227,-1,155,185,-7,-3,-1,46,170,-1,226,30,-5,-1,225,-1,14,
  259. 224,-1,93,213,-45,-25,-13,-7,-3,-1,124,199,-1,77,139,-1,212,-1,184,154,-7,-3,-1,169,108,-1,198,61,-1,211,210,
  260. -9,-5,-3,-1,45,13,29,-1,123,183,-5,-1,209,-1,92,208,-1,197,138,-17,-7,-3,-1,168,76,-1,196,107,-5,-1,182,-1,
  261. 153,12,-1,60,195,-9,-3,-1,122,167,-1,166,-1,192,11,-1,194,-1,44,91,-55,-29,-15,-7,-3,-1,181,28,-1,137,152,-3,
  262. -1,193,75,-1,180,106,-5,-3,-1,59,121,179,-3,-1,151,136,-1,43,90,-11,-5,-1,178,-1,165,27,-1,177,-1,176,105,-7,
  263. -3,-1,150,74,-1,164,120,-3,-1,135,58,163,-17,-7,-3,-1,89,149,-1,42,162,-3,-1,26,161,-3,-1,10,160,104,-7,-3,
  264. -1,134,73,-1,148,57,-5,-1,147,-1,119,9,-1,88,133,-53,-29,-13,-7,-3,-1,41,103,-1,118,146,-1,145,-1,25,144,-7,
  265. -3,-1,72,132,-1,87,117,-3,-1,56,131,-1,102,71,-7,-3,-1,40,130,-1,24,129,-7,-3,-1,116,8,-1,128,86,-3,-1,101,
  266. 55,-1,115,70,-17,-7,-3,-1,39,114,-1,100,23,-3,-1,85,113,-3,-1,7,112,54,-7,-3,-1,99,69,-1,84,38,-3,-1,98,22,
  267. -3,-1,6,96,53,-33,-19,-9,-5,-1,97,-1,83,68,-1,37,82,-3,-1,21,81,-3,-1,5,80,52,-7,-3,-1,67,36,-1,66,51,-1,
  268. 65,-1,20,4,-9,-3,-1,35,50,-3,-1,64,3,19,-3,-1,49,48,34,-9,-7,-3,-1,18,33,-1,2,32,17,-3,-1,1,16,0
  269. };
  270. const int16 huffmanTab16[] =
  271. {
  272. -509,-503,-461,-323,-103,-37,-27,-15,-7,-3,-1,239,254,-1,223,253,-3,-1,207,252,-1,191,251,-5,-1,175,-1,250,159,
  273. -3,-1,249,248,143,-7,-3,-1,127,247,-1,111,246,255,-9,-5,-3,-1,95,245,79,-1,244,243,-53,-1,240,-1,63,-29,-19,
  274. -13,-7,-5,-1,206,-1,236,221,222,-1,233,-1,234,217,-1,238,-1,237,235,-3,-1,190,205,-3,-1,220,219,174,-11,-5,
  275. -1,204,-1,173,218,-3,-1,126,172,202,-5,-3,-1,201,125,94,189,242,-93,-5,-3,-1,47,15,31,-1,241,-49,-25,-13,
  276. -5,-1,158,-1,188,203,-3,-1,142,232,-1,157,231,-7,-3,-1,187,141,-1,216,110,-1,230,156,-13,-7,-3,-1,171,186,
  277. -1,229,215,-1,78,-1,228,140,-3,-1,200,62,-1,109,-1,214,155,-19,-11,-5,-3,-1,185,170,225,-1,212,-1,184,169,
  278. -5,-1,123,-1,183,208,227,-7,-3,-1,14,224,-1,93,213,-3,-1,124,199,-1,77,139,-75,-45,-27,-13,-7,-3,-1,154,
  279. 108,-1,198,61,-3,-1,92,197,13,-7,-3,-1,138,168,-1,153,76,-3,-1,182,122,60,-11,-5,-3,-1,91,137,28,-1,192,-1,
  280. 152,121,-1,226,-1,46,30,-15,-7,-3,-1,211,45,-1,210,209,-5,-1,59,-1,151,136,29,-7,-3,-1,196,107,-1,195,167,-1,
  281. 44,-1,194,181,-23,-13,-7,-3,-1,193,12,-1,75,180,-3,-1,106,166,179,-5,-3,-1,90,165,43,-1,178,27,-13,-5,-1,177,
  282. -1,11,176,-3,-1,105,150,-1,74,164,-5,-3,-1,120,135,163,-3,-1,58,89,42,-97,-57,-33,-19,-11,-5,-3,-1,149,104,161,
  283. -3,-1,134,119,148,-5,-3,-1,73,87,103,162,-5,-1,26,-1,10,160,-3,-1,57,147,-1,88,133,-9,-3,-1,41,146,-3,-1,118,
  284. 9,25,-5,-1,145,-1,144,72,-3,-1,132,117,-1,56,131,-21,-11,-5,-3,-1,102,40,130,-3,-1,71,116,24,-3,-1,129,128,-3,
  285. -1,8,86,55,-9,-5,-1,115,-1,101,70,-1,39,114,-5,-3,-1,100,85,7,23,-23,-13,-5,-1,113,-1,112,54,-3,-1,99,69,-1,
  286. 84,38,-3,-1,98,22,-1,97,-1,6,96,-9,-5,-1,83,-1,53,68,-1,37,82,-1,81,-1,21,5,-33,-23,-13,-7,-3,-1,52,67,-1,80,
  287. 36,-3,-1,66,51,20,-5,-1,65,-1,4,64,-1,35,50,-3,-1,19,49,-3,-1,3,48,34,-3,-1,18,33,-1,2,32,-3,-1,17,1,16,0
  288. };
  289. const int16 huffmanTab24[] =
  290. {
  291. -451,-117,-43,-25,-15,-7,-3,-1,239,254,-1,223,253,-3,-1,207,252,-1,191,251,-5,-1,250,-1,175,159,-1,249,248,-9,
  292. -5,-3,-1,143,127,247,-1,111,246,-3,-1,95,245,-1,79,244,-71,-7,-3,-1,63,243,-1,47,242,-5,-1,241,-1,31,240,-25,-9,
  293. -1,15,-3,-1,238,222,-1,237,206,-7,-3,-1,236,221,-1,190,235,-3,-1,205,220,-1,174,234,-15,-7,-3,-1,189,219,-1,204,
  294. 158,-3,-1,233,173,-1,218,188,-7,-3,-1,203,142,-1,232,157,-3,-1,217,126,-1,231,172,255,-235,-143,-77,-45,-25,-15,
  295. -7,-3,-1,202,187,-1,141,216,-5,-3,-1,14,224,13,230,-5,-3,-1,110,156,201,-1,94,186,-9,-5,-1,229,-1,171,125,-1,215,
  296. 228,-3,-1,140,200,-3,-1,78,46,62,-15,-7,-3,-1,109,214,-1,227,155,-3,-1,185,170,-1,226,30,-7,-3,-1,225,93,-1,213,124,
  297. -3,-1,199,77,-1,139,184,-31,-15,-7,-3,-1,212,154,-1,169,108,-3,-1,198,61,-1,211,45,-7,-3,-1,210,29,-1,123,183,-3,-1,
  298. 209,92,-1,197,138,-17,-7,-3,-1,168,153,-1,76,196,-3,-1,107,182,-3,-1,208,12,60,-7,-3,-1,195,122,-1,167,44,-3,-1,194,
  299. 91,-1,181,28,-57,-35,-19,-7,-3,-1,137,152,-1,193,75,-5,-3,-1,192,11,59,-3,-1,176,10,26,-5,-1,180,-1,106,166,-3,-1,121,
  300. 151,-3,-1,160,9,144,-9,-3,-1,179,136,-3,-1,43,90,178,-7,-3,-1,165,27,-1,177,105,-1,150,164,-17,-9,-5,-3,-1,74,120,135,
  301. -1,58,163,-3,-1,89,149,-1,42,162,-7,-3,-1,161,104,-1,134,119,-3,-1,73,148,-1,57,147,-63,-31,-15,-7,-3,-1,88,133,-1,41,
  302. 103,-3,-1,118,146,-1,25,145,-7,-3,-1,72,132,-1,87,117,-3,-1,56,131,-1,102,40,-17,-7,-3,-1,130,24,-1,71,116,-5,-1,129,
  303. -1,8,128,-1,86,101,-7,-5,-1,23,-1,7,112,115,-3,-1,55,39,114,-15,-7,-3,-1,70,100,-1,85,113,-3,-1,54,99,-1,69,84,-7,-3,
  304. -1,38,98,-1,22,97,-5,-3,-1,6,96,53,-1,83,68,-51,-37,-23,-15,-9,-3,-1,37,82,-1,21,-1,5,80,-1,81,-1,52,67,-3,-1,36,66,
  305. -1,51,20,-9,-5,-1,65,-1,4,64,-1,35,50,-1,19,49,-7,-5,-3,-1,3,48,34,18,-1,33,-1,2,32,-3,-1,17,1,-1,16,0
  306. };
  307. struct BitsToTableMap
  308. {
  309. uint32 bits;
  310. const int16* table;
  311. };
  312. const BitsToTableMap huffmanTables1[] =
  313. {
  314. { 0, huffmanTab0 }, { 0, huffmanTab1 }, { 0, huffmanTab2 }, { 0, huffmanTab3 },
  315. { 0, huffmanTab0 }, { 0, huffmanTab5 }, { 0, huffmanTab6 }, { 0, huffmanTab7 },
  316. { 0, huffmanTab8 }, { 0, huffmanTab9 }, { 0, huffmanTab10 }, { 0, huffmanTab11 },
  317. { 0, huffmanTab12 }, { 0, huffmanTab13 }, { 0, huffmanTab0 }, { 0, huffmanTab15 },
  318. { 1, huffmanTab16 }, { 2, huffmanTab16 }, { 3, huffmanTab16 }, { 4, huffmanTab16 },
  319. { 6, huffmanTab16 }, { 8, huffmanTab16 }, { 10, huffmanTab16 }, { 13, huffmanTab16 },
  320. { 4, huffmanTab24 }, { 5, huffmanTab24 }, { 6, huffmanTab24 }, { 7, huffmanTab24 },
  321. { 8, huffmanTab24 }, { 9, huffmanTab24 }, { 11, huffmanTab24 }, { 13, huffmanTab24 }
  322. };
  323. const int16 huffmanTabC0[] = { -29,-21,-13,-7,-3,-1,11,15,-1,13,14,-3,-1,7,5,9,-3,-1,6,3,-1,10,12,-3,-1,2,1,-1,4,8,0 };
  324. const int16 huffmanTabC1[] = { -15,-7,-3,-1,15,14,-1,13,12,-3,-1,11,10,-1,9,8,-7,-3,-1,7,6,-1,5,4,-3,-1,3,2,-1,1,0 };
  325. const BitsToTableMap huffmanTables2[] = { { 0, huffmanTabC0 }, { 0, huffmanTabC1 } };
  326. //==============================================================================
  327. struct VBRTagData
  328. {
  329. bool read (const uint8* data) noexcept
  330. {
  331. flags = 0;
  332. const int layer = (data[1] >> 1) & 3;
  333. if (layer != 1)
  334. return false;
  335. const int type = (data[1] >> 3) & 1;
  336. const int sampleRateIndex = (data[2] >> 2) & 3;
  337. const int mode = (data[3] >> 6) & 3;
  338. static const int bitRates[3][16] =
  339. {
  340. { 0, 8, 16, 24, 32, 40, 48, 56, 64, 80, 96, 112, 128, 144, 160, -1 }, // MPEG2
  341. { 0, 32, 40, 48, 56, 64, 80, 96, 112, 128, 160, 192, 224, 256, 320, -1 }, // MPEG1
  342. { 0, 8, 16, 24, 32, 40, 48, 56, 64, -1, -1, -1, -1, -1, -1, -1 }, // MPEG 2.5
  343. };
  344. const int bitrate = bitRates[type][((data[2] >> 4) & 15)];
  345. const int sampleRates[3][4] =
  346. {
  347. { 22050, 24000, 16000, -1 }, // MPEG2
  348. { 44100, 48000, 32000, -1 }, // MPEG1
  349. { 11025, 12000, 8000, -1 }, // MPEG2.5
  350. };
  351. if ((data[1] >> 4) == 0xe)
  352. sampleRate = sampleRates[2][sampleRateIndex];
  353. else
  354. sampleRate = sampleRates[type][sampleRateIndex];
  355. data += type != 0 ? (mode != 3 ? (32 + 4) : (17 + 4))
  356. : (mode != 3 ? (17 + 4) : (9 + 4));
  357. if (! isVbrTag (data))
  358. return false;
  359. data += 4;
  360. flags = ByteOrder::bigEndianInt (data);
  361. data += 4;
  362. if (flags & 1)
  363. {
  364. frames = ByteOrder::bigEndianInt (data);
  365. data += 4;
  366. }
  367. if (flags & 2)
  368. {
  369. bytes = ByteOrder::bigEndianInt (data);
  370. data += 4;
  371. }
  372. if (flags & 4)
  373. {
  374. if (toc != nullptr)
  375. for (int i = 0; i < 100; ++i)
  376. toc[i] = data[i];
  377. data += 100;
  378. }
  379. vbrScale = -1;
  380. if (flags & 8)
  381. vbrScale = (int) ByteOrder::bigEndianInt (data);
  382. headersize = ((type + 1) * 72000 * bitrate) / sampleRate;
  383. return true;
  384. }
  385. uint8 toc[100];
  386. int sampleRate, vbrScale, headersize;
  387. unsigned int flags, frames, bytes;
  388. private:
  389. static bool isVbrTag (const uint8* const d) noexcept
  390. {
  391. return (d[0] == 'X' && d[1] == 'i' && d[2] == 'n' && d[3] == 'g')
  392. || (d[0] == 'I' && d[1] == 'n' && d[2] == 'f' && d[3] == 'o');
  393. }
  394. };
  395. //==============================================================================
  396. struct MP3Frame
  397. {
  398. MP3Frame()
  399. {
  400. zeromem (this, sizeof (MP3Frame));
  401. single = -1;
  402. }
  403. void selectLayer2Table()
  404. {
  405. static const int translate[3][2][16] =
  406. {
  407. { { 0, 2, 2, 2, 2, 2, 2, 0, 0, 0, 1, 1, 1, 1, 1, 0 }, { 0, 2, 2, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0 } },
  408. { { 0, 2, 2, 2, 2, 2, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0 }, { 0, 2, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 } },
  409. { { 0, 3, 3, 3, 3, 3, 3, 0, 0, 0, 1, 1, 1, 1, 1, 0 }, { 0, 3, 3, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0 } }
  410. };
  411. static const AllocationTable* const tables[] = { allocTable0, allocTable1, allocTable2, allocTable3, allocTable4 };
  412. static const int limits[] = { 27, 30, 8, 12, 30 };
  413. const int index = lsf ? 4 : translate[sampleRateIndex][2 - numChannels][bitrateIndex];
  414. layer2SubBandLimit = limits [index];
  415. allocationTable = tables [index];
  416. }
  417. int getFrequency() const noexcept
  418. {
  419. const int frequencies[] = { 44100, 48000, 32000, 22050, 24000, 16000, 11025, 12000, 8000 };
  420. return frequencies [sampleRateIndex];
  421. }
  422. void decodeHeader (const uint32 header)
  423. {
  424. jassert (((header >> 10) & 3) != 3);
  425. mpeg25 = (header & (1 << 20)) == 0;
  426. lsf = mpeg25 ? 1 : ((header & (1 << 19)) ? 0 : 1);
  427. layer = 4 - ((header >> 17) & 3);
  428. sampleRateIndex = mpeg25 ? (6 + ((header >> 10) & 3)) : ((int) ((header >> 10) & 3) + (lsf * 3));
  429. crc16FollowsHeader = ((header >> 16) & 1) == 0;
  430. bitrateIndex = (header >> 12) & 15;
  431. padding = (header >> 9) & 1;
  432. mode = (header >> 6) & 3;
  433. modeExt = (header >> 4) & 3;
  434. //extension = (header >> 8) & 1;
  435. //copyright = (header >> 3) & 1;
  436. //original = (header >> 2) & 1;
  437. //emphasis = header & 3;
  438. numChannels = (mode == 3) ? 1 : 2;
  439. static const int frameSizes [2][3][16] =
  440. {
  441. { { 0, 32, 64, 96, 128, 160, 192, 224, 256, 288, 320, 352, 384, 416, 448 },
  442. { 0, 32, 48, 56, 64, 80, 96, 112, 128, 160, 192, 224, 256, 320, 384 },
  443. { 0, 32, 40, 48, 56, 64, 80, 96, 112, 128, 160, 192, 224, 256, 320 } },
  444. { { 0, 32, 48, 56, 64, 80, 96, 112, 128, 144, 160, 176, 192, 224, 256 },
  445. { 0, 8, 16, 24, 32, 40, 48, 56, 64, 80, 96, 112, 128, 144, 160 },
  446. { 0, 8, 16, 24, 32, 40, 48, 56, 64, 80, 96, 112, 128, 144, 160 } }
  447. };
  448. switch (layer)
  449. {
  450. case 1: frameSize = (((frameSizes[lsf][0][bitrateIndex] * 12000) / getFrequency() + padding) * 4) - 4; break;
  451. case 2: frameSize = (frameSizes[lsf][1][bitrateIndex] * 144000) / getFrequency() + (padding - 4); break;
  452. case 3: frameSize = (bitrateIndex == 0) ? 0 : ((frameSizes[lsf][2][bitrateIndex] * 144000) / (getFrequency() << lsf) + (padding - 4)); break;
  453. default: break;
  454. }
  455. }
  456. int layer, frameSize, numChannels, single;
  457. int lsf; // 0 = mpeg-1, 1 = mpeg-2/LSF
  458. bool mpeg25; // true = mpeg-2.5, false = mpeg-1/2
  459. bool crc16FollowsHeader;
  460. int bitrateIndex, sampleRateIndex, padding;
  461. int mode, modeExt, layer2SubBandLimit;
  462. enum { downSampleLimit = 32 };
  463. const AllocationTable* allocationTable;
  464. };
  465. //==============================================================================
  466. struct Constants
  467. {
  468. Constants()
  469. {
  470. cosTables[0] = cos64; cosTables[1] = cos32; cosTables[2] = cos16; cosTables[3] = cos8; cosTables[4] = cos4;
  471. initDecodeTables();
  472. initLayer2Tables();
  473. initLayer3Tables();
  474. }
  475. const uint8* getGroupTable (const int16 d1, const uint32 index) const noexcept
  476. {
  477. switch (d1)
  478. {
  479. case 3: return &group3tab [3 * jmin (index, 3u * 3u * 3u)];
  480. case 5: return &group5tab [3 * jmin (index, 5u * 5u * 5u)];
  481. case 9: return &group9tab [3 * jmin (index, 9u * 9u * 9u)];
  482. default: break;
  483. }
  484. static const uint8 dummy[] = { 0, 0, 0 };
  485. return dummy;
  486. }
  487. float muls[27][64];
  488. float nToThe4Over3[8207];
  489. float antiAliasingCa[8], antiAliasingCs[8];
  490. float win[4][36];
  491. float win1[4][36];
  492. float powToGains[256 + 118 + 4];
  493. int longLimit[9][23];
  494. int shortLimit[9][14];
  495. float tan1_1[16], tan2_1[16], tan1_2[16], tan2_2[16];
  496. float pow1_1[2][16], pow2_1[2][16], pow1_2[2][16], pow2_2[2][16];
  497. int* map[9][3];
  498. int* mapEnd[9][3];
  499. uint32 nLength2[512];
  500. uint32 iLength2[256];
  501. float decodeWin[512 + 32];
  502. float* cosTables[5];
  503. private:
  504. int mapbuf0[9][152];
  505. int mapbuf1[9][156];
  506. int mapbuf2[9][44];
  507. float cos64[16], cos32[8], cos16[4], cos8[2], cos4[1];
  508. uint8 group3tab [32 * 3];
  509. uint8 group5tab [128 * 3];
  510. uint8 group9tab [1024 * 3];
  511. void initDecodeTables()
  512. {
  513. int i, j, scaleval = -1;
  514. float* table = decodeWin;
  515. for (i = 0; i < 5; ++i)
  516. {
  517. int kr = 0x10 >> i;
  518. int divv = 0x40 >> i;
  519. float* costab = cosTables[i];
  520. for (int k = 0; k < kr; ++k)
  521. costab[k] = (float) (1.0 / (2.0 * cos (double_Pi * (k * 2 + 1) / divv)));
  522. }
  523. for (i = 0, j = 0; i < 256; ++i, ++j, table += 32)
  524. {
  525. if (table < decodeWin + 512 + 16)
  526. table[16] = table[0] = (float) (decodeWindow[j] * scaleval);
  527. if (i % 32 == 31)
  528. table -= 1023;
  529. if (i % 64 == 63)
  530. scaleval = -scaleval;
  531. }
  532. for (; i < 512; ++i, --j, table += 32)
  533. {
  534. if (table < decodeWin + 512 + 16)
  535. table[16] = table[0] = (float) (decodeWindow[j] * scaleval);
  536. if (i % 32 == 31) table -= 1023;
  537. if (i % 64 == 63) scaleval = -scaleval;
  538. }
  539. }
  540. void initLayer2Tables()
  541. {
  542. static const uint8 base[3][9] =
  543. {
  544. { 1, 0, 2 },
  545. { 17, 18, 0, 19, 20 },
  546. { 21, 1, 22, 23, 0, 24, 25, 2, 26 }
  547. };
  548. static const int tableLengths[] = { 3, 5, 9 };
  549. static uint8* tables[] = { group3tab, group5tab, group9tab };
  550. for (int i = 0; i < 3; ++i)
  551. {
  552. uint8* table = tables[i];
  553. const int len = tableLengths[i];
  554. for (int j = 0; j < len; ++j)
  555. for (int k = 0; k < len; ++k)
  556. for (int l = 0; l < len; ++l)
  557. {
  558. *table++ = base[i][l];
  559. *table++ = base[i][k];
  560. *table++ = base[i][j];
  561. }
  562. }
  563. for (int k = 0; k < 27; ++k)
  564. {
  565. static const double multipliers[] =
  566. {
  567. 0, -2.0 / 3.0, 2.0 / 3.0, 2.0 / 7.0, 2.0 / 15.0, 2.0 / 31.0, 2.0 / 63.0, 2.0 / 127.0, 2.0 / 255.0,
  568. 2.0 / 511.0, 2.0 / 1023.0, 2.0 / 2047.0, 2.0 / 4095.0, 2.0 / 8191.0, 2.0 / 16383.0, 2.0 / 32767.0, 2.0 / 65535.0,
  569. -4.0 / 5.0, -2.0 / 5.0, 2.0 / 5.0, 4.0 / 5.0, -8.0 / 9.0, -4.0 / 9.0, -2.0 / 9.0, 2.0 / 9.0, 4.0 / 9.0, 8.0 / 9.0
  570. };
  571. float* table = muls[k];
  572. for (int j = 3, i = 0; i < 63; ++i, --j)
  573. *table++ = (float) (multipliers[k] * pow (2.0, j / 3.0));
  574. *table++ = 0;
  575. }
  576. }
  577. void initLayer3Tables()
  578. {
  579. int i, j;
  580. for (i = -256; i < 118 + 4; ++i)
  581. powToGains[i + 256] = (float) pow (2.0, -0.25 * (i + 210));
  582. for (i = 0; i < 8207; ++i)
  583. nToThe4Over3[i] = (float) pow ((double) i, 4.0 / 3.0);
  584. for (i = 0; i < 8; ++i)
  585. {
  586. static double Ci[] = { -0.6, -0.535, -0.33, -0.185, -0.095, -0.041, -0.0142, -0.0037 };
  587. const double sq = sqrt (1.0 + Ci[i] * Ci[i]);
  588. antiAliasingCs[i] = (float) (1.0 / sq);
  589. antiAliasingCa[i] = (float) (Ci[i] / sq);
  590. }
  591. for (i = 0; i < 18; ++i)
  592. {
  593. win[0][i] = win[1][i] = (float) (0.5 * sin (double_Pi / 72.0 * (2 * i + 1)) / cos (double_Pi * (2 * i + 19) / 72.0));
  594. win[0][i + 18] = win[3][i + 18] = (float) (0.5 * sin (double_Pi / 72.0 * (2 * (i + 18) + 1)) / cos (double_Pi * (2 * (i + 18) + 19) / 72.0));
  595. }
  596. const double piOver72 = double_Pi;
  597. for (i = 0; i < 6; ++i)
  598. {
  599. win[1][i + 18] = (float) (0.5 / cos (piOver72 * (2 * (i + 18) + 19)));
  600. win[3][i + 12] = (float) (0.5 / cos (piOver72 * (2 * (i + 12) + 19)));
  601. win[1][i + 24] = (float) (0.5 * sin (double_Pi / 24.0 * (2 * i + 13)) / cos (piOver72 * (2 * (i + 24) + 19)));
  602. win[1][i + 30] = win[3][i] = 0;
  603. win[3][i + 6] = (float) (0.5 * sin (double_Pi / 24.0 * (2 * i + 1)) / cos (piOver72 * (2 * (i + 6) + 19)));
  604. }
  605. for (i = 0; i < 12; ++i)
  606. win[2][i] = (float) (0.5 * sin (double_Pi / 24.0 * (2 * i + 1)) / cos (double_Pi * (2 * i + 7) / 24.0));
  607. for (j = 0; j < 4; ++j)
  608. {
  609. static const int len[4] = { 36, 36, 12, 36 };
  610. for (i = 0; i < len[j]; i += 2) win1[j][i] = win[j][i];
  611. for (i = 1; i < len[j]; i += 2) win1[j][i] = -win[j][i];
  612. }
  613. const double sqrt2 = 1.41421356237309504880168872420969808;
  614. for (i = 0; i < 16; ++i)
  615. {
  616. const double t = tan (i * double_Pi / 12.0);
  617. tan1_1[i] = (float) (t / (1.0 + t));
  618. tan2_1[i] = (float) (1.0 / (1.0 + t));
  619. tan1_2[i] = (float) (sqrt2 * t / (1.0 + t));
  620. tan2_2[i] = (float) (sqrt2 / (1.0 + t));
  621. for (j = 0; j < 2; ++j)
  622. {
  623. double p1 = 1.0, p2 = 1.0;
  624. if (i > 0)
  625. {
  626. const double base = pow (2.0, -0.25 * (j + 1));
  627. if (i & 1)
  628. p1 = pow (base, (i + 1) * 0.5);
  629. else
  630. p2 = pow (base, i * 0.5);
  631. }
  632. pow1_1[j][i] = (float) p1;
  633. pow2_1[j][i] = (float) p2;
  634. pow1_2[j][i] = (float) (sqrt2 * p1);
  635. pow2_2[j][i] = (float) (sqrt2 * p2);
  636. }
  637. }
  638. for (j = 0; j < 9; ++j)
  639. {
  640. const BandInfoStruct& bi = bandInfo[j];
  641. int cb;
  642. int* mp = map[j][0] = mapbuf0[j];
  643. const int16* bdf = bi.longDiff;
  644. for (i = 0, cb = 0; cb < 8; ++cb, i += *bdf++)
  645. {
  646. *mp++ = (*bdf) >> 1;
  647. *mp++ = i;
  648. *mp++ = 3;
  649. *mp++ = cb;
  650. }
  651. bdf = bi.shortDiff + 3;
  652. for (cb = 3; cb < 13; ++cb)
  653. {
  654. const int l = (*bdf++) >> 1;
  655. for (int lwin = 0; lwin < 3; ++lwin)
  656. {
  657. *mp++ = l;
  658. *mp++ = i + lwin;
  659. *mp++ = lwin;
  660. *mp++ = cb;
  661. }
  662. i += 6 * l;
  663. }
  664. mapEnd[j][0] = mp;
  665. mp = map[j][1] = mapbuf1[j];
  666. bdf = bi.shortDiff;
  667. for (i = 0, cb = 0; cb < 13; ++cb)
  668. {
  669. const int l = (*bdf++) >> 1;
  670. for (int lwin = 0; lwin < 3; ++lwin)
  671. {
  672. *mp++ = l;
  673. *mp++ = i + lwin;
  674. *mp++ = lwin;
  675. *mp++ = cb;
  676. }
  677. i += 6 * l;
  678. }
  679. mapEnd[j][1] = mp;
  680. mp = map[j][2] = mapbuf2[j];
  681. bdf = bi.longDiff;
  682. for (cb = 0; cb < 22; ++cb)
  683. {
  684. *mp++ = (*bdf++) >> 1;
  685. *mp++ = cb;
  686. }
  687. mapEnd[j][2] = mp;
  688. }
  689. for (j = 0; j < 9; ++j)
  690. {
  691. for (i = 0; i < 23; ++i) longLimit[j][i] = jmin (32, (bandInfo[j].longIndex[i] - 1 + 8) / 18 + 1);
  692. for (i = 0; i < 14; ++i) shortLimit[j][i] = jmin (32, (bandInfo[j].shortIndex[i] - 1) / 18 + 1);
  693. }
  694. for (i = 0; i < 5; ++i)
  695. for (j = 0; j < 6; ++j)
  696. for (int k = 0; k < 6; ++k)
  697. {
  698. const int n = k + j * 6 + i * 36;
  699. iLength2[n] = (unsigned int) (i | (j << 3) | (k << 6) | (3 << 12));
  700. }
  701. for (i = 0; i < 4; ++i)
  702. for (j = 0; j < 4; ++j)
  703. for (int k = 0; k < 4; ++k)
  704. {
  705. const int n = k + j * 4 + i * 16;
  706. iLength2[n + 180] = (unsigned int) (i | (j << 3) | (k << 6) | (4 << 12));
  707. }
  708. for (i = 0; i < 4; ++i)
  709. for (j = 0; j < 3; ++j)
  710. {
  711. const int n = j + i * 3;
  712. iLength2[n + 244] = (unsigned int) (i | (j << 3) | (5 << 12));
  713. nLength2[n + 500] = (unsigned int) (i | (j << 3) | (2 << 12) | (1 << 15));
  714. }
  715. for (i = 0; i < 5; ++i)
  716. for (j = 0; j < 5; ++j)
  717. for (int k = 0; k < 4; ++k)
  718. for (int l = 0; l < 4; ++l)
  719. {
  720. const int n = l + k * 4 + j * 16 + i * 80;
  721. nLength2[n] = (unsigned int) (i | (j << 3) | (k << 6) | (l << 9) | (0 << 12));
  722. }
  723. for (i = 0; i < 5; ++i)
  724. for (j = 0; j < 5; ++j)
  725. for (int k = 0; k < 4; ++k)
  726. {
  727. const int n = k + j * 4 + i * 20;
  728. nLength2[n + 400] = (unsigned int) (i | (j << 3) | (k << 6) | (1 << 12));
  729. }
  730. }
  731. };
  732. static const Constants constants;
  733. //==============================================================================
  734. struct Layer3SideInfo
  735. {
  736. struct Info
  737. {
  738. void doAntialias (float xr[32][18]) const noexcept
  739. {
  740. float* xr1 = xr[1];
  741. int sb;
  742. if (blockType == 2)
  743. {
  744. if (mixedBlockFlag == 0)
  745. return;
  746. sb = 1;
  747. }
  748. else
  749. sb = (int) maxb - 1;
  750. for (; sb != 0; --sb, xr1 += 10)
  751. {
  752. const float* cs = constants.antiAliasingCs;
  753. const float* ca = constants.antiAliasingCa;
  754. float* xr2 = xr1;
  755. for (int ss = 7; ss >= 0; --ss)
  756. {
  757. const float bu = *--xr2, bd = *xr1;
  758. *xr2 = (bu * *cs) - (bd * *ca);
  759. *xr1++ = (bd * *cs++) + (bu * *ca++);
  760. }
  761. }
  762. }
  763. void doIStereo (float xrBuffer[2][32][18], const int* const scaleFactors,
  764. const int sampleRate, const bool msStereo, const int lsf) const noexcept
  765. {
  766. float (*xr) [32 * 18] = (float (*) [32 * 18]) xrBuffer;
  767. const BandInfoStruct& bi = bandInfo[sampleRate];
  768. const float* tabl1, *tabl2;
  769. if (lsf != 0)
  770. {
  771. const int p = scaleFactorCompression & 1;
  772. if (msStereo)
  773. {
  774. tabl1 = constants.pow1_2[p];
  775. tabl2 = constants.pow2_2[p];
  776. }
  777. else
  778. {
  779. tabl1 = constants.pow1_1[p];
  780. tabl2 = constants.pow2_1[p];
  781. }
  782. }
  783. else
  784. {
  785. if (msStereo)
  786. {
  787. tabl1 = constants.tan1_2;
  788. tabl2 = constants.tan2_2;
  789. }
  790. else
  791. {
  792. tabl1 = constants.tan1_1;
  793. tabl2 = constants.tan2_1;
  794. }
  795. }
  796. if (blockType == 2)
  797. {
  798. bool doL = mixedBlockFlag != 0;
  799. for (uint32 lwin = 0; lwin < 3; ++lwin)
  800. {
  801. uint32 sfb = maxBand[lwin];
  802. doL = doL && (sfb <= 3);
  803. for (; sfb < 12; ++sfb)
  804. {
  805. const int p = scaleFactors[sfb * 3 + lwin - mixedBlockFlag];
  806. if (p != 7)
  807. {
  808. const float t1 = tabl1[p];
  809. const float t2 = tabl2[p];
  810. int sb = bi.shortDiff[sfb];
  811. uint32 index = (uint32) sb + lwin;
  812. for (; sb > 0; --sb, index += 3)
  813. {
  814. float v = xr[0][index];
  815. xr[0][index] = v * t1;
  816. xr[1][index] = v * t2;
  817. }
  818. }
  819. }
  820. const int p = scaleFactors[11 * 3 + lwin - mixedBlockFlag];
  821. if (p != 7)
  822. {
  823. const float t1 = tabl1[p];
  824. const float t2 = tabl2[p];
  825. int sb = bi.shortDiff[12];
  826. uint32 index = (uint32) sb + lwin;
  827. for (; sb > 0; --sb, index += 3)
  828. {
  829. float v = xr[0][index];
  830. xr[0][index] = v * t1;
  831. xr[1][index] = v * t2;
  832. }
  833. }
  834. }
  835. if (doL)
  836. {
  837. int index = bi.longIndex[maxBandl];
  838. for (uint32 sfb = maxBandl; sfb < 8; ++sfb)
  839. {
  840. int sb = bi.longDiff[sfb];
  841. const int p = scaleFactors[sfb];
  842. if (p != 7)
  843. {
  844. const float t1 = tabl1[p];
  845. const float t2 = tabl2[p];
  846. for (; sb > 0; --sb, ++index)
  847. {
  848. float v = xr[0][index];
  849. xr[0][index] = v * t1;
  850. xr[1][index] = v * t2;
  851. }
  852. }
  853. else
  854. index += sb;
  855. }
  856. }
  857. }
  858. else
  859. {
  860. int index = bi.longIndex[maxBandl];
  861. for (uint32 sfb = maxBandl; sfb < 21; ++sfb)
  862. {
  863. int sb = bi.longDiff[sfb];
  864. const int p = scaleFactors[sfb];
  865. if (p != 7)
  866. {
  867. const float t1 = tabl1[p];
  868. const float t2 = tabl2[p];
  869. for (; sb > 0; --sb, ++index)
  870. {
  871. const float v = xr[0][index];
  872. xr[0][index] = v * t1;
  873. xr[1][index] = v * t2;
  874. }
  875. }
  876. else
  877. index += sb;
  878. }
  879. const int p = scaleFactors[20];
  880. if (p != 7)
  881. {
  882. const float t1 = tabl1[p], t2 = tabl2[p];
  883. for (int sb = bi.longDiff[21]; sb > 0; --sb, ++index)
  884. {
  885. const float v = xr[0][index];
  886. xr[0][index] = v * t1;
  887. xr[1][index] = v * t2;
  888. }
  889. }
  890. }
  891. }
  892. int scfsi;
  893. uint32 part2_3Length, bigValues;
  894. uint32 scaleFactorCompression, blockType, mixedBlockFlag;
  895. uint32 tableSelect[3];
  896. uint32 maxBand[3];
  897. uint32 maxBandl, maxb, region1Start, region2Start;
  898. uint32 preflag, scaleFactorScale, count1TableSelect;
  899. const float* fullGain[3];
  900. const float* pow2gain;
  901. };
  902. struct InfoPair { Info gr[2]; };
  903. InfoPair ch[2];
  904. uint32 mainDataStart, privateBits;
  905. };
  906. //==============================================================================
  907. namespace DCT
  908. {
  909. enum { SBLIMIT = 32 };
  910. static const float cos6_1 = 0.866025388f;
  911. static const float cos6_2 = 0.5f;
  912. static const float cos9[] = { 1.0f, 0.98480773f, 0.939692616f, 0.866025388f, 0.766044438f, 0.642787635f, 0.5f, 0.342020154f, 0.173648179f };
  913. static const float cos36[] = { 0.501909912f, 0.517638087f, 0.551688969f, 0.610387266f, 0.707106769f, 0.871723413f, 1.18310082f, 1.93185163f, 5.73685646f };
  914. static const float cos12[] = { 0.517638087f, 0.707106769f, 1.93185163f };
  915. inline void dct36_0 (const int v, float* const ts, float* const out1, float* const out2,
  916. const float* const wintab, float sum0, const float sum1) noexcept
  917. {
  918. const float tmp = sum0 + sum1;
  919. out2[9 + v] = tmp * wintab[27 + v];
  920. out2[8 - v] = tmp * wintab[26 - v];
  921. sum0 -= sum1;
  922. ts[SBLIMIT * (8 - v)] = out1[8 - v] + sum0 * wintab[8 - v];
  923. ts[SBLIMIT * (9 + v)] = out1[9 + v] + sum0 * wintab[9 + v];
  924. }
  925. inline void dct36_1 (const int v, float* const ts, float* const out1, float* const out2, const float* const wintab,
  926. const float tmp1a, const float tmp1b, const float tmp2a, const float tmp2b) noexcept
  927. {
  928. dct36_0 (v, ts, out1, out2, wintab, tmp1a + tmp2a, (tmp1b + tmp2b) * cos36[v]);
  929. }
  930. inline void dct36_2 (const int v, float* const ts, float* const out1, float* const out2, const float* const wintab,
  931. const float tmp1a, const float tmp1b, const float tmp2a, const float tmp2b) noexcept
  932. {
  933. dct36_0 (v, ts, out1, out2, wintab, tmp2a - tmp1a, (tmp2b - tmp1b) * cos36[v]);
  934. }
  935. static void dct36 (float* const in, float* const out1, float* const out2, const float* const wintab, float* const ts) noexcept
  936. {
  937. in[17] += in[16]; in[16] += in[15]; in[15] += in[14]; in[14] += in[13]; in[13] += in[12];
  938. in[12] += in[11]; in[11] += in[10]; in[10] += in[9]; in[9] += in[8]; in[8] += in[7];
  939. in[7] += in[6]; in[6] += in[5]; in[5] += in[4]; in[4] += in[3]; in[3] += in[2];
  940. in[2] += in[1]; in[1] += in[0]; in[17] += in[15]; in[15] += in[13]; in[13] += in[11];
  941. in[11] += in[9]; in[9] += in[7]; in[7] += in[5]; in[5] += in[3]; in[3] += in[1];
  942. const float ta33 = in[6] * cos9[3];
  943. const float ta66 = in[12] * cos9[6];
  944. const float tb33 = in[7] * cos9[3];
  945. const float tb66 = in[13] * cos9[6];
  946. {
  947. const float tmp1a = in[2] * cos9[1] + ta33 + in[10] * cos9[5] + in[14] * cos9[7];
  948. const float tmp1b = in[3] * cos9[1] + tb33 + in[11] * cos9[5] + in[15] * cos9[7];
  949. const float tmp2a = in[0] + in[4] * cos9[2] + in[8] * cos9[4] + ta66 + in[16] * cos9[8];
  950. const float tmp2b = in[1] + in[5] * cos9[2] + in[9] * cos9[4] + tb66 + in[17] * cos9[8];
  951. dct36_1 (0, ts, out1, out2, wintab, tmp1a, tmp1b, tmp2a, tmp2b);
  952. dct36_2 (8, ts, out1, out2, wintab, tmp1a, tmp1b, tmp2a, tmp2b);
  953. }
  954. {
  955. const float tmp1a = (in[2] - in[10] - in[14]) * cos9[3];
  956. const float tmp1b = (in[3] - in[11] - in[15]) * cos9[3];
  957. const float tmp2a = (in[4] - in[8] - in[16]) * cos9[6] - in[12] + in[0];
  958. const float tmp2b = (in[5] - in[9] - in[17]) * cos9[6] - in[13] + in[1];
  959. dct36_1 (1, ts, out1, out2, wintab, tmp1a, tmp1b, tmp2a, tmp2b);
  960. dct36_2 (7, ts, out1, out2, wintab, tmp1a, tmp1b, tmp2a, tmp2b);
  961. }
  962. {
  963. const float tmp1a = in[2] * cos9[5] - ta33 - in[10] * cos9[7] + in[14] * cos9[1];
  964. const float tmp1b = in[3] * cos9[5] - tb33 - in[11] * cos9[7] + in[15] * cos9[1];
  965. const float tmp2a = in[0] - in[4] * cos9[8] - in[8] * cos9[2] + ta66 + in[16] * cos9[4];
  966. const float tmp2b = in[1] - in[5] * cos9[8] - in[9] * cos9[2] + tb66 + in[17] * cos9[4];
  967. dct36_1 (2, ts, out1, out2, wintab, tmp1a, tmp1b, tmp2a, tmp2b);
  968. dct36_2 (6, ts, out1, out2, wintab, tmp1a, tmp1b, tmp2a, tmp2b);
  969. }
  970. {
  971. const float tmp1a = in[2] * cos9[7] - ta33 + in[10] * cos9[1] - in[14] * cos9[5];
  972. const float tmp1b = in[3] * cos9[7] - tb33 + in[11] * cos9[1] - in[15] * cos9[5];
  973. const float tmp2a = in[0] - in[4] * cos9[4] + in[8] * cos9[8] + ta66 - in[16] * cos9[2];
  974. const float tmp2b = in[1] - in[5] * cos9[4] + in[9] * cos9[8] + tb66 - in[17] * cos9[2];
  975. dct36_1 (3, ts, out1, out2, wintab, tmp1a, tmp1b, tmp2a, tmp2b);
  976. dct36_2 (5, ts, out1, out2, wintab, tmp1a, tmp1b, tmp2a, tmp2b);
  977. }
  978. const float sum0 = in[0] - in[4] + in[8] - in[12] + in[16];
  979. const float sum1 = (in[1] - in[5] + in[9] - in[13] + in[17]) * cos36[4];
  980. dct36_0 (4, ts, out1, out2, wintab, sum0, sum1);
  981. }
  982. struct DCT12Inputs
  983. {
  984. float in0, in1, in2, in3, in4, in5;
  985. inline DCT12Inputs (const float* const in) noexcept
  986. {
  987. in5 = in[5*3] + (in4 = in[4*3]);
  988. in4 += (in3 = in[3*3]);
  989. in3 += (in2 = in[2*3]);
  990. in2 += (in1 = in[1*3]);
  991. in1 += (in0 = in[0*3]);
  992. in5 += in3; in3 += in1;
  993. in2 *= cos6_1;
  994. in3 *= cos6_1;
  995. }
  996. inline void process() noexcept
  997. {
  998. in0 += in4 * cos6_2;
  999. in4 = in0 + in2; in0 -= in2;
  1000. in1 += in5 * cos6_2;
  1001. in5 = (in1 + in3) * cos12[0];
  1002. in1 = (in1 - in3) * cos12[2];
  1003. in3 = in4 + in5; in4 -= in5;
  1004. in2 = in0 + in1; in0 -= in1;
  1005. }
  1006. };
  1007. static void dct12 (const float* in, float* const out1, float* const out2, const float* wi, float* ts) noexcept
  1008. {
  1009. {
  1010. ts[0] = out1[0];
  1011. ts[SBLIMIT * 1] = out1[1];
  1012. ts[SBLIMIT * 2] = out1[2];
  1013. ts[SBLIMIT * 3] = out1[3];
  1014. ts[SBLIMIT * 4] = out1[4];
  1015. ts[SBLIMIT * 5] = out1[5];
  1016. DCT12Inputs inputs (in);
  1017. {
  1018. float tmp1 = (inputs.in0 - inputs.in4);
  1019. const float tmp2 = (inputs.in1 - inputs.in5) * cos12[1];
  1020. const float tmp0 = tmp1 + tmp2;
  1021. tmp1 -= tmp2;
  1022. ts[16 * SBLIMIT] = out1[16] + tmp0 * wi[10];
  1023. ts[13 * SBLIMIT] = out1[13] + tmp0 * wi[7];
  1024. ts[7 * SBLIMIT] = out1[7] + tmp1 * wi[1];
  1025. ts[10 * SBLIMIT] = out1[10] + tmp1 * wi[4];
  1026. }
  1027. inputs.process();
  1028. ts[17 * SBLIMIT] = out1[17] + inputs.in2 * wi[11];
  1029. ts[12 * SBLIMIT] = out1[12] + inputs.in2 * wi[6];
  1030. ts[14 * SBLIMIT] = out1[14] + inputs.in3 * wi[8];
  1031. ts[15 * SBLIMIT] = out1[15] + inputs.in3 * wi[9];
  1032. ts[6 * SBLIMIT] = out1[6] + inputs.in0 * wi[0];
  1033. ts[11 * SBLIMIT] = out1[11] + inputs.in0 * wi[5];
  1034. ts[8 * SBLIMIT] = out1[8] + inputs.in4 * wi[2];
  1035. ts[9 * SBLIMIT] = out1[9] + inputs.in4 * wi[3];
  1036. }
  1037. {
  1038. DCT12Inputs inputs (++in);
  1039. float tmp1 = (inputs.in0 - inputs.in4);
  1040. const float tmp2 = (inputs.in1 - inputs.in5) * cos12[1];
  1041. const float tmp0 = tmp1 + tmp2;
  1042. tmp1 -= tmp2;
  1043. out2[4] = tmp0 * wi[10];
  1044. out2[1] = tmp0 * wi[7];
  1045. ts[13 * SBLIMIT] += tmp1 * wi[1];
  1046. ts[16 * SBLIMIT] += tmp1 * wi[4];
  1047. inputs.process();
  1048. out2[5] = inputs.in2 * wi[11];
  1049. out2[0] = inputs.in2 * wi[6];
  1050. out2[2] = inputs.in3 * wi[8];
  1051. out2[3] = inputs.in3 * wi[9];
  1052. ts[12 * SBLIMIT] += inputs.in0 * wi[0];
  1053. ts[17 * SBLIMIT] += inputs.in0 * wi[5];
  1054. ts[14 * SBLIMIT] += inputs.in4 * wi[2];
  1055. ts[15 * SBLIMIT] += inputs.in4 * wi[5 - 2];
  1056. }
  1057. {
  1058. DCT12Inputs inputs (++in);
  1059. out2[12] = out2[13] = out2[14] = out2[15] = out2[16] = out2[17] = 0;
  1060. float tmp1 = (inputs.in0 - inputs.in4);
  1061. const float tmp2 = (inputs.in1 - inputs.in5) * cos12[1];
  1062. const float tmp0 = tmp1 + tmp2;
  1063. tmp1 -= tmp2;
  1064. out2[10] = tmp0 * wi[10];
  1065. out2[7] = tmp0 * wi[7];
  1066. out2[1] += tmp1 * wi[1];
  1067. out2[4] += tmp1 * wi[4];
  1068. inputs.process();
  1069. out2[11] = inputs.in2 * wi[11];
  1070. out2[6] = inputs.in2 * wi[6];
  1071. out2[8] = inputs.in3 * wi[8];
  1072. out2[9] = inputs.in3 * wi[9];
  1073. out2[0] += inputs.in0 * wi[0];
  1074. out2[5] += inputs.in0 * wi[5];
  1075. out2[2] += inputs.in4 * wi[2];
  1076. out2[3] += inputs.in4 * wi[3];
  1077. }
  1078. }
  1079. static void dct64 (float* const out0, float* const out1, float* const b1, float* const b2, const float* const samples) noexcept
  1080. {
  1081. {
  1082. const float* const costab = constants.cosTables[0];
  1083. b1[0x00] = samples[0x00] + samples[0x1F]; b1[0x1F] = (samples[0x00] - samples[0x1F]) * costab[0x0];
  1084. b1[0x01] = samples[0x01] + samples[0x1E]; b1[0x1E] = (samples[0x01] - samples[0x1E]) * costab[0x1];
  1085. b1[0x02] = samples[0x02] + samples[0x1D]; b1[0x1D] = (samples[0x02] - samples[0x1D]) * costab[0x2];
  1086. b1[0x03] = samples[0x03] + samples[0x1C]; b1[0x1C] = (samples[0x03] - samples[0x1C]) * costab[0x3];
  1087. b1[0x04] = samples[0x04] + samples[0x1B]; b1[0x1B] = (samples[0x04] - samples[0x1B]) * costab[0x4];
  1088. b1[0x05] = samples[0x05] + samples[0x1A]; b1[0x1A] = (samples[0x05] - samples[0x1A]) * costab[0x5];
  1089. b1[0x06] = samples[0x06] + samples[0x19]; b1[0x19] = (samples[0x06] - samples[0x19]) * costab[0x6];
  1090. b1[0x07] = samples[0x07] + samples[0x18]; b1[0x18] = (samples[0x07] - samples[0x18]) * costab[0x7];
  1091. b1[0x08] = samples[0x08] + samples[0x17]; b1[0x17] = (samples[0x08] - samples[0x17]) * costab[0x8];
  1092. b1[0x09] = samples[0x09] + samples[0x16]; b1[0x16] = (samples[0x09] - samples[0x16]) * costab[0x9];
  1093. b1[0x0A] = samples[0x0A] + samples[0x15]; b1[0x15] = (samples[0x0A] - samples[0x15]) * costab[0xA];
  1094. b1[0x0B] = samples[0x0B] + samples[0x14]; b1[0x14] = (samples[0x0B] - samples[0x14]) * costab[0xB];
  1095. b1[0x0C] = samples[0x0C] + samples[0x13]; b1[0x13] = (samples[0x0C] - samples[0x13]) * costab[0xC];
  1096. b1[0x0D] = samples[0x0D] + samples[0x12]; b1[0x12] = (samples[0x0D] - samples[0x12]) * costab[0xD];
  1097. b1[0x0E] = samples[0x0E] + samples[0x11]; b1[0x11] = (samples[0x0E] - samples[0x11]) * costab[0xE];
  1098. b1[0x0F] = samples[0x0F] + samples[0x10]; b1[0x10] = (samples[0x0F] - samples[0x10]) * costab[0xF];
  1099. }
  1100. {
  1101. const float* const costab = constants.cosTables[1];
  1102. b2[0x00] = b1[0x00] + b1[0x0F]; b2[0x0F] = (b1[0x00] - b1[0x0F]) * costab[0];
  1103. b2[0x01] = b1[0x01] + b1[0x0E]; b2[0x0E] = (b1[0x01] - b1[0x0E]) * costab[1];
  1104. b2[0x02] = b1[0x02] + b1[0x0D]; b2[0x0D] = (b1[0x02] - b1[0x0D]) * costab[2];
  1105. b2[0x03] = b1[0x03] + b1[0x0C]; b2[0x0C] = (b1[0x03] - b1[0x0C]) * costab[3];
  1106. b2[0x04] = b1[0x04] + b1[0x0B]; b2[0x0B] = (b1[0x04] - b1[0x0B]) * costab[4];
  1107. b2[0x05] = b1[0x05] + b1[0x0A]; b2[0x0A] = (b1[0x05] - b1[0x0A]) * costab[5];
  1108. b2[0x06] = b1[0x06] + b1[0x09]; b2[0x09] = (b1[0x06] - b1[0x09]) * costab[6];
  1109. b2[0x07] = b1[0x07] + b1[0x08]; b2[0x08] = (b1[0x07] - b1[0x08]) * costab[7];
  1110. b2[0x10] = b1[0x10] + b1[0x1F]; b2[0x1F] = (b1[0x1F] - b1[0x10]) * costab[0];
  1111. b2[0x11] = b1[0x11] + b1[0x1E]; b2[0x1E] = (b1[0x1E] - b1[0x11]) * costab[1];
  1112. b2[0x12] = b1[0x12] + b1[0x1D]; b2[0x1D] = (b1[0x1D] - b1[0x12]) * costab[2];
  1113. b2[0x13] = b1[0x13] + b1[0x1C]; b2[0x1C] = (b1[0x1C] - b1[0x13]) * costab[3];
  1114. b2[0x14] = b1[0x14] + b1[0x1B]; b2[0x1B] = (b1[0x1B] - b1[0x14]) * costab[4];
  1115. b2[0x15] = b1[0x15] + b1[0x1A]; b2[0x1A] = (b1[0x1A] - b1[0x15]) * costab[5];
  1116. b2[0x16] = b1[0x16] + b1[0x19]; b2[0x19] = (b1[0x19] - b1[0x16]) * costab[6];
  1117. b2[0x17] = b1[0x17] + b1[0x18]; b2[0x18] = (b1[0x18] - b1[0x17]) * costab[7];
  1118. }
  1119. {
  1120. const float* const costab = constants.cosTables[2];
  1121. b1[0x00] = b2[0x00] + b2[0x07]; b1[0x07] = (b2[0x00] - b2[0x07]) * costab[0];
  1122. b1[0x01] = b2[0x01] + b2[0x06]; b1[0x06] = (b2[0x01] - b2[0x06]) * costab[1];
  1123. b1[0x02] = b2[0x02] + b2[0x05]; b1[0x05] = (b2[0x02] - b2[0x05]) * costab[2];
  1124. b1[0x03] = b2[0x03] + b2[0x04]; b1[0x04] = (b2[0x03] - b2[0x04]) * costab[3];
  1125. b1[0x08] = b2[0x08] + b2[0x0F]; b1[0x0F] = (b2[0x0F] - b2[0x08]) * costab[0];
  1126. b1[0x09] = b2[0x09] + b2[0x0E]; b1[0x0E] = (b2[0x0E] - b2[0x09]) * costab[1];
  1127. b1[0x0A] = b2[0x0A] + b2[0x0D]; b1[0x0D] = (b2[0x0D] - b2[0x0A]) * costab[2];
  1128. b1[0x0B] = b2[0x0B] + b2[0x0C]; b1[0x0C] = (b2[0x0C] - b2[0x0B]) * costab[3];
  1129. b1[0x10] = b2[0x10] + b2[0x17]; b1[0x17] = (b2[0x10] - b2[0x17]) * costab[0];
  1130. b1[0x11] = b2[0x11] + b2[0x16]; b1[0x16] = (b2[0x11] - b2[0x16]) * costab[1];
  1131. b1[0x12] = b2[0x12] + b2[0x15]; b1[0x15] = (b2[0x12] - b2[0x15]) * costab[2];
  1132. b1[0x13] = b2[0x13] + b2[0x14]; b1[0x14] = (b2[0x13] - b2[0x14]) * costab[3];
  1133. b1[0x18] = b2[0x18] + b2[0x1F]; b1[0x1F] = (b2[0x1F] - b2[0x18]) * costab[0];
  1134. b1[0x19] = b2[0x19] + b2[0x1E]; b1[0x1E] = (b2[0x1E] - b2[0x19]) * costab[1];
  1135. b1[0x1A] = b2[0x1A] + b2[0x1D]; b1[0x1D] = (b2[0x1D] - b2[0x1A]) * costab[2];
  1136. b1[0x1B] = b2[0x1B] + b2[0x1C]; b1[0x1C] = (b2[0x1C] - b2[0x1B]) * costab[3];
  1137. }
  1138. {
  1139. const float cos0 = constants.cosTables[3][0];
  1140. const float cos1 = constants.cosTables[3][1];
  1141. b2[0x00] = b1[0x00] + b1[0x03]; b2[0x03] = (b1[0x00] - b1[0x03]) * cos0;
  1142. b2[0x01] = b1[0x01] + b1[0x02]; b2[0x02] = (b1[0x01] - b1[0x02]) * cos1;
  1143. b2[0x04] = b1[0x04] + b1[0x07]; b2[0x07] = (b1[0x07] - b1[0x04]) * cos0;
  1144. b2[0x05] = b1[0x05] + b1[0x06]; b2[0x06] = (b1[0x06] - b1[0x05]) * cos1;
  1145. b2[0x08] = b1[0x08] + b1[0x0B]; b2[0x0B] = (b1[0x08] - b1[0x0B]) * cos0;
  1146. b2[0x09] = b1[0x09] + b1[0x0A]; b2[0x0A] = (b1[0x09] - b1[0x0A]) * cos1;
  1147. b2[0x0C] = b1[0x0C] + b1[0x0F]; b2[0x0F] = (b1[0x0F] - b1[0x0C]) * cos0;
  1148. b2[0x0D] = b1[0x0D] + b1[0x0E]; b2[0x0E] = (b1[0x0E] - b1[0x0D]) * cos1;
  1149. b2[0x10] = b1[0x10] + b1[0x13]; b2[0x13] = (b1[0x10] - b1[0x13]) * cos0;
  1150. b2[0x11] = b1[0x11] + b1[0x12]; b2[0x12] = (b1[0x11] - b1[0x12]) * cos1;
  1151. b2[0x14] = b1[0x14] + b1[0x17]; b2[0x17] = (b1[0x17] - b1[0x14]) * cos0;
  1152. b2[0x15] = b1[0x15] + b1[0x16]; b2[0x16] = (b1[0x16] - b1[0x15]) * cos1;
  1153. b2[0x18] = b1[0x18] + b1[0x1B]; b2[0x1B] = (b1[0x18] - b1[0x1B]) * cos0;
  1154. b2[0x19] = b1[0x19] + b1[0x1A]; b2[0x1A] = (b1[0x19] - b1[0x1A]) * cos1;
  1155. b2[0x1C] = b1[0x1C] + b1[0x1F]; b2[0x1F] = (b1[0x1F] - b1[0x1C]) * cos0;
  1156. b2[0x1D] = b1[0x1D] + b1[0x1E]; b2[0x1E] = (b1[0x1E] - b1[0x1D]) * cos1;
  1157. }
  1158. {
  1159. const float cos0 = constants.cosTables[4][0];
  1160. b1[0x00] = b2[0x00] + b2[0x01]; b1[0x01] = (b2[0x00] - b2[0x01]) * cos0;
  1161. b1[0x02] = b2[0x02] + b2[0x03]; b1[0x03] = (b2[0x03] - b2[0x02]) * cos0; b1[0x02] += b1[0x03];
  1162. b1[0x04] = b2[0x04] + b2[0x05]; b1[0x05] = (b2[0x04] - b2[0x05]) * cos0;
  1163. b1[0x06] = b2[0x06] + b2[0x07]; b1[0x07] = (b2[0x07] - b2[0x06]) * cos0;
  1164. b1[0x06] += b1[0x07]; b1[0x04] += b1[0x06]; b1[0x06] += b1[0x05]; b1[0x05] += b1[0x07];
  1165. b1[0x08] = b2[0x08] + b2[0x09]; b1[0x09] = (b2[0x08] - b2[0x09]) * cos0;
  1166. b1[0x0A] = b2[0x0A] + b2[0x0B]; b1[0x0B] = (b2[0x0B] - b2[0x0A]) * cos0; b1[0x0A] += b1[0x0B];
  1167. b1[0x0C] = b2[0x0C] + b2[0x0D]; b1[0x0D] = (b2[0x0C] - b2[0x0D]) * cos0;
  1168. b1[0x0E] = b2[0x0E] + b2[0x0F]; b1[0x0F] = (b2[0x0F] - b2[0x0E]) * cos0;
  1169. b1[0x0E] += b1[0x0F]; b1[0x0C] += b1[0x0E]; b1[0x0E] += b1[0x0D]; b1[0x0D] += b1[0x0F];
  1170. b1[0x10] = b2[0x10] + b2[0x11]; b1[0x11] = (b2[0x10] - b2[0x11]) * cos0;
  1171. b1[0x12] = b2[0x12] + b2[0x13]; b1[0x13] = (b2[0x13] - b2[0x12]) * cos0; b1[0x12] += b1[0x13];
  1172. b1[0x14] = b2[0x14] + b2[0x15]; b1[0x15] = (b2[0x14] - b2[0x15]) * cos0;
  1173. b1[0x16] = b2[0x16] + b2[0x17]; b1[0x17] = (b2[0x17] - b2[0x16]) * cos0;
  1174. b1[0x16] += b1[0x17]; b1[0x14] += b1[0x16]; b1[0x16] += b1[0x15]; b1[0x15] += b1[0x17];
  1175. b1[0x18] = b2[0x18] + b2[0x19]; b1[0x19] = (b2[0x18] - b2[0x19]) * cos0;
  1176. b1[0x1A] = b2[0x1A] + b2[0x1B]; b1[0x1B] = (b2[0x1B] - b2[0x1A]) * cos0; b1[0x1A] += b1[0x1B];
  1177. b1[0x1C] = b2[0x1C] + b2[0x1D]; b1[0x1D] = (b2[0x1C] - b2[0x1D]) * cos0;
  1178. b1[0x1E] = b2[0x1E] + b2[0x1F]; b1[0x1F] = (b2[0x1F] - b2[0x1E]) * cos0;
  1179. b1[0x1E] += b1[0x1F]; b1[0x1C] += b1[0x1E]; b1[0x1E] += b1[0x1D]; b1[0x1D] += b1[0x1F];
  1180. }
  1181. out0[0x10 * 16] = b1[0x00]; out0[0x10 * 12] = b1[0x04]; out0[0x10 * 8] = b1[0x02]; out0[0x10 * 4] = b1[0x06];
  1182. out0[0] = b1[0x01]; out1[0] = b1[0x01]; out1[0x10 * 4] = b1[0x05]; out1[0x10 * 8] = b1[0x03];
  1183. out1[0x10 * 12] = b1[0x07];
  1184. b1[0x08] += b1[0x0C]; out0[0x10 * 14] = b1[0x08]; b1[0x0C] += b1[0x0a]; out0[0x10 * 10] = b1[0x0C];
  1185. b1[0x0A] += b1[0x0E]; out0[0x10 * 6] = b1[0x0A]; b1[0x0E] += b1[0x09]; out0[0x10 * 2] = b1[0x0E];
  1186. b1[0x09] += b1[0x0D]; out1[0x10 * 2] = b1[0x09]; b1[0x0D] += b1[0x0B]; out1[0x10 * 6] = b1[0x0D];
  1187. b1[0x0B] += b1[0x0F]; out1[0x10 * 10] = b1[0x0B]; out1[0x10 * 14] = b1[0x0F];
  1188. b1[0x18] += b1[0x1C]; out0[0x10 * 15] = b1[0x10] + b1[0x18]; out0[0x10 * 13] = b1[0x18] + b1[0x14];
  1189. b1[0x1C] += b1[0x1a]; out0[0x10 * 11] = b1[0x14] + b1[0x1C]; out0[0x10 * 9] = b1[0x1C] + b1[0x12];
  1190. b1[0x1A] += b1[0x1E]; out0[0x10 * 7] = b1[0x12] + b1[0x1A]; out0[0x10 * 5] = b1[0x1A] + b1[0x16];
  1191. b1[0x1E] += b1[0x19]; out0[0x10 * 3] = b1[0x16] + b1[0x1E]; out0[0x10 * 1] = b1[0x1E] + b1[0x11];
  1192. b1[0x19] += b1[0x1D]; out1[0x10 * 1] = b1[0x11] + b1[0x19]; out1[0x10 * 3] = b1[0x19] + b1[0x15];
  1193. b1[0x1D] += b1[0x1B]; out1[0x10 * 5] = b1[0x15] + b1[0x1D]; out1[0x10 * 7] = b1[0x1D] + b1[0x13];
  1194. b1[0x1B] += b1[0x1F]; out1[0x10 * 9] = b1[0x13] + b1[0x1B]; out1[0x10 * 11] = b1[0x1B] + b1[0x17];
  1195. out1[0x10 * 13] = b1[0x17] + b1[0x1F]; out1[0x10 * 15] = b1[0x1F];
  1196. }
  1197. static void dct64 (float* const a, float* const b, const float* const c) noexcept
  1198. {
  1199. float temp[64];
  1200. dct64 (a, b, temp, temp + 32, c);
  1201. }
  1202. }
  1203. //==============================================================================
  1204. struct MP3Stream
  1205. {
  1206. MP3Stream (InputStream& source)
  1207. : stream (source, 8192),
  1208. numFrames (0), currentFrameIndex (0), vbrHeaderFound (false)
  1209. {
  1210. reset();
  1211. }
  1212. int decodeNextBlock (float* const out0, float* const out1, int& done)
  1213. {
  1214. if (! headerParsed)
  1215. {
  1216. int nextFrameOffset = scanForNextFrameHeader (false);
  1217. if (lastFrameSize == -1 || needToSyncBitStream)
  1218. {
  1219. needToSyncBitStream = false;
  1220. readVBRHeader();
  1221. if (vbrHeaderFound)
  1222. return 1;
  1223. }
  1224. if (nextFrameOffset < 0)
  1225. return -1;
  1226. if (nextFrameOffset > 0)
  1227. {
  1228. int size;
  1229. wasFreeFormat = false;
  1230. needToSyncBitStream = true;
  1231. size = (int) (bufferPointer - (bufferSpace[bufferSpaceIndex] + 512));
  1232. if (size > 2880)
  1233. {
  1234. size = 0;
  1235. bufferPointer = bufferSpace[bufferSpaceIndex] + 512;
  1236. }
  1237. const int toSkip = (size + nextFrameOffset) - 2880;
  1238. if (toSkip > 0)
  1239. {
  1240. stream.skipNextBytes (toSkip);
  1241. nextFrameOffset -= toSkip;
  1242. }
  1243. stream.read (bufferPointer, nextFrameOffset);
  1244. lastFrameSize += nextFrameOffset;
  1245. }
  1246. frame.decodeHeader ((uint32) stream.readIntBigEndian());
  1247. headerParsed = true;
  1248. frameSize = frame.frameSize;
  1249. isFreeFormat = (frameSize == 0);
  1250. sideInfoSize = frame.lsf != 0 ? ((frame.numChannels == 1) ? 9 : 17)
  1251. : ((frame.numChannels == 1) ? 17 : 32);
  1252. if (frame.crc16FollowsHeader)
  1253. sideInfoSize += 2;
  1254. bufferSpaceIndex = 1 - bufferSpaceIndex;
  1255. bufferPointer = bufferSpace[bufferSpaceIndex] + 512;
  1256. bitIndex = 0;
  1257. if (lastFrameSize == -1)
  1258. return 1;
  1259. }
  1260. if (! sideParsed)
  1261. {
  1262. if (frame.layer == 3)
  1263. {
  1264. stream.read (bufferPointer, sideInfoSize);
  1265. if (frame.crc16FollowsHeader)
  1266. getBits (16);
  1267. const int bits = jmax (0, decodeLayer3SideInfo());
  1268. dataSize = (bits + 7) / 8;
  1269. if (! isFreeFormat)
  1270. dataSize = jmin (dataSize, frame.frameSize - sideInfoSize);
  1271. }
  1272. else
  1273. {
  1274. dataSize = frame.frameSize;
  1275. sideInfoSize = 0;
  1276. }
  1277. sideParsed = true;
  1278. }
  1279. int result = 1;
  1280. if (! dataParsed)
  1281. {
  1282. stream.read (bufferPointer, dataSize);
  1283. if (out0 != nullptr)
  1284. {
  1285. if (frame.layer < 3 && frame.crc16FollowsHeader)
  1286. getBits (16);
  1287. switch (frame.layer)
  1288. {
  1289. case 1: decodeLayer1Frame (out0, out1, done); break;
  1290. case 2: decodeLayer2Frame (out0, out1, done); break;
  1291. case 3: decodeLayer3Frame (out0, out1, done); break;
  1292. default: break;
  1293. }
  1294. }
  1295. bufferPointer = bufferSpace[bufferSpaceIndex] + 512 + sideInfoSize + dataSize;
  1296. dataParsed = true;
  1297. result = 0;
  1298. }
  1299. if (isFreeFormat)
  1300. {
  1301. if (wasFreeFormat)
  1302. {
  1303. frameSize = lastFrameSizeNoPadding + frame.padding;
  1304. }
  1305. else
  1306. {
  1307. const int nextFrameOffset = scanForNextFrameHeader (true);
  1308. if (nextFrameOffset < 0)
  1309. return result;
  1310. frameSize = nextFrameOffset + sideInfoSize + dataSize;
  1311. lastFrameSizeNoPadding = frameSize - frame.padding;
  1312. }
  1313. }
  1314. if (result == 0)
  1315. return result;
  1316. int bytes = frameSize - (sideInfoSize + dataSize);
  1317. if (bytes > 0)
  1318. {
  1319. const int toSkip = bytes - 512;
  1320. if (toSkip > 0)
  1321. {
  1322. stream.skipNextBytes (toSkip);
  1323. bytes -= toSkip;
  1324. frameSize -= toSkip;
  1325. }
  1326. stream.read (bufferPointer, bytes);
  1327. bufferPointer += bytes;
  1328. }
  1329. lastFrameSize = frameSize;
  1330. wasFreeFormat = isFreeFormat;
  1331. frameSize = 0;
  1332. headerParsed = sideParsed = dataParsed = false;
  1333. return result;
  1334. }
  1335. bool seek (int frameIndex)
  1336. {
  1337. frameIndex = jmax (0, frameIndex);
  1338. while (frameIndex >= frameStreamPositions.size() * storedStartPosInterval)
  1339. {
  1340. int dummy = 0;
  1341. const int result = decodeNextBlock (nullptr, nullptr, dummy);
  1342. if (result < 0)
  1343. return false;
  1344. if (result > 0)
  1345. break;
  1346. }
  1347. frameIndex = jmin (frameIndex & ~(storedStartPosInterval - 1),
  1348. frameStreamPositions.size() * storedStartPosInterval - 1);
  1349. stream.setPosition (frameStreamPositions.getUnchecked (frameIndex / storedStartPosInterval));
  1350. currentFrameIndex = frameIndex;
  1351. reset();
  1352. return true;
  1353. }
  1354. MP3Frame frame;
  1355. VBRTagData vbrTagData;
  1356. BufferedInputStream stream;
  1357. int numFrames, currentFrameIndex;
  1358. bool vbrHeaderFound;
  1359. private:
  1360. bool headerParsed, sideParsed, dataParsed, needToSyncBitStream;
  1361. bool isFreeFormat, wasFreeFormat;
  1362. int sideInfoSize, dataSize;
  1363. int frameSize, lastFrameSize, lastFrameSizeNoPadding;
  1364. int bufferSpaceIndex;
  1365. Layer3SideInfo sideinfo;
  1366. uint8 bufferSpace[2][2880 + 1024];
  1367. uint8* bufferPointer;
  1368. int bitIndex, synthBo;
  1369. float hybridBlock[2][2][32 * 18];
  1370. int hybridBlockIndex[2];
  1371. float synthBuffers[2][2][0x110];
  1372. float hybridIn[2][32][18];
  1373. float hybridOut[2][18][32];
  1374. void reset() noexcept
  1375. {
  1376. headerParsed = sideParsed = dataParsed = isFreeFormat = wasFreeFormat = false;
  1377. lastFrameSize = -1;
  1378. needToSyncBitStream = true;
  1379. frameSize = sideInfoSize = dataSize = frameSize = bitIndex = 0;
  1380. lastFrameSizeNoPadding = bufferSpaceIndex = 0;
  1381. bufferPointer = bufferSpace[bufferSpaceIndex] + 512;
  1382. synthBo = 1;
  1383. zerostruct (sideinfo);
  1384. zeromem (bufferSpace, sizeof (bufferSpace));
  1385. zeromem (hybridBlock, sizeof (hybridBlock));
  1386. zeromem (hybridBlockIndex, sizeof (hybridBlockIndex));
  1387. zeromem (synthBuffers, sizeof (synthBuffers));
  1388. }
  1389. enum { storedStartPosInterval = 4 };
  1390. Array<int64> frameStreamPositions;
  1391. struct SideInfoLayer1
  1392. {
  1393. uint8 allocation[32][2];
  1394. uint8 scaleFactor[32][2];
  1395. };
  1396. struct SideInfoLayer2
  1397. {
  1398. uint8 allocation[32][2];
  1399. uint8 scaleFactor[32][2][3];
  1400. };
  1401. static bool isValidHeader (const uint32 header, const int oldLayer) noexcept
  1402. {
  1403. const int newLayer = 4 - ((header >> 17) & 3);
  1404. return (header & 0xffe00000) == 0xffe00000
  1405. && newLayer != 4
  1406. && (oldLayer <= 0 || newLayer == oldLayer)
  1407. && ((header >> 12) & 15) != 15
  1408. && ((header >> 10) & 3) != 3
  1409. && (header & 3) != 2;
  1410. }
  1411. bool rollBackBufferPointer (int backstep) noexcept
  1412. {
  1413. if (lastFrameSize < 0 && backstep > 0)
  1414. return false;
  1415. const uint8* oldBuffer = bufferSpace[1 - bufferSpaceIndex] + 512;
  1416. bufferPointer -= backstep;
  1417. if (backstep != 0)
  1418. memcpy (bufferPointer, oldBuffer + lastFrameSize - backstep, (size_t) backstep);
  1419. bitIndex = 0;
  1420. return true;
  1421. }
  1422. uint32 getBits (const int numBits) noexcept
  1423. {
  1424. if (numBits <= 0 || bufferPointer == nullptr)
  1425. return 0;
  1426. const uint32 result = ((((((bufferPointer[0] << 8) | bufferPointer[1]) << 8)
  1427. | bufferPointer[2]) << bitIndex) & 0xffffff) >> (24 - numBits);
  1428. bitIndex += numBits;
  1429. bufferPointer += (bitIndex >> 3);
  1430. bitIndex &= 7;
  1431. return result;
  1432. }
  1433. uint32 getOneBit() noexcept
  1434. {
  1435. const uint8 result = (uint8) (*bufferPointer << bitIndex);
  1436. ++bitIndex;
  1437. bufferPointer += (bitIndex >> 3);
  1438. bitIndex &= 7;
  1439. return result >> 7;
  1440. }
  1441. uint32 getBitsUnchecked (const int numBits) noexcept
  1442. {
  1443. const uint32 result = ((((bufferPointer[0] << 8) | bufferPointer[1]) << bitIndex) & 0xffff) >> (16 - numBits);
  1444. bitIndex += numBits;
  1445. bufferPointer += (bitIndex >> 3);
  1446. bitIndex &= 7;
  1447. return result;
  1448. }
  1449. inline uint8 getBitsUint8 (const int numBits) noexcept { return (uint8) getBitsUnchecked (numBits); }
  1450. inline uint16 getBitsUint16 (const int numBits) noexcept { return (uint16) getBitsUnchecked (numBits); }
  1451. int scanForNextFrameHeader (const bool checkTypeAgainstLastFrame) noexcept
  1452. {
  1453. const int64 oldPos = stream.getPosition();
  1454. int offset = -3;
  1455. uint32 header = 0;
  1456. for (;;)
  1457. {
  1458. if (stream.isExhausted() || stream.getPosition() > oldPos + 32768)
  1459. {
  1460. offset = -1;
  1461. break;
  1462. }
  1463. header = (header << 8) | (uint8) stream.readByte();
  1464. if (offset >= 0 && isValidHeader (header, frame.layer))
  1465. {
  1466. if (! checkTypeAgainstLastFrame)
  1467. break;
  1468. const bool mpeg25 = (header & (1 << 20)) == 0;
  1469. const uint32 lsf = mpeg25 ? 1 : ((header & (1 << 19)) ? 0 : 1);
  1470. const uint32 sampleRateIndex = mpeg25 ? (6 + ((header >> 10) & 3)) : (((header >> 10) & 3) + (lsf * 3));
  1471. const uint32 mode = (header >> 6) & 3;
  1472. const uint32 numChannels = (mode == 3) ? 1 : 2;
  1473. if (numChannels == (uint32) frame.numChannels && lsf == (uint32) frame.lsf
  1474. && mpeg25 == frame.mpeg25 && sampleRateIndex == (uint32) frame.sampleRateIndex)
  1475. break;
  1476. }
  1477. ++offset;
  1478. }
  1479. if (offset >= 0)
  1480. {
  1481. if ((currentFrameIndex & (storedStartPosInterval - 1)) == 0)
  1482. frameStreamPositions.set (currentFrameIndex / storedStartPosInterval, oldPos + offset);
  1483. ++currentFrameIndex;
  1484. }
  1485. stream.setPosition (oldPos);
  1486. return offset;
  1487. }
  1488. void readVBRHeader()
  1489. {
  1490. int64 oldPos = stream.getPosition();
  1491. uint8 xing[194];
  1492. stream.read (xing, sizeof (xing));
  1493. vbrHeaderFound = vbrTagData.read (xing);
  1494. if (vbrHeaderFound)
  1495. {
  1496. numFrames = (int) vbrTagData.frames;
  1497. oldPos += jmax (vbrTagData.headersize, 1);
  1498. }
  1499. stream.setPosition (oldPos);
  1500. }
  1501. void decodeLayer1Frame (float* const pcm0, float* const pcm1, int& samplesDone) noexcept
  1502. {
  1503. float fraction[2][32];
  1504. SideInfoLayer1 si;
  1505. layer1Step1 (si);
  1506. const int single = (frame.numChannels == 1 || frame.single == 3) ? 0 : frame.single;
  1507. if (single >= 0)
  1508. {
  1509. for (int i = 0; i < 12; ++i)
  1510. {
  1511. layer1Step2 (si, fraction);
  1512. synthesise (fraction[single], 0, pcm0, samplesDone);
  1513. }
  1514. }
  1515. else
  1516. {
  1517. for (int i = 0; i < 12; ++i)
  1518. {
  1519. layer1Step2 (si, fraction);
  1520. synthesiseStereo (fraction[0], fraction[1], pcm0, pcm1, samplesDone);
  1521. }
  1522. }
  1523. }
  1524. void decodeLayer2Frame (float* const pcm0, float* const pcm1, int& samplesDone)
  1525. {
  1526. float fraction[2][4][32];
  1527. frame.selectLayer2Table();
  1528. SideInfoLayer2 si;
  1529. layer2Step1 (si);
  1530. const int single = (frame.numChannels == 1 || frame.single == 3) ? 0 : frame.single;
  1531. if (single >= 0)
  1532. {
  1533. for (int i = 0; i < 12; ++i)
  1534. {
  1535. layer2Step2 (si, i >> 2, fraction);
  1536. for (int j = 0; j < 3; ++j)
  1537. synthesise (fraction[single][j], 0, pcm0, samplesDone);
  1538. }
  1539. }
  1540. else
  1541. {
  1542. for (int i = 0; i < 12; ++i)
  1543. {
  1544. layer2Step2 (si, i >> 2, fraction);
  1545. for (int j = 0; j < 3; ++j)
  1546. synthesiseStereo (fraction[0][j], fraction[1][j], pcm0, pcm1, samplesDone);
  1547. }
  1548. }
  1549. }
  1550. void decodeLayer3Frame (float* const pcm0, float* const pcm1, int& samplesDone) noexcept
  1551. {
  1552. if (! rollBackBufferPointer ((int) sideinfo.mainDataStart))
  1553. return;
  1554. const int single = frame.numChannels == 1 ? 0 : frame.single;
  1555. const int numChans = (frame.numChannels == 1 || single >= 0) ? 1 : 2;
  1556. const bool msStereo = (frame.mode == 1) && (frame.modeExt & 2) != 0;
  1557. const bool iStereo = (frame.mode == 1) && (frame.modeExt & 1) != 0;
  1558. const int granules = frame.lsf ? 1 : 2;
  1559. int scaleFactors[2][39];
  1560. for (int gr = 0; gr < granules; ++gr)
  1561. {
  1562. {
  1563. Layer3SideInfo::Info& granule = sideinfo.ch[0].gr[gr];
  1564. const int part2bits = frame.lsf ? getLayer3ScaleFactors2 (scaleFactors[0], granule, 0)
  1565. : getLayer3ScaleFactors1 (scaleFactors[0], granule);
  1566. if (layer3DequantizeSample (hybridIn[0], scaleFactors[0], granule, frame.sampleRateIndex, part2bits))
  1567. return;
  1568. }
  1569. if (frame.numChannels == 2)
  1570. {
  1571. Layer3SideInfo::Info& granule = sideinfo.ch[1].gr[gr];
  1572. const int part2bits = frame.lsf ? getLayer3ScaleFactors2 (scaleFactors[1], granule, iStereo)
  1573. : getLayer3ScaleFactors1 (scaleFactors[1], granule);
  1574. if (layer3DequantizeSample (hybridIn[1], scaleFactors[1], granule, frame.sampleRateIndex, part2bits))
  1575. return;
  1576. if (msStereo)
  1577. {
  1578. for (int i = 0; i < 32 * 18; ++i)
  1579. {
  1580. const float tmp0 = ((const float*) hybridIn[0]) [i];
  1581. const float tmp1 = ((const float*) hybridIn[1]) [i];
  1582. ((float*) hybridIn[1]) [i] = tmp0 - tmp1;
  1583. ((float*) hybridIn[0]) [i] = tmp0 + tmp1;
  1584. }
  1585. }
  1586. if (iStereo)
  1587. granule.doIStereo (hybridIn, scaleFactors[1], frame.sampleRateIndex, msStereo, frame.lsf);
  1588. if (msStereo || iStereo || single == 3)
  1589. {
  1590. if (granule.maxb > sideinfo.ch[0].gr[gr].maxb)
  1591. sideinfo.ch[0].gr[gr].maxb = granule.maxb;
  1592. else
  1593. granule.maxb = sideinfo.ch[0].gr[gr].maxb;
  1594. }
  1595. switch (single)
  1596. {
  1597. case 3:
  1598. {
  1599. float* in0 = (float*) hybridIn[0];
  1600. const float* in1 = (const float*) hybridIn[1];
  1601. for (int i = 0; i < (int) (18 * granule.maxb); ++i, ++in0)
  1602. *in0 = (*in0 + *in1++);
  1603. }
  1604. break;
  1605. case 1:
  1606. {
  1607. float* in0 = (float*) hybridIn[0];
  1608. const float* in1 = (const float*) hybridIn[1];
  1609. for (int i = 0; i < (int) (18 * granule.maxb); ++i)
  1610. *in0++ = *in1++;
  1611. }
  1612. break;
  1613. }
  1614. }
  1615. for (int ch = 0; ch < numChans; ++ch)
  1616. {
  1617. const Layer3SideInfo::Info& granule = sideinfo.ch[ch].gr[gr];
  1618. granule.doAntialias (hybridIn[ch]);
  1619. layer3Hybrid (hybridIn[ch], hybridOut[ch], ch, granule);
  1620. }
  1621. for (int ss = 0; ss < 18; ++ss)
  1622. {
  1623. if (single >= 0)
  1624. synthesise (hybridOut[0][ss], 0, pcm0, samplesDone);
  1625. else
  1626. synthesiseStereo (hybridOut[0][ss], hybridOut[1][ss], pcm0, pcm1, samplesDone);
  1627. }
  1628. }
  1629. }
  1630. int decodeLayer3SideInfo() noexcept
  1631. {
  1632. const int numChannels = frame.numChannels;
  1633. const int sampleRate = frame.sampleRateIndex;
  1634. const int single = (numChannels == 1) ? 0 : frame.single;
  1635. const bool msStereo = (frame.mode == 1) && (frame.modeExt & 2) != 0;
  1636. const int granules = frame.lsf ? 1 : 2;
  1637. if (frame.lsf == 0)
  1638. getLayer3SideInfo1 (numChannels, msStereo, sampleRate, single);
  1639. else
  1640. getLayer3SideInfo2 (numChannels, msStereo, sampleRate, single);
  1641. int databits = 0;
  1642. for (int gr = 0; gr < granules; ++gr)
  1643. for (int ch = 0; ch < numChannels; ++ch)
  1644. databits += sideinfo.ch[ch].gr[gr].part2_3Length;
  1645. return databits - 8 * (int) sideinfo.mainDataStart;
  1646. }
  1647. void layer1Step1 (SideInfoLayer1& si) noexcept
  1648. {
  1649. zerostruct (si);
  1650. int i, jsbound = (frame.mode == 1) ? (frame.modeExt << 2) + 4 : 32;
  1651. if (frame.numChannels == 2)
  1652. {
  1653. for (i = 0; i < jsbound; ++i)
  1654. {
  1655. si.allocation[i][0] = getBitsUint8 (4);
  1656. si.allocation[i][1] = getBitsUint8 (4);
  1657. }
  1658. for (i = jsbound; i < 32; ++i)
  1659. si.allocation[i][0] = si.allocation[i][1] = getBitsUint8 (4);
  1660. for (i = 0; i < 32; ++i)
  1661. {
  1662. si.scaleFactor[i][0] = si.allocation[i][0] ? getBitsUint8 (6) : 0;
  1663. si.scaleFactor[i][1] = si.allocation[i][1] ? getBitsUint8 (6) : 0;
  1664. }
  1665. }
  1666. else
  1667. {
  1668. for (i = 0; i < 32; ++i)
  1669. si.allocation[i][0] = getBitsUint8 (4);
  1670. for (i = 0; i < 32; ++i)
  1671. si.scaleFactor[i][0] = si.allocation[i][0] ? getBitsUint8 (6) : 0;
  1672. }
  1673. }
  1674. void layer1Step2 (SideInfoLayer1& si, float fraction[2][32]) noexcept
  1675. {
  1676. if (frame.numChannels == 2)
  1677. {
  1678. int i, jsbound = (frame.mode == 1) ? (frame.modeExt << 2) + 4 : 32;
  1679. for (i = 0; i < jsbound; ++i)
  1680. {
  1681. const uint8 n0 = si.allocation[i][0];
  1682. const uint8 n1 = si.allocation[i][1];
  1683. fraction[0][i] = n0 > 0 ? (float) (((-1 << n0) + getBitsUint16 (n0 + 1) + 1) * constants.muls[n0 + 1][si.scaleFactor[i][0]]) : 0;
  1684. fraction[1][i] = n1 > 0 ? (float) (((-1 << n1) + getBitsUint16 (n1 + 1) + 1) * constants.muls[n1 + 1][si.scaleFactor[i][1]]) : 0;
  1685. }
  1686. for (i = jsbound; i < 32; ++i)
  1687. {
  1688. const uint8 n = si.allocation[i][0];
  1689. if (n > 0)
  1690. {
  1691. const uint32 w = ((uint32) (-1 << n) + getBitsUint16 (n + 1) + 1);
  1692. fraction[0][i] = (float) (w * constants.muls[n + 1][si.scaleFactor[i][0]]);
  1693. fraction[1][i] = (float) (w * constants.muls[n + 1][si.scaleFactor[i][1]]);
  1694. }
  1695. else
  1696. fraction[0][i] = fraction[1][i] = 0;
  1697. }
  1698. }
  1699. else
  1700. {
  1701. for (int i = 0; i < 32; ++i)
  1702. {
  1703. const uint8 n = si.allocation[i][0];
  1704. const uint8 j = si.scaleFactor[i][0];
  1705. if (n > 0)
  1706. fraction[0][i] = (float) (((-1 << n) + getBitsUint16 (n + 1) + 1) * constants.muls[n + 1][j]);
  1707. else
  1708. fraction[0][i] = 0;
  1709. }
  1710. }
  1711. }
  1712. void layer2Step1 (SideInfoLayer2& si) noexcept
  1713. {
  1714. zerostruct (si);
  1715. const int sblimit = frame.layer2SubBandLimit;
  1716. const int jsbound = (frame.mode == 1) ? (frame.modeExt << 2) + 4 : frame.layer2SubBandLimit;
  1717. const AllocationTable* allocTable = frame.allocationTable;
  1718. uint8 scfsi[32][2];
  1719. if (frame.numChannels == 2)
  1720. {
  1721. for (int i = 0; i < jsbound; ++i)
  1722. {
  1723. const int16 step = allocTable->bits;
  1724. allocTable += (1 << step);
  1725. si.allocation[i][0] = getBitsUint8 (step);
  1726. si.allocation[i][1] = getBitsUint8 (step);
  1727. }
  1728. for (int i = jsbound; i < sblimit; ++i)
  1729. {
  1730. const int16 step = allocTable->bits;
  1731. const uint8 b0 = getBitsUint8 (step);
  1732. allocTable += (1 << step);
  1733. si.allocation[i][0] = b0;
  1734. si.allocation[i][1] = b0;
  1735. }
  1736. for (int i = 0; i < sblimit; ++i)
  1737. {
  1738. scfsi[i][0] = si.allocation[i][0] ? getBitsUint8 (2) : 0;
  1739. scfsi[i][1] = si.allocation[i][1] ? getBitsUint8 (2) : 0;
  1740. }
  1741. }
  1742. else
  1743. {
  1744. for (int i = 0; i < sblimit; ++i)
  1745. {
  1746. const int16 step = allocTable->bits;
  1747. allocTable += (1 << step);
  1748. si.allocation[i][0] = getBitsUint8 (step);
  1749. }
  1750. for (int i = 0; i < sblimit; ++i)
  1751. scfsi[i][0] = si.allocation[i][0] ? getBitsUint8 (2) : 0;
  1752. }
  1753. for (int i = 0; i < sblimit; ++i)
  1754. {
  1755. for (int ch = 0; ch < frame.numChannels; ++ch)
  1756. {
  1757. uint8 s0 = 0, s1 = 0, s2 = 0;
  1758. if (si.allocation[i][ch])
  1759. {
  1760. switch (scfsi[i][ch])
  1761. {
  1762. case 0:
  1763. s0 = getBitsUint8 (6);
  1764. s1 = getBitsUint8 (6);
  1765. s2 = getBitsUint8 (6);
  1766. break;
  1767. case 1:
  1768. s1 = s0 = getBitsUint8 (6);
  1769. s2 = getBitsUint8 (6);
  1770. break;
  1771. case 2:
  1772. s2 = s1 = s0 = getBitsUint8 (6);
  1773. break;
  1774. case 3:
  1775. s0 = getBitsUint8 (6);
  1776. s2 = s1 = getBitsUint8 (6);
  1777. break;
  1778. default:
  1779. break;
  1780. }
  1781. }
  1782. si.scaleFactor[i][ch][0] = s0;
  1783. si.scaleFactor[i][ch][1] = s1;
  1784. si.scaleFactor[i][ch][2] = s2;
  1785. }
  1786. }
  1787. }
  1788. void layer2Step2 (SideInfoLayer2& si, const int gr, float fraction[2][4][32]) noexcept
  1789. {
  1790. const AllocationTable* allocTable = frame.allocationTable;
  1791. const int jsbound = (frame.mode == 1) ? (frame.modeExt << 2) + 4 : frame.layer2SubBandLimit;
  1792. for (int i = 0; i < jsbound; ++i)
  1793. {
  1794. const int16 step = allocTable->bits;
  1795. for (int ch = 0; ch < frame.numChannels; ++ch)
  1796. {
  1797. const uint8 ba = si.allocation[i][ch];
  1798. if (ba != 0)
  1799. {
  1800. const uint8 x1 = jmin ((uint8) 63, si.scaleFactor[i][ch][gr]);
  1801. const AllocationTable* const alloc2 = allocTable + ba;
  1802. const int16 k = jmin ((int16) 16, alloc2->bits);
  1803. const int16 d1 = alloc2->d;
  1804. if (d1 < 0)
  1805. {
  1806. const double cm = constants.muls[k][x1];
  1807. fraction[ch][0][i] = (float) (((int) getBits (k) + d1) * cm);
  1808. fraction[ch][1][i] = (float) (((int) getBits (k) + d1) * cm);
  1809. fraction[ch][2][i] = (float) (((int) getBits (k) + d1) * cm);
  1810. }
  1811. else
  1812. {
  1813. const uint8* const tab = constants.getGroupTable (d1, getBits (k));
  1814. fraction[ch][0][i] = (float) constants.muls[tab[0]][x1];
  1815. fraction[ch][1][i] = (float) constants.muls[tab[1]][x1];
  1816. fraction[ch][2][i] = (float) constants.muls[tab[2]][x1];
  1817. }
  1818. }
  1819. else
  1820. {
  1821. fraction[ch][0][i] = fraction[ch][1][i] = fraction[ch][2][i] = 0;
  1822. }
  1823. }
  1824. allocTable += (1 << step);
  1825. }
  1826. for (int i = jsbound; i < frame.layer2SubBandLimit; ++i)
  1827. {
  1828. const int16 step = allocTable->bits;
  1829. const uint8 ba = si.allocation[i][0];
  1830. if (ba != 0)
  1831. {
  1832. const AllocationTable* const alloc2 = allocTable + ba;
  1833. int16 k = alloc2->bits;
  1834. int16 d1 = alloc2->d;
  1835. k = (k <= 16) ? k : 16;
  1836. if (d1 < 0)
  1837. {
  1838. const int v0 = (int) getBits (k);
  1839. const int v1 = (int) getBits (k);
  1840. const int v2 = (int) getBits (k);
  1841. for (int ch = 0; ch < frame.numChannels; ++ch)
  1842. {
  1843. const uint8 x1 = jmin ((uint8) 63, si.scaleFactor[i][ch][gr]);
  1844. const double cm = constants.muls[k][x1];
  1845. fraction[ch][0][i] = (float) ((v0 + d1) * cm);
  1846. fraction[ch][1][i] = (float) ((v1 + d1) * cm);
  1847. fraction[ch][2][i] = (float) ((v2 + d1) * cm);
  1848. }
  1849. }
  1850. else
  1851. {
  1852. const uint8* const tab = constants.getGroupTable (d1, getBits (k));
  1853. const uint8 k0 = tab[0];
  1854. const uint8 k1 = tab[1];
  1855. const uint8 k2 = tab[2];
  1856. for (int ch = 0; ch < frame.numChannels; ++ch)
  1857. {
  1858. const uint8 x1 = jmin ((uint8) 63, si.scaleFactor[i][ch][gr]);
  1859. fraction[ch][0][i] = (float) constants.muls[k0][x1];
  1860. fraction[ch][1][i] = (float) constants.muls[k1][x1];
  1861. fraction[ch][2][i] = (float) constants.muls[k2][x1];
  1862. }
  1863. }
  1864. }
  1865. else
  1866. {
  1867. fraction[0][0][i] = fraction[0][1][i] = fraction[0][2][i] = 0;
  1868. fraction[1][0][i] = fraction[1][1][i] = fraction[1][2][i] = 0;
  1869. }
  1870. allocTable += (1 << step);
  1871. }
  1872. for (int ch = 0; ch < frame.numChannels; ++ch)
  1873. for (int i = frame.layer2SubBandLimit; i < 32; ++i)
  1874. fraction[ch][0][i] = fraction[ch][1][i] = fraction[ch][2][i] = 0;
  1875. }
  1876. void getLayer3SideInfo1 (const int stereo, const bool msStereo, const int sampleRate, const int single) noexcept
  1877. {
  1878. const int powdiff = (single == 3) ? 4 : 0;
  1879. sideinfo.mainDataStart = getBits (9);
  1880. sideinfo.privateBits = getBitsUnchecked (stereo == 1 ? 5 : 3);
  1881. for (int ch = 0; ch < stereo; ++ch)
  1882. {
  1883. sideinfo.ch[ch].gr[0].scfsi = -1;
  1884. sideinfo.ch[ch].gr[1].scfsi = (int) getBitsUnchecked (4);
  1885. }
  1886. for (int gr = 0; gr < 2; ++gr)
  1887. {
  1888. for (int ch = 0; ch < stereo; ++ch)
  1889. {
  1890. Layer3SideInfo::Info& granule = sideinfo.ch[ch].gr[gr];
  1891. granule.part2_3Length = getBits (12);
  1892. granule.bigValues = jmin (288u, getBitsUnchecked (9));
  1893. const int qss = (int) getBitsUnchecked (8);
  1894. granule.pow2gain = constants.powToGains + 256 - qss + powdiff;
  1895. if (msStereo)
  1896. granule.pow2gain += 2;
  1897. granule.scaleFactorCompression = getBitsUnchecked (4);
  1898. if (getOneBit())
  1899. {
  1900. granule.blockType = getBitsUnchecked (2);
  1901. granule.mixedBlockFlag = getOneBit();
  1902. granule.tableSelect[0] = getBitsUnchecked (5);
  1903. granule.tableSelect[1] = getBitsUnchecked (5);
  1904. granule.tableSelect[2] = 0;
  1905. for (int i = 0; i < 3; ++i)
  1906. {
  1907. const uint32 sbg = (getBitsUnchecked (3) << 3);
  1908. granule.fullGain[i] = granule.pow2gain + sbg;
  1909. }
  1910. granule.region1Start = 36 >> 1;
  1911. granule.region2Start = 576 >> 1;
  1912. }
  1913. else
  1914. {
  1915. for (int i = 0; i < 3; ++i)
  1916. granule.tableSelect[i] = getBitsUnchecked (5);
  1917. const int r0c = (int) getBitsUnchecked (4);
  1918. const int r1c = (int) getBitsUnchecked (3);
  1919. const int region0index = jmin (22, r0c + 1);
  1920. const int region1index = jmin (22, r0c + 1 + r1c + 1);
  1921. granule.region1Start = (uint32) (bandInfo[sampleRate].longIndex[region0index] >> 1);
  1922. granule.region2Start = (uint32) (bandInfo[sampleRate].longIndex[region1index] >> 1);
  1923. granule.blockType = 0;
  1924. granule.mixedBlockFlag = 0;
  1925. }
  1926. granule.preflag = getOneBit();
  1927. granule.scaleFactorScale = getOneBit();
  1928. granule.count1TableSelect = getOneBit();
  1929. }
  1930. }
  1931. }
  1932. void getLayer3SideInfo2 (const int stereo, const bool msStereo, const int sampleRate, const int single) noexcept
  1933. {
  1934. const int powdiff = (single == 3) ? 4 : 0;
  1935. sideinfo.mainDataStart = getBits (8);
  1936. sideinfo.privateBits = stereo == 1 ? getOneBit() : getBitsUnchecked (2);
  1937. for (int ch = 0; ch < stereo; ++ch)
  1938. {
  1939. Layer3SideInfo::Info& granule = sideinfo.ch[ch].gr[0];
  1940. granule.part2_3Length = getBits (12);
  1941. granule.bigValues = jmin (288u, getBitsUnchecked (9));
  1942. const uint32 qss = getBitsUnchecked (8);
  1943. granule.pow2gain = constants.powToGains + 256 - qss + powdiff;
  1944. if (msStereo)
  1945. granule.pow2gain += 2;
  1946. granule.scaleFactorCompression = getBits (9);
  1947. if (getOneBit())
  1948. {
  1949. granule.blockType = getBitsUnchecked (2);
  1950. granule.mixedBlockFlag = getOneBit();
  1951. granule.tableSelect[0] = getBitsUnchecked (5);
  1952. granule.tableSelect[1] = getBitsUnchecked (5);
  1953. granule.tableSelect[2] = 0;
  1954. for (int i = 0; i < 3; ++i)
  1955. {
  1956. const uint32 sbg = (getBitsUnchecked (3) << 3);
  1957. granule.fullGain[i] = granule.pow2gain + sbg;
  1958. }
  1959. if (granule.blockType == 0)
  1960. {}
  1961. if (granule.blockType == 2)
  1962. granule.region1Start = sampleRate == 8 ? 36 : (36 >> 1);
  1963. else
  1964. granule.region1Start = sampleRate == 8 ? (108 >> 1) : (54 >> 1);
  1965. granule.region2Start = 576 >> 1;
  1966. }
  1967. else
  1968. {
  1969. for (int i = 0; i < 3; ++i)
  1970. granule.tableSelect[i] = getBitsUnchecked (5);
  1971. const int r0c = (int) getBitsUnchecked (4);
  1972. const int r1c = (int) getBitsUnchecked (3);
  1973. const int region0index = jmin (22, r0c + 1);
  1974. const int region1index = jmin (22, r0c + 1 + r1c + 1);
  1975. granule.region1Start = (uint32) (bandInfo[sampleRate].longIndex[region0index] >> 1);
  1976. granule.region2Start = (uint32) (bandInfo[sampleRate].longIndex[region1index] >> 1);
  1977. granule.blockType = 0;
  1978. granule.mixedBlockFlag = 0;
  1979. }
  1980. granule.scaleFactorScale = getOneBit();
  1981. granule.count1TableSelect = getOneBit();
  1982. }
  1983. }
  1984. int getLayer3ScaleFactors1 (int* scf, const Layer3SideInfo::Info& granule) noexcept
  1985. {
  1986. static const uint8 lengths[2][16] =
  1987. {
  1988. { 0, 0, 0, 0, 3, 1, 1, 1, 2, 2, 2, 3, 3, 3, 4, 4 },
  1989. { 0, 1, 2, 3, 0, 1, 2, 3, 1, 2, 3, 1, 2, 3, 2, 3 }
  1990. };
  1991. int numBits;
  1992. const int num0 = lengths[0][granule.scaleFactorCompression];
  1993. const int num1 = lengths[1][granule.scaleFactorCompression];
  1994. if (granule.blockType == 2)
  1995. {
  1996. int i = 18;
  1997. numBits = (num0 + num1) * 18;
  1998. if (granule.mixedBlockFlag)
  1999. {
  2000. for (int j = 8; --j >= 0;) *scf++ = (int) getBitsUnchecked (num0);
  2001. numBits -= num0;
  2002. i = 9;
  2003. }
  2004. for (; --i >= 0;) *scf++ = (int) getBitsUnchecked (num0);
  2005. for (i = 18; --i >= 0;) *scf++ = (int) getBitsUnchecked (num1);
  2006. *scf++ = 0;
  2007. *scf++ = 0;
  2008. *scf++ = 0;
  2009. }
  2010. else
  2011. {
  2012. const int scfsi = granule.scfsi;
  2013. if (scfsi < 0)
  2014. {
  2015. for (int i = 11; --i >= 0;) *scf++ = (int) getBitsUnchecked (num0);
  2016. for (int j = 10; --j >= 0;) *scf++ = (int) getBitsUnchecked (num1);
  2017. numBits = (num0 + num1) * 10 + num0;
  2018. }
  2019. else
  2020. {
  2021. numBits = 0;
  2022. if ((scfsi & 8) == 0)
  2023. {
  2024. for (int i = 6; --i >= 0;) *scf++ = (int) getBitsUnchecked (num0);
  2025. numBits += num0 * 6;
  2026. }
  2027. else
  2028. scf += 6;
  2029. if ((scfsi & 4) == 0)
  2030. {
  2031. for (int i = 5; --i >= 0;) *scf++ = (int) getBitsUnchecked (num0);
  2032. numBits += num0 * 5;
  2033. }
  2034. else
  2035. scf += 5;
  2036. if ((scfsi & 2) == 0)
  2037. {
  2038. for (int i = 5; --i >= 0;) *scf++ = (int) getBitsUnchecked (num1);
  2039. numBits += num1 * 5;
  2040. }
  2041. else
  2042. scf += 5;
  2043. if ((scfsi & 1) == 0)
  2044. {
  2045. for (int i = 5; --i >= 0;) *scf++ = (int) getBitsUnchecked (num1);
  2046. numBits += num1 * 5;
  2047. }
  2048. else
  2049. scf += 5;
  2050. }
  2051. *scf = 0;
  2052. }
  2053. return numBits;
  2054. }
  2055. int getLayer3ScaleFactors2 (int* scf, Layer3SideInfo::Info& granule, const bool iStereo) noexcept
  2056. {
  2057. static const uint8 scaleTable[3][6][4] =
  2058. {
  2059. { { 6, 5, 5, 5 }, { 6, 5, 7, 3 }, { 11, 10, 0, 0 }, { 7, 7, 7, 0 }, { 6, 6, 6, 3 }, { 8, 8, 5, 0 } },
  2060. { { 9, 9, 9, 9 }, { 9, 9, 12, 6 }, { 18, 18, 0, 0 }, { 12, 12, 12, 0 }, { 12, 9, 9, 6 }, { 15, 12, 9, 0 } },
  2061. { { 6, 9, 9, 9 }, { 6, 9, 12, 6 }, { 15, 18, 0, 0 }, { 6, 15, 12, 0 }, { 6, 12, 9, 6 }, { 6, 18, 9, 0 } }
  2062. };
  2063. uint32 len = iStereo ? constants.iLength2 [granule.scaleFactorCompression >> 1]
  2064. : constants.nLength2 [granule.scaleFactorCompression];
  2065. granule.preflag = (len >> 15) & 1;
  2066. int n = 0;
  2067. if (granule.blockType == 2)
  2068. {
  2069. ++n;
  2070. if (granule.mixedBlockFlag)
  2071. ++n;
  2072. }
  2073. const uint8* const data = scaleTable[n][(len >> 12) & 7];
  2074. int i, numBits = 0;
  2075. for (i = 0; i < 4; ++i)
  2076. {
  2077. int num = len & 7;
  2078. len >>= 3;
  2079. if (num)
  2080. {
  2081. for (int j = 0; j < (int) (data[i]); ++j)
  2082. *scf++ = (int) getBitsUnchecked (num);
  2083. numBits += data[i] * num;
  2084. }
  2085. else
  2086. {
  2087. for (int j = 0; j < (int) (data[i]); ++j)
  2088. *scf++ = 0;
  2089. }
  2090. }
  2091. n = (n << 1) + 1;
  2092. for (i = 0; i < n; ++i)
  2093. *scf++ = 0;
  2094. return numBits;
  2095. }
  2096. bool layer3DequantizeSample (float xr[32][18], int* scf, Layer3SideInfo::Info& granule, int sampleRate, int part2bits) noexcept
  2097. {
  2098. const uint32 shift = 1 + granule.scaleFactorScale;
  2099. float* xrpnt = (float*) xr;
  2100. int part2remain = (int) granule.part2_3Length - part2bits;
  2101. zeromem (xrpnt, sizeof (float) * (size_t) (&xr[32][0] - xrpnt));
  2102. const int bv = (int) granule.bigValues;
  2103. const int region1 = (int) granule.region1Start;
  2104. const int region2 = (int) granule.region2Start;
  2105. int l3 = ((576 >> 1) - bv) >> 1;
  2106. int l[3];
  2107. if (bv <= region1)
  2108. {
  2109. l[0] = bv;
  2110. l[1] = 0;
  2111. l[2] = 0;
  2112. }
  2113. else
  2114. {
  2115. l[0] = region1;
  2116. if (bv <= region2)
  2117. {
  2118. l[1] = bv - l[0];
  2119. l[2] = 0;
  2120. }
  2121. else
  2122. {
  2123. l[1] = region2 - l[0];
  2124. l[2] = bv - region2;
  2125. }
  2126. }
  2127. for (int i = 0; i < 3; ++i)
  2128. if (l[i] < 0)
  2129. l[i] = 0;
  2130. if (granule.blockType == 2)
  2131. {
  2132. int max[4];
  2133. int step = 0, lwin = 0, cb = 0, mc = 0;
  2134. float v = 0;
  2135. int* map;
  2136. int* mapEnd;
  2137. if (granule.mixedBlockFlag)
  2138. {
  2139. max[3] = -1;
  2140. max[0] = max[1] = max[2] = 2;
  2141. map = constants.map [sampleRate][0];
  2142. mapEnd = constants.mapEnd [sampleRate][0];
  2143. }
  2144. else
  2145. {
  2146. max[0] = max[1] = max[2] = max[3] = -1;
  2147. map = constants.map [sampleRate][1];
  2148. mapEnd = constants.mapEnd [sampleRate][1];
  2149. }
  2150. for (int i = 0; i < 2; ++i)
  2151. {
  2152. const BitsToTableMap* h = huffmanTables1 + granule.tableSelect[i];
  2153. for (int lp = l[i]; lp != 0; --lp, --mc)
  2154. {
  2155. int x, y;
  2156. if (mc == 0)
  2157. {
  2158. mc = *map++;
  2159. xrpnt = ((float*) xr) + (*map++);
  2160. lwin = *map++;
  2161. cb = *map++;
  2162. if (lwin == 3)
  2163. {
  2164. v = granule.pow2gain[ (*scf++) << shift];
  2165. step = 1;
  2166. }
  2167. else
  2168. {
  2169. v = granule.fullGain[lwin][ (*scf++) << shift];
  2170. step = 3;
  2171. }
  2172. }
  2173. const int16* val = h->table;
  2174. while ((y = *val++) < 0)
  2175. {
  2176. if (getOneBit())
  2177. val -= y;
  2178. --part2remain;
  2179. }
  2180. x = y >> 4;
  2181. y &= 15;
  2182. if (x == 15)
  2183. {
  2184. max[lwin] = cb;
  2185. part2remain -= h->bits + 1;
  2186. x += getBits ((int) h->bits);
  2187. *xrpnt = constants.nToThe4Over3[x] * (getOneBit() ? -v : v);
  2188. }
  2189. else if (x)
  2190. {
  2191. max[lwin] = cb;
  2192. *xrpnt = constants.nToThe4Over3[x] * (getOneBit() ? -v : v);
  2193. --part2remain;
  2194. }
  2195. else
  2196. *xrpnt = 0;
  2197. xrpnt += step;
  2198. if (y == 15)
  2199. {
  2200. max[lwin] = cb;
  2201. part2remain -= h->bits + 1;
  2202. y += getBits ((int) h->bits);
  2203. *xrpnt = constants.nToThe4Over3[y] * (getOneBit() ? -v : v);
  2204. }
  2205. else if (y)
  2206. {
  2207. max[lwin] = cb;
  2208. *xrpnt = constants.nToThe4Over3[y] * (getOneBit() ? -v : v);
  2209. --part2remain;
  2210. }
  2211. else
  2212. *xrpnt = 0;
  2213. xrpnt += step;
  2214. }
  2215. }
  2216. for (; l3 && (part2remain > 0); --l3)
  2217. {
  2218. const BitsToTableMap* h = huffmanTables2 + granule.count1TableSelect;
  2219. const int16* val = h->table;
  2220. int16 a;
  2221. while ((a = *val++) < 0)
  2222. {
  2223. if (part2remain <= 0)
  2224. {
  2225. a = 0;
  2226. break;
  2227. }
  2228. --part2remain;
  2229. if (getOneBit())
  2230. val -= a;
  2231. }
  2232. for (int i = 0; i < 4; ++i)
  2233. {
  2234. if ((i & 1) == 0)
  2235. {
  2236. if (mc == 0)
  2237. {
  2238. mc = *map++;
  2239. xrpnt = ((float*) xr) + (*map++);
  2240. lwin = *map++;
  2241. cb = *map++;
  2242. if (lwin == 3)
  2243. {
  2244. v = granule.pow2gain[ (*scf++) << shift];
  2245. step = 1;
  2246. }
  2247. else
  2248. {
  2249. v = granule.fullGain[lwin][ (*scf++) << shift];
  2250. step = 3;
  2251. }
  2252. }
  2253. --mc;
  2254. }
  2255. if ((a & (8 >> i)))
  2256. {
  2257. max[lwin] = cb;
  2258. if (part2remain == 0)
  2259. break;
  2260. --part2remain;
  2261. *xrpnt = getOneBit() ? -v : v;
  2262. }
  2263. else
  2264. *xrpnt = 0;
  2265. xrpnt += step;
  2266. }
  2267. }
  2268. while (map < mapEnd)
  2269. {
  2270. if (mc == 0)
  2271. {
  2272. mc = *map++;
  2273. xrpnt = ((float*) xr) + *map++;
  2274. step = (*map++ == 3) ? 1 : 3;
  2275. ++map;
  2276. }
  2277. --mc;
  2278. *xrpnt = 0; xrpnt += step;
  2279. *xrpnt = 0; xrpnt += step;
  2280. }
  2281. granule.maxBand[0] = (uint32) (max[0] + 1);
  2282. granule.maxBand[1] = (uint32) (max[1] + 1);
  2283. granule.maxBand[2] = (uint32) (max[2] + 1);
  2284. granule.maxBandl = (uint32) (max[3] + 1);
  2285. const int rmax = jmax (max[0], max[1], max[3]) + 1;
  2286. granule.maxb = rmax ? (uint32) constants.shortLimit[sampleRate][rmax]
  2287. : (uint32) constants.longLimit[sampleRate][max[3] + 1];
  2288. }
  2289. else
  2290. {
  2291. static const int pretab1[] = { 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 2, 2, 3, 3, 3, 2, 0 };
  2292. static const int pretab2[] = { 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 };
  2293. const int* pretab = (const int*) (granule.preflag ? pretab1 : pretab2);
  2294. int max = -1, cb = 0, mc = 0;
  2295. int* map = constants.map [sampleRate][2];
  2296. float v = 0;
  2297. for (int i = 0; i < 3; ++i)
  2298. {
  2299. const BitsToTableMap* h = huffmanTables1 + granule.tableSelect[i];
  2300. for (int lp = l[i]; lp != 0; --lp, --mc)
  2301. {
  2302. if (mc == 0)
  2303. {
  2304. mc = *map++;
  2305. v = granule.pow2gain [((*scf++) + (*pretab++)) << shift];
  2306. cb = *map++;
  2307. }
  2308. const int16* val = h->table;
  2309. int y;
  2310. while ((y = *val++) < 0)
  2311. {
  2312. if (getOneBit()) val -= y;
  2313. --part2remain;
  2314. }
  2315. int x = y >> 4;
  2316. y &= 15;
  2317. if (x == 15)
  2318. {
  2319. max = cb;
  2320. part2remain -= h->bits + 1;
  2321. x += getBits ((int) h->bits);
  2322. *xrpnt++ = constants.nToThe4Over3[x] * (getOneBit() ? -v : v);
  2323. }
  2324. else if (x)
  2325. {
  2326. max = cb;
  2327. *xrpnt++ = constants.nToThe4Over3[x] * (getOneBit() ? -v : v);
  2328. --part2remain;
  2329. }
  2330. else
  2331. *xrpnt++ = 0;
  2332. if (y == 15)
  2333. {
  2334. max = cb;
  2335. part2remain -= h->bits + 1;
  2336. y += getBits ((int) h->bits);
  2337. *xrpnt++ = constants.nToThe4Over3[y] * (getOneBit() ? -v : v);
  2338. }
  2339. else if (y)
  2340. {
  2341. max = cb;
  2342. *xrpnt++ = constants.nToThe4Over3[y] * (getOneBit() ? -v : v);
  2343. --part2remain;
  2344. }
  2345. else
  2346. *xrpnt++ = 0;
  2347. }
  2348. }
  2349. for (; l3 && part2remain > 0; --l3)
  2350. {
  2351. const BitsToTableMap* h = huffmanTables2 + granule.count1TableSelect;
  2352. const int16* values = h->table;
  2353. int16 a;
  2354. while ((a = *values++) < 0)
  2355. {
  2356. if (part2remain <= 0)
  2357. {
  2358. a = 0;
  2359. break;
  2360. }
  2361. --part2remain;
  2362. if (getOneBit())
  2363. values -= a;
  2364. }
  2365. for (int i = 0; i < 4; ++i)
  2366. {
  2367. if ((i & 1) == 0)
  2368. {
  2369. if (mc == 0)
  2370. {
  2371. mc = *map++;
  2372. cb = *map++;
  2373. v = granule.pow2gain [((*scf++) + (*pretab++)) << shift];
  2374. }
  2375. --mc;
  2376. }
  2377. if ((a & (0x8 >> i)))
  2378. {
  2379. max = cb;
  2380. if (part2remain <= 0)
  2381. break;
  2382. --part2remain;
  2383. *xrpnt++ = getOneBit() ? -v : v;
  2384. }
  2385. else
  2386. *xrpnt++ = 0;
  2387. }
  2388. }
  2389. zeromem (xrpnt, sizeof (float) * (size_t) (&xr[32][0] - xrpnt));
  2390. granule.maxBandl = (uint32) (max + 1);
  2391. granule.maxb = (uint32) constants.longLimit[sampleRate][granule.maxBandl];
  2392. }
  2393. while (part2remain > 16)
  2394. {
  2395. getBits (16);
  2396. part2remain -= 16;
  2397. }
  2398. if (part2remain > 0)
  2399. getBits (part2remain);
  2400. else if (part2remain < 0)
  2401. return true;
  2402. return false;
  2403. }
  2404. void layer3Hybrid (float fsIn[32][18], float tsOut[18][32], int ch, const Layer3SideInfo::Info& granule) noexcept
  2405. {
  2406. float* ts = (float*) tsOut;
  2407. float* rawout1, *rawout2;
  2408. int sb = 0;
  2409. {
  2410. int b = hybridBlockIndex[ch];
  2411. rawout1 = hybridBlock[b][ch];
  2412. b = 1 - b;
  2413. rawout2 = hybridBlock[b][ch];
  2414. hybridBlockIndex[ch] = b;
  2415. }
  2416. if (granule.mixedBlockFlag)
  2417. {
  2418. sb = 2;
  2419. DCT::dct36 (fsIn[0], rawout1, rawout2, constants.win[0], ts);
  2420. DCT::dct36 (fsIn[1], rawout1 + 18, rawout2 + 18, constants.win1[0], ts + 1);
  2421. rawout1 += 36;
  2422. rawout2 += 36;
  2423. ts += 2;
  2424. }
  2425. const uint32 bt = granule.blockType;
  2426. if (bt == 2)
  2427. {
  2428. for (; sb < (int) granule.maxb; sb += 2, ts += 2, rawout1 += 36, rawout2 += 36)
  2429. {
  2430. DCT::dct12 (fsIn[sb], rawout1, rawout2, constants.win[2], ts);
  2431. DCT::dct12 (fsIn[sb + 1], rawout1 + 18, rawout2 + 18, constants.win1[2], ts + 1);
  2432. }
  2433. }
  2434. else
  2435. {
  2436. for (; sb < (int) granule.maxb; sb += 2, ts += 2, rawout1 += 36, rawout2 += 36)
  2437. {
  2438. DCT::dct36 (fsIn[sb], rawout1, rawout2, constants.win[bt], ts);
  2439. DCT::dct36 (fsIn[sb + 1], rawout1 + 18, rawout2 + 18, constants.win1[bt], ts + 1);
  2440. }
  2441. }
  2442. for (; sb < 32; ++sb, ++ts)
  2443. {
  2444. for (int i = 0; i < 18; ++i)
  2445. {
  2446. ts[i * 32] = *rawout1++;
  2447. *rawout2++ = 0;
  2448. }
  2449. }
  2450. }
  2451. void synthesiseStereo (const float* bandPtr0, const float* bandPtr1, float* out0, float* out1, int& samplesDone) noexcept
  2452. {
  2453. int dummy = samplesDone;
  2454. synthesise (bandPtr0, 0, out0, dummy);
  2455. synthesise (bandPtr1, 1, out1, samplesDone);
  2456. }
  2457. void synthesise (const float* bandPtr, const int channel, float* out, int& samplesDone)
  2458. {
  2459. out += samplesDone;
  2460. const int bo = channel == 0 ? ((synthBo - 1) & 15) : synthBo;
  2461. float (*buf)[0x110] = synthBuffers[channel];
  2462. float* b0;
  2463. int j, bo1 = bo;
  2464. if (bo & 1)
  2465. {
  2466. b0 = buf[0];
  2467. DCT::dct64 (buf[1] + ((bo + 1) & 15), buf[0] + bo, bandPtr);
  2468. }
  2469. else
  2470. {
  2471. ++bo1;
  2472. b0 = buf[1];
  2473. DCT::dct64 (buf[0] + bo, buf[1] + bo1, bandPtr);
  2474. }
  2475. synthBo = bo;
  2476. const float* window = constants.decodeWin + 16 - bo1;
  2477. for (j = 16; j != 0; --j, b0 += 16, window += 32)
  2478. {
  2479. float sum = window[0] * b0[0]; sum -= window[1] * b0[1];
  2480. sum += window[2] * b0[2]; sum -= window[3] * b0[3];
  2481. sum += window[4] * b0[4]; sum -= window[5] * b0[5];
  2482. sum += window[6] * b0[6]; sum -= window[7] * b0[7];
  2483. sum += window[8] * b0[8]; sum -= window[9] * b0[9];
  2484. sum += window[10] * b0[10]; sum -= window[11] * b0[11];
  2485. sum += window[12] * b0[12]; sum -= window[13] * b0[13];
  2486. sum += window[14] * b0[14]; sum -= window[15] * b0[15];
  2487. *out++ = sum;
  2488. }
  2489. {
  2490. float sum = window[0] * b0[0]; sum += window[2] * b0[2];
  2491. sum += window[4] * b0[4]; sum += window[6] * b0[6];
  2492. sum += window[8] * b0[8]; sum += window[10] * b0[10];
  2493. sum += window[12] * b0[12]; sum += window[14] * b0[14];
  2494. *out++ = sum;
  2495. b0 -= 16; window -= 32;
  2496. window += bo1 << 1;
  2497. }
  2498. for (j = 15; j != 0; --j, b0 -= 16, window -= 32)
  2499. {
  2500. float sum = -window[-1] * b0[0]; sum -= window[-2] * b0[1];
  2501. sum -= window[-3] * b0[2]; sum -= window[-4] * b0[3];
  2502. sum -= window[-5] * b0[4]; sum -= window[-6] * b0[5];
  2503. sum -= window[-7] * b0[6]; sum -= window[-8] * b0[7];
  2504. sum -= window[-9] * b0[8]; sum -= window[-10] * b0[9];
  2505. sum -= window[-11] * b0[10]; sum -= window[-12] * b0[11];
  2506. sum -= window[-13] * b0[12]; sum -= window[-14] * b0[13];
  2507. sum -= window[-15] * b0[14]; sum -= window[0] * b0[15];
  2508. *out++ = sum;
  2509. }
  2510. samplesDone += 32;
  2511. }
  2512. JUCE_DECLARE_NON_COPYABLE_WITH_LEAK_DETECTOR (MP3Stream)
  2513. };
  2514. //==============================================================================
  2515. static const char* const mp3FormatName = "MP3 file";
  2516. static const char* const mp3Extensions[] = { ".mp3", nullptr };
  2517. //==============================================================================
  2518. class MP3Reader : public AudioFormatReader
  2519. {
  2520. public:
  2521. MP3Reader (InputStream* const in)
  2522. : AudioFormatReader (in, TRANS (mp3FormatName)),
  2523. stream (*in), currentPosition (0),
  2524. decodedStart (0), decodedEnd (0)
  2525. {
  2526. skipID3();
  2527. const int64 streamPos = stream.stream.getPosition();
  2528. if (readNextBlock())
  2529. {
  2530. bitsPerSample = 32;
  2531. usesFloatingPointData = true;
  2532. sampleRate = stream.frame.getFrequency();
  2533. numChannels = (unsigned int) stream.frame.numChannels;
  2534. lengthInSamples = findLength (streamPos);
  2535. }
  2536. }
  2537. //==============================================================================
  2538. bool readSamples (int** destSamples, int numDestChannels, int startOffsetInDestBuffer,
  2539. int64 startSampleInFile, int numSamples) override
  2540. {
  2541. jassert (destSamples != nullptr);
  2542. if (currentPosition != startSampleInFile)
  2543. {
  2544. if (! stream.seek ((int) (startSampleInFile / 1152 - 1)))
  2545. {
  2546. currentPosition = -1;
  2547. createEmptyDecodedData();
  2548. }
  2549. else
  2550. {
  2551. decodedStart = decodedEnd = 0;
  2552. const int64 streamPos = stream.currentFrameIndex * 1152;
  2553. int toSkip = (int) (startSampleInFile - streamPos);
  2554. jassert (toSkip >= 0);
  2555. while (toSkip > 0)
  2556. {
  2557. if (! readNextBlock())
  2558. {
  2559. createEmptyDecodedData();
  2560. break;
  2561. }
  2562. const int numReady = decodedEnd - decodedStart;
  2563. if (numReady > toSkip)
  2564. {
  2565. decodedStart += toSkip;
  2566. break;
  2567. }
  2568. toSkip -= numReady;
  2569. }
  2570. currentPosition = startSampleInFile;
  2571. }
  2572. }
  2573. while (numSamples > 0)
  2574. {
  2575. if (decodedEnd <= decodedStart && ! readNextBlock())
  2576. {
  2577. for (int i = numDestChannels; --i >= 0;)
  2578. if (destSamples[i] != nullptr)
  2579. zeromem (destSamples[i] + startOffsetInDestBuffer, sizeof (float) * (size_t) numSamples);
  2580. return false;
  2581. }
  2582. const int numToCopy = jmin (decodedEnd - decodedStart, numSamples);
  2583. float* const* const dst = reinterpret_cast <float**> (destSamples);
  2584. memcpy (dst[0] + startOffsetInDestBuffer, decoded0 + decodedStart, sizeof (float) * (size_t) numToCopy);
  2585. if (numDestChannels > 1 && dst[1] != nullptr)
  2586. memcpy (dst[1] + startOffsetInDestBuffer, (numChannels < 2 ? decoded0 : decoded1) + decodedStart, sizeof (float) * (size_t) numToCopy);
  2587. startOffsetInDestBuffer += numToCopy;
  2588. decodedStart += numToCopy;
  2589. currentPosition += numToCopy;
  2590. numSamples -= numToCopy;
  2591. }
  2592. return true;
  2593. }
  2594. private:
  2595. MP3Stream stream;
  2596. int64 currentPosition;
  2597. enum { decodedDataSize = 1152 };
  2598. float decoded0 [decodedDataSize], decoded1 [decodedDataSize];
  2599. int decodedStart, decodedEnd;
  2600. void createEmptyDecodedData() noexcept
  2601. {
  2602. zeromem (decoded0, sizeof (decoded0));
  2603. zeromem (decoded1, sizeof (decoded1));
  2604. decodedStart = 0;
  2605. decodedEnd = decodedDataSize;
  2606. }
  2607. bool readNextBlock()
  2608. {
  2609. for (int attempts = 10; --attempts >= 0;)
  2610. {
  2611. int samplesDone = 0;
  2612. const int result = stream.decodeNextBlock (decoded0, decoded1, samplesDone);
  2613. if (result > 0 && stream.stream.isExhausted())
  2614. {
  2615. createEmptyDecodedData();
  2616. return true;
  2617. }
  2618. if (result <= 0)
  2619. {
  2620. decodedStart = 0;
  2621. decodedEnd = samplesDone;
  2622. return result == 0;
  2623. }
  2624. }
  2625. return false;
  2626. }
  2627. void skipID3()
  2628. {
  2629. const int64 originalPosition = stream.stream.getPosition();
  2630. const uint32 firstWord = (uint32) stream.stream.readInt();
  2631. if ((firstWord & 0xffffff) == 0x334449)
  2632. {
  2633. uint8 buffer[6];
  2634. if (stream.stream.read (buffer, 6) == 6
  2635. && buffer[0] != 0xff
  2636. && ((buffer[2] | buffer[3] | buffer[4] | buffer[5]) & 0x80) == 0)
  2637. {
  2638. const uint32 length = (((uint32) buffer[2]) << 21)
  2639. | (((uint32) buffer[3]) << 14)
  2640. | (((uint32) buffer[4]) << 7)
  2641. | ((uint32) buffer[5]);
  2642. stream.stream.skipNextBytes (length);
  2643. return;
  2644. }
  2645. }
  2646. stream.stream.setPosition (originalPosition);
  2647. }
  2648. int64 findLength (int64 streamStartPos)
  2649. {
  2650. int64 numFrames = stream.numFrames;
  2651. if (numFrames <= 0)
  2652. {
  2653. const int64 streamSize = stream.stream.getTotalLength();
  2654. if (streamSize > 0)
  2655. numFrames = (streamSize - streamStartPos) / (stream.frame.frameSize);
  2656. }
  2657. return numFrames * 1152;
  2658. }
  2659. JUCE_DECLARE_NON_COPYABLE_WITH_LEAK_DETECTOR (MP3Reader)
  2660. };
  2661. }
  2662. //==============================================================================
  2663. MP3AudioFormat::MP3AudioFormat()
  2664. : AudioFormat (MP3Decoder::mp3FormatName, StringArray (MP3Decoder::mp3Extensions))
  2665. {}
  2666. MP3AudioFormat::~MP3AudioFormat() {}
  2667. Array<int> MP3AudioFormat::getPossibleSampleRates() { return Array<int>(); }
  2668. Array<int> MP3AudioFormat::getPossibleBitDepths() { return Array<int>(); }
  2669. bool MP3AudioFormat::canDoStereo() { return true; }
  2670. bool MP3AudioFormat::canDoMono() { return true; }
  2671. bool MP3AudioFormat::isCompressed() { return true; }
  2672. StringArray MP3AudioFormat::getQualityOptions() { return StringArray(); }
  2673. AudioFormatReader* MP3AudioFormat::createReaderFor (InputStream* sourceStream, const bool deleteStreamIfOpeningFails)
  2674. {
  2675. ScopedPointer<MP3Decoder::MP3Reader> r (new MP3Decoder::MP3Reader (sourceStream));
  2676. if (r->lengthInSamples > 0)
  2677. return r.release();
  2678. if (! deleteStreamIfOpeningFails)
  2679. r->input = nullptr;
  2680. return nullptr;
  2681. }
  2682. AudioFormatWriter* MP3AudioFormat::createWriterFor (OutputStream*, double /*sampleRateToUse*/,
  2683. unsigned int /*numberOfChannels*/, int /*bitsPerSample*/,
  2684. const StringPairArray& /*metadataValues*/, int /*qualityOptionIndex*/)
  2685. {
  2686. return nullptr;
  2687. }
  2688. #endif