/*
 *  mpadec - MPEG audio decoder
 *  Copyright (C) 2002-2004 Dmitriy Startsev (dstartsev@rambler.ru)
 *
 *  This library is free software; you can redistribute it and/or
 *  modify it under the terms of the GNU Lesser General Public
 *  License as published by the Free Software Foundation; either
 *  version 2.1 of the License, or (at your option) any later version.
 *
 *  This library is distributed in the hope that it will be useful,
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 *  Lesser General Public License for more details.
 *
 *  You should have received a copy of the GNU Lesser General Public
 *  License along with this library; if not, write to the Free Software
 *  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
 */

/* $Id: synth.c,v 1.3 2004/08/03 05:22:22 metal_man Exp $ */

#include "mpadec_internal.h"

#define ROUND(x) (floor((x) + 0.5))
#define LROUND(x) ((int32_t)(ROUND(x)))
#define LLROUND(x) ((int64_t)(ROUND(x)))

static const FLOAT costab[32] = {
  0.50060299823519630134550410676638, 0.50547095989754365998444458560696,
  0.51544730992262454697495130564925, 0.53104259108978417447573393235719,
  0.55310389603444452782938083813705, 0.58293496820613387367383070125262,
  0.62250412303566481615725615676281, 0.67480834145500574602596871104104,
  0.74453627100229844977698119197295, 0.83934964541552703873926374662543,
  0.97256823786196069369768941405256, 1.16943993343288495515577028404220,
  1.48416461631416627724332693742810, 2.05778100995341155085655447971040,
  3.40760841846871878570119133345910, 10.1900081235480568112121092010360,
  0.50241928618815570551167011928012, 0.52249861493968888062857531905669,
  0.56694403481635770368053791515488, 0.64682178335999012954836011165200,
  0.78815462345125022473398248719742, 1.06067768599034747134045174723310,
  1.72244709823833392781591536415660, 5.10114861868916385810624549234540,
  0.50979557910415916894193980398784, 0.60134488693504528054372182390922,
  0.89997622313641570463850954094189, 2.56291544774150617879608629617770,
  0.54119610014619698439972320536639, 1.30656296487637652785664317342720,
  0.70710678118654752440084436210485, 0.0
};

