winb.c

著名的 helix realplayer 基于手机 symbian 系统的 播放器全套源代码

/* ***** BEGIN LICENSE BLOCK ***** 
 * Version: RCSL 1.0/RPSL 1.0 
 *  
 * Portions Copyright (c) 1995-2002 RealNetworks, Inc. All Rights Reserved. 
 *      
 * The contents of this file, and the files included with this file, are 
 * subject to the current version of the RealNetworks Public Source License 
 * Version 1.0 (the "RPSL") available at 
 * http://www.helixcommunity.org/content/rpsl unless you have licensed 
 * the file under the RealNetworks Community Source License Version 1.0 
 * (the "RCSL") available at http://www.helixcommunity.org/content/rcsl, 
 * in which case the RCSL will apply. You may also obtain the license terms 
 * directly from RealNetworks.  You may not use this file except in 
 * compliance with the RPSL or, if you have a valid RCSL with RealNetworks 
 * applicable to this file, the RCSL.  Please see the applicable RPSL or 
 * RCSL for the rights, obligations and limitations governing use of the 
 * contents of the file.  
 *  
 * This file is part of the Helix DNA Technology. RealNetworks is the 
 * developer of the Original Code and owns the copyrights in the portions 
 * it created. 
 *  
 * This file, and the files included with this file, is distributed and made 
 * available on an 'AS IS' basis, WITHOUT WARRANTY OF ANY KIND, EITHER 
 * EXPRESS OR IMPLIED, AND REALNETWORKS HEREBY DISCLAIMS ALL SUCH WARRANTIES, 
 * INCLUDING WITHOUT LIMITATION, ANY WARRANTIES OF MERCHANTABILITY, FITNESS 
 * FOR A PARTICULAR PURPOSE, QUIET ENJOYMENT OR NON-INFRINGEMENT. 
 * 
 * Technology Compatibility Kit Test Suite(s) Location: 
 *    http://www.helixcommunity.org/content/tck 
 * 
 * Contributor(s): 
 *  
 * ***** END LICENSE BLOCK ***** */ 

