git.oblomov.eu Git - wine/blob - dlls/winemp3.acm/layer3.c

   1 /*
   2  * Mpeg Layer-3 audio decoder
   3  * --------------------------
   4  * copyright (c) 1995,1996,1997 by Michael Hipp.
   5  * All rights reserved. See also 'README'
   6  * This file has been copied from mpglib.
   7  *
   8  * This library is free software; you can redistribute it and/or
   9  * modify it under the terms of the GNU Lesser General Public
  10  * License as published by the Free Software Foundation; either
  11  * version 2.1 of the License, or (at your option) any later version.
  12  *
  13  * This library is distributed in the hope that it will be useful,
  14  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  15  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  16  * Lesser General Public License for more details.
  17  *
  18  * You should have received a copy of the GNU Lesser General Public
  19  * License along with this library; if not, write to the Free Software
  20  * Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301, USA
  21  */
  22
  23 #include <stdlib.h>
  24 #include "mpg123.h"
  25 #include "mpglib.h"
  26 #include "huffman.h"
  27
  28 extern struct mpstr *gmp;
  29
  30 #define MPEG1
  31
  32
  33 static real ispow[8207];
  34 static real aa_ca[8],aa_cs[8];
  35 static real COS1[12][6];
  36 static real win[4][36];
  37 static real win1[4][36];
  38 static real gainpow2[256+118+4];
  39 static real COS9[9];
  40 static real COS6_1,COS6_2;
  41 static real tfcos36[9];
  42 static real tfcos12[3];
  43
  44 struct bandInfoStruct {
  45   short longIdx[23];
  46   short longDiff[22];
  47   short shortIdx[14];
  48   short shortDiff[13];
  49 };
  50
  51 static int longLimit[9][23];
  52 static int shortLimit[9][14];
  53
  54 static const struct bandInfoStruct bandInfo[9] = {
  55
  56 /* MPEG 1.0 */
  57  { {0,4,8,12,16,20,24,30,36,44,52,62,74, 90,110,134,162,196,238,288,342,418,576},
  58    {4,4,4,4,4,4,6,6,8, 8,10,12,16,20,24,28,34,42,50,54, 76,158},
  59    {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},
  60    {4,4,4,4,6,8,10,12,14,18,22,30,56} } ,
  61
  62  { {0,4,8,12,16,20,24,30,36,42,50,60,72, 88,106,128,156,190,230,276,330,384,576},
  63    {4,4,4,4,4,4,6,6,6, 8,10,12,16,18,22,28,34,40,46,54, 54,192},
  64    {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},
  65    {4,4,4,4,6,6,10,12,14,16,20,26,66} } ,
  66
  67  { {0,4,8,12,16,20,24,30,36,44,54,66,82,102,126,156,194,240,296,364,448,550,576} ,
  68    {4,4,4,4,4,4,6,6,8,10,12,16,20,24,30,38,46,56,68,84,102, 26} ,
  69    {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} ,
  70    {4,4,4,4,6,8,12,16,20,26,34,42,12} }  ,
  71
  72 /* MPEG 2.0 */
  73  { {0,6,12,18,24,30,36,44,54,66,80,96,116,140,168,200,238,284,336,396,464,522,576},
  74    {6,6,6,6,6,6,8,10,12,14,16,20,24,28,32,38,46,52,60,68,58,54 } ,
  75    {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} ,
  76    {4,4,4,6,6,8,10,14,18,26,32,42,18 } } ,
  77
  78  { {0,6,12,18,24,30,36,44,54,66,80,96,114,136,162,194,232,278,330,394,464,540,576},
  79    {6,6,6,6,6,6,8,10,12,14,16,18,22,26,32,38,46,52,64,70,76,36 } ,
  80    {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} ,
  81    {4,4,4,6,8,10,12,14,18,24,32,44,12 } } ,
  82
  83  { {0,6,12,18,24,30,36,44,54,66,80,96,116,140,168,200,238,284,336,396,464,522,576},
  84    {6,6,6,6,6,6,8,10,12,14,16,20,24,28,32,38,46,52,60,68,58,54 },
  85    {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},
  86    {4,4,4,6,8,10,12,14,18,24,30,40,18 } } ,
  87 /* MPEG 2.5 */
  88  { {0,6,12,18,24,30,36,44,54,66,80,96,116,140,168,200,238,284,336,396,464,522,576} ,
  89    {6,6,6,6,6,6,8,10,12,14,16,20,24,28,32,38,46,52,60,68,58,54},
  90    {0,12,24,36,54,78,108,144,186,240,312,402,522,576},
  91    {4,4,4,6,8,10,12,14,18,24,30,40,18} },
  92  { {0,6,12,18,24,30,36,44,54,66,80,96,116,140,168,200,238,284,336,396,464,522,576} ,
  93    {6,6,6,6,6,6,8,10,12,14,16,20,24,28,32,38,46,52,60,68,58,54},
  94    {0,12,24,36,54,78,108,144,186,240,312,402,522,576},
  95    {4,4,4,6,8,10,12,14,18,24,30,40,18} },
  96  { {0,12,24,36,48,60,72,88,108,132,160,192,232,280,336,400,476,566,568,570,572,574,576},
  97    {12,12,12,12,12,12,16,20,24,28,32,40,48,56,64,76,90,2,2,2,2,2},
  98    {0, 24, 48, 72,108,156,216,288,372,480,486,492,498,576},
  99    {8,8,8,12,16,20,24,28,36,2,2,2,26} } ,
 100 };
 101
 102 static int mapbuf0[9][152];
 103 static int mapbuf1[9][156];
 104 static int mapbuf2[9][44];
 105 static int *map[9][3];
 106 static int *mapend[9][3];
 107
 108 static unsigned int n_slen2[512]; /* MPEG 2.0 slen for 'normal' mode */
 109 static unsigned int i_slen2[256]; /* MPEG 2.0 slen for intensity stereo */
 110
 111 static real tan1_1[16],tan2_1[16],tan1_2[16],tan2_2[16];
 112 static real pow1_1[2][16],pow2_1[2][16],pow1_2[2][16],pow2_2[2][16];
 113
 114 /*
 115  * init tables for layer-3
 116  */
 117 void init_layer3(int down_sample_sblimit)
 118 {
 119   int i,j,k,l;
 120
 121   for(i=-256;i<118+4;i++)
 122     gainpow2[i+256] = pow((double)2.0,-0.25 * (double) (i+210) );
 123
 124   for(i=0;i<8207;i++)
 125     ispow[i] = pow((double)i,(double)4.0/3.0);
 126
 127   for (i=0;i<8;i++)
 128   {
 129     static const double Ci[8]={-0.6,-0.535,-0.33,-0.185,-0.095,-0.041,-0.0142,-0.0037};
 130     double sq=sqrt(1.0+Ci[i]*Ci[i]);
 131     aa_cs[i] = 1.0/sq;
 132     aa_ca[i] = Ci[i]/sq;
 133   }
 134
 135   for(i=0;i<18;i++)
 136   {
 137     win[0][i]    = win[1][i]    = 0.5 * sin( M_PI / 72.0 * (double) (2*(i+0) +1) ) / cos ( M_PI * (double) (2*(i+0) +19) / 72.0 );
 138     win[0][i+18] = win[3][i+18] = 0.5 * sin( M_PI / 72.0 * (double) (2*(i+18)+1) ) / cos ( M_PI * (double) (2*(i+18)+19) / 72.0 );
 139   }
 140   for(i=0;i<6;i++)
 141   {
 142     win[1][i+18] = 0.5 / cos ( M_PI * (double) (2*(i+18)+19) / 72.0 );
 143     win[3][i+12] = 0.5 / cos ( M_PI * (double) (2*(i+12)+19) / 72.0 );
 144     win[1][i+24] = 0.5 * sin( M_PI / 24.0 * (double) (2*i+13) ) / cos ( M_PI * (double) (2*(i+24)+19) / 72.0 );
 145     win[1][i+30] = win[3][i] = 0.0;
 146     win[3][i+6 ] = 0.5 * sin( M_PI / 24.0 * (double) (2*i+1) )  / cos ( M_PI * (double) (2*(i+6 )+19) / 72.0 );
 147   }
 148
 149   for(i=0;i<9;i++)
 150     COS9[i] = cos( M_PI / 18.0 * (double) i);
 151
 152   for(i=0;i<9;i++)
 153     tfcos36[i] = 0.5 / cos ( M_PI * (double) (i*2+1) / 36.0 );
 154   for(i=0;i<3;i++)
 155     tfcos12[i] = 0.5 / cos ( M_PI * (double) (i*2+1) / 12.0 );
 156
 157   COS6_1 = cos( M_PI / 6.0 * (double) 1);
 158   COS6_2 = cos( M_PI / 6.0 * (double) 2);
 159
 160   for(i=0;i<12;i++)
 161   {
 162     win[2][i]  = 0.5 * sin( M_PI / 24.0 * (double) (2*i+1) ) / cos ( M_PI * (double) (2*i+7) / 24.0 );
 163     for(j=0;j<6;j++)
 164       COS1[i][j] = cos( M_PI / 24.0 * (double) ((2*i+7)*(2*j+1)) );
 165   }
 166
 167   for(j=0;j<4;j++) {
 168     static const int len[4] = { 36,36,12,36 };
 169     for(i=0;i<len[j];i+=2)
 170       win1[j][i] = + win[j][i];
 171     for(i=1;i<len[j];i+=2)
 172       win1[j][i] = - win[j][i];
 173   }
 174
 175   for(i=0;i<16;i++)
 176   {
 177     double t = tan( (double) i * M_PI / 12.0 );
 178     tan1_1[i] = t / (1.0+t);
 179     tan2_1[i] = 1.0 / (1.0 + t);
 180     tan1_2[i] = M_SQRT2 * t / (1.0+t);
 181     tan2_2[i] = M_SQRT2 / (1.0 + t);
 182
 183     for(j=0;j<2;j++) {
 184       double base = pow(2.0,-0.25*(j+1.0));
 185       double p1=1.0,p2=1.0;
 186       if(i > 0) {
 187         if( i & 1 )
 188           p1 = pow(base,(i+1.0)*0.5);
 189         else
 190           p2 = pow(base,i*0.5);
 191       }
 192       pow1_1[j][i] = p1;
 193       pow2_1[j][i] = p2;
 194       pow1_2[j][i] = M_SQRT2 * p1;
 195       pow2_2[j][i] = M_SQRT2 * p2;
 196     }
 197   }
 198
 199   for(j=0;j<9;j++)
 200   {
 201    const struct bandInfoStruct *bi = &bandInfo[j];
 202    int *mp;
 203    int cb,lwin;
 204    const short *bdf;
 205
 206    mp = map[j][0] = mapbuf0[j];
 207    bdf = bi->longDiff;
 208    for(i=0,cb = 0; cb < 8 ; cb++,i+=*bdf++) {
 209      *mp++ = (*bdf) >> 1;
 210      *mp++ = i;
 211      *mp++ = 3;
 212      *mp++ = cb;
 213    }
 214    bdf = bi->shortDiff+3;
 215    for(cb=3;cb<13;cb++) {
 216      int l = (*bdf++) >> 1;
 217      for(lwin=0;lwin<3;lwin++) {
 218        *mp++ = l;
 219        *mp++ = i + lwin;
 220        *mp++ = lwin;
 221        *mp++ = cb;
 222      }
 223      i += 6*l;
 224    }
 225    mapend[j][0] = mp;
 226
 227    mp = map[j][1] = mapbuf1[j];
 228    bdf = bi->shortDiff+0;
 229    for(i=0,cb=0;cb<13;cb++) {
 230      int l = (*bdf++) >> 1;
 231      for(lwin=0;lwin<3;lwin++) {
 232        *mp++ = l;
 233        *mp++ = i + lwin;
 234        *mp++ = lwin;
 235        *mp++ = cb;
 236      }
 237      i += 6*l;
 238    }
 239    mapend[j][1] = mp;
 240
 241    mp = map[j][2] = mapbuf2[j];
 242    bdf = bi->longDiff;
 243    for(cb = 0; cb < 22 ; cb++) {
 244      *mp++ = (*bdf++) >> 1;
 245      *mp++ = cb;
 246    }
 247    mapend[j][2] = mp;
 248
 249   }
 250
 251   for(j=0;j<9;j++) {
 252     for(i=0;i<23;i++) {
 253       longLimit[j][i] = (bandInfo[j].longIdx[i] - 1 + 8) / 18 + 1;
 254       if(longLimit[j][i] > (down_sample_sblimit) )
 255         longLimit[j][i] = down_sample_sblimit;
 256     }
 257     for(i=0;i<14;i++) {
 258       shortLimit[j][i] = (bandInfo[j].shortIdx[i] - 1) / 18 + 1;
 259       if(shortLimit[j][i] > (down_sample_sblimit) )
 260         shortLimit[j][i] = down_sample_sblimit;
 261     }
 262   }
 263
 264   for(i=0;i<5;i++) {
 265     for(j=0;j<6;j++) {
 266       for(k=0;k<6;k++) {
 267         int n = k + j * 6 + i * 36;
 268         i_slen2[n] = i|(j<<3)|(k<<6)|(3<<12);
 269       }
 270     }
 271   }
 272   for(i=0;i<4;i++) {
 273     for(j=0;j<4;j++) {
 274       for(k=0;k<4;k++) {
 275         int n = k + j * 4 + i * 16;
 276         i_slen2[n+180] = i|(j<<3)|(k<<6)|(4<<12);
 277       }
 278     }
 279   }
 280   for(i=0;i<4;i++) {
 281     for(j=0;j<3;j++) {
 282       int n = j + i * 3;
 283       i_slen2[n+244] = i|(j<<3) | (5<<12);
 284       n_slen2[n+500] = i|(j<<3) | (2<<12) | (1<<15);
 285     }
 286   }
 287
 288   for(i=0;i<5;i++) {
 289     for(j=0;j<5;j++) {
 290       for(k=0;k<4;k++) {
 291         for(l=0;l<4;l++) {
 292           int n = l + k * 4 + j * 16 + i * 80;
 293           n_slen2[n] = i|(j<<3)|(k<<6)|(l<<9)|(0<<12);
 294         }
 295       }
 296     }
 297   }
 298   for(i=0;i<5;i++) {
 299     for(j=0;j<5;j++) {
 300       for(k=0;k<4;k++) {
 301         int n = k + j * 4 + i * 20;
 302         n_slen2[n+400] = i|(j<<3)|(k<<6)|(1<<12);
 303       }
 304     }
 305   }
 306 }
 307
 308 /*
 309  * read additional side information
 310  */
 311 #ifdef MPEG1
 312 static int III_get_side_info_1(struct III_sideinfo *si,int stereo,
 313  int ms_stereo,long sfreq,int single)
 314 {
 315    int ch, gr;
 316    int powdiff = (single == 3) ? 4 : 0;
 317
 318    si->main_data_begin = getbits(9);
 319    if (stereo == 1)
 320      si->private_bits = getbits_fast(5);
 321    else
 322      si->private_bits = getbits_fast(3);
 323
 324    for (ch=0; ch<stereo; ch++) {
 325        si->ch[ch].gr[0].scfsi = -1;
 326        si->ch[ch].gr[1].scfsi = getbits_fast(4);
 327    }
 328
 329    for (gr=0; gr<2; gr++)
 330    {
 331      for (ch=0; ch<stereo; ch++)
 332      {
 333        register struct gr_info_s *gr_info = &(si->ch[ch].gr[gr]);
 334
 335        gr_info->part2_3_length = getbits(12);
 336        gr_info->big_values = getbits_fast(9);
 337        if(gr_info->big_values > 288) {
 338           fprintf(stderr,"big_values too large!\n");
 339           gr_info->big_values = 288;
 340        }
 341        gr_info->pow2gain = gainpow2+256 - getbits_fast(8) + powdiff;
 342        if(ms_stereo)
 343          gr_info->pow2gain += 2;
 344        gr_info->scalefac_compress = getbits_fast(4);
 345 /* window-switching flag == 1 for block_Type != 0 .. and block-type == 0 -> win-sw-flag = 0 */
 346        if(get1bit())
 347        {
 348          int i;
 349          gr_info->block_type = getbits_fast(2);
 350          gr_info->mixed_block_flag = get1bit();
 351          gr_info->table_select[0] = getbits_fast(5);
 352          gr_info->table_select[1] = getbits_fast(5);
 353          /*
 354           * table_select[2] not needed, because there is no region2,
 355           * but to satisfy some verifications tools we set it either.
 356           */
 357          gr_info->table_select[2] = 0;
 358          for(i=0;i<3;i++)
 359            gr_info->full_gain[i] = gr_info->pow2gain + (getbits_fast(3)<<3);
 360
 361          if(gr_info->block_type == 0) {
 362            fprintf(stderr,"Blocktype == 0 and window-switching == 1 not allowed.\n");
 363            return 0;
 364          }
 365          /* region_count/start parameters are implicit in this case. */
 366          gr_info->region1start = 36>>1;
 367          gr_info->region2start = 576>>1;
 368        }
 369        else
 370        {
 371          int i,r0c,r1c;
 372          for (i=0; i<3; i++)
 373            gr_info->table_select[i] = getbits_fast(5);
 374          r0c = getbits_fast(4);
 375          r1c = getbits_fast(3);
 376          gr_info->region1start = bandInfo[sfreq].longIdx[r0c+1] >> 1 ;
 377          gr_info->region2start = bandInfo[sfreq].longIdx[r0c+1+r1c+1] >> 1;
 378          gr_info->block_type = 0;
 379          gr_info->mixed_block_flag = 0;
 380        }
 381        gr_info->preflag = get1bit();
 382        gr_info->scalefac_scale = get1bit();
 383        gr_info->count1table_select = get1bit();
 384      }
 385    }
 386    return !0;
 387 }
 388 #endif
 389
 390 /*
 391  * Side Info for MPEG 2.0 / LSF
 392  */
 393 static int III_get_side_info_2(struct III_sideinfo *si,int stereo,
 394  int ms_stereo,long sfreq,int single)
 395 {
 396    int ch;
 397    int powdiff = (single == 3) ? 4 : 0;
 398
 399    si->main_data_begin = getbits(8);
 400    if (stereo == 1)
 401      si->private_bits = get1bit();
 402    else
 403      si->private_bits = getbits_fast(2);
 404
 405    for (ch=0; ch<stereo; ch++)
 406    {
 407        register struct gr_info_s *gr_info = &(si->ch[ch].gr[0]);
 408
 409        gr_info->part2_3_length = getbits(12);
 410        gr_info->big_values = getbits_fast(9);
 411        if(gr_info->big_values > 288) {
 412          fprintf(stderr,"big_values too large!\n");
 413          gr_info->big_values = 288;
 414        }
 415        gr_info->pow2gain = gainpow2+256 - getbits_fast(8) + powdiff;
 416        if(ms_stereo)
 417          gr_info->pow2gain += 2;
 418        gr_info->scalefac_compress = getbits(9);
 419 /* window-switching flag == 1 for block_Type != 0 .. and block-type == 0 -> win-sw-flag = 0 */
 420        if(get1bit())
 421        {
 422          int i;
 423          gr_info->block_type = getbits_fast(2);
 424          gr_info->mixed_block_flag = get1bit();
 425          gr_info->table_select[0] = getbits_fast(5);
 426          gr_info->table_select[1] = getbits_fast(5);
 427          /*
 428           * table_select[2] not needed, because there is no region2,
 429           * but to satisfy some verifications tools we set it either.
 430           */
 431          gr_info->table_select[2] = 0;
 432          for(i=0;i<3;i++)
 433            gr_info->full_gain[i] = gr_info->pow2gain + (getbits_fast(3)<<3);
 434
 435          if(gr_info->block_type == 0) {
 436            fprintf(stderr,"Blocktype == 0 and window-switching == 1 not allowed.\n");
 437            return 0;
 438          }
 439          /* region_count/start parameters are implicit in this case. */
 440 /* check this again! */
 441          if(gr_info->block_type == 2)
 442            gr_info->region1start = 36>>1;
 443          else if(sfreq == 8)
 444 /* check this for 2.5 and sfreq=8 */
 445            gr_info->region1start = 108>>1;
 446          else
 447            gr_info->region1start = 54>>1;
 448          gr_info->region2start = 576>>1;
 449        }
 450        else
 451        {
 452          int i,r0c,r1c;
 453          for (i=0; i<3; i++)
 454            gr_info->table_select[i] = getbits_fast(5);
 455          r0c = getbits_fast(4);
 456          r1c = getbits_fast(3);
 457          gr_info->region1start = bandInfo[sfreq].longIdx[r0c+1] >> 1 ;
 458          gr_info->region2start = bandInfo[sfreq].longIdx[r0c+1+r1c+1] >> 1;
 459          gr_info->block_type = 0;
 460          gr_info->mixed_block_flag = 0;
 461        }
 462        gr_info->scalefac_scale = get1bit();
 463        gr_info->count1table_select = get1bit();
 464    }
 465    return !0;
 466 }
 467
 468 /*
 469  * read scalefactors
 470  */
 471 #ifdef MPEG1
 472 static int III_get_scale_factors_1(int *scf,struct gr_info_s *gr_info)
 473 {
 474    static const unsigned char slen[2][16] = {
 475      {0, 0, 0, 0, 3, 1, 1, 1, 2, 2, 2, 3, 3, 3, 4, 4},
 476      {0, 1, 2, 3, 0, 1, 2, 3, 1, 2, 3, 1, 2, 3, 2, 3}
 477    };
 478    int numbits;
 479    int num0 = slen[0][gr_info->scalefac_compress];
 480    int num1 = slen[1][gr_info->scalefac_compress];
 481
 482     if (gr_info->block_type == 2) {
 483       int i=18;
 484       numbits = (num0 + num1) * 18;
 485
 486       if (gr_info->mixed_block_flag) {
 487          for (i=8;i;i--)
 488            *scf++ = getbits_fast(num0);
 489          i = 9;
 490          numbits -= num0; /* num0 * 17 + num1 * 18 */
 491       }
 492
 493       for (;i;i--)
 494         *scf++ = getbits_fast(num0);
 495       for (i = 18; i; i--)
 496         *scf++ = getbits_fast(num1);
 497       *scf++ = 0; *scf++ = 0; *scf++ = 0; /* short[13][0..2] = 0 */
 498     }
 499     else {
 500       int i;
 501       int scfsi = gr_info->scfsi;
 502
 503       if(scfsi < 0) { /* scfsi < 0 => granule == 0 */
 504          for(i=11;i;i--)
 505            *scf++ = getbits_fast(num0);
 506          for(i=10;i;i--)
 507            *scf++ = getbits_fast(num1);
 508          numbits = (num0 + num1) * 10 + num0;
 509          *scf++ = 0;
 510       }
 511       else {
 512         numbits = 0;
 513         if(!(scfsi & 0x8)) {
 514           for (i=0;i<6;i++)
 515             *scf++ = getbits_fast(num0);
 516           numbits += num0 * 6;
 517         }
 518         else {
 519           scf += 6;
 520         }
 521
 522         if(!(scfsi & 0x4)) {
 523           for (i=0;i<5;i++)
 524             *scf++ = getbits_fast(num0);
 525           numbits += num0 * 5;
 526         }
 527         else {
 528           scf += 5;
 529         }
 530
 531         if(!(scfsi & 0x2)) {
 532           for(i=0;i<5;i++)
 533             *scf++ = getbits_fast(num1);
 534           numbits += num1 * 5;
 535         }
 536         else {
 537           scf += 5;
 538         }
 539
 540         if(!(scfsi & 0x1)) {
 541           for (i=0;i<5;i++)
 542             *scf++ = getbits_fast(num1);
 543           numbits += num1 * 5;
 544         }
 545         else {
 546            scf += 5;
 547         }
 548         *scf++ = 0;  /* no l[21] in original sources */
 549       }
 550     }
 551     return numbits;
 552 }
 553 #endif
 554
 555
 556 static int III_get_scale_factors_2(int *scf,struct gr_info_s *gr_info,int i_stereo)
 557 {
 558   const unsigned char *pnt;
 559   int i,j;
 560   unsigned int slen;
 561   int n = 0;
 562   int numbits = 0;
 563
 564   static const unsigned char stab[3][6][4] = {
 565    { { 6, 5, 5,5 } , { 6, 5, 7,3 } , { 11,10,0,0} ,
 566      { 7, 7, 7,0 } , { 6, 6, 6,3 } , {  8, 8,5,0} } ,
 567    { { 9, 9, 9,9 } , { 9, 9,12,6 } , { 18,18,0,0} ,
 568      {12,12,12,0 } , {12, 9, 9,6 } , { 15,12,9,0} } ,
 569    { { 6, 9, 9,9 } , { 6, 9,12,6 } , { 15,18,0,0} ,
 570      { 6,15,12,0 } , { 6,12, 9,6 } , {  6,18,9,0} } };
 571
 572   if(i_stereo) /* i_stereo AND second channel -> do_layer3() checks this */
 573     slen = i_slen2[gr_info->scalefac_compress>>1];
 574   else
 575     slen = n_slen2[gr_info->scalefac_compress];
 576
 577   gr_info->preflag = (slen>>15) & 0x1;
 578
 579   n = 0;
 580   if( gr_info->block_type == 2 ) {
 581     n++;
 582     if(gr_info->mixed_block_flag)
 583       n++;
 584   }
 585
 586   pnt = stab[n][(slen>>12)&0x7];
 587
 588   for(i=0;i<4;i++) {
 589     int num = slen & 0x7;
 590     slen >>= 3;
 591     if(num) {
 592       for(j=0;j<(int)(pnt[i]);j++)
 593         *scf++ = getbits_fast(num);
 594       numbits += pnt[i] * num;
 595     }
 596     else {
 597       for(j=0;j<(int)(pnt[i]);j++)
 598         *scf++ = 0;
 599     }
 600   }
 601
 602   n = (n << 1) + 1;
 603   for(i=0;i<n;i++)
 604     *scf++ = 0;
 605
 606   return numbits;
 607 }
 608
 609 static const int pretab1[22] = {0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,2,2,3,3,3,2,0};
 610 static const int pretab2[22] = {0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0};
 611
 612 /*
 613  * don't forget to apply the same changes to III_dequantize_sample_ms() !!!
 614  */
 615 static int III_dequantize_sample(real xr[SBLIMIT][SSLIMIT],int *scf,
 616    struct gr_info_s *gr_info,int sfreq,int part2bits)
 617 {
 618   int shift = 1 + gr_info->scalefac_scale;
 619   real *xrpnt = (real *) xr;
 620   int l[3],l3;
 621   int part2remain = gr_info->part2_3_length - part2bits;
 622   int *me;
 623
 624   {
 625     int bv       = gr_info->big_values;
 626     int region1  = gr_info->region1start;
 627     int region2  = gr_info->region2start;
 628
 629     l3 = ((576>>1)-bv)>>1;
 630 /*
 631  * we may lose the 'odd' bit here !!
 632  * check this later again
 633  */
 634     if(bv <= region1) {
 635       l[0] = bv; l[1] = 0; l[2] = 0;
 636     }
 637     else {
 638       l[0] = region1;
 639       if(bv <= region2) {
 640         l[1] = bv - l[0];  l[2] = 0;
 641       }
 642       else {
 643         l[1] = region2 - l[0]; l[2] = bv - region2;
 644       }
 645     }
 646   }
 647
 648   if(gr_info->block_type == 2) {
 649     /*
 650      * decoding with short or mixed mode BandIndex table
 651      */
 652     int i,max[4];
 653     int step=0,lwin=0,cb=0;
 654     register real v = 0.0;
 655     register int *m,mc;
 656
 657     if(gr_info->mixed_block_flag) {
 658       max[3] = -1;
 659       max[0] = max[1] = max[2] = 2;
 660       m = map[sfreq][0];
 661       me = mapend[sfreq][0];
 662     }
 663     else {
 664       max[0] = max[1] = max[2] = max[3] = -1;
 665       /* max[3] not really needed in this case */
 666       m = map[sfreq][1];
 667       me = mapend[sfreq][1];
 668     }
 669
 670     mc = 0;
 671     for(i=0;i<2;i++) {
 672       int lp = l[i];
 673       const struct newhuff *h = ht+gr_info->table_select[i];
 674       for(;lp;lp--,mc--) {
 675         register int x,y;
 676         if( (!mc) ) {
 677           mc = *m++;
 678           xrpnt = ((real *) xr) + (*m++);
 679           lwin = *m++;
 680           cb = *m++;
 681           if(lwin == 3) {
 682             v = gr_info->pow2gain[(*scf++) << shift];
 683             step = 1;
 684           }
 685           else {
 686             v = gr_info->full_gain[lwin][(*scf++) << shift];
 687             step = 3;
 688           }
 689         }
 690         {
 691           register const short *val = h->table;
 692           while((y=*val++)<0) {
 693             if (get1bit())
 694               val -= y;
 695             part2remain--;
 696           }
 697           x = y >> 4;
 698           y &= 0xf;
 699         }
 700         if(x == 15) {
 701           max[lwin] = cb;
 702           part2remain -= h->linbits+1;
 703           x += getbits(h->linbits);
 704           if(get1bit())
 705             *xrpnt = -ispow[x] * v;
 706           else
 707             *xrpnt =  ispow[x] * v;
 708         }
 709         else if(x) {
 710           max[lwin] = cb;
 711           if(get1bit())
 712             *xrpnt = -ispow[x] * v;
 713           else
 714             *xrpnt =  ispow[x] * v;
 715           part2remain--;
 716         }
 717         else
 718           *xrpnt = 0.0;
 719         xrpnt += step;
 720         if(y == 15) {
 721           max[lwin] = cb;
 722           part2remain -= h->linbits+1;
 723           y += getbits(h->linbits);
 724           if(get1bit())
 725             *xrpnt = -ispow[y] * v;
 726           else
 727             *xrpnt =  ispow[y] * v;
 728         }
 729         else if(y) {
 730           max[lwin] = cb;
 731           if(get1bit())
 732             *xrpnt = -ispow[y] * v;
 733           else
 734             *xrpnt =  ispow[y] * v;
 735           part2remain--;
 736         }
 737         else
 738           *xrpnt = 0.0;
 739         xrpnt += step;
 740       }
 741     }
 742     for(;l3 && (part2remain > 0);l3--) {
 743       const struct newhuff *h = htc+gr_info->count1table_select;
 744       const short *val = h->table;
 745       short a;
 746
 747       while((a=*val++)<0) {
 748         part2remain--;
 749         if(part2remain < 0) {
 750           part2remain++;
 751           a = 0;
 752           break;
 753         }
 754         if (get1bit())
 755           val -= a;
 756       }
 757
 758       for(i=0;i<4;i++) {
 759         if(!(i & 1)) {
 760           if(!mc) {
 761             mc = *m++;
 762             xrpnt = ((real *) xr) + (*m++);
 763             lwin = *m++;
 764             cb = *m++;
 765             if(lwin == 3) {
 766               v = gr_info->pow2gain[(*scf++) << shift];
 767               step = 1;
 768             }
 769             else {
 770               v = gr_info->full_gain[lwin][(*scf++) << shift];
 771               step = 3;
 772             }
 773           }
 774           mc--;
 775         }
 776         if( (a & (0x8>>i)) ) {
 777           max[lwin] = cb;
 778           part2remain--;
 779           if(part2remain < 0) {
 780             part2remain++;
 781             break;
 782           }
 783           if(get1bit())
 784             *xrpnt = -v;
 785           else
 786             *xrpnt = v;
 787         }
 788         else
 789           *xrpnt = 0.0;
 790         xrpnt += step;
 791       }
 792     }
 793
 794     while( m < me ) {
 795       if(!mc) {
 796         mc = *m++;
 797         xrpnt = ((real *) xr) + *m++;
 798         if( (*m++) == 3)
 799           step = 1;
 800         else
 801           step = 3;
 802         m++; /* cb */
 803       }
 804       mc--;
 805       *xrpnt = 0.0;
 806       xrpnt += step;
 807       *xrpnt = 0.0;
 808       xrpnt += step;
 809 /* we could add a little opt. here:
 810  * if we finished a band for window 3 or a long band
 811  * further bands could copied in a simple loop without a
 812  * special 'map' decoding
 813  */
 814     }
 815
 816     gr_info->maxband[0] = max[0]+1;
 817     gr_info->maxband[1] = max[1]+1;
 818     gr_info->maxband[2] = max[2]+1;
 819     gr_info->maxbandl = max[3]+1;
 820
 821     {
 822       int rmax = max[0] > max[1] ? max[0] : max[1];
 823       rmax = (rmax > max[2] ? rmax : max[2]) + 1;
 824       gr_info->maxb = rmax ? shortLimit[sfreq][rmax] : longLimit[sfreq][max[3]+1];
 825     }
 826
 827   }
 828   else {
 829         /*
 830      * decoding with 'long' BandIndex table (block_type != 2)
 831      */
 832     const int *pretab = gr_info->preflag ? pretab1 : pretab2;
 833     int i,max = -1;
 834     int cb = 0;
 835     register int *m = map[sfreq][2];
 836     register real v = 0.0;
 837     register int mc = 0;
 838 #if 0
 839     me = mapend[sfreq][2];
 840 #endif
 841
 842         /*
 843      * long hash table values
 844      */
 845     for(i=0;i<3;i++) {
 846       int lp = l[i];
 847       const struct newhuff *h = ht+gr_info->table_select[i];
 848
 849       for(;lp;lp--,mc--) {
 850         int x,y;
 851
 852         if(!mc) {
 853           mc = *m++;
 854           v = gr_info->pow2gain[((*scf++) + (*pretab++)) << shift];
 855           cb = *m++;
 856         }
 857         {
 858           register const short *val = h->table;
 859           while((y=*val++)<0) {
 860             if (get1bit())
 861               val -= y;
 862             part2remain--;
 863           }
 864           x = y >> 4;
 865           y &= 0xf;
 866         }
 867         if (x == 15) {
 868           max = cb;
 869           part2remain -= h->linbits+1;
 870           x += getbits(h->linbits);
 871           if(get1bit())
 872             *xrpnt++ = -ispow[x] * v;
 873           else
 874             *xrpnt++ =  ispow[x] * v;
 875         }
 876         else if(x) {
 877           max = cb;
 878           if(get1bit())
 879             *xrpnt++ = -ispow[x] * v;
 880           else
 881             *xrpnt++ =  ispow[x] * v;
 882           part2remain--;
 883         }
 884         else
 885           *xrpnt++ = 0.0;
 886
 887         if (y == 15) {
 888           max = cb;
 889           part2remain -= h->linbits+1;
 890           y += getbits(h->linbits);
 891           if(get1bit())
 892             *xrpnt++ = -ispow[y] * v;
 893           else
 894             *xrpnt++ =  ispow[y] * v;
 895         }
 896         else if(y) {
 897           max = cb;
 898           if(get1bit())
 899             *xrpnt++ = -ispow[y] * v;
 900           else
 901             *xrpnt++ =  ispow[y] * v;
 902           part2remain--;
 903         }
 904         else
 905           *xrpnt++ = 0.0;
 906       }
 907     }
 908
 909         /*
 910      * short (count1table) values
 911      */
 912     for(;l3 && (part2remain > 0);l3--) {
 913       const struct newhuff *h = htc+gr_info->count1table_select;
 914       const short *val = h->table;
 915       short a;
 916
 917       while((a=*val++)<0) {
 918         part2remain--;
 919         if(part2remain < 0) {
 920           part2remain++;
 921           a = 0;
 922           break;
 923         }
 924         if (get1bit())
 925           val -= a;
 926       }
 927
 928       for(i=0;i<4;i++) {
 929         if(!(i & 1)) {
 930           if(!mc) {
 931             mc = *m++;
 932             cb = *m++;
 933             v = gr_info->pow2gain[((*scf++) + (*pretab++)) << shift];
 934           }
 935           mc--;
 936         }
 937         if ( (a & (0x8>>i)) ) {
 938           max = cb;
 939           part2remain--;
 940           if(part2remain < 0) {
 941             part2remain++;
 942             break;
 943           }
 944           if(get1bit())
 945             *xrpnt++ = -v;
 946           else
 947             *xrpnt++ = v;
 948         }
 949         else
 950           *xrpnt++ = 0.0;
 951       }
 952     }
 953
 954         /*
 955      * zero part
 956      */
 957     for(i=(&xr[SBLIMIT][0]-xrpnt)>>1;i;i--) {
 958       *xrpnt++ = 0.0;
 959       *xrpnt++ = 0.0;
 960     }
 961
 962     gr_info->maxbandl = max+1;
 963     gr_info->maxb = longLimit[sfreq][gr_info->maxbandl];
 964   }
 965
 966   while( part2remain > 16 ) {
 967     getbits(16); /* Dismiss stuffing Bits */
 968     part2remain -= 16;
 969   }
 970   if(part2remain > 0)
 971     getbits(part2remain);
 972   else if(part2remain < 0) {
 973     fprintf(stderr,"mpg123: Can't rewind stream by %d bits!\n",-part2remain);
 974     return 1; /* -> error */
 975   }
 976   return 0;
 977 }
 978
 979 #if 0
 980 static int III_dequantize_sample_ms(real xr[2][SBLIMIT][SSLIMIT],int *scf,
 981    struct gr_info_s *gr_info,int sfreq,int part2bits)
 982 {
 983   int shift = 1 + gr_info->scalefac_scale;
 984   real *xrpnt = (real *) xr[1];
 985   real *xr0pnt = (real *) xr[0];
 986   int l[3],l3;
 987   int part2remain = gr_info->part2_3_length - part2bits;
 988   int *me;
 989
 990   {
 991     int bv       = gr_info->big_values;
 992     int region1  = gr_info->region1start;
 993     int region2  = gr_info->region2start;
 994
 995     l3 = ((576>>1)-bv)>>1;
 996 /*
 997  * we may lose the 'odd' bit here !!
 998  * check this later gain
 999  */
1000     if(bv <= region1) {
1001       l[0] = bv; l[1] = 0; l[2] = 0;
1002     }
1003     else {
1004       l[0] = region1;
1005       if(bv <= region2) {
1006         l[1] = bv - l[0];  l[2] = 0;
1007       }
1008       else {
1009         l[1] = region2 - l[0]; l[2] = bv - region2;
1010       }
1011     }
1012   }
1013
1014   if(gr_info->block_type == 2) {
1015     int i,max[4];
1016     int step=0,lwin=0,cb=0;
1017     register real v = 0.0;
1018     register int *m,mc = 0;
1019
1020     if(gr_info->mixed_block_flag) {
1021       max[3] = -1;
1022       max[0] = max[1] = max[2] = 2;
1023       m = map[sfreq][0];
1024       me = mapend[sfreq][0];
1025     }
1026     else {
1027       max[0] = max[1] = max[2] = max[3] = -1;
1028       /* max[3] not really needed in this case */
1029       m = map[sfreq][1];
1030       me = mapend[sfreq][1];
1031     }
1032
1033     for(i=0;i<2;i++) {
1034       int lp = l[i];
1035       struct newhuff *h = ht+gr_info->table_select[i];
1036       for(;lp;lp--,mc--) {
1037         int x,y;
1038
1039         if(!mc) {
1040           mc = *m++;
1041           xrpnt = ((real *) xr[1]) + *m;
1042           xr0pnt = ((real *) xr[0]) + *m++;
1043           lwin = *m++;
1044           cb = *m++;
1045           if(lwin == 3) {
1046             v = gr_info->pow2gain[(*scf++) << shift];
1047             step = 1;
1048           }
1049           else {
1050             v = gr_info->full_gain[lwin][(*scf++) << shift];
1051             step = 3;
1052           }
1053         }
1054         {
1055           register short *val = h->table;
1056           while((y=*val++)<0) {
1057             if (get1bit())
1058               val -= y;
1059             part2remain--;
1060           }
1061           x = y >> 4;
1062           y &= 0xf;
1063         }
1064         if(x == 15) {
1065           max[lwin] = cb;
1066           part2remain -= h->linbits+1;
1067           x += getbits(h->linbits);
1068           if(get1bit()) {
1069             real a = ispow[x] * v;
1070             *xrpnt = *xr0pnt + a;
1071             *xr0pnt -= a;
1072           }
1073           else {
1074             real a = ispow[x] * v;
1075             *xrpnt = *xr0pnt - a;
1076             *xr0pnt += a;
1077           }
1078         }
1079         else if(x) {
1080           max[lwin] = cb;
1081           if(get1bit()) {
1082             real a = ispow[x] * v;
1083             *xrpnt = *xr0pnt + a;
1084             *xr0pnt -= a;
1085           }
1086           else {
1087             real a = ispow[x] * v;
1088             *xrpnt = *xr0pnt - a;
1089             *xr0pnt += a;
1090           }
1091           part2remain--;
1092         }
1093         else
1094           *xrpnt = *xr0pnt;
1095         xrpnt += step;
1096         xr0pnt += step;
1097
1098         if(y == 15) {
1099           max[lwin] = cb;
1100           part2remain -= h->linbits+1;
1101           y += getbits(h->linbits);
1102           if(get1bit()) {
1103             real a = ispow[y] * v;
1104             *xrpnt = *xr0pnt + a;
1105             *xr0pnt -= a;
1106           }
1107           else {
1108             real a = ispow[y] * v;
1109             *xrpnt = *xr0pnt - a;
1110             *xr0pnt += a;
1111           }
1112         }
1113         else if(y) {
1114           max[lwin] = cb;
1115           if(get1bit()) {
1116             real a = ispow[y] * v;
1117             *xrpnt = *xr0pnt + a;
1118             *xr0pnt -= a;
1119           }
1120           else {
1121             real a = ispow[y] * v;
1122             *xrpnt = *xr0pnt - a;
1123             *xr0pnt += a;
1124           }
1125           part2remain--;
1126         }
1127         else
1128           *xrpnt = *xr0pnt;
1129         xrpnt += step;
1130         xr0pnt += step;
1131       }
1132     }
1133
1134     for(;l3 && (part2remain > 0);l3--) {
1135       struct newhuff *h = htc+gr_info->count1table_select;
1136       register short *val = h->table,a;
1137
1138       while((a=*val++)<0) {
1139         part2remain--;
1140         if(part2remain < 0) {
1141           part2remain++;
1142           a = 0;
1143           break;
1144         }
1145         if (get1bit())
1146           val -= a;
1147       }
1148
1149       for(i=0;i<4;i++) {
1150         if(!(i & 1)) {
1151           if(!mc) {
1152             mc = *m++;
1153             xrpnt = ((real *) xr[1]) + *m;
1154             xr0pnt = ((real *) xr[0]) + *m++;
1155             lwin = *m++;
1156             cb = *m++;
1157             if(lwin == 3) {
1158               v = gr_info->pow2gain[(*scf++) << shift];
1159               step = 1;
1160             }
1161             else {
1162               v = gr_info->full_gain[lwin][(*scf++) << shift];
1163               step = 3;
1164             }
1165           }
1166           mc--;
1167         }
1168         if( (a & (0x8>>i)) ) {
1169           max[lwin] = cb;
1170           part2remain--;
1171           if(part2remain < 0) {
1172             part2remain++;
1173             break;
1174           }
1175           if(get1bit()) {
1176             *xrpnt = *xr0pnt + v;
1177             *xr0pnt -= v;
1178           }
1179           else {
1180             *xrpnt = *xr0pnt - v;
1181             *xr0pnt += v;
1182           }
1183         }
1184         else
1185           *xrpnt = *xr0pnt;
1186         xrpnt += step;
1187         xr0pnt += step;
1188       }
1189     }
1190
1191     while( m < me ) {
1192       if(!mc) {
1193         mc = *m++;
1194         xrpnt = ((real *) xr[1]) + *m;
1195         xr0pnt = ((real *) xr[0]) + *m++;
1196         if(*m++ == 3)
1197           step = 1;
1198         else
1199           step = 3;
1200         m++; /* cb */
1201       }
1202       mc--;
1203       *xrpnt = *xr0pnt;
1204       xrpnt += step;
1205       xr0pnt += step;
1206       *xrpnt = *xr0pnt;
1207       xrpnt += step;
1208       xr0pnt += step;
1209 /* we could add a little opt. here:
1210  * if we finished a band for window 3 or a long band
1211  * further bands could copied in a simple loop without a
1212  * special 'map' decoding
1213  */
1214     }
1215
1216     gr_info->maxband[0] = max[0]+1;
1217     gr_info->maxband[1] = max[1]+1;
1218     gr_info->maxband[2] = max[2]+1;
1219     gr_info->maxbandl = max[3]+1;
1220
1221     {
1222       int rmax = max[0] > max[1] ? max[0] : max[1];
1223       rmax = (rmax > max[2] ? rmax : max[2]) + 1;
1224       gr_info->maxb = rmax ? shortLimit[sfreq][rmax] : longLimit[sfreq][max[3]+1];
1225     }
1226   }
1227   else {
1228     int *pretab = gr_info->preflag ? pretab1 : pretab2;
1229     int i,max = -1;
1230     int cb = 0;
1231     register int mc=0,*m = map[sfreq][2];
1232     register real v = 0.0;
1233 #if 0
1234     me = mapend[sfreq][2];
1235 #endif
1236
1237     for(i=0;i<3;i++) {
1238       int lp = l[i];
1239       struct newhuff *h = ht+gr_info->table_select[i];
1240
1241       for(;lp;lp--,mc--) {
1242         int x,y;
1243         if(!mc) {
1244           mc = *m++;
1245           cb = *m++;
1246           v = gr_info->pow2gain[((*scf++) + (*pretab++)) << shift];
1247         }
1248         {
1249           register short *val = h->table;
1250           while((y=*val++)<0) {
1251             if (get1bit())
1252               val -= y;
1253             part2remain--;
1254           }
1255           x = y >> 4;
1256           y &= 0xf;
1257         }
1258         if (x == 15) {
1259           max = cb;
1260           part2remain -= h->linbits+1;
1261           x += getbits(h->linbits);
1262           if(get1bit()) {
1263             real a = ispow[x] * v;
1264             *xrpnt++ = *xr0pnt + a;
1265             *xr0pnt++ -= a;
1266           }
1267           else {
1268             real a = ispow[x] * v;
1269             *xrpnt++ = *xr0pnt - a;
1270             *xr0pnt++ += a;
1271           }
1272         }
1273         else if(x) {
1274           max = cb;
1275           if(get1bit()) {
1276             real a = ispow[x] * v;
1277             *xrpnt++ = *xr0pnt + a;
1278             *xr0pnt++ -= a;
1279           }
1280           else {
1281             real a = ispow[x] * v;
1282             *xrpnt++ = *xr0pnt - a;
1283             *xr0pnt++ += a;
1284           }
1285           part2remain--;
1286         }
1287         else
1288           *xrpnt++ = *xr0pnt++;
1289
1290         if (y == 15) {
1291           max = cb;
1292           part2remain -= h->linbits+1;
1293           y += getbits(h->linbits);
1294           if(get1bit()) {
1295             real a = ispow[y] * v;
1296             *xrpnt++ = *xr0pnt + a;
1297             *xr0pnt++ -= a;
1298           }
1299           else {
1300             real a = ispow[y] * v;
1301             *xrpnt++ = *xr0pnt - a;
1302             *xr0pnt++ += a;
1303           }
1304         }
1305         else if(y) {
1306           max = cb;
1307           if(get1bit()) {
1308             real a = ispow[y] * v;
1309             *xrpnt++ = *xr0pnt + a;
1310             *xr0pnt++ -= a;
1311           }
1312           else {
1313             real a = ispow[y] * v;
1314             *xrpnt++ = *xr0pnt - a;
1315             *xr0pnt++ += a;
1316           }
1317           part2remain--;
1318         }
1319         else
1320           *xrpnt++ = *xr0pnt++;
1321       }
1322     }
1323
1324     for(;l3 && (part2remain > 0);l3--) {
1325       struct newhuff *h = htc+gr_info->count1table_select;
1326       register short *val = h->table,a;
1327
1328       while((a=*val++)<0) {
1329         part2remain--;
1330         if(part2remain < 0) {
1331           part2remain++;
1332           a = 0;
1333           break;
1334         }
1335         if (get1bit())
1336           val -= a;
1337       }
1338
1339       for(i=0;i<4;i++) {
1340         if(!(i & 1)) {
1341           if(!mc) {
1342             mc = *m++;
1343             cb = *m++;
1344             v = gr_info->pow2gain[((*scf++) + (*pretab++)) << shift];
1345           }
1346           mc--;
1347         }
1348         if ( (a & (0x8>>i)) ) {
1349           max = cb;
1350           part2remain--;
1351           if(part2remain <= 0) {
1352             part2remain++;
1353             break;
1354           }
1355           if(get1bit()) {
1356             *xrpnt++ = *xr0pnt + v;
1357             *xr0pnt++ -= v;
1358           }
1359           else {
1360             *xrpnt++ = *xr0pnt - v;
1361             *xr0pnt++ += v;
1362           }
1363         }
1364         else
1365           *xrpnt++ = *xr0pnt++;
1366       }
1367     }
1368     for(i=(&xr[1][SBLIMIT][0]-xrpnt)>>1;i;i--) {
1369       *xrpnt++ = *xr0pnt++;
1370       *xrpnt++ = *xr0pnt++;
1371     }
1372
1373     gr_info->maxbandl = max+1;
1374     gr_info->maxb = longLimit[sfreq][gr_info->maxbandl];
1375   }
1376
1377   while ( part2remain > 16 ) {
1378     getbits(16); /* Dismiss stuffing Bits */
1379     part2remain -= 16;
1380   }
1381   if(part2remain > 0 )
1382     getbits(part2remain);
1383   else if(part2remain < 0) {
1384     fprintf(stderr,"mpg123_ms: Can't rewind stream by %d bits!\n",-part2remain);
1385     return 1; /* -> error */
1386   }
1387   return 0;
1388 }
1389 #endif
1390
1391 /*
1392  * III_stereo: calculate real channel values for Joint-I-Stereo-mode
1393  */
1394 static void III_i_stereo(real xr_buf[2][SBLIMIT][SSLIMIT],int *scalefac,
1395    struct gr_info_s *gr_info,int sfreq,int ms_stereo,int lsf)
1396 {
1397       real (*xr)[SBLIMIT*SSLIMIT] = (real (*)[SBLIMIT*SSLIMIT] ) xr_buf;
1398       const struct bandInfoStruct *bi = &bandInfo[sfreq];
1399       real *tab1,*tab2;
1400
1401       if(lsf) {
1402         int p = gr_info->scalefac_compress & 0x1;
1403             if(ms_stereo) {
1404           tab1 = pow1_2[p]; tab2 = pow2_2[p];
1405         }
1406         else {
1407           tab1 = pow1_1[p]; tab2 = pow2_1[p];
1408         }
1409       }
1410       else {
1411         if(ms_stereo) {
1412           tab1 = tan1_2; tab2 = tan2_2;
1413         }
1414         else {
1415           tab1 = tan1_1; tab2 = tan2_1;
1416         }
1417       }
1418
1419       if (gr_info->block_type == 2)
1420       {
1421          int lwin,do_l = 0;
1422          if( gr_info->mixed_block_flag )
1423            do_l = 1;
1424
1425          for (lwin=0;lwin<3;lwin++) /* process each window */
1426          {
1427              /* get first band with zero values */
1428            int is_p,sb,idx,sfb = gr_info->maxband[lwin];  /* sfb is minimal 3 for mixed mode */
1429            if(sfb > 3)
1430              do_l = 0;
1431
1432            for(;sfb<12;sfb++)
1433            {
1434              is_p = scalefac[sfb*3+lwin-gr_info->mixed_block_flag]; /* scale: 0-15 */
1435              if(is_p != 7) {
1436                real t1,t2;
1437                sb = bi->shortDiff[sfb];
1438                idx = bi->shortIdx[sfb] + lwin;
1439                t1 = tab1[is_p]; t2 = tab2[is_p];
1440                for (; sb > 0; sb--,idx+=3)
1441                {
1442                  real v = xr[0][idx];
1443                  xr[0][idx] = v * t1;
1444                  xr[1][idx] = v * t2;
1445                }
1446              }
1447            }
1448
1449 #if 1
1450 /* in the original: copy 10 to 11 , here: copy 11 to 12
1451 maybe still wrong??? (copy 12 to 13?) */
1452            is_p = scalefac[11*3+lwin-gr_info->mixed_block_flag]; /* scale: 0-15 */
1453            sb = bi->shortDiff[12];
1454            idx = bi->shortIdx[12] + lwin;
1455 #else
1456            is_p = scalefac[10*3+lwin-gr_info->mixed_block_flag]; /* scale: 0-15 */
1457            sb = bi->shortDiff[11];
1458            idx = bi->shortIdx[11] + lwin;
1459 #endif
1460            if(is_p != 7)
1461            {
1462              real t1,t2;
1463              t1 = tab1[is_p]; t2 = tab2[is_p];
1464              for ( ; sb > 0; sb--,idx+=3 )
1465              {
1466                real v = xr[0][idx];
1467                xr[0][idx] = v * t1;
1468                xr[1][idx] = v * t2;
1469              }
1470            }
1471          } /* end for(lwin; .. ; . ) */
1472
1473          if (do_l)
1474          {
1475 /* also check l-part, if ALL bands in the three windows are 'empty'
1476  * and mode = mixed_mode
1477  */
1478            int sfb = gr_info->maxbandl;
1479            int idx = bi->longIdx[sfb];
1480
1481            for ( ; sfb<8; sfb++ )
1482            {
1483              int sb = bi->longDiff[sfb];
1484              int is_p = scalefac[sfb]; /* scale: 0-15 */
1485              if(is_p != 7) {
1486                real t1,t2;
1487                t1 = tab1[is_p]; t2 = tab2[is_p];
1488                for ( ; sb > 0; sb--,idx++)
1489                {
1490                  real v = xr[0][idx];
1491                  xr[0][idx] = v * t1;
1492                  xr[1][idx] = v * t2;
1493                }
1494              }
1495              else
1496                idx += sb;
1497            }
1498          }
1499       }
1500       else /* ((gr_info->block_type != 2)) */
1501       {
1502         int sfb = gr_info->maxbandl;
1503         int is_p,idx = bi->longIdx[sfb];
1504         for ( ; sfb<21; sfb++)
1505         {
1506           int sb = bi->longDiff[sfb];
1507           is_p = scalefac[sfb]; /* scale: 0-15 */
1508           if(is_p != 7) {
1509             real t1,t2;
1510             t1 = tab1[is_p]; t2 = tab2[is_p];
1511             for ( ; sb > 0; sb--,idx++)
1512             {
1513                real v = xr[0][idx];
1514                xr[0][idx] = v * t1;
1515                xr[1][idx] = v * t2;
1516             }
1517           }
1518           else
1519             idx += sb;
1520         }
1521
1522         is_p = scalefac[20]; /* copy l-band 20 to l-band 21 */
1523         if(is_p != 7)
1524         {
1525           int sb;
1526           real t1 = tab1[is_p],t2 = tab2[is_p];
1527
1528           for ( sb = bi->longDiff[21]; sb > 0; sb--,idx++ )
1529           {
1530             real v = xr[0][idx];
1531             xr[0][idx] = v * t1;
1532             xr[1][idx] = v * t2;
1533           }
1534         }
1535       } /* ... */
1536 }
1537
1538 static void III_antialias(real xr[SBLIMIT][SSLIMIT],struct gr_info_s *gr_info)
1539 {
1540    int sblim;
1541
1542    if(gr_info->block_type == 2)
1543    {
1544       if(!gr_info->mixed_block_flag)
1545         return;
1546       sblim = 1;
1547    }
1548    else {
1549      sblim = gr_info->maxb-1;
1550    }
1551
1552    /* 31 alias-reduction operations between each pair of sub-bands */
1553    /* with 8 butterflies between each pair                         */
1554
1555    {
1556      int sb;
1557      real *xr1=(real *) xr[1];
1558
1559      for(sb=sblim;sb;sb--,xr1+=10)
1560      {
1561        int ss;
1562        real *cs=aa_cs,*ca=aa_ca;
1563        real *xr2 = xr1;
1564
1565        for(ss=7;ss>=0;ss--)
1566        {       /* upper and lower butterfly inputs */
1567          register real bu = *--xr2,bd = *xr1;
1568          *xr2   = (bu * (*cs)   ) - (bd * (*ca)   );
1569          *xr1++ = (bd * (*cs++) ) + (bu * (*ca++) );
1570        }
1571      }
1572   }
1573 }
1574
1575 /*
1576  DCT insipired by Jeff Tsay's DCT from the maplay package
1577  this is an optimized version with manual unroll.
1578
1579  References:
1580  [1] S. Winograd: "On Computing the Discrete Fourier Transform",
1581      Mathematics of Computation, Volume 32, Number 141, January 1978,
1582      Pages 175-199
1583 */
1584
1585 static void dct36(real *inbuf,real *o1,real *o2,real *wintab,real *tsbuf)
1586 {
1587   {
1588     register real *in = inbuf;
1589
1590     in[17]+=in[16]; in[16]+=in[15]; in[15]+=in[14];
1591     in[14]+=in[13]; in[13]+=in[12]; in[12]+=in[11];
1592     in[11]+=in[10]; in[10]+=in[9];  in[9] +=in[8];
1593     in[8] +=in[7];  in[7] +=in[6];  in[6] +=in[5];
1594     in[5] +=in[4];  in[4] +=in[3];  in[3] +=in[2];
1595     in[2] +=in[1];  in[1] +=in[0];
1596
1597     in[17]+=in[15]; in[15]+=in[13]; in[13]+=in[11]; in[11]+=in[9];
1598     in[9] +=in[7];  in[7] +=in[5];  in[5] +=in[3];  in[3] +=in[1];
1599
1600
1601   {
1602
1603 #define MACRO0(v) { \
1604     real tmp; \
1605     out2[9+(v)] = (tmp = sum0 + sum1) * w[27+(v)]; \
1606     out2[8-(v)] = tmp * w[26-(v)];  } \
1607     sum0 -= sum1; \
1608     ts[SBLIMIT*(8-(v))] = out1[8-(v)] + sum0 * w[8-(v)]; \
1609     ts[SBLIMIT*(9+(v))] = out1[9+(v)] + sum0 * w[9+(v)];
1610 #define MACRO1(v) { \
1611         real sum0,sum1; \
1612     sum0 = tmp1a + tmp2a; \
1613         sum1 = (tmp1b + tmp2b) * tfcos36[(v)]; \
1614         MACRO0(v); }
1615 #define MACRO2(v) { \
1616     real sum0,sum1; \
1617     sum0 = tmp2a - tmp1a; \
1618     sum1 = (tmp2b - tmp1b) * tfcos36[(v)]; \
1619         MACRO0(v); }
1620
1621     register const real *c = COS9;
1622     register real *out2 = o2;
1623         register real *w = wintab;
1624         register real *out1 = o1;
1625         register real *ts = tsbuf;
1626
1627     real ta33,ta66,tb33,tb66;
1628
1629     ta33 = in[2*3+0] * c[3];
1630     ta66 = in[2*6+0] * c[6];
1631     tb33 = in[2*3+1] * c[3];
1632     tb66 = in[2*6+1] * c[6];
1633
1634     {
1635       real tmp1a,tmp2a,tmp1b,tmp2b;
1636       tmp1a =             in[2*1+0] * c[1] + ta33 + in[2*5+0] * c[5] + in[2*7+0] * c[7];
1637       tmp1b =             in[2*1+1] * c[1] + tb33 + in[2*5+1] * c[5] + in[2*7+1] * c[7];
1638       tmp2a = in[2*0+0] + in[2*2+0] * c[2] + in[2*4+0] * c[4] + ta66 + in[2*8+0] * c[8];
1639       tmp2b = in[2*0+1] + in[2*2+1] * c[2] + in[2*4+1] * c[4] + tb66 + in[2*8+1] * c[8];
1640
1641       MACRO1(0);
1642       MACRO2(8);
1643     }
1644
1645     {
1646       real tmp1a,tmp2a,tmp1b,tmp2b;
1647       tmp1a = ( in[2*1+0] - in[2*5+0] - in[2*7+0] ) * c[3];
1648       tmp1b = ( in[2*1+1] - in[2*5+1] - in[2*7+1] ) * c[3];
1649       tmp2a = ( in[2*2+0] - in[2*4+0] - in[2*8+0] ) * c[6] - in[2*6+0] + in[2*0+0];
1650       tmp2b = ( in[2*2+1] - in[2*4+1] - in[2*8+1] ) * c[6] - in[2*6+1] + in[2*0+1];
1651
1652       MACRO1(1);
1653       MACRO2(7);
1654     }
1655
1656     {
1657       real tmp1a,tmp2a,tmp1b,tmp2b;
1658       tmp1a =             in[2*1+0] * c[5] - ta33 - in[2*5+0] * c[7] + in[2*7+0] * c[1];
1659       tmp1b =             in[2*1+1] * c[5] - tb33 - in[2*5+1] * c[7] + in[2*7+1] * c[1];
1660       tmp2a = in[2*0+0] - in[2*2+0] * c[8] - in[2*4+0] * c[2] + ta66 + in[2*8+0] * c[4];
1661       tmp2b = in[2*0+1] - in[2*2+1] * c[8] - in[2*4+1] * c[2] + tb66 + in[2*8+1] * c[4];
1662
1663       MACRO1(2);
1664       MACRO2(6);
1665     }
1666
1667     {
1668       real tmp1a,tmp2a,tmp1b,tmp2b;
1669       tmp1a =             in[2*1+0] * c[7] - ta33 + in[2*5+0] * c[1] - in[2*7+0] * c[5];
1670       tmp1b =             in[2*1+1] * c[7] - tb33 + in[2*5+1] * c[1] - in[2*7+1] * c[5];
1671       tmp2a = in[2*0+0] - in[2*2+0] * c[4] + in[2*4+0] * c[8] + ta66 - in[2*8+0] * c[2];
1672       tmp2b = in[2*0+1] - in[2*2+1] * c[4] + in[2*4+1] * c[8] + tb66 - in[2*8+1] * c[2];
1673
1674       MACRO1(3);
1675       MACRO2(5);
1676     }
1677
1678         {
1679                 real sum0,sum1;
1680         sum0 =  in[2*0+0] - in[2*2+0] + in[2*4+0] - in[2*6+0] + in[2*8+0];
1681         sum1 = (in[2*0+1] - in[2*2+1] + in[2*4+1] - in[2*6+1] + in[2*8+1] ) * tfcos36[4];
1682                 MACRO0(4);
1683         }
1684   }
1685
1686   }
1687 }
1688
1689 /*
1690  * new DCT12
1691  */
1692 static void dct12(real *in,real *rawout1,real *rawout2,register real *wi,register real *ts)
1693 {
1694 #define DCT12_PART1 \
1695              in5 = in[5*3];  \
1696      in5 += (in4 = in[4*3]); \
1697      in4 += (in3 = in[3*3]); \
1698      in3 += (in2 = in[2*3]); \
1699      in2 += (in1 = in[1*3]); \
1700      in1 += (in0 = in[0*3]); \
1701                              \
1702      in5 += in3; in3 += in1; \
1703                              \
1704      in2 *= COS6_1; \
1705      in3 *= COS6_1; \
1706
1707 #define DCT12_PART2 \
1708      in0 += in4 * COS6_2; \
1709                           \
1710      in4 = in0 + in2;     \
1711      in0 -= in2;          \
1712                           \
1713      in1 += in5 * COS6_2; \
1714                           \
1715      in5 = (in1 + in3) * tfcos12[0]; \
1716      in1 = (in1 - in3) * tfcos12[2]; \
1717                          \
1718      in3 = in4 + in5;    \
1719      in4 -= in5;         \
1720                          \
1721      in2 = in0 + in1;    \
1722      in0 -= in1;
1723
1724
1725    {
1726      real in0,in1,in2,in3,in4,in5;
1727      register real *out1 = rawout1;
1728      ts[SBLIMIT*0] = out1[0]; ts[SBLIMIT*1] = out1[1]; ts[SBLIMIT*2] = out1[2];
1729      ts[SBLIMIT*3] = out1[3]; ts[SBLIMIT*4] = out1[4]; ts[SBLIMIT*5] = out1[5];
1730
1731      DCT12_PART1
1732
1733      {
1734        real tmp0,tmp1 = (in0 - in4);
1735        {
1736          real tmp2 = (in1 - in5) * tfcos12[1];
1737          tmp0 = tmp1 + tmp2;
1738          tmp1 -= tmp2;
1739        }
1740        ts[(17-1)*SBLIMIT] = out1[17-1] + tmp0 * wi[11-1];
1741        ts[(12+1)*SBLIMIT] = out1[12+1] + tmp0 * wi[6+1];
1742        ts[(6 +1)*SBLIMIT] = out1[6 +1] + tmp1 * wi[1];
1743        ts[(11-1)*SBLIMIT] = out1[11-1] + tmp1 * wi[5-1];
1744      }
1745
1746      DCT12_PART2
1747
1748      ts[(17-0)*SBLIMIT] = out1[17-0] + in2 * wi[11-0];
1749      ts[(12+0)*SBLIMIT] = out1[12+0] + in2 * wi[6+0];
1750      ts[(12+2)*SBLIMIT] = out1[12+2] + in3 * wi[6+2];
1751      ts[(17-2)*SBLIMIT] = out1[17-2] + in3 * wi[11-2];
1752
1753      ts[(6+0)*SBLIMIT]  = out1[6+0] + in0 * wi[0];
1754      ts[(11-0)*SBLIMIT] = out1[11-0] + in0 * wi[5-0];
1755      ts[(6+2)*SBLIMIT]  = out1[6+2] + in4 * wi[2];
1756      ts[(11-2)*SBLIMIT] = out1[11-2] + in4 * wi[5-2];
1757   }
1758
1759   in++;
1760
1761   {
1762      real in0,in1,in2,in3,in4,in5;
1763      register real *out2 = rawout2;
1764
1765      DCT12_PART1
1766
1767      {
1768        real tmp0,tmp1 = (in0 - in4);
1769        {
1770          real tmp2 = (in1 - in5) * tfcos12[1];
1771          tmp0 = tmp1 + tmp2;
1772          tmp1 -= tmp2;
1773        }
1774        out2[5-1] = tmp0 * wi[11-1];
1775        out2[0+1] = tmp0 * wi[6+1];
1776        ts[(12+1)*SBLIMIT] += tmp1 * wi[1];
1777        ts[(17-1)*SBLIMIT] += tmp1 * wi[5-1];
1778      }
1779
1780      DCT12_PART2
1781
1782      out2[5-0] = in2 * wi[11-0];
1783      out2[0+0] = in2 * wi[6+0];
1784      out2[0+2] = in3 * wi[6+2];
1785      out2[5-2] = in3 * wi[11-2];
1786
1787      ts[(12+0)*SBLIMIT] += in0 * wi[0];
1788      ts[(17-0)*SBLIMIT] += in0 * wi[5-0];
1789      ts[(12+2)*SBLIMIT] += in4 * wi[2];
1790      ts[(17-2)*SBLIMIT] += in4 * wi[5-2];
1791   }
1792
1793   in++;
1794
1795   {
1796      real in0,in1,in2,in3,in4,in5;
1797      register real *out2 = rawout2;
1798      out2[12]=out2[13]=out2[14]=out2[15]=out2[16]=out2[17]=0.0;
1799
1800      DCT12_PART1
1801
1802      {
1803        real tmp0,tmp1 = (in0 - in4);
1804        {
1805          real tmp2 = (in1 - in5) * tfcos12[1];
1806          tmp0 = tmp1 + tmp2;
1807          tmp1 -= tmp2;
1808        }
1809        out2[11-1] = tmp0 * wi[11-1];
1810        out2[6 +1] = tmp0 * wi[6+1];
1811        out2[0+1] += tmp1 * wi[1];
1812        out2[5-1] += tmp1 * wi[5-1];
1813      }
1814
1815      DCT12_PART2
1816
1817      out2[11-0] = in2 * wi[11-0];
1818      out2[6 +0] = in2 * wi[6+0];
1819      out2[6 +2] = in3 * wi[6+2];
1820      out2[11-2] = in3 * wi[11-2];
1821
1822      out2[0+0] += in0 * wi[0];
1823      out2[5-0] += in0 * wi[5-0];
1824      out2[0+2] += in4 * wi[2];
1825      out2[5-2] += in4 * wi[5-2];
1826   }
1827 }
1828
1829 /*
1830  * III_hybrid
1831  */
1832 static void III_hybrid(real fsIn[SBLIMIT][SSLIMIT],real tsOut[SSLIMIT][SBLIMIT],
1833    int ch,struct gr_info_s *gr_info)
1834 {
1835    real *tspnt = (real *) tsOut;
1836    real (*block)[2][SBLIMIT*SSLIMIT] = gmp->hybrid_block;
1837    int *blc = gmp->hybrid_blc;
1838    real *rawout1,*rawout2;
1839    int bt;
1840    int sb = 0;
1841
1842    {
1843      int b = blc[ch];
1844      rawout1=block[b][ch];
1845      b=-b+1;
1846      rawout2=block[b][ch];
1847      blc[ch] = b;
1848    }
1849
1850
1851    if(gr_info->mixed_block_flag) {
1852      sb = 2;
1853      dct36(fsIn[0],rawout1,rawout2,win[0],tspnt);
1854      dct36(fsIn[1],rawout1+18,rawout2+18,win1[0],tspnt+1);
1855      rawout1 += 36; rawout2 += 36; tspnt += 2;
1856    }
1857
1858    bt = gr_info->block_type;
1859    if(bt == 2) {
1860      for (; sb<gr_info->maxb; sb+=2,tspnt+=2,rawout1+=36,rawout2+=36) {
1861        dct12(fsIn[sb],rawout1,rawout2,win[2],tspnt);
1862        dct12(fsIn[sb+1],rawout1+18,rawout2+18,win1[2],tspnt+1);
1863      }
1864    }
1865    else {
1866      for (; sb<gr_info->maxb; sb+=2,tspnt+=2,rawout1+=36,rawout2+=36) {
1867        dct36(fsIn[sb],rawout1,rawout2,win[bt],tspnt);
1868        dct36(fsIn[sb+1],rawout1+18,rawout2+18,win1[bt],tspnt+1);
1869      }
1870    }
1871
1872    for(;sb<SBLIMIT;sb++,tspnt++) {
1873      int i;
1874      for(i=0;i<SSLIMIT;i++) {
1875        tspnt[i*SBLIMIT] = *rawout1++;
1876        *rawout2++ = 0.0;
1877      }
1878    }
1879 }
1880
1881 /*
1882  * main layer3 handler
1883  */
1884 int do_layer3(struct frame *fr,unsigned char *pcm_sample,int *pcm_point)
1885 {
1886   int gr, ch, ss,clip=0;
1887   int scalefacs[2][39]; /* max 39 for short[13][3] mode, mixed: 38, long: 22 */
1888   struct III_sideinfo sideinfo;
1889   int stereo = fr->stereo;
1890   int single = fr->single;
1891   int ms_stereo,i_stereo;
1892   int sfreq = fr->sampling_frequency;
1893   int stereo1,granules;
1894
1895   if(stereo == 1) { /* stream is mono */
1896     stereo1 = 1;
1897     single = 0;
1898   }
1899   else if(single >= 0) /* stream is stereo, but force to mono */
1900     stereo1 = 1;
1901   else
1902     stereo1 = 2;
1903
1904   if(fr->mode == MPG_MD_JOINT_STEREO) {
1905     ms_stereo = fr->mode_ext & 0x2;
1906     i_stereo  = fr->mode_ext & 0x1;
1907   }
1908   else
1909     ms_stereo = i_stereo = 0;
1910
1911   if(fr->lsf) {
1912     granules = 1;
1913     if(!III_get_side_info_2(&sideinfo,stereo,ms_stereo,sfreq,single))
1914       return -1;
1915   }
1916   else {
1917     granules = 2;
1918 #ifdef MPEG1
1919     if(!III_get_side_info_1(&sideinfo,stereo,ms_stereo,sfreq,single))
1920       return -1;
1921 #else
1922     fprintf(stderr,"Not supported\n");
1923 #endif
1924   }
1925
1926   if(set_pointer(sideinfo.main_data_begin) == MP3_ERR)
1927     return -1;
1928
1929   for (gr=0;gr<granules;gr++)
1930   {
1931     real hybridIn[2][SBLIMIT][SSLIMIT];
1932     real hybridOut[2][SSLIMIT][SBLIMIT];
1933
1934     {
1935       struct gr_info_s *gr_info = &(sideinfo.ch[0].gr[gr]);
1936       long part2bits;
1937       if(fr->lsf)
1938         part2bits = III_get_scale_factors_2(scalefacs[0],gr_info,0);
1939       else {
1940 #ifdef MPEG1
1941         part2bits = III_get_scale_factors_1(scalefacs[0],gr_info);
1942 #else
1943         fprintf(stderr,"Not supported\n");
1944 #endif
1945       }
1946       if(III_dequantize_sample(hybridIn[0], scalefacs[0],gr_info,sfreq,part2bits))
1947         return clip;
1948     }
1949     if(stereo == 2) {
1950       struct gr_info_s *gr_info = &(sideinfo.ch[1].gr[gr]);
1951       long part2bits;
1952       if(fr->lsf)
1953         part2bits = III_get_scale_factors_2(scalefacs[1],gr_info,i_stereo);
1954       else {
1955 #ifdef MPEG1
1956         part2bits = III_get_scale_factors_1(scalefacs[1],gr_info);
1957 #else
1958         fprintf(stderr,"Not supported\n");
1959 #endif
1960       }
1961
1962       if(III_dequantize_sample(hybridIn[1],scalefacs[1],gr_info,sfreq,part2bits))
1963           return clip;
1964
1965       if(ms_stereo) {
1966         int i;
1967         for(i=0;i<SBLIMIT*SSLIMIT;i++) {
1968           real tmp0,tmp1;
1969           tmp0 = ((real *) hybridIn[0])[i];
1970           tmp1 = ((real *) hybridIn[1])[i];
1971           ((real *) hybridIn[0])[i] = tmp0 + tmp1;
1972           ((real *) hybridIn[1])[i] = tmp0 - tmp1;
1973         }
1974       }
1975
1976       if(i_stereo)
1977         III_i_stereo(hybridIn,scalefacs[1],gr_info,sfreq,ms_stereo,fr->lsf);
1978
1979       if(ms_stereo || i_stereo || (single == 3) ) {
1980         if(gr_info->maxb > sideinfo.ch[0].gr[gr].maxb)
1981           sideinfo.ch[0].gr[gr].maxb = gr_info->maxb;
1982         else
1983           gr_info->maxb = sideinfo.ch[0].gr[gr].maxb;
1984       }
1985
1986       switch(single) {
1987         case 3:
1988           {
1989             register int i;
1990             register real *in0 = (real *) hybridIn[0],*in1 = (real *) hybridIn[1];
1991             for(i=0;i<SSLIMIT*gr_info->maxb;i++,in0++)
1992               *in0 = (*in0 + *in1++); /* *0.5 done by pow-scale */
1993           }
1994           break;
1995         case 1:
1996           {
1997             register int i;
1998             register real *in0 = (real *) hybridIn[0],*in1 = (real *) hybridIn[1];
1999             for(i=0;i<SSLIMIT*gr_info->maxb;i++)
2000               *in0++ = *in1++;
2001           }
2002           break;
2003       }
2004     }
2005
2006     for(ch=0;ch<stereo1;ch++) {
2007       struct gr_info_s *gr_info = &(sideinfo.ch[ch].gr[gr]);
2008       III_antialias(hybridIn[ch],gr_info);
2009       III_hybrid(hybridIn[ch], hybridOut[ch], ch,gr_info);
2010     }
2011
2012     for(ss=0;ss<SSLIMIT;ss++) {
2013       if(single >= 0) {
2014         clip += synth_1to1_mono(hybridOut[0][ss],pcm_sample,pcm_point);
2015       }
2016       else {
2017         int p1 = *pcm_point;
2018         clip += synth_1to1(hybridOut[0][ss],0,pcm_sample,&p1);
2019         clip += synth_1to1(hybridOut[1][ss],1,pcm_sample,pcm_point);
2020       }
2021     }
2022   }
2023
2024   return clip;
2025 }