static void dct64(FLOAT *outptr0, FLOAT *outptr1, FLOAT *samples)
{
  FLOAT tmp1[32], tmp2[32];

  {
    register FLOAT *in = samples;

    tmp1[0] = in[0] + in[31];
    tmp1[1] = in[1] + in[30];
    tmp1[31] = (in[0] - in[31])*costab[0];
    tmp1[30] = (in[1] - in[30])*costab[1];

    tmp1[2] = in[2] + in[29];
    tmp1[3] = in[3] + in[28];
    tmp1[29] = (in[2] - in[29])*costab[2];
    tmp1[28] = (in[3] - in[28])*costab[3];

    tmp1[4] = in[4] + in[27];
    tmp1[5] = in[5] + in[26];
    tmp1[27] = (in[4] - in[27])*costab[4];
    tmp1[26] = (in[5] - in[26])*costab[5];

    tmp1[6] = in[6] + in[25];
    tmp1[7] = in[7] + in[24];
    tmp1[25] = (in[6] - in[25])*costab[6];
    tmp1[24] = (in[7] - in[24])*costab[7];

    tmp1[8] = in[8] + in[23];
    tmp1[9] = in[9] + in[22];
    tmp1[23] = (in[8] - in[23])*costab[8];
    tmp1[22] = (in[9] - in[22])*costab[9];

    tmp1[10] = in[10] + in[21];
    tmp1[11] = in[11] + in[20];
    tmp1[21] = (in[10] - in[21])*costab[10];
    tmp1[20] = (in[11] - in[20])*costab[11];

    tmp1[12] = in[12] + in[19];
    tmp1[13] = in[13] + in[18];
    tmp1[19] = (in[12] - in[19])*costab[12];
    tmp1[18] = (in[13] - in[18])*costab[13];

    tmp1[14] = in[14] + in[17];
    tmp1[15] = in[15] + in[16];
    tmp1[17] = (in[14] - in[17])*costab[14];
    tmp1[16] = (in[15] - in[16])*costab[15];
  }
  {
    tmp2[0] = tmp1[0] + tmp1[15];
    tmp2[1] = tmp1[1] + tmp1[14];
    tmp2[15] = (tmp1[0] - tmp1[15])*costab[16 + 0];
    tmp2[14] = (tmp1[1] - tmp1[14])*costab[16 + 1];

    tmp2[2] = tmp1[2] + tmp1[13];
    tmp2[3] = tmp1[3] + tmp1[12];
    tmp2[13] = (tmp1[2] - tmp1[13])*costab[16 + 2];
    tmp2[12] = (tmp1[3] - tmp1[12])*costab[16 + 3];

    tmp2[4] = tmp1[4] + tmp1[11];
    tmp2[5] = tmp1[5] + tmp1[10];
    tmp2[11] = (tmp1[4] - tmp1[11])*costab[16 + 4];
    tmp2[10] = (tmp1[5] - tmp1[10])*costab[16 + 5];

    tmp2[6] = tmp1[6] + tmp1[9];
    tmp2[7] = tmp1[7] + tmp1[8];
    tmp2[9] = (tmp1[6] - tmp1[9])*costab[16 + 6];
    tmp2[8] = (tmp1[7] - tmp1[8])*costab[16 + 7];

    tmp2[16] = tmp1[16] + tmp1[31];
    tmp2[17] = tmp1[17] + tmp1[30];
    tmp2[31] = (tmp1[31] - tmp1[16])*costab[16 + 0];
    tmp2[30] = (tmp1[30] - tmp1[17])*costab[16 + 1];

    tmp2[18] = tmp1[18] + tmp1[29];
    tmp2[19] = tmp1[19] + tmp1[28];
    tmp2[29] = (tmp1[29] - tmp1[18])*costab[16 + 2];
    tmp2[28] = (tmp1[28] - tmp1[19])*costab[16 + 3];

    tmp2[20] = tmp1[20] + tmp1[27];
    tmp2[21] = tmp1[21] + tmp1[26];
    tmp2[27] = (tmp1[27] - tmp1[20])*costab[16 + 4];
    tmp2[26] = (tmp1[26] - tmp1[21])*costab[16 + 5];

    tmp2[22] = tmp1[22] + tmp1[25];
    tmp2[23] = tmp1[23] + tmp1[24];
    tmp2[25] = (tmp1[25] - tmp1[22])*costab[16 + 6];
    tmp2[24] = (tmp1[24] - tmp1[23])*costab[16 + 7];
  }
  {
    tmp1[0] = tmp2[0] + tmp2[7];
    tmp1[7] = (tmp2[0] - tmp2[7])*costab[16 + 8 + 0];
    tmp1[1] = tmp2[1] + tmp2[6];
    tmp1[6] = (tmp2[1] - tmp2[6])*costab[16 + 8 + 1];
    tmp1[2] = tmp2[2] + tmp2[5];
    tmp1[5] = (tmp2[2] - tmp2[5])*costab[16 + 8 + 2];
    tmp1[3] = tmp2[3] + tmp2[4];
    tmp1[4] = (tmp2[3] - tmp2[4])*costab[16 + 8 + 3];

    tmp1[8] = tmp2[8] + tmp2[15];
    tmp1[15] = (tmp2[15] - tmp2[8])*costab[16 + 8 + 0];
    tmp1[9] = tmp2[9] + tmp2[14];
    tmp1[14] = (tmp2[14] - tmp2[9])*costab[16 + 8 + 1];
    tmp1[10] = tmp2[10] + tmp2[13];
    tmp1[13] = (tmp2[13] - tmp2[10])*costab[16 + 8 + 2];
    tmp1[11] = tmp2[11] + tmp2[12];
    tmp1[12] = (tmp2[12] - tmp2[11])*costab[16 + 8 + 3];

    tmp1[16] = tmp2[16] + tmp2[23];
    tmp1[23] = (tmp2[16] - tmp2[23])*costab[16 + 8 + 0];
    tmp1[17] = tmp2[17] + tmp2[22];
    tmp1[22] = (tmp2[17] - tmp2[22])*costab[16 + 8 + 1];
    tmp1[18] = tmp2[18] + tmp2[21];
    tmp1[21] = (tmp2[18] - tmp2[21])*costab[16 + 8 + 2];
    tmp1[19] = tmp2[19] + tmp2[20];
    tmp1[20] = (tmp2[19] - tmp2[20])*costab[16 + 8 + 3];

    tmp1[24] = tmp2[24] + tmp2[31];
    tmp1[31] = (tmp2[31] - tmp2[24])*costab[16 + 8 + 0];
    tmp1[25] = tmp2[25] + tmp2[30];
    tmp1[30] = (tmp2[30] - tmp2[25])*costab[16 + 8 + 1];
    tmp1[26] = tmp2[26] + tmp2[29];
    tmp1[29] = (tmp2[29] - tmp2[26])*costab[16 + 8 + 2];
    tmp1[27] = tmp2[27] + tmp2[28];
    tmp1[28] = (tmp2[28] - tmp2[27])*costab[16 + 8 + 3];
  }
  {
    tmp2[0] = tmp1[0] + tmp1[3];
    tmp2[3] = (tmp1[0] - tmp1[3])*costab[16 + 8 + 4 + 0];
    tmp2[1] = tmp1[1] + tmp1[2];
    tmp2[2] = (tmp1[1] - tmp1[2])*costab[16 + 8 + 4 + 1];

    tmp2[4] = tmp1[4] + tmp1[7];
    tmp2[7] = (tmp1[7] - tmp1[4])*costab[16 + 8 + 4 + 0];
    tmp2[5] = tmp1[5] + tmp1[6];
    tmp2[6] = (tmp1[6] - tmp1[5])*costab[16 + 8 + 4 + 1];

    tmp2[8] = tmp1[8] + tmp1[11];
    tmp2[11] = (tmp1[8] - tmp1[11])*costab[16 + 8 + 4 + 0];
    tmp2[9] = tmp1[9] + tmp1[10];
    tmp2[10] = (tmp1[9] - tmp1[10])*costab[16 + 8 + 4 + 1];

    tmp2[12] = tmp1[12] + tmp1[15];
    tmp2[15] = (tmp1[15] - tmp1[12])*costab[16 + 8 + 4 + 0];
    tmp2[13] = tmp1[13] + tmp1[14];
    tmp2[14] = (tmp1[14] - tmp1[13])*costab[16 + 8 + 4 + 1];

    tmp2[16] = tmp1[16] + tmp1[19];
    tmp2[19] = (tmp1[16] - tmp1[19])*costab[16 + 8 + 4 + 0];
    tmp2[17] = tmp1[17] + tmp1[18];
    tmp2[18] = (tmp1[17] - tmp1[18])*costab[16 + 8 + 4 + 1];

    tmp2[20] = tmp1[20] + tmp1[23];
    tmp2[23] = (tmp1[23] - tmp1[20])*costab[16 + 8 + 4 + 0];
    tmp2[21] = tmp1[21] + tmp1[22];
    tmp2[22] = (tmp1[22] - tmp1[21])*costab[16 + 8 + 4 + 1];

    tmp2[24] = tmp1[24] + tmp1[27];
    tmp2[27] = (tmp1[24] - tmp1[27])*costab[16 + 8 + 4 + 0];
    tmp2[25] = tmp1[25] + tmp1[26];
    tmp2[26] = (tmp1[25] - tmp1[26])*costab[16 + 8 + 4 + 1];

    tmp2[28] = tmp1[28] + tmp1[31];
    tmp2[31] = (tmp1[31] - tmp1[28])*costab[16 + 8 + 4 + 0];
    tmp2[29] = tmp1[29] + tmp1[30];
    tmp2[30] = (tmp1[30] - tmp1[29])*costab[16 + 8 + 4 + 1];
  }
  {
    tmp1[0] = tmp2[0] + tmp2[1];
    tmp1[1] = (tmp2[0] - tmp2[1])*costab[16 + 8 + 4 + 2];
    tmp1[2] = tmp2[2] + tmp2[3];
    tmp1[3] = (tmp2[3] - tmp2[2])*costab[16 + 8 + 4 + 2];
    tmp1[2] += tmp1[3];

    tmp1[4] = tmp2[4] + tmp2[5];
    tmp1[5] = (tmp2[4] - tmp2[5])*costab[16 + 8 + 4 + 2];
    tmp1[6] = tmp2[6] + tmp2[7];
    tmp1[7] = (tmp2[7] - tmp2[6])*costab[16 + 8 + 4 + 2];
    tmp1[6] += tmp1[7];
    tmp1[4] += tmp1[6];
    tmp1[6] += tmp1[5];
    tmp1[5] += tmp1[7];

    tmp1[8] = tmp2[8] + tmp2[9];
    tmp1[9] = (tmp2[8] - tmp2[9])*costab[16 + 8 + 4 + 2];
    tmp1[10] = tmp2[10] + tmp2[11];
    tmp1[11] = (tmp2[11] - tmp2[10])*costab[16 + 8 + 4 + 2];
    tmp1[10] += tmp1[11];

    tmp1[12] = tmp2[12] + tmp2[13];
    tmp1[13] = (tmp2[12] - tmp2[13])*costab[16 + 8 + 4 + 2];
    tmp1[14] = tmp2[14] + tmp2[15];
    tmp1[15] = (tmp2[15] - tmp2[14])*costab[16 + 8 + 4 + 2];
    tmp1[14] += tmp1[15];
    tmp1[12] += tmp1[14];
    tmp1[14] += tmp1[13];
    tmp1[13] += tmp1[15];

    tmp1[16] = tmp2[16] + tmp2[17];
    tmp1[17] = (tmp2[16] - tmp2[17])*costab[16 + 8 + 4 + 2];
    tmp1[18] = tmp2[18] + tmp2[19];
    tmp1[19] = (tmp2[19] - tmp2[18])*costab[16 + 8 + 4 + 2];
    tmp1[18] += tmp1[19];

    tmp1[20] = tmp2[20] + tmp2[21];
    tmp1[21] = (tmp2[20] - tmp2[21])*costab[16 + 8 + 4 + 2];
    tmp1[22] = tmp2[22] + tmp2[23];
    tmp1[23] = (tmp2[23] - tmp2[22])*costab[16 + 8 + 4 + 2];
    tmp1[22] += tmp1[23];
    tmp1[20] += tmp1[22];
    tmp1[22] += tmp1[21];
    tmp1[21] += tmp1[23];

    tmp1[24] = tmp2[24] + tmp2[25];
    tmp1[25] = (tmp2[24] - tmp2[25])*costab[16 + 8 + 4 + 2];
    tmp1[26] = tmp2[26] + tmp2[27];
    tmp1[27] = (tmp2[27] - tmp2[26])*costab[16 + 8 + 4 + 2];
    tmp1[26] += tmp1[27];

    tmp1[28] = tmp2[28] + tmp2[29];
    tmp1[29] = (tmp2[28] - tmp2[29])*costab[16 + 8 + 4 + 2];
    tmp1[30] = tmp2[30] + tmp2[31];
    tmp1[31] = (tmp2[31] - tmp2[30])*costab[16 + 8 + 4 + 2];
    tmp1[30] += tmp1[31];
    tmp1[28] += tmp1[30];
    tmp1[30] += tmp1[29];
    tmp1[29] += tmp1[31];
  }
  {
    register FLOAT tmp, *out0 = outptr0, *out1 = outptr1;

    out0[16*16] = tmp1[0];
    out0[12*16] = tmp1[4];
    out0[8*16] = tmp1[2];
    out0[4*16] = tmp1[6];
    out0[0*16] = tmp1[1];
    out1[0*16] = tmp1[1];
    out1[4*16] = tmp1[5];
    out1[8*16] = tmp1[3];
    out1[12*16] = tmp1[7];

    out0[14*16] = tmp1[8] + tmp1[12];
    out0[10*16] = tmp1[12] + tmp1[10];
    out0[6*16] = tmp1[10] + tmp1[14];
    out0[2*16] = tmp1[14] + tmp1[9];
    out1[2*16] = tmp1[9] + tmp1[13];
    out1[6*16] = tmp1[13] + tmp1[11];
    out1[10*16] = tmp1[11] + tmp1[15];
    out1[14*16] = tmp1[15];

    tmp = tmp1[24] + tmp1[28];
    out0[15*16] = tmp + tmp1[16];
    out0[13*16] = tmp + tmp1[20];
    tmp = tmp1[28] + tmp1[26];
    out0[11*16] = tmp + tmp1[20];
    out0[9*16] = tmp + tmp1[18];
    tmp = tmp1[26] + tmp1[30];
    out0[7*16] = tmp + tmp1[18];
    out0[5*16] = tmp + tmp1[22];
    tmp = tmp1[30] + tmp1[25];
    out0[3*16] = tmp + tmp1[22];
    out0[1*16] = tmp + tmp1[17];
    tmp = tmp1[25] + tmp1[29];
    out1[1*16] = tmp + tmp1[17];
    out1[3*16] = tmp + tmp1[21];
    tmp = tmp1[29] + tmp1[27];
    out1[5*16] = tmp + tmp1[21];
    out1[7*16] = tmp + tmp1[19];
    tmp = tmp1[27] + tmp1[31];
    out1[9*16] = tmp + tmp1[19];
    out1[11*16] = tmp + tmp1[23];
    out1[13*16] = tmp1[23] + tmp1[31];
    out1[15*16] = tmp1[31];
  }
}

