rpcrt4: Use RegDeleteTree to properly unregister interfaces.
[wine] / 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 "config.h"
24 #include <stdlib.h>
25 #include "wine/debug.h"
26 #include "mpg123.h"
27 #include "mpglib.h"
28 #include "huffman.h"
29
30 #define MPEG1
31 WINE_DEFAULT_DEBUG_CHANNEL(mpeg3);
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           FIXME("big_values (%d) too large!\n", gr_info->big_values);
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            FIXME("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          FIXME("big_values(%d) too large!\n", gr_info->big_values);
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            FIXME("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     FIXME("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     FIXME("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 inspired 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,struct mpstr *mp)
1834 {
1835    real *tspnt = (real *) tsOut;
1836    real (*block)[2][SBLIMIT*SSLIMIT] = mp->hybrid_block;
1837    int *blc = mp->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     FIXME("Not supported\n");
1923 #endif
1924   }
1925
1926   if(set_pointer(fr->mp,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         FIXME("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         FIXME("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, fr->mp);
2010     }
2011
2012     for(ss=0;ss<SSLIMIT;ss++) {
2013       if(single >= 0) {
2014         clip += synth_1to1_mono(fr->mp,hybridOut[0][ss],pcm_sample,pcm_point);
2015       }
2016       else {
2017         int p1 = *pcm_point;
2018         clip += synth_1to1(fr->mp,hybridOut[0][ss],0,pcm_sample,&p1);
2019         clip += synth_1to1(fr->mp,hybridOut[1][ss],1,pcm_sample,pcm_point);
2020       }
2021     }
2022   }
2023
2024   return clip;
2025 }