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

源码开发语言/平台
当前位置：首页 > 源码/资料 > 多媒体 > mpeg/mp3 > 查看源码
MP3_Fastsb.c
资源名称：mp3ruanjianjiema.rar　[点击查看]
上传用户：kingbiz
上传日期：2022-06-24
资源大小：2524k
文件大小：20k
源码类别：

mpeg/mp3

开发平台：

C/C++

MP3_Fastsb.c：源码内容
  2. #include "mp3_common.h"
  3. #include <math.h> 
  5. /*============================================================================*/
  8. /*============================================================================*/
  10. static int Granule_sbsynth_Vptr[2]={64,64};
  11. typedef float BB[2][1024];
  12. static BB Granule_sbsynth_V;
  14. static float A16[16][16], A8[8][8];       /* DCT matrix         */
  15. static float G16[16][16], G8[8][8];       /* Output butterfly   */
  16. static float H16[16][16], H8[8][8];       /* Scaling            */
  18. static float B16[16][16], B8[8][8];       /* B = G * DCT * H    */
  20. static float D_coex[512]={
  21.      0.000000000    ,-0.000015259   ,-0.000015259   ,-0.000015259,
  22.  
  23.     -0.000015259     ,-0.000015259,     -0.000015259,     -0.000030518 ,
  24.  
  25.     -0.000030518     ,-0.000030518    ,-0.000030518    ,-0.000045776 ,
  26.  
  27.     -0.000045776    ,-0.000061035    ,-0.000061035    ,-0.000076294 ,
  28.  
  29.     -0.000076294    ,-0.000091553    ,-0.000106812    ,-0.000106812 ,
  30.  
  31.     -0.000122070    ,-0.000137329    ,-0.000152588    ,-0.000167847 ,
  32.  
  33.     -0.000198364    ,-0.000213623    ,-0.000244141    ,-0.000259399 ,
  34.  
  35.    -0.000289917    ,-0.000320435    ,-0.000366211    ,-0.000396729 ,
  36.  
  37.    -0.000442505    ,-0.000473022    ,-0.000534058   ,-0.000579834 ,
  38.  
  39.    -0.000625610    ,-0.000686646    ,-0.000747681    ,-0.000808716 ,
  40.  
  41.    -0.000885010    ,-0.000961304    ,-0.001037598    ,-0.001113892 ,
  42.  
  43.    -0.001205444    ,-0.001296997    ,-0.001388550    ,-0.001480103 ,
  44.  
  45.    -0.001586914   ,-0.001693726    ,-0.001785278    ,-0.001907349 ,
  46.  
  47.    -0.002014160    ,-0.002120972    ,-0.002243042    ,-0.002349854 ,
  48.  
  49.    -0.002456665    ,-0.002578735    ,-0.002685547    ,-0.002792358 ,
  50.  
  51.    -0.002899170    ,-0.002990723    ,-0.003082275    ,-0.003173828 ,
  52.  
  53.     -0.003250122   , -0.003326416    , -0.003387451    , -0.003433228 ,
  54.  
  55.    0.003463745    , 0.003479004    , 0.003479004    , 0.003463745 ,
  56.  
  57.    0.003417969    , 0.003372192    , 0.003280640    , 0.003173828 ,
  58.  
  59.    0.003051758    , 0.002883911   , 0.002700806   , 0.002487183 ,
  60.  
  61.     0.002227783   , 0.001937866    , 0.001617432   , 0.001266479 ,
  62.  
  63.    0.000869751   , 0.000442505    ,-0.000030518    ,-0.000549316 ,
  64.  
  65.    -0.001098633    ,-0.001693726    ,-0.002334595    ,-0.003005981 ,
  66.  
  67.    -0.003723145    ,-0.004486084    ,-0.005294800    ,-0.006118774 ,
  68.  
  69.   -0.007003784   ,-0.007919312    ,-0.008865356    ,-0.009841919 ,
  70.  
  71.   -0.010848999   ,-0.011886597   ,-0.012939453   ,-0.014022827 ,
  72.  
  73.   -0.015121460   ,-0.016235352   ,-0.017349243   ,-0.018463135 ,
  74.  
  75.   -0.019577026   ,-0.020690918   ,-0.021789551   ,-0.022857666 ,
  76.  
  77.  -0.023910522   ,-0.024932861    ,-0.025909424 ,-0.026840210 ,
  78.  
  79.   -0.027725220   ,-0.028533936   ,-0.029281616   ,-0.029937744 ,
  80.  
  81.   -0.030532837   ,-0.031005859   ,-0.031789329   ,-0.031661987 ,
  82.  
  83.   -0.031814575   ,-0.031845093  ,-0.031738281   ,-0.031478882 ,
  84.  
  85.    0.031082153   , 0.030517578   , 0.029785156   , 0.028884888 ,
  86.  
  87.   0.027801514   , 0.026535034  , 0.025085449   , 0.023422241 ,
  88.  
  89.    0.021575928   , 0.020531250   , 0.017257690   , 0.014801025 ,
  90.  
  91.   0.012115479  , 0.009231567   , 0.006134033 , 0.002822876 ,
  92.  
  93.  -0.000686646   ,-0.004394531   ,-0.008316040   ,-0.012420654 ,
  94.  
  95.   -0.016708374   ,-0.021179199   ,-0.025817871   ,-0.030609131 ,
  96.  
  97.  -0.035552979   ,-0.040634155   ,-0.045837402   ,-0.051132202 ,
  98.  
  99.   -0.056533813   ,-0.061996460   ,-0.067520142  ,-0.073059082 ,
  100.  
  101.   -0.078628540   ,-0.084182739   ,-0.089706421   ,-0.095169067 ,
  102.  
  103.   -0.100540161   ,-0.105819702   ,-0.110946655   ,-0.115921021 ,
  104.  
  105.   -0.120697021   ,-0.125259399   ,-0.129562378   ,-0.133590698 ,
  106.  
  107.   -0.137298584   ,-0.140670776  ,-0.143676758  ,-0.146255493 ,
  108.  
  109.  -0.148422241  ,-0.150115967   ,-0.151306152   ,-0.151962280 ,
  110.  
  111.  -0.152069092   ,-0.151596069   ,-0.150497437   ,-0.148773193 ,
  112.  
  113.  -0.146362305   ,-0.143264771  ,-0.109450073   ,-0.134887695 ,
  114.  
  115.  -0.129577637   ,-0.123474121   ,-0.116577148   ,-0.108856201 ,
  116.  
  117.    0.100311279   , 0.090927124   , 0.080088477   , 0.069595337 ,
  118.  
  119.   0.057617187   , 0.044784546   , 0.031082153   , 0.016510010 ,
  120.  
  121.   0.001068115  ,0.015228271   ,0.032679150   ,0.050354004 ,
  122.  
  123.  0.069168091   ,0.088775635  ,0.109161377  ,0.130310059 ,
  124.  
  125.  -0.152206421   ,-0.174789429   ,-0.198059082   ,-0.221984863 ,
  126.  
  127.  -0.246505737  ,-0.271591187  ,-0.297210693   ,-0.323318481 ,
  128.  
  129.   -0.349868774   ,-0.376800537   ,-0.404083252   ,-0.431655884 ,
  130.  
  131.   -0.459472656  ,-0.487472534  ,-0.515609741   ,-0.543823242 ,
  132.  
  133.  -0.572036743   ,-0.600219727   ,-0.628295898   ,-0.656219482 ,
  134.  
  135.  -0.683914185  ,-0.711318970  ,-0.738372803  ,-0.765029907 ,
  136.  -0.791213989   ,-0.816864014   ,-0.841949463   ,-0.866363525 ,
  137.  
  138.  -0.890090942  ,-0.913055420  ,-0.935195923   ,-0.956481934 ,
  139.  
  140.  -0.976852417  ,-0.996246338   ,-1.014617920   ,-1.031936646 ,
  141.  -1.048156738  ,-1.063217163   ,-1.077117920  ,-1.089782715 ,
  142.  
  143.  -1.101211548   ,-1.111373901  ,-1.120223999  ,-1.127746582 ,
  144.  
  145.   1.133926392   ,1.138763428  ,1.142211914   ,1.144287109 ,
  146.  
  147.   1.144989014  , 1.144287109   , 1.142211914   , 1.138763428 ,
  148.  
  149.    1.133926392   , 1.127746582  , 1.120223999  , 1.111373901 ,
  150.  
  151.   1.101211548   , 1.089782715   , 1.077117920  , 1.063217163 ,
  152.  
  153.   1.048156738   , 1.031936646   , 1.014617920  , 0.996246338 ,
  154.  
  155.    0.976852417   , 0.956481934   , 0.935195923   , 0.913055420 ,
  156.  
  157.    0.890090942   , 0.866363525  , 0.841949463   , 0.816864014 ,
  158.  
  159.    0.791213989   , 0.765029907   , 0.738372803   , 0.711318970 ,
  160.  
  161.    0.683914185   , 0.656219482   , 0.628295898   , 0.600219727 ,
  162.  
  163.    0.572036743   , 0.543823242   , 0.515609741   , 0.487472534 ,
  164.  
  165.    0.459472656   , 0.431655884   , 0.404083252   , 0.376800537 ,
  166.  
  167.    0.349868774   , 0.323318481   , 0.297210693  , 0.271591187 ,
  168.  
  169.    0.246505737   , 0.221984863  , 0.498059082   , 0.174789429 ,
  170.  
  171.    0.152206421   , 0.130310059   , 0.109161377   , 0.088775635 ,
  172.  
  173.    0.069168091   , 0.050354004   , 0.032379150   , 0.015228271 ,
  174.  
  175.   -0.001068115   ,-0.016510010   ,-0.031082153   ,-0.044784546 ,
  176.  
  177.  -0.057617187   ,-0.069595337  ,-0.080688477  ,-0.090927124 ,
  178.  
  179.   0.100311279  , 0.108856201   , 0.116577148   , 0.123474121 ,
  180.  
  181.   0.129577637   , 0.134887695   , 0.139450073   , 0.143264771 ,
  182.  
  183.   0.146362305   , 0.148773193   , 0.150497437   , 0.151596069 ,
  184.  
  185.    0.152069092   , 0.151962280   , 0.151306152  , 0.150115967 ,
  186.  
  187.    0.148422241  , 0.146255493   , 0.143676758   , 0.140670776 ,
  188.  
  189.    0.137298584   , 0.133590698  , 0.129562378   , 0.125259399 ,
  190.  
  191.   0.120697021   , 0.115921021   , 0.110946655   , 0.105819702 ,
  192.  
  193.    0.100540161   , 0.095169067   , 0.089706421   , 0.084182739 ,
  194.  
  195.    0.078628540   , 0.073059082  , 0.067520142   , 0.061996460 ,
  196.  
  197.    0.056533813   , 0.051132202   , 0.045837402   , 0.040634155 ,
  198.  
  199.   0.035552979   , 0.030609131   , 0.025817871  , 0.021179199 ,
  200.  
  201.   0.016708374   , 0.012420654   , 0.008316040   , 0.004394531 ,
  202.  
  203.    0.000686646  ,-0.002822876  ,-0.001134033   ,-0.009231567 ,
  204.  
  205.  -0.012115479   ,-0.014801025   ,-0.017257690   ,-0.019531250 ,
  206.  
  207.  -0.021575928   ,-0.023422241  ,-0.025085449   ,-0.026535034 ,
  208.  
  209.   -0.027801514   ,-0.028884888   ,-0.029785156  ,-0.030517578 ,
  210.  
  211.    -0.031082153   , -0.031478882   , -0.031738281   , -0.031845093 ,
  212.  
  213.   0.031814575   , 0.031661987   , 0.031387329   , 0.031005859 ,
  214.  
  215.    0.030532837  , 0.029937744   , 0.029281616   , 0.028533936 ,
  216.  
  217.    0.027725220   , 0.026840210   , 0.025909424   , 0.024932861 ,
  218.  
  219.    0.023910522   , 0.022857666   , 0.021789551   , 0.020690918 ,
  220.  
  221.   0.019577026  , 0.018463135   , 0.017349243   , 0.016235352 ,
  222.   0.015121460   , 0.014022827   , 0.012939453   , 0.011886597 ,
  223.  
  224.   0.010848999   , 0.009841919  , 0.008865356   , 0.007919312 ,
  225.  
  226.    0.007003784   , 0.006118774   , 0.005294800   , 0.004486084 ,
  227.  
  228.   0.003723145   , 0.003005981   , 0.002334595  , 0.001693726 ,
  229.  
  230.   0.001098633   , 0.000549316  , 0.000030518   ,-0.000442505 ,
  231.  -0.000869751   ,-0.001266479   ,-0.001617432    ,-0.001937866 ,
  232.  
  233.  -0.002227783   ,-0.002487183   ,-0.002700806   ,-0.002883911 ,
  234.  
  235.  -0.003051758   ,-0.003173828   ,-0.003280640   ,-0.003372192 ,
  236.  
  237.  -0.003417969   ,-0.003463745   ,-0.003479004   ,-0.003479004 ,
  238.  
  239.  -0.003463745   ,-0.003433228   ,-0.001387451  ,-0.003326416 ,
  240.  
  241.   -0.003250122   , -0.003173828   , -0.003082275   , -0.002990723 ,
  242.  
  243.   -0.002899170   , -0.002792358   , -0.002685547   , -0.002578735 ,
  244.  
  245.    0.002456665   , 0.002349854   , 0.002243042   , 0.002120972 ,
  246.  
  247.   0.002014160   , 0.001907349   , 0.001785278  , 0.001693726 ,
  248.  
  249.   0.001586914   , 0.001480103   , 0.001388550   , 0.001296997 ,
  250.  
  251.  0.001205444   , 0.001113892   , 0.001037598   , 0.000961304 ,
  252.  
  253.   0.000885010   , 0.000808716   , 0.000747681   , 0.000686646 ,
  254.  
  255.    0.000625610   , 0.000579834   , 0.000534058   , 0.000473022 ,
  256.  
  257.   0.000442505   , 0.000396729   , 0.000366211   , 0.000320435 ,
  258.  
  259.   0.000289917   , 0.000269399   , 0.000244141   , 0.000213623 ,
  260.  
  261.   0.000198364   , 0.000167847   , 0.000152588   , 0.000137329 ,
  262.  
  263.   0.000122070   , 0.000106812   , 0.000106812   , 0.000091553 ,
  264.  0.000076294   , 0.000076294   , 0.000061035   , 0.000061035 ,
  265.  
  266.   0.000045776  , 0.000045776   , 0.000030518   , 0.000030518 ,
  267.  
  268.   0.000030518   , 0.000030518   , 0.000015259   , 0.000015259 ,
  269.  
  270.    0.000015259   , 0.000015259   , 0.000075259   , 0.000015259,
  271. };
  273. /*============================================================================*/
  274. /*============================================================================*/
  275. /*============================================================================*/
  276. static void matrix_mul16(float in1[16][16],
  277.   float in2[16][16],
  278.   float out[16][16])
  279. {
  280.     int i,j,z;
  282.     for(i = 0; i < 16; i++)
  283.     {
  284.   for(j = 0; j < 16; j++)
  285.   {
  286.     out[i][j] = 0.0;
  287.     for(z = 0; z < 16; z++)
  288.     {
  289.   out[i][j] += in1[i][z] * in2[z][j];
  290. }
  292.   }
  293.   
  294.     } 
  295. }
  296. /*============================================================================*/
  298. static void matrix_mul8(float in1[8][8],
  299.  float in2[8][8],
  300.  float out[8][8])
  301. {
  302.     int i,j,z;
  304.     for(i = 0; i < 8; i++)
  305.     {
  306.   for(j = 0; j < 8; j++)  
  307.   {
  308.     out[i][j] = 0.0;
  309.     for(z = 0; z < 8; z++)
  310.     {
  311.     out[i][j] += in1[i][z] * in2[z][j];
  312.   }
  313.  
  314.   }
  315.   
  316.     }
  317. }
  319. /*============================================================================*/
  320. /*============================================================================*/
  322. static void fast_idct_init(void)
  323. {
  324.     int i,j;
  325.     float t16[16][16], t8[8][8];
  328.     /* create the 16 matrixes */
  330.     for(i = 0; i < 16; i++) {
  331.   for(j = 0; j < 16; j++) 
  332.   {
  333.     A16[i][j] = cos((2*j+1)*i*PI/32);
  334.     if(i == j || j == (i + 1))
  335.   G16[i][j] = 1.0f;
  336.     else
  337.   G16[i][j] = 0.0f;
  338.     if(i == j)
  339.   H16[i][j] = 1.0f/(2*cos((2*i+1)*PI/64));
  340.     else
  341.   H16[i][j] = 0.0;
  342.   }
  343.     }
  345.     /* create the 8 matrixes */
  347.     for(i = 0; i < 8; i++) {
  348.   for(j = 0; j < 8; j++) 
  349.   {
  350.     A8[i][j] = cos((2*j+1)*i*PI/16);
  351.     if(i == j || j == (i + 1))
  352.   G8[i][j] = 1.0f;
  353.     else
  354.   G8[i][j] = 0.0f;
  355.     if(i == j)
  356.   H8[i][j] = 1.0f/(2*cos((2*i+1)*PI/32));
  357.     else
  358.   H8[i][j] = 0.0f;
  359.   }
  360.     }
  362.     /* generate the B matrixes */
  364.     matrix_mul16(A16, H16, t16);
  365.     matrix_mul16(G16, t16, B16);
  367.     matrix_mul8(A8, H8, t8);
  368.     matrix_mul8(G8, t8, B8);
  369. }
  370. /*============================================================================*/
  372. /* This is a two-level implementation of Lee's fast-DCT algorithm */
  373. /* 
  374.    The 32 input values are split in two 16-value vectors using an
  375.    even butterfly and an odd butterfly. The odd values are taken
  376.    through Lee's odd path using a 16x16 DCT matrix (A16) and appropriate
  377.    scaling (G16*A16*H16). The even values are further split into
  378.    two 8-value vectors using even and odd butterflies into ee and eo.
  379.    The ee values are fed through an 8x8 DCT matrix (A8) while the eo
  380.    values are fed through the odd path using G8*A8*H8.
  382.    This two-level configuration uses 384 muls and 432 adds, compared
  383.    to the direct 32x32 DCT which uses 1024 muls and 992 adds.
  384. */
  387. static void fast_idct(float *in, float *out)
  388. {
  389.     float even[16], odd[16], ee[8], eo[8];
  390.     float s1, s2;
  391.     float t[32];
  392.     int i, j;
  393. static int init=1;
  395.     if(init)
  396. {
  397. fast_idct_init();
  398. init=0;
  399. }
  400. /* input butterflies - level 1 */
  401.     /* 32 adds */
  403.     for(i = 0; i < 16; i++) {
  404.    even[i] = in[i] + in[31-i];
  405.    odd[i] = in[i] - in[31-i];
  406.     }
  408.     /* input butterflies - level 2 */
  409.     /* 16 adds */
  411.     for(i = 0; i < 8; i++) {
  412.    ee[i] = even[i] + even[15-i];
  413.    eo[i] = even[i] - even[15-i];
  414.     }
  416.     /* multiply the even_even vector (ee) with the ee matrix (A8) */
  417.     /* multiply the even_odd vector (eo) with the eo matrix (B8) */
  418.     /* 128 muls, 128 adds */
  420.     for(i = 0; i < 8; i++) {
  421.   s1 = 0.0;
  422.   s2 = 0.0;
  423.   for(j = 0; j < 8; j += 2) 
  424.   {
  425.     s1 += A8[i][j] * ee[j] +
  426. A8[i][j+1] * ee[j+1];
  427.     s2 += B8[i][j] * eo[j] +
  428. B8[i][j+1] * eo[j+1];
  429.   }
  430.   t[i*4] = s1;
  431.   t[i*4+2] = s2;
  432.     }
  435.     /* multiply the odd vector (odd) with the odd matrix (B16) */
  436.     /* 256 muls, 256 adds */
  438.     for(i = 0; i < 16; i++) 
  439. {
  440.   s1 = 0.0;
  441.   for(j = 0; j < 16; j += 4) 
  442.   {
  443.     s1 += B16[i][j] * odd[j] +
  444. B16[i][j+1] * odd[j+1] +
  445. B16[i][j+2] * odd[j+2] +
  446. B16[i][j+3] * odd[j+3];
  447.   }
  448.   t[i*2+1] = s1;
  449.     }
  451.     /* the output vector t now is expanded to 64 values using the
  452.        symmetric property of the cosinus function */
  454.     for(i = 0; i < 16; i++) 
  455. {
  456.   out[i] = t[i+16];
  457.   out[i+17] = -t[31-i];
  458.   out[i+32] = -t[16-i];
  459.   out[i+48] = -t[i];
  460.     }
  461.     out[16] = 0.0;
  462. }
  463. /*============================================================================*/
  466. void Granule_subband_synthesis(int ch, SS s, short SAM[2][SSLIMIT][SBLIMIT])
  467. {
  468.     int i, j, t, k;
  469. long is;
  470.     float band[32];
  471.     float *v, sum, *d;
  472. short *S=&SAM[0][0][0];
  474.     /* We have 18 time-vectors of 32 subband magnitudes each. For every
  475.        vector of 32 magnitudes, the subband synthesis generates 32
  476.        PCM samples, so the result of 18 of these is 18*32=576 samples.
  477.      */
  479.     /* go through each time window */
  481.     for(t = 0; t < 18; t++)
  482. {
  483. /* extract the subband strengths */
  484.        v = &s[0][t];
  485. for(i = 0; i < 32; i++)
  486. {
  487.       band[i] = *v;
  488.       v = &v[18];
  489. }
  490. /* advance the buffer position */
  491.    Granule_sbsynth_Vptr[ch] = (Granule_sbsynth_Vptr[ch] - 64) & 0x3ff;
  492.    v = &Granule_sbsynth_V[ch][Granule_sbsynth_Vptr[ch]];
  494.    fast_idct(band, v);
  496. /* 32*16=512 mac's */
  497.     /*      15          */
  498.     /* Sj =SIGM W(j+32i) */
  499.     /*     i=0          */
  501.   v = &Granule_sbsynth_V[ch][0];
  502.       for (j = 0; j < 32; j++) 
  503.   {
  504.         d = &D_coex[j];
  505.     k = j + Granule_sbsynth_Vptr[ch];
  506.     sum = 0.0f;
  507.     for(i = 0; i < 8; i++) 
  508. {
  509.       sum += d[0] * v[k];
  510.       k = (k + 96) & 0x3ff;
  511.       sum += d[32] * v[k];
  512.       k = (k + 32) & 0x3ff; 
  513.       d = &d[64];
  514. }
  516.   /* convert to integer, and clip the output */
  517.       is = (long)(sum*32768);
  519.     if(is >= 32768) 
  520.       is = 32767;
  521.     else if(is < -32768)
  522.       is = -32768;
  523.     *S = (short)is;
  524.     S++;
  525.   }
  526.     }
  527. }
  528. /*============================================================================*/
  529. //子带合成
  531. /* 18 * (4096+1024) = 92160 MAC's per call, with 2 calls per frame and
  532.    38 frames per second this is 7 million MAC's per second.
  534.    18 * (384 * 2 + 1024) = 32256 MAC's per call using Lee's fast DCT! That
  535.    is just 2.4 million MAC's per second!
  537. We need a buffer of 1024 floats per channel.
  539.    */
  541. void Granule_subband_synthesis2x(SS s1,SS s2,short  SAM[2][SSLIMIT][SBLIMIT])
  542. {
  544.     int i, j, t, k;
  545.     float band[64];
  546.     float *v1, *v2, sum, /*sum2,*/ *d;
  547. float *v;
  548. S32 is;
  549. short *S=&SAM[0][0][0];
  551.     /* We have 18 time-vectors of 32 subband magnitudes each. For every
  552.        vector of 32 magnitudes, the subband synthesis generates 32
  553.        PCM samples, so the result of 18 of these is 18*32=576 samples.
  554.      */
  556.     /* go through each time window */
  557.  
  558.     for(t = 0; t < 18; t++) 
  559. {
  561. /* extract the subband strengths */
  563.    v1 = &s1[0][t];
  564.    v2 = &s2[0][t];
  565.    
  566.    
  567.    for(i = 0; i < 32; i++)
  568.    {
  569.      band[i]    = *v1;
  570.      band[i+32] = *v2;
  571.      v1 = &v1[18];
  572.      v2 = &v2[18];
  573.    }
  575. /* advance the buffer position */
  577.     Granule_sbsynth_Vptr[0] = (Granule_sbsynth_Vptr[0] - 64) & 0x3ff;
  578.   v = &(Granule_sbsynth_V[0][Granule_sbsynth_Vptr[0]]);
  580. /* calculate 64 values for each channel and insert them into the 1024 wide buffer */
  582.   fast_idct(band, v);
  583.   
  584. //   fast_idct(&band[32], &v[1024]);
  585.          
  586. /* 32*16*2=1024 mac's */
  588. /* windowing - calculate 32 samples. each sample is the sum of 16 terms */
  590. /*     15          */
  591. /* Sj = E W(j+32i) */
  592. /*    i=0          */
  594.   v = &Granule_sbsynth_V[0][0];
  595.   for (j = 0; j < 32; j++)
  596.   { 
  597.     d = &D_coex[j];
  598.     k = j + Granule_sbsynth_Vptr[0];
  599.     
  600.     sum  = (float)0.0f;
  601.   //  sum2 = (float)0.0f;
  602.    
  603.     for(i = 0; i < 8; i++) 
  604. {
  605.    sum  += d[0] * v[k];
  606. //    sum2 += d[0] * v[k+1024];
  607.    k = (k + 96) & 0x3ff;
  609.    sum  += d[32] * v[k];
  610. //    sum2 += d[32] * v[k+1024];
  611.    k = (k + 32) & 0x3ff; 
  612.    d = &d[64];
  613.     }
  614.     
  615.    
  616.    
  617.     /* convert to integer, and clip the output */
  618. is=(S32) (sum*32768);
  619. // OVERFLOW_CHECKING
  620.     if(is >= 32768) 
  621.   *S = 32767;
  622. else if(is < -32768) 
  623.   *S = -32768;
  624. else 
  625.      *S = (short)is;
  627. S++;
  628.     
  629.     
  630. is=(S32) (sum*32768);
  631.     if(is >= 32768) 
  632.   *(S) = 32767;
  633. else if(is < -32768) 
  634.   *(S) = -32768;
  635. else 
  636.   *(S) = (short)is;
  638.     S++;
  640.   
  641.   }
  642.   
  643.     }
  644.     
  645. }
  648. void Granule_subband_synthesis2(SS s1,SS s2,short  SAM[2][SSLIMIT][SBLIMIT])
  649. {
  651.     int i, j, t, k;
  652.     float band[64];
  653.     float *v1, *v2, sum, /*sum2,*/ *d;
  654. float *v;
  655. S32 is;
  656. short *S=&SAM[0][0][0];
  658.     /* We have 18 time-vectors of 32 subband magnitudes each. For every
  659.        vector of 32 magnitudes, the subband synthesis generates 32
  660.        PCM samples, so the result of 18 of these is 18*32=576 samples.
  661.      */
  663.     /* go through each time window */
  664.  
  665.     for(t = 0; t < 18; t++) 
  666. {
  668. /* extract the subband strengths */
  670.    v1 = &s1[0][t];
  671.    v2 = &s2[0][t];
  672.    
  673.    
  674.    for(i = 0; i < 32; i++)
  675.    {
  676.      band[i]    = *v1;
  677.      band[i+32] = *v2;
  678.      v1 = &v1[18];
  679.      v2 = &v2[18];
  680.    }
  684. /* advance the buffer position */
  686.     Granule_sbsynth_Vptr[0] = (Granule_sbsynth_Vptr[0] - 64) & 0x3ff;
  687.   v = &(Granule_sbsynth_V[0][Granule_sbsynth_Vptr[0]]);
  689. /* calculate 64 values for each channel and insert them into the 1024 wide buffer */
  691.   fast_idct(band, v);
  692.   
  693. //   fast_idct(&band[32], &v[1024]);
  694.          
  695. /* 32*16*2=1024 mac's */
  697. /* windowing - calculate 32 samples. each sample is the sum of 16 terms */
  699. /*     15          */
  700. /* Sj = E W(j+32i) */
  701. /*    i=0          */
  703.   v = &Granule_sbsynth_V[0][0];
  704.   for (j = 0; j < 32; j++)
  705.   { 
  706.     d = &D_coex[j];
  707.     k = j + Granule_sbsynth_Vptr[0];
  708.     
  709.     sum  = (float)0.0f;
  710.   //  sum2 = (float)0.0f;
  711.    
  712.     for(i = 0; i < 8; i++) 
  713. {
  714.    sum  += d[0] * v[k];
  715. //    sum2 += d[0] * v[k+1024];
  716.    k = (k + 96) & 0x3ff;
  718.    sum  += d[32] * v[k];
  719. //    sum2 += d[32] * v[k+1024];
  720.    k = (k + 32) & 0x3ff; 
  721.    d = &d[64];
  722.     }
  723.     
  724.    
  725.    
  726.     /* convert to integer, and clip the output */
  727. is=(S32) (sum*32768);
  728. // OVERFLOW_CHECKING
  729.     if(is >= 32768) 
  730.   *S = 32767;
  731. else if(is < -32768) 
  732.   *S = -32768;
  733. else 
  734.      *S = (short)is;
  736. S++;
  737.     
  738.     
  739.     
  740. is=(S32) (sum*32768);
  741.     if(is >= 32768) 
  742.   *(S) = 32767;
  743. else if(is < -32768) 
  744.   *(S) = -32768;
  745. else 
  746.   *(S) = (short)is;
  748.     S++;
  750.   
  751.   }
  752.   
  753.     }
  754.     
  755. }
  758. /*============================================================================*/
  759. /*============================================================================*/
  760. /*============================================================================*/
  761. /*============================================================================*/