static void synth_full(mpadec_t mpadec, FLOAT *bandptr, int channel, FLOAT *buffer)
{
  register struct mpadec_t *mpa = (struct mpadec_t *)mpadec;
  unsigned bo;
  FLOAT *b0, (*buf)[0x110];

  if (!channel) {
    mpa->synth_bufoffs--;
    mpa->synth_bufoffs &= 0x0F;
    buf = mpa->synth_buffers[0];
  } else buf = mpa->synth_buffers[1];
  if (mpa->synth_bufoffs & 1) {
    b0 = buf[0];
    bo = mpa->synth_bufoffs;
    dct64(buf[1] + ((mpa->synth_bufoffs + 1) & 0x0F), buf[0] + mpa->synth_bufoffs, bandptr);
  } else {
    b0 = buf[1];
    bo = mpa->synth_bufoffs + 1;
    dct64(buf[0] + mpa->synth_bufoffs, buf[1] + (mpa->synth_bufoffs + 1), bandptr);
  }
  {
    register int i;
    register FLOAT *out = buffer;
    register FLOAT *win = mpa->tables.decwin + (16 - bo);

    for (i = 16; i; i--, win += 32, b0 += 16) {
      register FLOAT sum = win[0]*b0[0];
      sum -= win[1]*b0[1];
      sum += win[2]*b0[2];
      sum -= win[3]*b0[3];
      sum += win[4]*b0[4];
      sum -= win[5]*b0[5];
      sum += win[6]*b0[6];
      sum -= win[7]*b0[7];
      sum += win[8]*b0[8];
      sum -= win[9]*b0[9];
      sum += win[10]*b0[10];
      sum -= win[11]*b0[11];
      sum += win[12]*b0[12];
      sum -= win[13]*b0[13];
      sum += win[14]*b0[14];
      sum -= win[15]*b0[15];
      *out++ = sum;
    }
    {
      register FLOAT sum = win[0]*b0[0];
      sum += win[2]*b0[2];
      sum += win[4]*b0[4];
      sum += win[6]*b0[6];
      sum += win[8]*b0[8];
      sum += win[10]*b0[10];
      sum += win[12]*b0[12];
      sum += win[14]*b0[14];
      *out++ = sum;
      win -= 32; b0 -= 16;
    }
    win += (bo << 1);
    for (i = 15; i; i--, win -= 32, b0 -= 16)
    {
      register FLOAT sum = -win[-1]*b0[0];
      sum -= win[-2]*b0[1];
      sum -= win[-3]*b0[2];
      sum -= win[-4]*b0[3];
      sum -= win[-5]*b0[4];
      sum -= win[-6]*b0[5];
      sum -= win[-7]*b0[6];
      sum -= win[-8]*b0[7];
      sum -= win[-9]*b0[8];
      sum -= win[-10]*b0[9];
      sum -= win[-11]*b0[10];
      sum -= win[-12]*b0[11];
      sum -= win[-13]*b0[12];
      sum -= win[-14]*b0[13];
      sum -= win[-15]*b0[14];
      sum -= win[-0]*b0[15];
      *out++ = sum;
    }
  }
}

static void synth_half(mpadec_t mpadec, FLOAT *bandptr, int channel, FLOAT *buffer)
{
  register struct mpadec_t *mpa = (struct mpadec_t *)mpadec;
  unsigned bo;
  FLOAT *b0, (*buf)[0x110];

  if (!channel) {
    mpa->synth_bufoffs--;
    mpa->synth_bufoffs &= 0x0F;
    buf = mpa->synth_buffers[0];
  } else buf = mpa->synth_buffers[1];
  if (mpa->synth_bufoffs & 1) {
    b0 = buf[0];
    bo = mpa->synth_bufoffs;
    dct64(buf[1] + ((mpa->synth_bufoffs + 1) & 0x0F), buf[0] + mpa->synth_bufoffs, bandptr);
  } else {
    b0 = buf[1];
    bo = mpa->synth_bufoffs + 1;
    dct64(buf[0] + mpa->synth_bufoffs, buf[1] + (mpa->synth_bufoffs + 1), bandptr);
  }
  {
    register int i;
    register FLOAT *out = buffer;
    register FLOAT *win = mpa->tables.decwin + (16 - bo);

    for (i = 8; i; i--, win += 64, b0 += 32) {
      register FLOAT sum = win[0]*b0[0];
      sum -= win[1]*b0[1];
      sum += win[2]*b0[2];
      sum -= win[3]*b0[3];
      sum += win[4]*b0[4];
      sum -= win[5]*b0[5];
      sum += win[6]*b0[6];
      sum -= win[7]*b0[7];
      sum += win[8]*b0[8];
      sum -= win[9]*b0[9];
      sum += win[10]*b0[10];
      sum -= win[11]*b0[11];
      sum += win[12]*b0[12];
      sum -= win[13]*b0[13];
      sum += win[14]*b0[14];
      sum -= win[15]*b0[15];
      *out++ = sum;
    }
    {
      register FLOAT sum = win[0]*b0[0];
      sum += win[2]*b0[2];
      sum += win[4]*b0[4];
      sum += win[6]*b0[6];
      sum += win[8]*b0[8];
      sum += win[10]*b0[10];
      sum += win[12]*b0[12];
      sum += win[14]*b0[14];
      *out++ = sum;
      win -= 64; b0 -= 32;
    }
    win += (bo << 1);
    for (i = 7; i; i--, win -= 64, b0 -= 32)
    {
      register FLOAT sum = -win[-1]*b0[0];
      sum -= win[-2]*b0[1];
      sum -= win[-3]*b0[2];
      sum -= win[-4]*b0[3];
      sum -= win[-5]*b0[4];
      sum -= win[-6]*b0[5];
      sum -= win[-7]*b0[6];
      sum -= win[-8]*b0[7];
      sum -= win[-9]*b0[8];
      sum -= win[-10]*b0[9];
      sum -= win[-11]*b0[10];
      sum -= win[-12]*b0[11];
      sum -= win[-13]*b0[12];
      sum -= win[-14]*b0[13];
      sum -= win[-15]*b0[14];
      sum -= win[-0]*b0[15];
      *out++ = sum;
    }
  }
}

