开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 通讯/手机编程 > Symbian > 查看源码
winb.c
资源名称：helix.src.0812.rar　[点击查看]
上传用户：zhongxx05
上传日期：2007-06-06
资源大小：33641k
文件大小：10k
源码类别：

Symbian

开发平台：

C/C++

winb.c：源码内容
  1. /* ***** BEGIN LICENSE BLOCK ***** 
  2.  * Version: RCSL 1.0/RPSL 1.0 
  3.  *  
  4.  * Portions Copyright (c) 1995-2002 RealNetworks, Inc. All Rights Reserved. 
  5.  *      
  6.  * The contents of this file, and the files included with this file, are 
  7.  * subject to the current version of the RealNetworks Public Source License 
  8.  * Version 1.0 (the "RPSL") available at 
  9.  * http://www.helixcommunity.org/content/rpsl unless you have licensed 
  10.  * the file under the RealNetworks Community Source License Version 1.0 
  11.  * (the "RCSL") available at http://www.helixcommunity.org/content/rcsl, 
  12.  * in which case the RCSL will apply. You may also obtain the license terms 
  13.  * directly from RealNetworks.  You may not use this file except in 
  14.  * compliance with the RPSL or, if you have a valid RCSL with RealNetworks 
  15.  * applicable to this file, the RCSL.  Please see the applicable RPSL or 
  16.  * RCSL for the rights, obligations and limitations governing use of the 
  17.  * contents of the file.  
  18.  *  
  19.  * This file is part of the Helix DNA Technology. RealNetworks is the 
  20.  * developer of the Original Code and owns the copyrights in the portions 
  21.  * it created. 
  22.  *  
  23.  * This file, and the files included with this file, is distributed and made 
  24.  * available on an 'AS IS' basis, WITHOUT WARRANTY OF ANY KIND, EITHER 
  25.  * EXPRESS OR IMPLIED, AND REALNETWORKS HEREBY DISCLAIMS ALL SUCH WARRANTIES, 
  26.  * INCLUDING WITHOUT LIMITATION, ANY WARRANTIES OF MERCHANTABILITY, FITNESS 
  27.  * FOR A PARTICULAR PURPOSE, QUIET ENJOYMENT OR NON-INFRINGEMENT. 
  28.  * 
  29.  * Technology Compatibility Kit Test Suite(s) Location: 
  30.  *    http://www.helixcommunity.org/content/tck 
  31.  * 
  32.  * Contributor(s): 
  33.  *  
  34.  * ***** END LICENSE BLOCK ***** */ 
  36. void windowB(float *vbuf, int vb_ptr, unsigned char *pcm)
  37. {
  38. int i, j;
  39. int si, bx;
  40. float *coef;
  41. float sum;
  42. long tmp;
  44. si = vb_ptr + 16;
  45. bx = (si+32) & 511;
  46. coef = wincoef;
  48. /*-- first 16 --*/
  49. for(i=0;i<16;i++) {
  50.    sum = 0.0F;
  51.    for(j=0;j<8;j++) {
  52.       sum += (*coef++)*vbuf[si];
  53.       si = (si+64) & 511;
  54.       sum -= (*coef++)*vbuf[bx];
  55.       bx = (bx+64) & 511;
  56.    }
  57.    si++;
  58.    bx--;
  59.    tmp = (long)sum;
  60.    if( tmp > 32767 )        tmp =  32767;
  61.    else if( tmp < -32768 )  tmp = -32768;
  62.    *pcm++ = ((unsigned char)(tmp>>8))^0x80;
  63. }
  64. /*--  special case --*/
  65. sum = 0.0F;
  66. for(j=0;j<8;j++) {
  67.    sum += (*coef++)*vbuf[bx];
  68.    bx = (bx+64) & 511;
  69. }
  70. tmp = (long)sum;
  71. if( tmp > 32767 )        tmp =  32767;
  72. else if( tmp < -32768 )  tmp = -32768;
  73. *pcm++ = ((unsigned char)(tmp>>8))^0x80;
  74. /*-- last 15 --*/
  75. coef = wincoef + 255;    /* back pass through coefs */
  76. for(i=0;i<15;i++) {
  77.    si--;
  78.    bx++;
  79.    sum = 0.0F;
  80.    for(j=0;j<8;j++) {
  81.       sum += (*coef--)*vbuf[si];
  82.       si = (si+64) & 511;
  83.       sum += (*coef--)*vbuf[bx];
  84.       bx = (bx+64) & 511;
  85.    }
  86.    tmp = (long)sum;
  87.    if( tmp > 32767 )        tmp =  32767;
  88.    else if( tmp < -32768 )  tmp = -32768;
  89.    *pcm++ = ((unsigned char)(tmp>>8))^0x80;
  90. }
  91. }
  92. /*------------------------------------------------------------*/
  93. void windowB_dual(float *vbuf, int vb_ptr, unsigned char *pcm)
  94. {
  95. int i, j;       /* dual window interleaves output */
  96. int si, bx;
  97. float *coef;
  98. float sum;
  99. long tmp;
  101. si = vb_ptr + 16;
  102. bx = (si+32) & 511;
  103. coef = wincoef;
  105. /*-- first 16 --*/
  106. for(i=0;i<16;i++) {
  107.    sum = 0.0F;
  108.    for(j=0;j<8;j++) {
  109.       sum += (*coef++)*vbuf[si];
  110.       si = (si+64) & 511;
  111.       sum -= (*coef++)*vbuf[bx];
  112.       bx = (bx+64) & 511;
  113.    }
  114.    si++;
  115.    bx--;
  116.    tmp = (long)sum;
  117.    if( tmp > 32767 )        tmp =  32767;
  118.    else if( tmp < -32768 )  tmp = -32768;
  119.    *pcm = ((unsigned char)(tmp>>8))^0x80;
  120.    pcm += 2;
  121. }
  122. /*--  special case --*/
  123. sum = 0.0F;
  124. for(j=0;j<8;j++) {
  125.    sum += (*coef++)*vbuf[bx];
  126.    bx = (bx+64) & 511;
  127. }
  128. tmp = (long)sum;
  129. if( tmp > 32767 )        tmp =  32767;
  130. else if( tmp < -32768 )  tmp = -32768;
  131. *pcm = ((unsigned char)(tmp>>8))^0x80;
  132. pcm += 2;
  133. /*-- last 15 --*/
  134. coef = wincoef + 255;    /* back pass through coefs */
  135. for(i=0;i<15;i++) {
  136.    si--;
  137.    bx++;
  138.    sum = 0.0F;
  139.    for(j=0;j<8;j++) {
  140.       sum += (*coef--)*vbuf[si];
  141.       si = (si+64) & 511;
  142.       sum += (*coef--)*vbuf[bx];
  143.       bx = (bx+64) & 511;
  144.    }
  145.    tmp = (long)sum;
  146.    if( tmp > 32767 )        tmp =  32767;
  147.    else if( tmp < -32768 )  tmp = -32768;
  148.    *pcm = ((unsigned char)(tmp>>8))^0x80;
  149.    pcm += 2;
  150. }
  151. }
  152. /*------------------------------------------------------------*/
  153. #ifdef REDUCTION
  154. /*------------------------------------------------------------*/
  155. /*------------------- 16 pt window ------------------------------*/
  156. void windowB16(float *vbuf, int vb_ptr, unsigned char *pcm)
  157. {
  158. int i, j;
  159. unsigned char si, bx;
  160. float *coef;
  161. float sum;
  162. long tmp;
  164. si = vb_ptr + 8;
  165. bx = si+16;
  166. coef = wincoef;
  168. /*-- first 8 --*/
  169. for(i=0;i<8;i++) {
  170.    sum = 0.0F;
  171.    for(j=0;j<8;j++) {
  172.       sum += (*coef++)*vbuf[si];
  173.       si += 32;
  174.       sum -= (*coef++)*vbuf[bx];
  175.       bx += 32;
  176.    }
  177.    si++;
  178.    bx--;
  179.    coef+=16;
  180.    tmp = (long)sum;
  181.    if( tmp > 32767 )        tmp =  32767;
  182.    else if( tmp < -32768 )  tmp = -32768;
  183.    *pcm++ = ((unsigned char)(tmp>>8))^0x80;
  184. }
  185. /*--  special case --*/
  186. sum = 0.0F;
  187. for(j=0;j<8;j++) {
  188.    sum += (*coef++)*vbuf[bx];
  189.    bx += 32;
  190. }
  191. tmp = (long)sum;
  192. if( tmp > 32767 )        tmp =  32767;
  193. else if( tmp < -32768 )  tmp = -32768;
  194. *pcm++ = ((unsigned char)(tmp>>8))^0x80;
  195. /*-- last 7 --*/
  196. coef = wincoef + 255;    /* back pass through coefs */
  197. for(i=0;i<7;i++) {
  198.    coef-=16;
  199.    si--;
  200.    bx++;
  201.    sum = 0.0F;
  202.    for(j=0;j<8;j++) {
  203.       sum += (*coef--)*vbuf[si];
  204.       si +=32;
  205.       sum += (*coef--)*vbuf[bx];
  206.       bx +=32;
  207.    }
  208.    tmp = (long)sum;
  209.    if( tmp > 32767 )        tmp =  32767;
  210.    else if( tmp < -32768 )  tmp = -32768;
  211.    *pcm++ = ((unsigned char)(tmp>>8))^0x80;
  212. }
  213. }
  214. /*--------------- 16 pt dual window (interleaved output) -----------------*/
  215. void windowB16_dual(float *vbuf, int vb_ptr, unsigned char *pcm)
  216. {
  217. int i, j;
  218. unsigned char si, bx;
  219. float *coef;
  220. float sum;
  221. long tmp;
  223. si = vb_ptr + 8;
  224. bx = si+16;
  225. coef = wincoef;
  227. /*-- first 8 --*/
  228. for(i=0;i<8;i++) {
  229.    sum = 0.0F;
  230.    for(j=0;j<8;j++) {
  231.       sum += (*coef++)*vbuf[si];
  232.       si += 32;
  233.       sum -= (*coef++)*vbuf[bx];
  234.       bx += 32;
  235.    }
  236.    si++;
  237.    bx--;
  238.    coef+=16;
  239.    tmp = (long)sum;
  240.    if( tmp > 32767 )        tmp =  32767;
  241.    else if( tmp < -32768 )  tmp = -32768;
  242.    *pcm = ((unsigned char)(tmp>>8))^0x80;
  243.    pcm+=2;
  244. }
  245. /*--  special case --*/
  246. sum = 0.0F;
  247. for(j=0;j<8;j++) {
  248.    sum += (*coef++)*vbuf[bx];
  249.    bx += 32;
  250. }
  251. tmp = (long)sum;
  252. if( tmp > 32767 )        tmp =  32767;
  253. else if( tmp < -32768 )  tmp = -32768;
  254. *pcm = ((unsigned char)(tmp>>8))^0x80;
  255. pcm+=2;
  256. /*-- last 7 --*/
  257. coef = wincoef + 255;    /* back pass through coefs */
  258. for(i=0;i<7;i++) {
  259.    coef-=16;
  260.    si--;
  261.    bx++;
  262.    sum = 0.0F;
  263.    for(j=0;j<8;j++) {
  264.       sum += (*coef--)*vbuf[si];
  265.       si +=32;
  266.       sum += (*coef--)*vbuf[bx];
  267.       bx +=32;
  268.    }
  269.    tmp = (long)sum;
  270.    if( tmp > 32767 )        tmp =  32767;
  271.    else if( tmp < -32768 )  tmp = -32768;
  272.    *pcm = ((unsigned char)(tmp>>8))^0x80;
  273.    pcm+=2;
  274. }
  275. }
  276. /*------------------- 8 pt window ------------------------------*/
  277. void windowB8(float *vbuf, int vb_ptr, unsigned char *pcm)
  278. {
  279. int i, j;
  280. int si, bx;
  281. float *coef;
  282. float sum;
  283. long tmp;
  285. si = vb_ptr + 4;
  286. bx = (si+8)&127;
  287. coef = wincoef;
  289. /*-- first 4 --*/
  290. for(i=0;i<4;i++) {
  291.    sum = 0.0F;
  292.    for(j=0;j<8;j++) {
  293.       sum += (*coef++)*vbuf[si];
  294.       si = (si+16) & 127;
  295.       sum -= (*coef++)*vbuf[bx];
  296.       bx = (bx+16) & 127;
  297.    }
  298.    si++;
  299.    bx--;
  300.    coef+=48;
  301.    tmp = (long)sum;
  302.    if( tmp > 32767 )        tmp =  32767;
  303.    else if( tmp < -32768 )  tmp = -32768;
  304.    *pcm++ = ((unsigned char)(tmp>>8))^0x80;
  305. }
  306. /*--  special case --*/
  307. sum = 0.0F;
  308. for(j=0;j<8;j++) {
  309.    sum += (*coef++)*vbuf[bx];
  310.    bx = (bx+16) & 127;
  311. }
  312. tmp = (long)sum;
  313. if( tmp > 32767 )        tmp =  32767;
  314. else if( tmp < -32768 )  tmp = -32768;
  315. *pcm++ = ((unsigned char)(tmp>>8))^0x80;
  316. /*-- last 3 --*/
  317. coef = wincoef + 255;    /* back pass through coefs */
  318. for(i=0;i<3;i++) {
  319.    coef-=48;
  320.    si--;
  321.    bx++;
  322.    sum = 0.0F;
  323.    for(j=0;j<8;j++) {
  324.       sum += (*coef--)*vbuf[si];
  325.       si = (si+16) & 127;
  326.       sum += (*coef--)*vbuf[bx];
  327.       bx = (bx+16) & 127;
  328.    }
  329.    tmp = (long)sum;
  330.    if( tmp > 32767 )        tmp =  32767;
  331.    else if( tmp < -32768 )  tmp = -32768;
  332.    *pcm++ = ((unsigned char)(tmp>>8))^0x80;
  333. }
  334. }
  335. /*--------------- 8 pt dual window (interleaved output) -----------------*/
  336. void windowB8_dual(float *vbuf, int vb_ptr, unsigned char *pcm)
  337. {
  338. int i, j;
  339. int si, bx;
  340. float *coef;
  341. float sum;
  342. long tmp;
  344. si = vb_ptr + 4;
  345. bx = (si+8) & 127;
  346. coef = wincoef;
  348. /*-- first 4 --*/
  349. for(i=0;i<4;i++) {
  350.    sum = 0.0F;
  351.    for(j=0;j<8;j++) {
  352.       sum += (*coef++)*vbuf[si];
  353.       si = (si+16) & 127;
  354.       sum -= (*coef++)*vbuf[bx];
  355.       bx = (bx+16) & 127;
  356.    }
  357.    si++;
  358.    bx--;
  359.    coef+=48;
  360.    tmp = (long)sum;
  361.    if( tmp > 32767 )        tmp =  32767;
  362.    else if( tmp < -32768 )  tmp = -32768;
  363.    *pcm = ((unsigned char)(tmp>>8))^0x80;
  364.    pcm+=2;
  365. }
  366. /*--  special case --*/
  367. sum = 0.0F;
  368. for(j=0;j<8;j++) {
  369.    sum += (*coef++)*vbuf[bx];
  370.    bx = (bx+16) & 127;
  371. }
  372. tmp = (long)sum;
  373. if( tmp > 32767 )        tmp =  32767;
  374. else if( tmp < -32768 )  tmp = -32768;
  375. *pcm = ((unsigned char)(tmp>>8))^0x80;
  376. pcm+=2;
  377. /*-- last 3 --*/
  378. coef = wincoef + 255;    /* back pass through coefs */
  379. for(i=0;i<3;i++) {
  380.    coef-=48;
  381.    si--;
  382.    bx++;
  383.    sum = 0.0F;
  384.    for(j=0;j<8;j++) {
  385.       sum += (*coef--)*vbuf[si];
  386.       si = (si+16) & 127;
  387.       sum += (*coef--)*vbuf[bx];
  388.       bx = (bx+16) & 127;
  389.    }
  390.    tmp = (long)sum;
  391.    if( tmp > 32767 )        tmp =  32767;
  392.    else if( tmp < -32768 )  tmp = -32768;
  393.    *pcm = ((unsigned char)(tmp>>8))^0x80;
  394.    pcm+=2;
  395. }
  396. }
  397. /*------------------------------------------------------------*/
  398. #endif      // end reduction