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

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

Symbian

开发平台：

C/C++

hwin.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. #include "statname.h"
  38. #ifdef __cplusplus
  39. extern "C" {
  40. #endif
  41. int hybrid(float xin[], float xprev[], float y[18][32],
  42.             int btype, int nlong, int ntot, int nprev, int band_limit_nsb);
  43. int hybrid_sum(float xin[], float xin_left[], float y[18][32],
  44.             int btype, int nlong, int ntot);
  45. void sum_f_bands( float a[], float b[], int n);
  46. void FreqInvert(float y[18][32], int n);
  47. #ifdef __cplusplus
  48. }
  49. #endif
  51. typedef float ARRAY36[36];
  52. /*-- windows by block type --*
  54. pi = 4.0*atan(1.0);
  56. ** type 0 **
  57. for(i=0;i<36;i++)  win[0][i] = (float)sin( pi/36 *(i+0.5) );
  59. ** type 1**
  60. for(i=0;i<18;i++)  win[1][i] = (float)sin( pi/36 *(i+0.5) );
  61. for(i=18;i<24;i++) win[1][i] = 1.0F;
  62. for(i=24;i<30;i++) win[1][i] = (float)sin( pi/12 *(i+0.5-18) );
  63. for(i=30;i<36;i++) win[1][i] = 0.0F;
  65. ** type 3**
  66. for(i=0;i<6;i++)   win[3][i] = 0.0F;
  67. for(i=6;i<12;i++)  win[3][i] = (float)sin( pi/12 *(i+0.5-6) );
  68. for(i=12;i<18;i++) win[3][i] = 1.0F;
  69. for(i=18;i<36;i++) win[3][i] = (float)sin( pi/36*(i+0.5) );
  71. ** type 2**
  72. for(i=0;i<12;i++)  win[2][i] = (float)sin( pi/12*(i+0.5) ) ;
  73. for(i=12;i<36;i++) win[2][i] = 0.0F;
  75. **--- invert signs by region to match mdct 18pt --> 36pt mapping **
  76. for(j=0;j<4;j++) {
  77.     if( j == 2 ) continue;
  78.     for(i=9;i<36;i++ ) win[j][i] = -win[j][i];
  79. }
  81. **-- invert signs for short blocks --**
  82. for(i=3;i<12;i++ ) win[2][i] = -win[2][i];
  83. */
  84. static const float win[4][36] = {
  85. {
  86. 4.3619386852e-002f,  1.3052618504e-001f,  2.1643961966e-001f,  3.0070579052e-001f, 
  87.  3.8268342614e-001f,  4.6174860001e-001f,  5.3729963303e-001f,  6.0876142979e-001f, 
  88.  6.7559021711e-001f, -7.3727732897e-001f, -7.9335331917e-001f, -8.4339141846e-001f, 
  89. -8.8701081276e-001f, -9.2387950420e-001f, -9.5371693373e-001f, -9.7629600763e-001f, 
  90. -9.9144488573e-001f, -9.9904823303e-001f, -9.9904823303e-001f, -9.9144488573e-001f, 
  91. -9.7629600763e-001f, -9.5371693373e-001f, -9.2387950420e-001f, -8.8701081276e-001f, 
  92. -8.4339141846e-001f, -7.9335331917e-001f, -7.3727732897e-001f, -6.7559021711e-001f, 
  93. -6.0876142979e-001f, -5.3729963303e-001f, -4.6174860001e-001f, -3.8268342614e-001f, 
  94. -3.0070579052e-001f, -2.1643961966e-001f, -1.3052618504e-001f, -4.3619386852e-002f, 
  95. },
  96. {
  97.  4.3619386852e-002f,  1.3052618504e-001f,  2.1643961966e-001f,  3.0070579052e-001f, 
  98.  3.8268342614e-001f,  4.6174860001e-001f,  5.3729963303e-001f,  6.0876142979e-001f, 
  99.  6.7559021711e-001f, -7.3727732897e-001f, -7.9335331917e-001f, -8.4339141846e-001f, 
  100. -8.8701081276e-001f, -9.2387950420e-001f, -9.5371693373e-001f, -9.7629600763e-001f, 
  101. -9.9144488573e-001f, -9.9904823303e-001f, -1.0000000000e+000f, -1.0000000000e+000f, 
  102. -1.0000000000e+000f, -1.0000000000e+000f, -1.0000000000e+000f, -1.0000000000e+000f, 
  103. -9.9144488573e-001f, -9.2387950420e-001f, -7.9335331917e-001f, -6.0876142979e-001f, 
  104. -3.8268342614e-001f, -1.3052618504e-001f,  0.0000000000e+000f,  0.0000000000e+000f, 
  105.  0.0000000000e+000f,  0.0000000000e+000f,  0.0000000000e+000f,  0.0000000000e+000f, 
  106. },
  107. {
  108.  1.3052618504e-001f,  3.8268342614e-001f,  6.0876142979e-001f, -7.9335331917e-001f, 
  109. -9.2387950420e-001f, -9.9144488573e-001f, -9.9144488573e-001f, -9.2387950420e-001f, 
  110. -7.9335331917e-001f, -6.0876142979e-001f, -3.8268342614e-001f, -1.3052618504e-001f, 
  111.  0.0000000000e+000f,  0.0000000000e+000f,  0.0000000000e+000f,  0.0000000000e+000f, 
  112.  0.0000000000e+000f,  0.0000000000e+000f,  0.0000000000e+000f,  0.0000000000e+000f, 
  113.  0.0000000000e+000f,  0.0000000000e+000f,  0.0000000000e+000f,  0.0000000000e+000f, 
  114.  0.0000000000e+000f,  0.0000000000e+000f,  0.0000000000e+000f,  0.0000000000e+000f, 
  115.  0.0000000000e+000f,  0.0000000000e+000f,  0.0000000000e+000f,  0.0000000000e+000f, 
  116.  0.0000000000e+000f,  0.0000000000e+000f,  0.0000000000e+000f,  0.0000000000e+000f, 
  117.  },
  118.  {
  119.  0.0000000000e+000f,  0.0000000000e+000f,  0.0000000000e+000f,  0.0000000000e+000f, 
  120.  0.0000000000e+000f,  0.0000000000e+000f,  1.3052618504e-001f,  3.8268342614e-001f, 
  121.  6.0876142979e-001f, -7.9335331917e-001f, -9.2387950420e-001f, -9.9144488573e-001f, 
  122. -1.0000000000e+000f, -1.0000000000e+000f, -1.0000000000e+000f, -1.0000000000e+000f, 
  123. -1.0000000000e+000f, -1.0000000000e+000f, -9.9904823303e-001f, -9.9144488573e-001f, 
  124. -9.7629600763e-001f, -9.5371693373e-001f, -9.2387950420e-001f, -8.8701081276e-001f, 
  125. -8.4339141846e-001f, -7.9335331917e-001f, -7.3727732897e-001f, -6.7559021711e-001f, 
  126. -6.0876142979e-001f, -5.3729963303e-001f, -4.6174860001e-001f, -3.8268342614e-001f, 
  127. -3.0070579052e-001f, -2.1643961966e-001f, -1.3052618504e-001f, -4.3619386852e-002f, 
  128. }
  129. };
  132. /*====================================================================*/
  133. void imdct18(float f[]);  /* 18 point */
  134. void imdct6_3(float f[]);  /* 6 point */
  135. /*====================================================================*/
  136. int hybrid(float xin[], float xprev[], float y[18][32],
  137.             int btype, int nlong, int ntot, int nprev, int band_limit_nsb)
  138. {
  139. int i, j;
  140. float *x, *x0;
  141. float xa, xb;
  142. int n;
  143. int nout;
  147. if( btype == 2 ) btype = 0;
  148. x = xin;
  149. x0 = xprev;
  151. /*-- do long blocks (if any) --*/
  152. n = (nlong+17)/18;       /* number of dct's to do */
  153. for(i=0;i<n;i++) {
  154.     imdct18(x);
  155.     for(j=0;j<9;j++) {
  156.         y[j][i]    = x0[j]   + win[btype][j]*x[9+j];
  157.         y[9+j][i]  = x0[9+j] + win[btype][9+j]*x[17-j];
  158.     }
  159.     /* window x for next time x0 */
  160.     for(j=0;j<4;j++) {
  161.         xa = x[j];
  162.         xb = x[8-j];
  163.         x[j]    = win[btype][18+j]*xb;
  164.         x[8-j]  = win[btype][(18+8)-j]*xa;
  165.         x[9+j]  = win[btype][(18+9)+j]*xa;
  166.         x[17-j] = win[btype][(18+17)-j]*xb;
  167.     }
  168.     xa = x[j];
  169.     x[j]    = win[btype][18+j]*xa;
  170.     x[9+j]  = win[btype][(18+9)+j]*xa;
  171.     
  172.     x+=18;
  173.     x0+=18;
  174. }
  176. /*-- do short blocks (if any) --*/
  177. n = (ntot+17)/18;       /* number of 6 pt dct's triples to do */
  178. for(   ;i<n;i++) {
  179.     imdct6_3(x);
  180.     for(j=0;j<3;j++) {
  181.         y[j][i]    = x0[j];
  182.         y[3+j][i]  = x0[3+j];
  184.         y[6+j][i]  = x0[6+j]  +  win[2][j]  *x[3+j];
  185.         y[9+j][i]  = x0[9+j]  +  win[2][3+j]*x[5-j];
  187.         y[12+j][i] = x0[12+j] + win[2][6+j]*x[2-j] + win[2][j]*x[(6+3)+j];
  188.         y[15+j][i] = x0[15+j] + win[2][9+j]*x[j]   + win[2][3+j]*x[(6+5)-j];
  189.     }
  190.     /* window x for next time x0 */
  191.     for(j=0;j<3;j++) {
  192.         x[j]    = win[2][6+j]*x[(6+2)-j] + win[2][j]  *x[(12+3)+j];
  193.         x[3+j]  = win[2][9+j]*x[6+j]     + win[2][3+j]*x[(12+5)-j];
  194.     }
  195.     for(j=0;j<3;j++) {
  196.         x[6+j]  = win[2][6+j]*x[(12+2)-j];
  197.         x[9+j]  = win[2][9+j]*x[12+j];
  198.     }
  199.     for(j=0;j<3;j++) {
  200.         x[12+j] = 0.0f;  
  201.         x[15+j] = 0.0f;  
  202.     }
  203.     x+=18;
  204.     x0+=18;
  205. }
  207. /*--- overlap prev if prev longer that current --*/
  208. n = (nprev+17)/18;
  209. for(;i<n;i++) {
  210.     for(j=0;j<18;j++) y[j][i]  = x0[j];
  211.     x0+=18;
  212. }
  213. nout = 18*i;
  215. /*--- clear remaining only to band limit --*/
  216. for(;i<band_limit_nsb;i++) {
  217.     for(j=0;j<18;j++)  y[j][i]  = 0.0f;
  218. }
  220. return nout;
  221. }
  222. /*====================================================================*/
  223. /*--------------------------------------------------------------------*/
  224. /*--------------------------------------------------------------------*/
  225. /*-- convert to mono, add curr result to y,
  226.     window and add next time to current left */
  227. int hybrid_sum(float xin[], float xin_left[], float y[18][32],
  228.             int btype, int nlong, int ntot)
  229. {
  230. int i, j;
  231. float *x, *x0;
  232. float xa, xb;
  233. int n;
  234. int nout;
  238. if( btype == 2 ) btype = 0;
  239. x = xin;
  240. x0 = xin_left;
  242. /*-- do long blocks (if any) --*/
  243. n = (nlong+17)/18;       /* number of dct's to do */
  244. for(i=0;i<n;i++) {
  245.     imdct18(x);
  246.     for(j=0;j<9;j++) {
  247.         y[j][i]    += win[btype][j]*x[9+j];
  248.         y[9+j][i]  += win[btype][9+j]*x[17-j];
  249.     }
  250.     /* window x for next time x0 */
  251.     for(j=0;j<4;j++) {
  252.         xa = x[j];
  253.         xb = x[8-j];
  254.         x0[j]    += win[btype][18+j]*xb;
  255.         x0[8-j]  += win[btype][(18+8)-j]*xa;
  256.         x0[9+j]  += win[btype][(18+9)+j]*xa;
  257.         x0[17-j] += win[btype][(18+17)-j]*xb;
  258.     }
  259.     xa = x[j];
  260.     x0[j]    += win[btype][18+j]*xa;
  261.     x0[9+j]  += win[btype][(18+9)+j]*xa;
  262.     
  263.     x+=18;
  264.     x0+=18;
  265. }
  267. /*-- do short blocks (if any) --*/
  268. n = (ntot+17)/18;       /* number of 6 pt dct's triples to do */
  269. for(   ;i<n;i++) {
  270.     imdct6_3(x);
  271.     for(j=0;j<3;j++) {
  272.         y[6+j][i]  +=  win[2][j]  *x[3+j];
  273.         y[9+j][i]  +=  win[2][3+j]*x[5-j];
  275.         y[12+j][i] +=  win[2][6+j]*x[2-j] + win[2][j]*x[(6+3)+j];
  276.         y[15+j][i] +=  win[2][9+j]*x[j]   + win[2][3+j]*x[(6+5)-j];
  277.     }
  278.     /* window x for next time */
  279.     for(j=0;j<3;j++) {
  280.         x0[j]    += win[2][6+j]*x[(6+2)-j] + win[2][j]  *x[(12+3)+j];
  281.         x0[3+j]  += win[2][9+j]*x[6+j]     + win[2][3+j]*x[(12+5)-j];
  282.     }
  283.     for(j=0;j<3;j++) {
  284.         x0[6+j]  += win[2][6+j]*x[(12+2)-j];
  285.         x0[9+j]  += win[2][9+j]*x[12+j];
  286.     }
  287.     x+=18;
  288.     x0+=18;
  289. }
  291. nout = 18*i;
  293. return nout;
  294. }
  295. /*--------------------------------------------------------------------*/
  296. void sum_f_bands( float a[], float b[], int n)
  297. {
  298. int i;
  299.     for(i=0;i<n;i++) a[i] += b[i];
  300. }
  301. /*--------------------------------------------------------------------*/
  302. /*--------------------------------------------------------------------*/
  303. void FreqInvert(float y[18][32], int n)
  304. {
  305. int i, j;
  307. n = (n+17)/18;
  308. for(j=0;j<18;j+=2) {
  309. for(i=0;i<n; i+=2)  {
  310.         y[1+j][1+i] = -y[1+j][1+i];
  311. }
  312. }
  316. }
  317. /*--------------------------------------------------------------------*/