/* Full quality */

/* 16 bit, little-endian */

static void synth_full16lmm(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register int16_t *out = (int16_t *)buffer;
  FLOAT buf[SBLIMIT];

  synth_full(mpadec, bandptr, channel, buf);
  for (i = 0; i < SBLIMIT; i++, out++) {
    register int32_t tmp = LROUND(buf[i]);
    if (tmp > 32767) tmp = 32767; else if (tmp < -32768) tmp = -32768;
    ((uint8_t *)out)[0] = (uint8_t)tmp;
    ((int8_t *)out)[1] = (int8_t)(tmp >> 8);
  }
}

#define synth_full16lsm synth_full16lmm

static void synth_full16lms(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register int16_t *out = (int16_t *)buffer;
  FLOAT buf[SBLIMIT];

  synth_full(mpadec, bandptr, channel, buf);
  for (i = 0; i < SBLIMIT; i++, out += 2) {
    register int32_t tmp = LROUND(buf[i]);
    if (tmp > 32767) tmp = 32767; else if (tmp < -32768) tmp = -32768;
    ((uint8_t *)out)[0] = ((uint8_t *)out)[2] = (uint8_t)tmp;
    ((int8_t *)out)[1] = ((int8_t *)out)[3] = (int8_t)(tmp >> 8);
  }
}

static void synth_full16lss(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register int16_t *out = (int16_t *)buffer;
  FLOAT buf[SBLIMIT];

  synth_full(mpadec, bandptr, channel, buf);
  if (channel) out++;
  for (i = 0; i < SBLIMIT; i++, out += 2) {
    register int32_t tmp = LROUND(buf[i]);
    if (tmp > 32767) tmp = 32767; else if (tmp < -32768) tmp = -32768;
    ((uint8_t *)out)[0] = (uint8_t)tmp;
    ((int8_t *)out)[1] = (int8_t)(tmp >> 8);
  }
}

/* 16 bit, big-endian */

static void synth_full16bmm(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register int16_t *out = (int16_t *)buffer;
  FLOAT buf[SBLIMIT];

  synth_full(mpadec, bandptr, channel, buf);
  for (i = 0; i < SBLIMIT; i++, out++) {
    register int32_t tmp = LROUND(buf[i]);
    if (tmp > 32767) tmp = 32767; else if (tmp < -32768) tmp = -32768;
    ((uint8_t *)out)[1] = (uint8_t)tmp;
    ((int8_t *)out)[0] = (int8_t)(tmp >> 8);
  }
}

#define synth_full16bsm synth_full16bmm

static void synth_full16bms(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register int16_t *out = (int16_t *)buffer;
  FLOAT buf[SBLIMIT];

  synth_full(mpadec, bandptr, channel, buf);
  for (i = 0; i < SBLIMIT; i++, out += 2) {
    register int32_t tmp = LROUND(buf[i]);
    if (tmp > 32767) tmp = 32767; else if (tmp < -32768) tmp = -32768;
    ((uint8_t *)out)[1] = ((uint8_t *)out)[3] = (uint8_t)tmp;
    ((int8_t *)out)[0] = ((int8_t *)out)[2] = (int8_t)(tmp >> 8);
  }
}

static void synth_full16bss(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register int16_t *out = (int16_t *)buffer;
  FLOAT buf[SBLIMIT];

  synth_full(mpadec, bandptr, channel, buf);
  if (channel) out++;
  for (i = 0; i < SBLIMIT; i++, out += 2) {
    register int32_t tmp = LROUND(buf[i]);
    if (tmp > 32767) tmp = 32767; else if (tmp < -32768) tmp = -32768;
    ((uint8_t *)out)[1] = (uint8_t)tmp;
    ((int8_t *)out)[0] = (int8_t)(tmp >> 8);
  }
}

/* 24 bit, little-endian */

static void synth_full24lmm(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register uint8_t *out = (uint8_t *)buffer;
  FLOAT buf[SBLIMIT];

  synth_full(mpadec, bandptr, channel, buf);
  for (i = 0; i < SBLIMIT; i++, out += 3) {
    register int32_t tmp = LROUND(buf[i]);
    if (tmp > 0x7FFFFF) tmp = 0x7FFFFF; else if (tmp < -0x800000) tmp = -0x800000;
    ((uint8_t *)out)[0] = (uint8_t)tmp;
    ((uint8_t *)out)[1] = (uint8_t)(tmp >> 8);
    ((int8_t *)out)[2] = (int8_t)(tmp >> 16);
  }
}

#define synth_full24lsm synth_full24lmm

static void synth_full24lms(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register uint8_t *out = (uint8_t *)buffer;
  FLOAT buf[SBLIMIT];

  synth_full(mpadec, bandptr, channel, buf);
  for (i = 0; i < SBLIMIT; i++, out += 6) {
    register int32_t tmp = LROUND(buf[i]);
    if (tmp > 0x7FFFFF) tmp = 0x7FFFFF; else if (tmp < -0x800000) tmp = -0x800000;
    ((uint8_t *)out)[0] = ((uint8_t *)out)[3] = (uint8_t)tmp;
    ((uint8_t *)out)[1] = ((uint8_t *)out)[4] = (uint8_t)(tmp >> 8);
    ((int8_t *)out)[2] = ((int8_t *)out)[5] = (int8_t)(tmp >> 16);
  }
}

static void synth_full24lss(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register uint8_t *out = (uint8_t *)buffer;
  FLOAT buf[SBLIMIT];

  synth_full(mpadec, bandptr, channel, buf);
  if (channel) out += 3;
  for (i = 0; i < SBLIMIT; i++, out += 6) {
    register int32_t tmp = LROUND(buf[i]);
    if (tmp > 0x7FFFFF) tmp = 0x7FFFFF; else if (tmp < -0x800000) tmp = -0x800000;
    ((uint8_t *)out)[0] = (uint8_t)tmp;
    ((uint8_t *)out)[1] = (uint8_t)(tmp >> 8);
    ((int8_t *)out)[2] = (int8_t)(tmp >> 16);
  }
}

/* 24 bit, big-endian */

static void synth_full24bmm(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register uint8_t *out = (uint8_t *)buffer;
  FLOAT buf[SBLIMIT];

  synth_full(mpadec, bandptr, channel, buf);
  for (i = 0; i < SBLIMIT; i++, out += 3) {
    register int32_t tmp = LROUND(buf[i]);
    if (tmp > 0x7FFFFF) tmp = 0x7FFFFF; else if (tmp < -0x800000) tmp = -0x800000;
    ((uint8_t *)out)[2] = (uint8_t)tmp;
    ((uint8_t *)out)[1] = (uint8_t)(tmp >> 8);
    ((int8_t *)out)[0] = (int8_t)(tmp >> 16);
  }
}

#define synth_full24bsm synth_full24bmm

static void synth_full24bms(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register uint8_t *out = (uint8_t *)buffer;
  FLOAT buf[SBLIMIT];

  synth_full(mpadec, bandptr, channel, buf);
  for (i = 0; i < SBLIMIT; i++, out += 6) {
    register int32_t tmp = LROUND(buf[i]);
    if (tmp > 0x7FFFFF) tmp = 0x7FFFFF; else if (tmp < -0x800000) tmp = -0x800000;
    ((uint8_t *)out)[2] = ((uint8_t *)out)[5] = (uint8_t)tmp;
    ((uint8_t *)out)[1] = ((uint8_t *)out)[4] = (uint8_t)(tmp >> 8);
    ((int8_t *)out)[0] = ((int8_t *)out)[3] = (int8_t)(tmp >> 16);
  }
}