void windowB(float *vbuf, int vb_ptr, unsigned char *pcm)
{
int i, j;
int si, bx;
float *coef;
float sum;
long tmp;

si = vb_ptr + 16;
bx = (si+32) & 511;
coef = wincoef;

/*-- first 16 --*/
for(i=0;i<16;i++) {
   sum = 0.0F;
   for(j=0;j<8;j++) {
      sum += (*coef++)*vbuf[si];
      si = (si+64) & 511;
      sum -= (*coef++)*vbuf[bx];
      bx = (bx+64) & 511;
   }
   si++;
   bx--;
   tmp = (long)sum;
   if( tmp > 32767 )        tmp =  32767;
   else if( tmp < -32768 )  tmp = -32768;
   *pcm++ = ((unsigned char)(tmp>>8))^0x80;
}
/*--  special case --*/
sum = 0.0F;
for(j=0;j<8;j++) {
   sum += (*coef++)*vbuf[bx];
   bx = (bx+64) & 511;
}
tmp = (long)sum;
if( tmp > 32767 )        tmp =  32767;
else if( tmp < -32768 )  tmp = -32768;
*pcm++ = ((unsigned char)(tmp>>8))^0x80;
/*-- last 15 --*/
coef = wincoef + 255;    /* back pass through coefs */
for(i=0;i<15;i++) {
   si--;
   bx++;
   sum = 0.0F;
   for(j=0;j<8;j++) {
      sum += (*coef--)*vbuf[si];
      si = (si+64) & 511;
      sum += (*coef--)*vbuf[bx];
      bx = (bx+64) & 511;
   }
   tmp = (long)sum;
   if( tmp > 32767 )        tmp =  32767;
   else if( tmp < -32768 )  tmp = -32768;
   *pcm++ = ((unsigned char)(tmp>>8))^0x80;
}
}
/*------------------------------------------------------------*/
void windowB_dual(float *vbuf, int vb_ptr, unsigned char *pcm)
{
int i, j;       /* dual window interleaves output */
int si, bx;
float *coef;
float sum;
long tmp;

si = vb_ptr + 16;
bx = (si+32) & 511;
coef = wincoef;

/*-- first 16 --*/
for(i=0;i<16;i++) {
   sum = 0.0F;
   for(j=0;j<8;j++) {
      sum += (*coef++)*vbuf[si];
      si = (si+64) & 511;
      sum -= (*coef++)*vbuf[bx];
      bx = (bx+64) & 511;
   }
   si++;
   bx--;
   tmp = (long)sum;
   if( tmp > 32767 )        tmp =  32767;
   else if( tmp < -32768 )  tmp = -32768;
   *pcm = ((unsigned char)(tmp>>8))^0x80;
   pcm += 2;
}
/*--  special case --*/
sum = 0.0F;
for(j=0;j<8;j++) {
   sum += (*coef++)*vbuf[bx];
   bx = (bx+64) & 511;
}
tmp = (long)sum;
if( tmp > 32767 )        tmp =  32767;
else if( tmp < -32768 )  tmp = -32768;
*pcm = ((unsigned char)(tmp>>8))^0x80;
pcm += 2;
/*-- last 15 --*/
coef = wincoef + 255;    /* back pass through coefs */
for(i=0;i<15;i++) {
   si--;
   bx++;
   sum = 0.0F;
   for(j=0;j<8;j++) {
      sum += (*coef--)*vbuf[si];
      si = (si+64) & 511;
      sum += (*coef--)*vbuf[bx];
      bx = (bx+64) & 511;
   }
   tmp = (long)sum;
   if( tmp > 32767 )        tmp =  32767;
   else if( tmp < -32768 )  tmp = -32768;
   *pcm = ((unsigned char)(tmp>>8))^0x80;
   pcm += 2;
}
}
/*------------------------------------------------------------*/
#ifdef REDUCTION
/*------------------------------------------------------------*/
/*------------------- 16 pt window ------------------------------*/
void windowB16(float *vbuf, int vb_ptr, unsigned char *pcm)
{
int i, j;
unsigned char si, bx;
float *coef;
float sum;
long tmp;

si = vb_ptr + 8;
bx = si+16;
coef = wincoef;

/*-- first 8 --*/
for(i=0;i<8;i++) {
   sum = 0.0F;
   for(j=0;j<8;j++) {
      sum += (*coef++)*vbuf[si];
      si += 32;
      sum -= (*coef++)*vbuf[bx];
      bx += 32;
   }
   si++;
   bx--;
   coef+=16;
   tmp = (long)sum;
   if( tmp > 32767 )        tmp =  32767;
   else if( tmp < -32768 )  tmp = -32768;
   *pcm++ = ((unsigned char)(tmp>>8))^0x80;
}
/*--  special case --*/
sum = 0.0F;
for(j=0;j<8;j++) {
   sum += (*coef++)*vbuf[bx];
   bx += 32;
}
tmp = (long)sum;
if( tmp > 32767 )        tmp =  32767;
else if( tmp < -32768 )  tmp = -32768;
*pcm++ = ((unsigned char)(tmp>>8))^0x80;
/*-- last 7 --*/
coef = wincoef + 255;    /* back pass through coefs */
for(i=0;i<7;i++) {
   coef-=16;
   si--;
   bx++;
   sum = 0.0F;
   for(j=0;j<8;j++) {
      sum += (*coef--)*vbuf[si];
      si +=32;
      sum += (*coef--)*vbuf[bx];
      bx +=32;
   }
   tmp = (long)sum;
   if( tmp > 32767 )        tmp =  32767;
   else if( tmp < -32768 )  tmp = -32768;
   *pcm++ = ((unsigned char)(tmp>>8))^0x80;
}
}
/*--------------- 16 pt dual window (interleaved output) -----------------*/
void windowB16_dual(float *vbuf, int vb_ptr, unsigned char *pcm)
{
int i, j;
unsigned char si, bx;
float *coef;
float sum;
long tmp;

si = vb_ptr + 8;
bx = si+16;
coef = wincoef;

/*-- first 8 --*/
for(i=0;i<8;i++) {
   sum = 0.0F;
   for(j=0;j<8;j++) {
      sum += (*coef++)*vbuf[si];
      si += 32;
      sum -= (*coef++)*vbuf[bx];
      bx += 32;
   }
   si++;
   bx--;
   coef+=16;
   tmp = (long)sum;
   if( tmp > 32767 )        tmp =  32767;
   else if( tmp < -32768 )  tmp = -32768;
   *pcm = ((unsigned char)(tmp>>8))^0x80;
   pcm+=2;
}
/*--  special case --*/
sum = 0.0F;
for(j=0;j<8;j++) {
   sum += (*coef++)*vbuf[bx];
   bx += 32;
}
tmp = (long)sum;
if( tmp > 32767 )        tmp =  32767;
else if( tmp < -32768 )  tmp = -32768;
*pcm = ((unsigned char)(tmp>>8))^0x80;
pcm+=2;
/*-- last 7 --*/
coef = wincoef + 255;    /* back pass through coefs */
for(i=0;i<7;i++) {
   coef-=16;
   si--;
   bx++;
   sum = 0.0F;
   for(j=0;j<8;j++) {
      sum += (*coef--)*vbuf[si];
      si +=32;
      sum += (*coef--)*vbuf[bx];
      bx +=32;
   }
   tmp = (long)sum;
   if( tmp > 32767 )        tmp =  32767;
   else if( tmp < -32768 )  tmp = -32768;
   *pcm = ((unsigned char)(tmp>>8))^0x80;
   pcm+=2;
}
}
/*------------------- 8 pt window ------------------------------*/
void windowB8(float *vbuf, int vb_ptr, unsigned char *pcm)
{
int i, j;
int si, bx;
float *coef;
float sum;
long tmp;

si = vb_ptr + 4;
bx = (si+8)&127;
coef = wincoef;

/*-- first 4 --*/
for(i=0;i<4;i++) {
   sum = 0.0F;
   for(j=0;j<8;j++) {
      sum += (*coef++)*vbuf[si];
      si = (si+16) & 127;
      sum -= (*coef++)*vbuf[bx];
      bx = (bx+16) & 127;
   }
   si++;
   bx--;
   coef+=48;
   tmp = (long)sum;
   if( tmp > 32767 )        tmp =  32767;
   else if( tmp < -32768 )  tmp = -32768;
   *pcm++ = ((unsigned char)(tmp>>8))^0x80;
}
/*--  special case --*/
sum = 0.0F;
for(j=0;j<8;j++) {
   sum += (*coef++)*vbuf[bx];
   bx = (bx+16) & 127;
}
tmp = (long)sum;
if( tmp > 32767 )        tmp =  32767;
else if( tmp < -32768 )  tmp = -32768;
*pcm++ = ((unsigned char)(tmp>>8))^0x80;
/*-- last 3 --*/
coef = wincoef + 255;    /* back pass through coefs */
for(i=0;i<3;i++) {
   coef-=48;
   si--;
   bx++;
   sum = 0.0F;
   for(j=0;j<8;j++) {
      sum += (*coef--)*vbuf[si];
      si = (si+16) & 127;
      sum += (*coef--)*vbuf[bx];
      bx = (bx+16) & 127;
   }
   tmp = (long)sum;
   if( tmp > 32767 )        tmp =  32767;
   else if( tmp < -32768 )  tmp = -32768;
   *pcm++ = ((unsigned char)(tmp>>8))^0x80;
}
}
/*--------------- 8 pt dual window (interleaved output) -----------------*/
void windowB8_dual(float *vbuf, int vb_ptr, unsigned char *pcm)
{
int i, j;
int si, bx;
float *coef;
float sum;
long tmp;

si = vb_ptr + 4;
bx = (si+8) & 127;
coef = wincoef;

/*-- first 4 --*/
for(i=0;i<4;i++) {
   sum = 0.0F;
   for(j=0;j<8;j++) {
      sum += (*coef++)*vbuf[si];
      si = (si+16) & 127;
      sum -= (*coef++)*vbuf[bx];
      bx = (bx+16) & 127;
   }
   si++;
   bx--;
   coef+=48;
   tmp = (long)sum;
   if( tmp > 32767 )        tmp =  32767;
   else if( tmp < -32768 )  tmp = -32768;
   *pcm = ((unsigned char)(tmp>>8))^0x80;
   pcm+=2;
}
/*--  special case --*/
sum = 0.0F;
for(j=0;j<8;j++) {
   sum += (*coef++)*vbuf[bx];
   bx = (bx+16) & 127;
}
tmp = (long)sum;
if( tmp > 32767 )        tmp =  32767;
else if( tmp < -32768 )  tmp = -32768;
*pcm = ((unsigned char)(tmp>>8))^0x80;
pcm+=2;
/*-- last 3 --*/
coef = wincoef + 255;    /* back pass through coefs */
for(i=0;i<3;i++) {
   coef-=48;
   si--;
   bx++;
   sum = 0.0F;
   for(j=0;j<8;j++) {
      sum += (*coef--)*vbuf[si];
      si = (si+16) & 127;
      sum += (*coef--)*vbuf[bx];
      bx = (bx+16) & 127;
   }
   tmp = (long)sum;
   if( tmp > 32767 )        tmp =  32767;
   else if( tmp < -32768 )  tmp = -32768;
   *pcm = ((unsigned char)(tmp>>8))^0x80;
   pcm+=2;
}
}
/*------------------------------------------------------------*/
#endif      // end reduction