OpenAL32/sample_cvt.cpp

   1
   2 #include "config.h"
   3
   4 #include "sample_cvt.h"
   5
   6 #include "AL/al.h"
   7 #include "alu.h"
   8 #include "alBuffer.h"
   9
  10
  11 /* A quick'n'dirty lookup table to decode a muLaw-encoded byte sample into a
  12  * signed 16-bit sample */
  13 const ALshort muLawDecompressionTable[256] = {
  14     -32124,-31100,-30076,-29052,-28028,-27004,-25980,-24956,
  15     -23932,-22908,-21884,-20860,-19836,-18812,-17788,-16764,
  16     -15996,-15484,-14972,-14460,-13948,-13436,-12924,-12412,
  17     -11900,-11388,-10876,-10364, -9852, -9340, -8828, -8316,
  18      -7932, -7676, -7420, -7164, -6908, -6652, -6396, -6140,
  19      -5884, -5628, -5372, -5116, -4860, -4604, -4348, -4092,
  20      -3900, -3772, -3644, -3516, -3388, -3260, -3132, -3004,
  21      -2876, -2748, -2620, -2492, -2364, -2236, -2108, -1980,
  22      -1884, -1820, -1756, -1692, -1628, -1564, -1500, -1436,
  23      -1372, -1308, -1244, -1180, -1116, -1052,  -988,  -924,
  24       -876,  -844,  -812,  -780,  -748,  -716,  -684,  -652,
  25       -620,  -588,  -556,  -524,  -492,  -460,  -428,  -396,
  26       -372,  -356,  -340,  -324,  -308,  -292,  -276,  -260,
  27       -244,  -228,  -212,  -196,  -180,  -164,  -148,  -132,
  28       -120,  -112,  -104,   -96,   -88,   -80,   -72,   -64,
  29        -56,   -48,   -40,   -32,   -24,   -16,    -8,     0,
  30      32124, 31100, 30076, 29052, 28028, 27004, 25980, 24956,
  31      23932, 22908, 21884, 20860, 19836, 18812, 17788, 16764,
  32      15996, 15484, 14972, 14460, 13948, 13436, 12924, 12412,
  33      11900, 11388, 10876, 10364,  9852,  9340,  8828,  8316,
  34       7932,  7676,  7420,  7164,  6908,  6652,  6396,  6140,
  35       5884,  5628,  5372,  5116,  4860,  4604,  4348,  4092,
  36       3900,  3772,  3644,  3516,  3388,  3260,  3132,  3004,
  37       2876,  2748,  2620,  2492,  2364,  2236,  2108,  1980,
  38       1884,  1820,  1756,  1692,  1628,  1564,  1500,  1436,
  39       1372,  1308,  1244,  1180,  1116,  1052,   988,   924,
  40        876,   844,   812,   780,   748,   716,   684,   652,
  41        620,   588,   556,   524,   492,   460,   428,   396,
  42        372,   356,   340,   324,   308,   292,   276,   260,
  43        244,   228,   212,   196,   180,   164,   148,   132,
  44        120,   112,   104,    96,    88,    80,    72,    64,
  45         56,    48,    40,    32,    24,    16,     8,     0
  46 };
  47
  48 /* A quick'n'dirty lookup table to decode an aLaw-encoded byte sample into a
  49  * signed 16-bit sample */
  50 const ALshort aLawDecompressionTable[256] = {
  51      -5504, -5248, -6016, -5760, -4480, -4224, -4992, -4736,
  52      -7552, -7296, -8064, -7808, -6528, -6272, -7040, -6784,
  53      -2752, -2624, -3008, -2880, -2240, -2112, -2496, -2368,
  54      -3776, -3648, -4032, -3904, -3264, -3136, -3520, -3392,
  55     -22016,-20992,-24064,-23040,-17920,-16896,-19968,-18944,
  56     -30208,-29184,-32256,-31232,-26112,-25088,-28160,-27136,
  57     -11008,-10496,-12032,-11520, -8960, -8448, -9984, -9472,
  58     -15104,-14592,-16128,-15616,-13056,-12544,-14080,-13568,
  59       -344,  -328,  -376,  -360,  -280,  -264,  -312,  -296,
  60       -472,  -456,  -504,  -488,  -408,  -392,  -440,  -424,
  61        -88,   -72,  -120,  -104,   -24,    -8,   -56,   -40,
  62       -216,  -200,  -248,  -232,  -152,  -136,  -184,  -168,
  63      -1376, -1312, -1504, -1440, -1120, -1056, -1248, -1184,
  64      -1888, -1824, -2016, -1952, -1632, -1568, -1760, -1696,
  65       -688,  -656,  -752,  -720,  -560,  -528,  -624,  -592,
  66       -944,  -912, -1008,  -976,  -816,  -784,  -880,  -848,
  67       5504,  5248,  6016,  5760,  4480,  4224,  4992,  4736,
  68       7552,  7296,  8064,  7808,  6528,  6272,  7040,  6784,
  69       2752,  2624,  3008,  2880,  2240,  2112,  2496,  2368,
  70       3776,  3648,  4032,  3904,  3264,  3136,  3520,  3392,
  71      22016, 20992, 24064, 23040, 17920, 16896, 19968, 18944,
  72      30208, 29184, 32256, 31232, 26112, 25088, 28160, 27136,
  73      11008, 10496, 12032, 11520,  8960,  8448,  9984,  9472,
  74      15104, 14592, 16128, 15616, 13056, 12544, 14080, 13568,
  75        344,   328,   376,   360,   280,   264,   312,   296,
  76        472,   456,   504,   488,   408,   392,   440,   424,
  77         88,    72,   120,   104,    24,     8,    56,    40,
  78        216,   200,   248,   232,   152,   136,   184,   168,
  79       1376,  1312,  1504,  1440,  1120,  1056,  1248,  1184,
  80       1888,  1824,  2016,  1952,  1632,  1568,  1760,  1696,
  81        688,   656,   752,   720,   560,   528,   624,   592,
  82        944,   912,  1008,   976,   816,   784,   880,   848
  83 };
  84
  85 namespace {
  86
  87 /* IMA ADPCM Stepsize table */
  88 constexpr int IMAStep_size[89] = {
  89        7,    8,    9,   10,   11,   12,   13,   14,   16,   17,   19,
  90       21,   23,   25,   28,   31,   34,   37,   41,   45,   50,   55,
  91       60,   66,   73,   80,   88,   97,  107,  118,  130,  143,  157,
  92      173,  190,  209,  230,  253,  279,  307,  337,  371,  408,  449,
  93      494,  544,  598,  658,  724,  796,  876,  963, 1060, 1166, 1282,
  94     1411, 1552, 1707, 1878, 2066, 2272, 2499, 2749, 3024, 3327, 3660,
  95     4026, 4428, 4871, 5358, 5894, 6484, 7132, 7845, 8630, 9493,10442,
  96    11487,12635,13899,15289,16818,18500,20350,22358,24633,27086,29794,
  97    32767
  98 };
  99
 100 /* IMA4 ADPCM Codeword decode table */
 101 constexpr int IMA4Codeword[16] = {
 102     1, 3, 5, 7, 9, 11, 13, 15,
 103    -1,-3,-5,-7,-9,-11,-13,-15,
 104 };
 105
 106 /* IMA4 ADPCM Step index adjust decode table */
 107 constexpr int IMA4Index_adjust[16] = {
 108    -1,-1,-1,-1, 2, 4, 6, 8,
 109    -1,-1,-1,-1, 2, 4, 6, 8
 110 };
 111
 112
 113 /* MSADPCM Adaption table */
 114 constexpr int MSADPCMAdaption[16] = {
 115     230, 230, 230, 230, 307, 409, 512, 614,
 116     768, 614, 512, 409, 307, 230, 230, 230
 117 };
 118
 119 /* MSADPCM Adaption Coefficient tables */
 120 constexpr int MSADPCMAdaptionCoeff[7][2] = {
 121     { 256,    0 },
 122     { 512, -256 },
 123     {   0,    0 },
 124     { 192,   64 },
 125     { 240,    0 },
 126     { 460, -208 },
 127     { 392, -232 }
 128 };
 129
 130 void DecodeIMA4Block(ALshort *dst, const ALubyte *src, ALint numchans, ALsizei align)
 131 {
 132     ALint sample[MAX_INPUT_CHANNELS]{};
 133     ALint index[MAX_INPUT_CHANNELS]{};
 134     ALuint code[MAX_INPUT_CHANNELS]{};
 135
 136     for(int c{0};c < numchans;c++)
 137     {
 138         sample[c]  = *(src++);
 139         sample[c] |= *(src++) << 8;
 140         sample[c]  = (sample[c]^0x8000) - 32768;
 141         index[c]  = *(src++);
 142         index[c] |= *(src++) << 8;
 143         index[c]  = (index[c]^0x8000) - 32768;
 144
 145         index[c] = clampi(index[c], 0, 88);
 146
 147         dst[c] = sample[c];
 148     }
 149
 150     for(int i{1};i < align;i++)
 151     {
 152         if((i&7) == 1)
 153         {
 154             for(int c{0};c < numchans;c++)
 155             {
 156                 code[c]  = *(src++);
 157                 code[c] |= *(src++) << 8;
 158                 code[c] |= *(src++) << 16;
 159                 code[c] |= *(src++) << 24;
 160             }
 161         }
 162
 163         for(int c{0};c < numchans;c++)
 164         {
 165             int nibble = code[c]&0xf;
 166             code[c] >>= 4;
 167
 168             sample[c] += IMA4Codeword[nibble] * IMAStep_size[index[c]] / 8;
 169             sample[c] = clampi(sample[c], -32768, 32767);
 170
 171             index[c] += IMA4Index_adjust[nibble];
 172             index[c] = clampi(index[c], 0, 88);
 173
 174             *(dst++) = sample[c];
 175         }
 176     }
 177 }
 178
 179 void DecodeMSADPCMBlock(ALshort *dst, const ALubyte *src, ALint numchans, ALsizei align)
 180 {
 181     ALubyte blockpred[MAX_INPUT_CHANNELS]{};
 182     ALint delta[MAX_INPUT_CHANNELS]{};
 183     ALshort samples[MAX_INPUT_CHANNELS][2]{};
 184
 185     for(int c{0};c < numchans;c++)
 186     {
 187         blockpred[c] = *(src++);
 188         blockpred[c] = minu(blockpred[c], 6);
 189     }
 190     for(int c{0};c < numchans;c++)
 191     {
 192         delta[c]  = *(src++);
 193         delta[c] |= *(src++) << 8;
 194         delta[c]  = (delta[c]^0x8000) - 32768;
 195     }
 196     for(int c{0};c < numchans;c++)
 197     {
 198         samples[c][0]  = *(src++);
 199         samples[c][0] |= *(src++) << 8;
 200         samples[c][0]  = (samples[c][0]^0x8000) - 32768;
 201     }
 202     for(int c{0};c < numchans;c++)
 203     {
 204         samples[c][1]  = *(src++);
 205         samples[c][1] |= *(src++) << 8;
 206         samples[c][1]  = (samples[c][1]^0x8000) - 0x8000;
 207     }
 208
 209     /* Second sample is written first. */
 210     for(int c{0};c < numchans;c++)
 211         *(dst++) = samples[c][1];
 212     for(int c{0};c < numchans;c++)
 213         *(dst++) = samples[c][0];
 214
 215     for(int i{2};i < align;i++)
 216     {
 217         for(int c{0};c < numchans;c++)
 218         {
 219             const ALint num = (i*numchans) + c;
 220             ALint nibble, pred;
 221
 222             /* Read the nibble (first is in the upper bits). */
 223             if(!(num&1))
 224                 nibble = (*src>>4)&0x0f;
 225             else
 226                 nibble = (*(src++))&0x0f;
 227
 228             pred  = (samples[c][0]*MSADPCMAdaptionCoeff[blockpred[c]][0] +
 229                      samples[c][1]*MSADPCMAdaptionCoeff[blockpred[c]][1]) / 256;
 230             pred += ((nibble^0x08) - 0x08) * delta[c];
 231             pred  = clampi(pred, -32768, 32767);
 232
 233             samples[c][1] = samples[c][0];
 234             samples[c][0] = pred;
 235
 236             delta[c] = (MSADPCMAdaption[nibble] * delta[c]) / 256;
 237             delta[c] = maxi(16, delta[c]);
 238
 239             *(dst++) = pred;
 240         }
 241     }
 242 }
 243
 244 } // namespace
 245
 246 void Convert_ALshort_ALima4(ALshort *dst, const ALubyte *src, ALsizei numchans, ALsizei len,
 247                             ALsizei align)
 248 {
 249     ALsizei byte_align = ((align-1)/2 + 4) * numchans;
 250
 251     assert(align > 0 && (len%align) == 0);
 252     len /= align;
 253     while(len--)
 254     {
 255         DecodeIMA4Block(dst, src, numchans, align);
 256         src += byte_align;
 257         dst += align*numchans;
 258     }
 259 }
 260
 261 void Convert_ALshort_ALmsadpcm(ALshort *dst, const ALubyte *src, ALsizei numchans, ALsizei len,
 262                                ALsizei align)
 263 {
 264     const ALsizei byte_align = ((align-2)/2 + 7) * numchans;
 265
 266     assert(align > 1 && (len%align) == 0);
 267     len /= align;
 268     while(len--)
 269     {
 270         DecodeMSADPCMBlock(dst, src, numchans, align);
 271         src += byte_align;
 272         dst += align*numchans;
 273     }
 274 }