static void synth_full24bss(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register uint8_t *out = (uint8_t *)buffer;
  FLOAT buf[SBLIMIT];

  synth_full(mpadec, bandptr, channel, buf);
  if (channel) out += 3;
  for (i = 0; i < SBLIMIT; i++, out += 6) {
    register int32_t tmp = LROUND(buf[i]);
    if (tmp > 0x7FFFFF) tmp = 0x7FFFFF; else if (tmp < -0x800000) tmp = -0x800000;
    ((uint8_t *)out)[2] = (uint8_t)tmp;
    ((uint8_t *)out)[1] = (uint8_t)(tmp >> 8);
    ((int8_t *)out)[0] = (int8_t)(tmp >> 16);
  }
}

/* 32 bit , little-endian */

static void synth_full32lmm(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register int32_t *out = (int32_t *)buffer;
  FLOAT buf[SBLIMIT];

  synth_full(mpadec, bandptr, channel, buf);
  for (i = 0; i < SBLIMIT; i++, out++) {
    register int64_t tmp = LLROUND(buf[i]);
    if (tmp > 0x7FFFFFFF) tmp = 0x7FFFFFFF; else if (tmp < (-0x7FFFFFFF - 1)) tmp = (-0x7FFFFFFF - 1);
    ((uint8_t *)out)[0] = (uint8_t)tmp;
    ((uint8_t *)out)[1] = (uint8_t)(tmp >> 8);
    ((uint8_t *)out)[2] = (uint8_t)(tmp >> 16);
    ((int8_t *)out)[3] = (int8_t)(tmp >> 24);
  }
}

#define synth_full32lsm synth_full32lmm

static void synth_full32lms(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register int32_t *out = (int32_t *)buffer;
  FLOAT buf[SBLIMIT];

  synth_full(mpadec, bandptr, channel, buf);
  for (i = 0; i < SBLIMIT; i++, out += 2) {
    register int64_t tmp = LLROUND(buf[i]);
    if (tmp > 0x7FFFFFFF) tmp = 0x7FFFFFFF; else if (tmp < (-0x7FFFFFFF - 1)) tmp = (-0x7FFFFFFF - 1);
    ((uint8_t *)out)[0] = ((uint8_t *)out)[4] = (uint8_t)tmp;
    ((uint8_t *)out)[1] = ((uint8_t *)out)[5] = (uint8_t)(tmp >> 8);
    ((uint8_t *)out)[2] = ((uint8_t *)out)[6] = (uint8_t)(tmp >> 16);
    ((int8_t *)out)[3] = ((int8_t *)out)[7] = (int8_t)(tmp >> 24);
  }
}

static void synth_full32lss(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register int32_t *out = (int32_t *)buffer;
  FLOAT buf[SBLIMIT];

  synth_full(mpadec, bandptr, channel, buf);
  if (channel) out++;
  for (i = 0; i < SBLIMIT; i++, out += 2) {
    register int64_t tmp = LLROUND(buf[i]);
    if (tmp > 0x7FFFFFFF) tmp = 0x7FFFFFFF; else if (tmp < (-0x7FFFFFFF - 1)) tmp = (-0x7FFFFFFF - 1);
    ((uint8_t *)out)[0] = (uint8_t)tmp;
    ((uint8_t *)out)[1] = (uint8_t)(tmp >> 8);
    ((uint8_t *)out)[2] = (uint8_t)(tmp >> 16);
    ((int8_t *)out)[3] = (int8_t)(tmp >> 24);
  }
}

/* 32 bit, big-endian */

static void synth_full32bmm(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register int32_t *out = (int32_t *)buffer;
  FLOAT buf[SBLIMIT];

  synth_full(mpadec, bandptr, channel, buf);
  for (i = 0; i < SBLIMIT; i++, out++) {
    register int64_t tmp = LLROUND(buf[i]);
    if (tmp > 0x7FFFFFFF) tmp = 0x7FFFFFFF; else if (tmp < (-0x7FFFFFFF - 1)) tmp = (-0x7FFFFFFF - 1);
    ((uint8_t *)out)[3] = (uint8_t)tmp;
    ((uint8_t *)out)[2] = (uint8_t)(tmp >> 8);
    ((uint8_t *)out)[1] = (uint8_t)(tmp >> 16);
    ((int8_t *)out)[0] = (int8_t)(tmp >> 24);
  }
}

#define synth_full32bsm synth_full32bmm

static void synth_full32bms(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register int32_t *out = (int32_t *)buffer;
  FLOAT buf[SBLIMIT];

  synth_full(mpadec, bandptr, channel, buf);
  for (i = 0; i < SBLIMIT; i++, out += 2) {
    register int64_t tmp = LLROUND(buf[i]);
    if (tmp > 0x7FFFFFFF) tmp = 0x7FFFFFFF; else if (tmp < (-0x7FFFFFFF - 1)) tmp = (-0x7FFFFFFF - 1);
    ((uint8_t *)out)[3] = ((uint8_t *)out)[7] = (uint8_t)tmp;
    ((uint8_t *)out)[2] = ((uint8_t *)out)[6] = (uint8_t)(tmp >> 8);
    ((uint8_t *)out)[1] = ((uint8_t *)out)[5] = (uint8_t)(tmp >> 16);
    ((int8_t *)out)[0] = ((int8_t *)out)[4] = (int8_t)(tmp >> 24);
  }
}

static void synth_full32bss(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register int32_t *out = (int32_t *)buffer;
  FLOAT buf[SBLIMIT];

  synth_full(mpadec, bandptr, channel, buf);
  if (channel) out++;
  for (i = 0; i < SBLIMIT; i++, out += 2) {
    register int64_t tmp = LLROUND(buf[i]);
    if (tmp > 0x7FFFFFFF) tmp = 0x7FFFFFFF; else if (tmp < (-0x7FFFFFFF - 1)) tmp = (-0x7FFFFFFF - 1);
    ((uint8_t *)out)[3] = (uint8_t)tmp;
    ((uint8_t *)out)[2] = (uint8_t)(tmp >> 8);
    ((uint8_t *)out)[1] = (uint8_t)(tmp >> 16);
    ((int8_t *)out)[0] = (int8_t)(tmp >> 24);
  }
}

/* 32 bit floating-point, little-endian */

static void synth_full32flmm(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register float *out = (float *)buffer;
  FLOAT buf[SBLIMIT];

  synth_full(mpadec, bandptr, channel, buf);
  for (i = 0; i < SBLIMIT; i++, out++) {
    int32_t tmp;
    *((float *)(&tmp)) = (float)buf[i];
    ((uint8_t *)out)[0] = (uint8_t)tmp;
    ((uint8_t *)out)[1] = (uint8_t)(tmp >> 8);
    ((uint8_t *)out)[2] = (uint8_t)(tmp >> 16);
    ((int8_t *)out)[3] = (int8_t)(tmp >> 24);
  }
}

#define synth_full32flsm synth_full32flmm

static void synth_full32flms(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register float *out = (float *)buffer;
  FLOAT buf[SBLIMIT];

  synth_full(mpadec, bandptr, channel, buf);
  for (i = 0; i < SBLIMIT; i++, out += 2) {
    int32_t tmp;
    *((float *)(&tmp)) = (float)buf[i];
    ((uint8_t *)out)[0] = ((uint8_t *)out)[4] = (uint8_t)tmp;
    ((uint8_t *)out)[1] = ((uint8_t *)out)[5] = (uint8_t)(tmp >> 8);
    ((uint8_t *)out)[2] = ((uint8_t *)out)[6] = (uint8_t)(tmp >> 16);
    ((int8_t *)out)[3] = ((int8_t *)out)[7] = (int8_t)(tmp >> 24);
  }
}

