hwin.c
上传用户:dangjiwu
上传日期:2013-07-19
资源大小:42019k
文件大小:11k
源码类别:

Symbian

开发平台:

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