static void synth_full32flss(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register float *out = (float *)buffer;
  FLOAT buf[SBLIMIT];

  synth_full(mpadec, bandptr, channel, buf);
  if (channel) out++;
  for (i = 0; i < SBLIMIT; i++, out += 2) {
    int32_t tmp;
    *((float *)(&tmp)) = (float)buf[i];
    ((uint8_t *)out)[0] = (uint8_t)tmp;
    ((uint8_t *)out)[1] = (uint8_t)(tmp >> 8);
    ((uint8_t *)out)[2] = (uint8_t)(tmp >> 16);
    ((int8_t *)out)[3] = (int8_t)(tmp >> 24);
  }
}

/* 32 bit floating-point, big-endian */

static void synth_full32fbmm(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register float *out = (float *)buffer;
  FLOAT buf[SBLIMIT];

  synth_full(mpadec, bandptr, channel, buf);
  for (i = 0; i < SBLIMIT; i++, out++) {
    int32_t tmp;
    *((float *)(&tmp)) = (float)buf[i];
    ((uint8_t *)out)[3] = (uint8_t)tmp;
    ((uint8_t *)out)[2] = (uint8_t)(tmp >> 8);
    ((uint8_t *)out)[1] = (uint8_t)(tmp >> 16);
    ((int8_t *)out)[0] = (int8_t)(tmp >> 24);
  }
}

#define synth_full32fbsm synth_full32fbmm

static void synth_full32fbms(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register float *out = (float *)buffer;
  FLOAT buf[SBLIMIT];

  synth_full(mpadec, bandptr, channel, buf);
  for (i = 0; i < SBLIMIT; i++, out += 2) {
    int32_t tmp;
    *((float *)(&tmp)) = (float)buf[i];
    ((uint8_t *)out)[3] = ((uint8_t *)out)[7] = (uint8_t)tmp;
    ((uint8_t *)out)[2] = ((uint8_t *)out)[6] = (uint8_t)(tmp >> 8);
    ((uint8_t *)out)[1] = ((uint8_t *)out)[5] = (uint8_t)(tmp >> 16);
    ((int8_t *)out)[0] = ((int8_t *)out)[4] = (int8_t)(tmp >> 24);
  }
}

static void synth_full32fbss(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register float *out = (float *)buffer;
  FLOAT buf[SBLIMIT];

  synth_full(mpadec, bandptr, channel, buf);
  if (channel) out++;
  for (i = 0; i < SBLIMIT; i++, out += 2) {
    int32_t tmp;
    *((float *)(&tmp)) = (float)buf[i];
    ((uint8_t *)out)[3] = (uint8_t)tmp;
    ((uint8_t *)out)[2] = (uint8_t)(tmp >> 8);
    ((uint8_t *)out)[1] = (uint8_t)(tmp >> 16);
    ((int8_t *)out)[0] = (int8_t)(tmp >> 24);
  }
}

/* Half quality */

/* 16 bit, little-endian */

static void synth_half16lmm(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register int16_t *out = (int16_t *)buffer;
  FLOAT buf[SBLIMIT/2];

  synth_half(mpadec, bandptr, channel, buf);
  for (i = 0; i < SBLIMIT/2; i++, out++) {
    register int32_t tmp = LROUND(buf[i]);
    if (tmp > 32767) tmp = 32767; else if (tmp < -32768) tmp = -32768;
    ((uint8_t *)out)[0] = (uint8_t)tmp;
    ((int8_t *)out)[1] = (int8_t)(tmp >> 8);
  }
}

#define synth_half16lsm synth_half16lmm

static void synth_half16lms(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register int16_t *out = (int16_t *)buffer;
  FLOAT buf[SBLIMIT/2];

  synth_half(mpadec, bandptr, channel, buf);
  for (i = 0; i < SBLIMIT/2; i++, out += 2) {
    register int32_t tmp = LROUND(buf[i]);
    if (tmp > 32767) tmp = 32767; else if (tmp < -32768) tmp = -32768;
    ((uint8_t *)out)[0] = ((uint8_t *)out)[2] = (uint8_t)tmp;
    ((int8_t *)out)[1] = ((int8_t *)out)[3] = (int8_t)(tmp >> 8);
  }
}

static void synth_half16lss(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register int16_t *out = (int16_t *)buffer;
  FLOAT buf[SBLIMIT/2];

  synth_half(mpadec, bandptr, channel, buf);
  if (channel) out++;
  for (i = 0; i < SBLIMIT/2; i++, out += 2) {
    register int32_t tmp = LROUND(buf[i]);
    if (tmp > 32767) tmp = 32767; else if (tmp < -32768) tmp = -32768;
    ((uint8_t *)out)[0] = (uint8_t)tmp;
    ((int8_t *)out)[1] = (int8_t)(tmp >> 8);
  }
}

/* 16 bit, big-endian */

static void synth_half16bmm(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register int16_t *out = (int16_t *)buffer;
  FLOAT buf[SBLIMIT/2];

  synth_half(mpadec, bandptr, channel, buf);
  for (i = 0; i < SBLIMIT/2; i++, out++) {
    register int32_t tmp = LROUND(buf[i]);
    if (tmp > 32767) tmp = 32767; else if (tmp < -32768) tmp = -32768;
    ((uint8_t *)out)[1] = (uint8_t)tmp;
    ((int8_t *)out)[0] = (int8_t)(tmp >> 8);
  }
}

#define synth_half16bsm synth_half16bmm

static void synth_half16bms(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register int16_t *out = (int16_t *)buffer;
  FLOAT buf[SBLIMIT/2];

  synth_half(mpadec, bandptr, channel, buf);
  for (i = 0; i < SBLIMIT/2; i++, out += 2) {
    register int32_t tmp = LROUND(buf[i]);
    if (tmp > 32767) tmp = 32767; else if (tmp < -32768) tmp = -32768;
    ((uint8_t *)out)[1] = ((uint8_t *)out)[3] = (uint8_t)tmp;
    ((int8_t *)out)[0] = ((int8_t *)out)[2] = (int8_t)(tmp >> 8);
  }
}

static void synth_half16bss(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register int16_t *out = (int16_t *)buffer;
  FLOAT buf[SBLIMIT/2];

  synth_half(mpadec, bandptr, channel, buf);
  if (channel) out++;
  for (i = 0; i < SBLIMIT/2; i++, out += 2) {
    register int32_t tmp = LROUND(buf[i]);
    if (tmp > 32767) tmp = 32767; else if (tmp < -32768) tmp = -32768;
    ((uint8_t *)out)[1] = (uint8_t)tmp;
    ((int8_t *)out)[0] = (int8_t)(tmp >> 8);
  }
}

/* 24 bit, little-endian */

static void synth_half24lmm(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register uint8_t *out = (uint8_t *)buffer;
  FLOAT buf[SBLIMIT/2];

  synth_half(mpadec, bandptr, channel, buf);
  for (i = 0; i < SBLIMIT/2; i++, out += 3) {
    register int32_t tmp = LROUND(buf[i]);
    if (tmp > 0x7FFFFF) tmp = 0x7FFFFF; else if (tmp < -0x800000) tmp = -0x800000;
    ((uint8_t *)out)[0] = (uint8_t)tmp;
    ((uint8_t *)out)[1] = (uint8_t)(tmp >> 8);
    ((int8_t *)out)[2] = (int8_t)(tmp >> 16);
  }
}

#define synth_half24lsm synth_half24lmm

static void synth_half24lms(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register uint8_t *out = (uint8_t *)buffer;
  FLOAT buf[SBLIMIT/2];

  synth_half(mpadec, bandptr, channel, buf);
  for (i = 0; i < SBLIMIT/2; i++, out += 6) {
    register int32_t tmp = LROUND(buf[i]);
    if (tmp > 0x7FFFFF) tmp = 0x7FFFFF; else if (tmp < -0x800000) tmp = -0x800000;
    ((uint8_t *)out)[0] = ((uint8_t *)out)[3] = (uint8_t)tmp;
    ((uint8_t *)out)[1] = ((uint8_t *)out)[4] = (uint8_t)(tmp >> 8);
    ((int8_t *)out)[2] = ((int8_t *)out)[5] = (int8_t)(tmp >> 16);
  }
}

static void synth_half24lss(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register uint8_t *out = (uint8_t *)buffer;
  FLOAT buf[SBLIMIT/2];

  synth_half(mpadec, bandptr, channel, buf);
  if (channel) out += 3;
  for (i = 0; i < SBLIMIT/2; i++, out += 6) {
    register int32_t tmp = LROUND(buf[i]);
    if (tmp > 0x7FFFFF) tmp = 0x7FFFFF; else if (tmp < -0x800000) tmp = -0x800000;
    ((uint8_t *)out)[0] = (uint8_t)tmp;
    ((uint8_t *)out)[1] = (uint8_t)(tmp >> 8);
    ((int8_t *)out)[2] = (int8_t)(tmp >> 16);
  }
}

/* 24 bit, big-endian */

static void synth_half24bmm(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register uint8_t *out = (uint8_t *)buffer;
  FLOAT buf[SBLIMIT/2];

  synth_half(mpadec, bandptr, channel, buf);
  for (i = 0; i < SBLIMIT/2; i++, out += 3) {
    register int32_t tmp = LROUND(buf[i]);
    if (tmp > 0x7FFFFF) tmp = 0x7FFFFF; else if (tmp < -0x800000) tmp = -0x800000;
    ((uint8_t *)out)[2] = (uint8_t)tmp;
    ((uint8_t *)out)[1] = (uint8_t)(tmp >> 8);
    ((int8_t *)out)[0] = (int8_t)(tmp >> 16);
  }
}

#define synth_half24bsm synth_half24bmm

static void synth_half24bms(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register uint8_t *out = (uint8_t *)buffer;
  FLOAT buf[SBLIMIT/2];

  synth_half(mpadec, bandptr, channel, buf);
  for (i = 0; i < SBLIMIT/2; i++, out += 6) {
    register int32_t tmp = LROUND(buf[i]);
    if (tmp > 0x7FFFFF) tmp = 0x7FFFFF; else if (tmp < -0x800000) tmp = -0x800000;
    ((uint8_t *)out)[2] = ((uint8_t *)out)[5] = (uint8_t)tmp;
    ((uint8_t *)out)[1] = ((uint8_t *)out)[4] = (uint8_t)(tmp >> 8);
    ((int8_t *)out)[0] = ((int8_t *)out)[3] = (int8_t)(tmp >> 16);
  }
}

static void synth_half24bss(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register uint8_t *out = (uint8_t *)buffer;
  FLOAT buf[SBLIMIT/2];

  synth_half(mpadec, bandptr, channel, buf);
  if (channel) out += 3;
  for (i = 0; i < SBLIMIT/2; i++, out += 6) {
    register int32_t tmp = LROUND(buf[i]);
    if (tmp > 0x7FFFFF) tmp = 0x7FFFFF; else if (tmp < -0x800000) tmp = -0x800000;
    ((uint8_t *)out)[2] = (uint8_t)tmp;
    ((uint8_t *)out)[1] = (uint8_t)(tmp >> 8);
    ((int8_t *)out)[0] = (int8_t)(tmp >> 16);
  }
}

/* 32 bit, little-endian */

static void synth_half32lmm(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register int32_t *out = (int32_t *)buffer;
  FLOAT buf[SBLIMIT/2];

  synth_half(mpadec, bandptr, channel, buf);
  for (i = 0; i < SBLIMIT/2; i++, out++) {
    register int64_t tmp = LLROUND(buf[i]);
    if (tmp > 0x7FFFFFFF) tmp = 0x7FFFFFFF; else if (tmp < (-0x7FFFFFFF - 1)) tmp = (-0x7FFFFFFF - 1);
    ((uint8_t *)out)[0] = (uint8_t)tmp;
    ((uint8_t *)out)[1] = (uint8_t)(tmp >> 8);
    ((uint8_t *)out)[2] = (uint8_t)(tmp >> 16);
    ((int8_t *)out)[3] = (int8_t)(tmp >> 24);
  }
}

#define synth_half32lsm synth_half32lmm

static void synth_half32lms(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register int32_t *out = (int32_t *)buffer;
  FLOAT buf[SBLIMIT/2];

  synth_half(mpadec, bandptr, channel, buf);
  for (i = 0; i < SBLIMIT/2; i++, out += 2) {
    register int64_t tmp = LLROUND(buf[i]);
    if (tmp > 0x7FFFFFFF) tmp = 0x7FFFFFFF; else if (tmp < (-0x7FFFFFFF - 1)) tmp = (-0x7FFFFFFF - 1);
    ((uint8_t *)out)[0] = ((uint8_t *)out)[4] = (uint8_t)tmp;
    ((uint8_t *)out)[1] = ((uint8_t *)out)[5] = (uint8_t)(tmp >> 8);
    ((uint8_t *)out)[2] = ((uint8_t *)out)[6] = (uint8_t)(tmp >> 16);
    ((int8_t *)out)[3] = ((int8_t *)out)[7] = (int8_t)(tmp >> 24);
  }
}

static void synth_half32lss(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register int32_t *out = (int32_t *)buffer;
  FLOAT buf[SBLIMIT/2];

  synth_half(mpadec, bandptr, channel, buf);
  if (channel) out++;
  for (i = 0; i < SBLIMIT/2; i++, out += 2) {
    register int64_t tmp = LLROUND(buf[i]);
    if (tmp > 0x7FFFFFFF) tmp = 0x7FFFFFFF; else if (tmp < (-0x7FFFFFFF - 1)) tmp = (-0x7FFFFFFF - 1);
    ((uint8_t *)out)[0] = (uint8_t)tmp;
    ((uint8_t *)out)[1] = (uint8_t)(tmp >> 8);
    ((uint8_t *)out)[2] = (uint8_t)(tmp >> 16);
    ((int8_t *)out)[3] = (int8_t)(tmp >> 24);
  }
}

/* 32 bit, big-endian */

static void synth_half32bmm(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register int32_t *out = (int32_t *)buffer;
  FLOAT buf[SBLIMIT/2];

  synth_half(mpadec, bandptr, channel, buf);
  for (i = 0; i < SBLIMIT/2; i++, out++) {
    register int64_t tmp = LLROUND(buf[i]);
    if (tmp > 0x7FFFFFFF) tmp = 0x7FFFFFFF; else if (tmp < (-0x7FFFFFFF - 1)) tmp = (-0x7FFFFFFF - 1);
    ((uint8_t *)out)[3] = (uint8_t)tmp;
    ((uint8_t *)out)[2] = (uint8_t)(tmp >> 8);
    ((uint8_t *)out)[1] = (uint8_t)(tmp >> 16);
    ((int8_t *)out)[0] = (int8_t)(tmp >> 24);
  }
}

#define synth_half32bsm synth_half32bmm

static void synth_half32bms(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register int32_t *out = (int32_t *)buffer;
  FLOAT buf[SBLIMIT/2];

  synth_half(mpadec, bandptr, channel, buf);
  for (i = 0; i < SBLIMIT/2; i++, out += 2) {
    register int64_t tmp = LLROUND(buf[i]);
    if (tmp > 0x7FFFFFFF) tmp = 0x7FFFFFFF; else if (tmp < (-0x7FFFFFFF - 1)) tmp = (-0x7FFFFFFF - 1);
    ((uint8_t *)out)[3] = ((uint8_t *)out)[7] = (uint8_t)tmp;
    ((uint8_t *)out)[2] = ((uint8_t *)out)[6] = (uint8_t)(tmp >> 8);
    ((uint8_t *)out)[1] = ((uint8_t *)out)[5] = (uint8_t)(tmp >> 16);
    ((int8_t *)out)[0] = ((int8_t *)out)[4] = (int8_t)(tmp >> 24);
  }
}

static void synth_half32bss(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register int32_t *out = (int32_t *)buffer;
  FLOAT buf[SBLIMIT/2];

  synth_half(mpadec, bandptr, channel, buf);
  if (channel) out++;
  for (i = 0; i < SBLIMIT/2; i++, out += 2) {
    register int64_t tmp = LLROUND(buf[i]);
    if (tmp > 0x7FFFFFFF) tmp = 0x7FFFFFFF; else if (tmp < (-0x7FFFFFFF - 1)) tmp = (-0x7FFFFFFF - 1);
    ((uint8_t *)out)[3] = (uint8_t)tmp;
    ((uint8_t *)out)[2] = (uint8_t)(tmp >> 8);
    ((uint8_t *)out)[1] = (uint8_t)(tmp >> 16);
    ((int8_t *)out)[0] = (int8_t)(tmp >> 24);
  }
}

/* 32 bit floating-point, little-endian */

static void synth_half32flmm(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register float *out = (float *)buffer;
  FLOAT buf[SBLIMIT/2];

  synth_half(mpadec, bandptr, channel, buf);
  for (i = 0; i < SBLIMIT/2; i++, out++) {
    int32_t tmp;
    *((float *)(&tmp)) = (float)buf[i];
    ((uint8_t *)out)[0] = (uint8_t)tmp;
    ((uint8_t *)out)[1] = (uint8_t)(tmp >> 8);
    ((uint8_t *)out)[2] = (uint8_t)(tmp >> 16);
    ((int8_t *)out)[3] = (int8_t)(tmp >> 24);
  }
}

#define synth_half32flsm synth_half32flmm

static void synth_half32flms(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register float *out = (float *)buffer;
  FLOAT buf[SBLIMIT/2];

  synth_half(mpadec, bandptr, channel, buf);
  for (i = 0; i < SBLIMIT/2; i++, out += 2) {
    int32_t tmp;
    *((float *)(&tmp)) = (float)buf[i];
    ((uint8_t *)out)[0] = ((uint8_t *)out)[4] = (uint8_t)tmp;
    ((uint8_t *)out)[1] = ((uint8_t *)out)[5] = (uint8_t)(tmp >> 8);
    ((uint8_t *)out)[2] = ((uint8_t *)out)[6] = (uint8_t)(tmp >> 16);
    ((int8_t *)out)[3] = ((int8_t *)out)[7] = (int8_t)(tmp >> 24);
  }
}

static void synth_half32flss(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register float *out = (float *)buffer;
  FLOAT buf[SBLIMIT/2];

  synth_half(mpadec, bandptr, channel, buf);
  if (channel) out++;
  for (i = 0; i < SBLIMIT/2; i++, out += 2) {
    int32_t tmp;
    *((float *)(&tmp)) = (float)buf[i];
    ((uint8_t *)out)[0] = (uint8_t)tmp;
    ((uint8_t *)out)[1] = (uint8_t)(tmp >> 8);
    ((uint8_t *)out)[2] = (uint8_t)(tmp >> 16);
    ((int8_t *)out)[3] = (int8_t)(tmp >> 24);
  }
}

/* 32 bit floating-point, big-endian */

static void synth_half32fbmm(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register float *out = (float *)buffer;
  FLOAT buf[SBLIMIT/2];

  synth_half(mpadec, bandptr, channel, buf);
  for (i = 0; i < SBLIMIT/2; i++, out++) {
    int32_t tmp;
    *((float *)(&tmp)) = (float)buf[i];
    ((uint8_t *)out)[3] = (uint8_t)tmp;
    ((uint8_t *)out)[2] = (uint8_t)(tmp >> 8);
    ((uint8_t *)out)[1] = (uint8_t)(tmp >> 16);
    ((int8_t *)out)[0] = (int8_t)(tmp >> 24);
  }
}

#define synth_half32fbsm synth_half32fbmm

static void synth_half32fbms(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register float *out = (float *)buffer;
  FLOAT buf[SBLIMIT/2];

  synth_half(mpadec, bandptr, channel, buf);
  for (i = 0; i < SBLIMIT/2; i++, out += 2) {
    int32_t tmp;
    *((float *)(&tmp)) = (float)buf[i];
    ((uint8_t *)out)[3] = ((uint8_t *)out)[7] = (uint8_t)tmp;
    ((uint8_t *)out)[2] = ((uint8_t *)out)[6] = (uint8_t)(tmp >> 8);
    ((uint8_t *)out)[1] = ((uint8_t *)out)[5] = (uint8_t)(tmp >> 16);
    ((int8_t *)out)[0] = ((int8_t *)out)[4] = (int8_t)(tmp >> 24);
  }
}

static void synth_half32fbss(mpadec_t mpadec, FLOAT *bandptr, int channel, uint8_t *buffer)
{
  register int i;
  register float *out = (float *)buffer;
  FLOAT buf[SBLIMIT/2];

  synth_half(mpadec, bandptr, channel, buf);
  if (channel) out++;
  for (i = 0; i < SBLIMIT/2; i++, out += 2) {
    int32_t tmp;
    *((float *)(&tmp)) = (float)buf[i];
    ((uint8_t *)out)[3] = (uint8_t)tmp;
    ((uint8_t *)out)[2] = (uint8_t)(tmp >> 8);
    ((uint8_t *)out)[1] = (uint8_t)(tmp >> 16);
    ((int8_t *)out)[0] = (int8_t)(tmp >> 24);
  }
}

void *synth_table[2][2][4][4] = {
  { { { synth_full16lmm,  synth_full16lms,  synth_full16lsm,  synth_full16lss  },
      { synth_full24lmm,  synth_full24lms,  synth_full24lsm,  synth_full24lss  },
      { synth_full32lmm,  synth_full32lms,  synth_full32lsm,  synth_full32lss  },
      { synth_full32flmm, synth_full32flms, synth_full32flsm, synth_full32flss } },
    { { synth_full16bmm,  synth_full16bms,  synth_full16bsm,  synth_full16bss  },
      { synth_full24bmm,  synth_full24bms,  synth_full24bsm,  synth_full24bss  },
      { synth_full32bmm,  synth_full32bms,  synth_full32bsm,  synth_full32bss  },
      { synth_full32fbmm, synth_full32fbms, synth_full32fbsm, synth_full32fbss } } },
  { { { synth_half16lmm,  synth_half16lms,  synth_half16lsm,  synth_half16lss  },
      { synth_half24lmm,  synth_half24lms,  synth_half24lsm,  synth_half24lss  },
      { synth_half32lmm,  synth_half32lms,  synth_half32lsm,  synth_half32lss  },
      { synth_half32flmm, synth_half32flms, synth_half32flsm, synth_half32flss } },
    { { synth_half16bmm,  synth_half16bms,  synth_half16bsm,  synth_half16bss  },
      { synth_half24bmm,  synth_half24bms,  synth_half24bsm,  synth_half24bss  },
      { synth_half32bmm,  synth_half32bms,  synth_half32bsm,  synth_half32bss  },
      { synth_half32fbmm, synth_half32fbms, synth_half32fbsm, synth_half32fbss } } }
};