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

源码开发语言/平台
当前位置：首页 > 源码/资料 > Windows编程 > 多媒体编程 > 查看源码
synfilt.cpp
资源名称：Visual C++视频音频开发实用工程案例精选.rar　[点击查看]
上传用户：tuheem
上传日期：2007-05-01
资源大小：21889k
文件大小：33k
源码类别：

多媒体编程

开发平台：

Visual C++

synfilt.cpp：源码内容
  1. /* synfilt.cpp
  3. Synthesis Filter implementation */
  5. /* -- 03/20/97 --
  6.  *  compute_new_v() -- reoptimized with the assumption that constant offsets
  7.  *    to memory are free.  Common subexpression were redone for better
  8.  *    optimization.
  9.  *  compute_pcm_samples() -- reoptimized with constant offsets.
  10.  *
  11.  *  -- Conrad Wei-Li Song (conradsong@mail.utexas.edu)
  12.  */
  14. /*
  15.  *  @(#) synthesis_filter.cc 1.14, last edit: 6/21/94 11:22:20
  16.  *  @(#) Copyright (C) 1993, 1994 Tobias Bading (bading@cs.tu-berlin.de)
  17.  *  @(#) Berlin University of Technology
  18.  *
  19.  *  This program is free software; you can redistribute it and/or modify
  20.  *  it under the terms of the GNU General Public License as published by
  21.  *  the Free Software Foundation; either version 2 of the License, or
  22.  *  (at your option) any later version.
  23.  *
  24.  *  This program is distributed in the hope that it will be useful,
  25.  *  but WITHOUT ANY WARRANTY; without even the implied warranty of
  26.  *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  27.  *  GNU General Public License for more details.
  28.  *
  29.  *  You should have received a copy of the GNU General Public License
  30.  *  along with this program; if not, write to the Free Software
  31.  *  Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
  32.  */
  34. /*
  35.  *  Changes from version 1.1 to 1.2:
  36.  *    - compute_new_v() uses a 32 point fast cosine transform as described by
  37.  *      Byeong Gi Lee in IEEE Transactions ASSP-32 Part 2, August 1984,
  38.  *      "A New Algorithm to Compute the Discrete Cosine Transform"
  39.  *      instead of the matrix-vector multiplication in V1.1
  40.  *    - loop unrolling done in compute_pcm_samples()
  41.  *    - if ULAW is defined, the synthesis filter does a downsampling
  42.  *      to 8 kHz by dropping samples and ignoring subbands above 4 kHz
  43.  */
  45. #include <math.h>
  46. #include "all.h"
  47. #include "synfilt.h"
  49. #define MY_PI 3.14159265358979323846
  51. #include "scalfact.h"
  53. #ifdef __WIN32__
  54. #pragma warning (disable: 4244 4305)
  55. #endif
  57. static const real cos1_64  = (real)(1.0 / (2.0 * cos(MY_PI / 64.0)));
  58. static const real cos3_64  = (real)(.0 / (2.0 * cos(MY_PI * 3.0  / 64.0)));
  59. static const real cos5_64  = (real)(1.0 / (2.0 * cos(MY_PI * 5.0  / 64.0)));
  60. static const real cos7_64  = (real)(1.0 / (2.0 * cos(MY_PI * 7.0  / 64.0)));
  61. static const real cos9_64  = (real)(1.0 / (2.0 * cos(MY_PI * 9.0  / 64.0)));
  62. static const real cos11_64 = (real)(1.0 / (2.0 * cos(MY_PI * 11.0 / 64.0)));
  63. static const real cos13_64 = (real)(1.0 / (2.0 * cos(MY_PI * 13.0 / 64.0)));
  64. static const real cos15_64 = (real)(1.0 / (2.0 * cos(MY_PI * 15.0 / 64.0)));
  65. static const real cos17_64 = (real)(1.0 / (2.0 * cos(MY_PI * 17.0 / 64.0)));
  66. static const real cos19_64 = (real)(1.0 / (2.0 * cos(MY_PI * 19.0 / 64.0)));
  67. static const real cos21_64 = (real)(1.0 / (2.0 * cos(MY_PI * 21.0 / 64.0)));
  68. static const real cos23_64 = (real)(1.0 / (2.0 * cos(MY_PI * 23.0 / 64.0)));
  69. static const real cos25_64 = (real)(1.0 / (2.0 * cos(MY_PI * 25.0 / 64.0)));
  70. static const real cos27_64 = (real)(1.0 / (2.0 * cos(MY_PI * 27.0 / 64.0)));
  71. static const real cos29_64 = (real)(1.0 / (2.0 * cos(MY_PI * 29.0 / 64.0)));
  72. static const real cos31_64 = (real)(1.0 / (2.0 * cos(MY_PI * 31.0 / 64.0)));
  73. static const real cos1_32  = (real)(1.0 / (2.0 * cos(MY_PI        / 32.0)));
  74. static const real cos3_32  = (real)(1.0 / (2.0 * cos(MY_PI * 3.0  / 32.0)));
  75. static const real cos5_32  = (real)(1.0 / (2.0 * cos(MY_PI * 5.0  / 32.0)));
  76. static const real cos7_32  = (real)(1.0 / (2.0 * cos(MY_PI * 7.0  / 32.0)));
  77. static const real cos9_32  = (real)(1.0 / (2.0 * cos(MY_PI * 9.0  / 32.0)));
  78. static const real cos11_32 = (real)(1.0 / (2.0 * cos(MY_PI * 11.0 / 32.0)));
  79. static const real cos13_32 = (real)(1.0 / (2.0 * cos(MY_PI * 13.0 / 32.0)));
  80. static const real cos15_32 = (real)(1.0 / (2.0 * cos(MY_PI * 15.0 / 32.0)));
  81. static const real cos1_16  = (real)(1.0 / (2.0 * cos(MY_PI        / 16.0)));
  82. static const real cos3_16  = (real)(1.0 / (2.0 * cos(MY_PI * 3.0  / 16.0)));
  83. static const real cos5_16  = (real)(1.0 / (2.0 * cos(MY_PI * 5.0  / 16.0)));
  84. static const real cos7_16  = (real)(1.0 / (2.0 * cos(MY_PI * 7.0  / 16.0)));
  85. static const real cos1_8   = (real)(1.0 / (2.0 * cos(MY_PI        / 8.0)));
  86. static const real cos3_8   = (real)(1.0 / (2.0 * cos(MY_PI * 3.0  / 8.0)));
  87. static const real cos1_4   = (real)(1.0 / (2.0 * cos(MY_PI / 4.0)));
  90. const real SynthesisFilter::d[512] = {
  91. // Note: These values are not in the same order
  92. // as in Annex 3-B.3 of the ISO/IEC DIS 11172-3
  93. 0.000000000f, -0.000442505f,  0.003250122f, -0.007003784f,
  94. 0.031082153f, -0.078628540f,  0.100311279f, -0.572036743f,
  95. 1.144989014f,  0.572036743f,  0.100311279f,  0.078628540f,
  96. 0.031082153f,  0.007003784f,  0.003250122f,  0.000442505f,
  97. -0.000015259f, -0.000473022f,  0.003326416f, -0.007919312f,
  98. 0.030517578f, -0.084182739f,  0.090927124f, -0.600219727f,
  99. 1.144287109f,  0.543823242f,  0.108856201f,  0.073059082f,
  100. 0.031478882f,  0.006118774f,  0.003173828f,  0.000396729f,
  101.   -0.000015259f, -0.000534058f,  0.003387451f, -0.008865356f,
  102. 0.029785156f, -0.089706421f,  0.080688477f, -0.628295898f,
  103. 1.142211914f,  0.515609741f,  0.116577148f,  0.067520142f,
  104.  0.031738281f,  0.005294800f,  0.003082275f,  0.000366211f,
  105. -0.000015259f, -0.000579834f,  0.003433228f, -0.009841919f,
  106.  0.028884888f, -0.095169067f,  0.069595337f, -0.656219482f,
  107. 1.138763428f,  0.487472534f,  0.123474121f,  0.061996460f,
  108.    0.031845093f,  0.004486084f,  0.002990723f,  0.000320435f,
  109.   -0.000015259f, -0.000625610f,  0.003463745f, -0.010848999f,
  110.  0.027801514f, -0.100540161f,  0.057617188f, -0.683914185f,
  111. 1.133926392f,  0.459472656f,  0.129577637f,  0.056533813f,
  112. 0.031814575f,  0.003723145f,  0.002899170f,  0.000289917f,
  113.   -0.000015259f, -0.000686646f,  0.003479004f, -0.011886597f,
  114. 0.026535034f, -0.105819702f,  0.044784546f, -0.711318970f,
  115. 1.127746582f,  0.431655884f,  0.134887695f,  0.051132202f,
  116. 0.031661987f,  0.003005981f,  0.002792358f,  0.000259399f,
  117. -0.000015259f, -0.000747681f,  0.003479004f, -0.012939453f,
  118. 0.025085449f, -0.110946655f,  0.031082153f, -0.738372803f,
  119.    1.120223999f,  0.404083252f,  0.139450073f,  0.045837402f,
  120.    0.031387329f,  0.002334595f,  0.002685547f,  0.000244141f,
  121.   -0.000030518f, -0.000808716f,  0.003463745f, -0.014022827f,
  122.    0.023422241f, -0.115921021f,  0.016510010f, -0.765029907f,
  123. 1.111373901f,  0.376800537f,  0.143264771f,  0.040634155f,
  124.    0.031005859f,  0.001693726f,  0.002578735f,  0.000213623f,
  125. -0.000030518f, -0.000885010f,  0.003417969f, -0.015121460f,
  126. 0.021575928f, -0.120697021f,  0.001068115f, -0.791213989f,
  127.    1.101211548f,  0.349868774f,  0.146362305f,  0.035552979f,
  128. 0.030532837f,  0.001098633f,  0.002456665f,  0.000198364f,
  129. -0.000030518f, -0.000961304f,  0.003372192f, -0.016235352f,
  130.    0.019531250f, -0.125259399f, -0.015228271f, -0.816864014f,
  131.  1.089782715f,  0.323318481f,  0.148773193f,  0.030609131f,
  132. 0.029937744f,  0.000549316f,  0.002349854f,  0.000167847f,
  133. -0.000030518f, -0.001037598f,  0.003280640f, -0.017349243f,
  134. 0.017257690f, -0.129562378f, -0.032379150f, -0.841949463f,
  135.    1.077117920f,  0.297210693f,  0.150497437f,  0.025817871f,
  136.    0.029281616f,  0.000030518f,  0.002243042f,  0.000152588f,
  137.   -0.000045776f, -0.001113892f,  0.003173828f, -0.018463135f,
  138. 0.014801025f, -0.133590698f, -0.050354004f, -0.866363525f,
  139. 1.063217163f,  0.271591187f,  0.151596069f,  0.021179199f,
  140. 0.028533936f, -0.000442505f,  0.002120972f,  0.000137329f,
  141. -0.000045776f, -0.001205444f,  0.003051758f, -0.019577026f,
  142. 0.012115479f, -0.137298584f, -0.069168091f, -0.890090942f,
  143.  1.048156738f,  0.246505737f,  0.152069092f,  0.016708374f,
  144. 0.027725220f, -0.000869751f,  0.002014160f,  0.000122070f,
  145.   -0.000061035f, -0.001296997f,  0.002883911f, -0.020690918f,
  146.    0.009231567f, -0.140670776f, -0.088775635f, -0.913055420f,
  147.  1.031936646f,  0.221984863f,  0.151962280f,  0.012420654f,
  148.    0.026840210f, -0.001266479f,  0.001907349f,  0.000106812f,
  149.   -0.000061035f, -0.001388550f,  0.002700806f, -0.021789551f,
  150. 0.006134033f, -0.143676758f, -0.109161377f, -0.935195923f,
  151.    1.014617920f,  0.198059082f,  0.151306152f,  0.008316040f,
  152. 0.025909424f, -0.001617432f,  0.001785278f,  0.000106812f,
  153. -0.000076294f, -0.001480103f,  0.002487183f, -0.022857666f,
  154.  0.002822876f, -0.146255493f, -0.130310059f, -0.956481934f,
  155. 0.996246338f,  0.174789429f,  0.150115967f,  0.004394531f,
  156.    0.024932861f, -0.001937866f,  0.001693726f,  0.000091553f,
  157.   -0.000076294f, -0.001586914f,  0.002227783f, -0.023910522f,
  158. -0.000686646f, -0.148422241f, -0.152206421f, -0.976852417f,
  159.    0.976852417f,  0.152206421f,  0.148422241f,  0.000686646f,
  160. 0.023910522f, -0.002227783f,  0.001586914f,  0.000076294f,
  161. -0.000091553f, -0.001693726f,  0.001937866f, -0.024932861f,
  162. -0.004394531f, -0.150115967f, -0.174789429f, -0.996246338f,
  163.    0.956481934f,  0.130310059f,  0.146255493f, -0.002822876f,
  164.    0.022857666f, -0.002487183f,  0.001480103f,  0.000076294f,
  165. -0.000106812f, -0.001785278f,  0.001617432f, -0.025909424f,
  166.   -0.008316040f, -0.151306152f, -0.198059082f, -1.014617920f,
  167.    0.935195923f,  0.109161377f,  0.143676758f, -0.006134033f,
  168.    0.021789551f, -0.002700806f,  0.001388550f,  0.000061035f,
  169.   -0.000106812f, -0.001907349f,  0.001266479f, -0.026840210f,
  170. -0.012420654f, -0.151962280f, -0.221984863f, -1.031936646f,
  171. 0.913055420f,  0.088775635f,  0.140670776f, -0.009231567f,
  172.  0.020690918f, -0.002883911f,  0.001296997f,  0.000061035f,
  173. -0.000122070f, -0.002014160f,  0.000869751f, -0.027725220f,
  174.   -0.016708374f, -0.152069092f, -0.246505737f, -1.048156738f,
  175.    0.890090942f,  0.069168091f,  0.137298584f, -0.012115479f,
  176. 0.019577026f, -0.003051758f,  0.001205444f,  0.000045776f,
  177. -0.000137329f, -0.002120972f,  0.000442505f, -0.028533936f,
  178.   -0.021179199f, -0.151596069f, -0.271591187f, -1.063217163f,
  179.    0.866363525f,  0.050354004f,  0.133590698f, -0.014801025f,
  180.    0.018463135f, -0.003173828f,  0.001113892f,  0.000045776f,
  181.   -0.000152588f, -0.002243042f, -0.000030518f, -0.029281616f,
  182.   -0.025817871f, -0.150497437f, -0.297210693f, -1.077117920f,
  183.  0.841949463f,  0.032379150f,  0.129562378f, -0.017257690f,
  184.  0.017349243f, -0.003280640f,  0.001037598f,  0.000030518f,
  185.   -0.000167847f, -0.002349854f, -0.000549316f, -0.029937744f,
  186.   -0.030609131f, -0.148773193f, -0.323318481f, -1.089782715f,
  187. 0.816864014f,  0.015228271f,  0.125259399f, -0.019531250f,
  188.    0.016235352f, -0.003372192f,  0.000961304f,  0.000030518f,
  189. -0.000198364f, -0.002456665f, -0.001098633f, -0.030532837f,
  190. -0.035552979f, -0.146362305f, -0.349868774f, -1.101211548f,
  191.  0.791213989f, -0.001068115f,  0.120697021f, -0.021575928f,
  192. 0.015121460f, -0.003417969f,  0.000885010f,  0.000030518f,
  193.   -0.000213623f, -0.002578735f, -0.001693726f, -0.031005859f,
  194. -0.040634155f, -0.143264771f, -0.376800537f, -1.111373901f,
  195.    0.765029907f, -0.016510010f,  0.115921021f, -0.023422241f,
  196.    0.014022827f, -0.003463745f,  0.000808716f,  0.000030518f,
  197.   -0.000244141f, -0.002685547f, -0.002334595f, -0.031387329f,
  198.   -0.045837402f, -0.139450073f, -0.404083252f, -1.120223999f,
  199.    0.738372803f, -0.031082153f,  0.110946655f, -0.025085449f,
  200. 0.012939453f, -0.003479004f,  0.000747681f,  0.000015259f,
  201. -0.000259399f, -0.002792358f, -0.003005981f, -0.031661987f,
  202.   -0.051132202f, -0.134887695f, -0.431655884f, -1.127746582f,
  203. 0.711318970f, -0.044784546f,  0.105819702f, -0.026535034f,
  204.    0.011886597f, -0.003479004f,  0.000686646f,  0.000015259f,
  205. -0.000289917f, -0.002899170f, -0.003723145f, -0.031814575f,
  206. -0.056533813f, -0.129577637f, -0.459472656f, -1.133926392f,
  207.    0.683914185f, -0.057617188f,  0.100540161f, -0.027801514f,
  208. 0.010848999f, -0.003463745f,  0.000625610f,  0.000015259f,
  209.   -0.000320435f, -0.002990723f, -0.004486084f, -0.031845093f,
  210.   -0.061996460f, -0.123474121f, -0.487472534f, -1.138763428f,
  211. 0.656219482f, -0.069595337f,  0.095169067f, -0.028884888f,
  212. 0.009841919f, -0.003433228f,  0.000579834f,  0.000015259f,
  213. -0.000366211f, -0.003082275f, -0.005294800f, -0.031738281f,
  214.   -0.067520142f, -0.116577148f, -0.515609741f, -1.142211914f,
  215. 0.628295898f, -0.080688477f,  0.089706421f, -0.029785156f,
  216. 0.008865356f, -0.003387451f,  0.000534058f,  0.000015259f,
  217.   -0.000396729f, -0.003173828f, -0.006118774f, -0.031478882f,
  218.   -0.073059082f, -0.108856201f, -0.543823242f, -1.144287109f,
  219. 0.600219727f, -0.090927124f,  0.084182739f, -0.030517578f,
  220. 0.007919312f, -0.003326416f,  0.000473022f,  0.000015259f
  221. };
  223. SynthesisFilter::SynthesisFilter(uint32 channelnumber, real factor)
  224. {
  225.   channel = channelnumber;
  226.   scalefactor = factor;
  228.   reset();
  229. }
  231. void SynthesisFilter::reset()
  232. {
  233.   register real *floatp, *floatp2;
  235.   // initialize v1[] and v2[]:
  236.   for (floatp = v1 + 512, floatp2 = v2 + 512; floatp > v1; )
  237.  *--floatp = *--floatp2 = 0.0;
  239.   // initialize samples[]:
  240.   for (floatp = samples + 32; floatp > samples; )
  241.  *--floatp = 0.0;
  243.   actual_v = v1;
  244.   actual_write_pos = 15;
  245. }
  247. void SynthesisFilter::compute_new_v()
  248. {
  249. real new_v[32];         // new V[0-15] and V[33-48] of Figure 3-A.2 in ISO DIS 11172-3
  250. real p[16];
  251. real pp[16];
  253. // compute new values via a fast cosine transform:
  254. {
  255. register real *x1 = samples;
  257. p[0] = x1[0] + x1[31];
  258. p[1] = x1[1] + x1[30];
  259. p[2] = x1[2] + x1[29];
  260. p[3] = x1[3] + x1[28];
  261. p[4] = x1[4] + x1[27];
  262. p[5] = x1[5] + x1[26];
  263. p[6] = x1[6] + x1[25];
  264. p[7] = x1[7] + x1[24];
  265. p[8] = x1[8] + x1[23];
  266. p[9] = x1[9] + x1[22];
  267. p[10] = x1[10] + x1[21];
  268. p[11] = x1[11] + x1[20];
  269. p[12] = x1[12] + x1[19];
  270. p[13] = x1[13] + x1[18];
  271. p[14] = x1[14] + x1[17];
  272. p[15] = x1[15] + x1[16];
  273. }
  275. {
  276. pp[0] = p[0] + p[15];
  277. pp[1] = p[1] + p[14];
  278. pp[2] = p[2] + p[13];
  279. pp[3] = p[3] + p[12];
  280. pp[4] = p[4] + p[11];
  281. pp[5] = p[5] + p[10];
  282. pp[6] = p[6] + p[9];
  283. pp[7] = p[7] + p[8];
  284. pp[8] = (p[0] - p[15]) * cos1_32;
  285. pp[9] = (p[1] - p[14]) * cos3_32;
  286. pp[10] = (p[2] - p[13]) * cos5_32;
  287. pp[11] = (p[3] - p[12]) * cos7_32;
  288. pp[12] = (p[4] - p[11]) * cos9_32;
  289. pp[13] = (p[5] - p[10]) * cos11_32;
  290. pp[14] = (p[6] - p[9]) * cos13_32;
  291. pp[15] = (p[7] - p[8]) * cos15_32;
  292. }
  294. {
  295. p[0] = pp[0] + pp[7];
  296. p[1] = pp[1] + pp[6];
  297. p[2] = pp[2] + pp[5];
  298. p[3] = pp[3] + pp[4];
  299. p[4] = (pp[0] - pp[7]) * cos1_16;
  300. p[5] = (pp[1] - pp[6]) * cos3_16;
  301. p[6] = (pp[2] - pp[5]) * cos5_16;
  302. p[7] = (pp[3] - pp[4]) * cos7_16;
  303. p[8] = pp[8] + pp[15];
  304. p[9] = pp[9] + pp[14];
  305. p[10] = pp[10] + pp[13];
  306. p[11] = pp[11] + pp[12];
  307. p[12] = (pp[8] - pp[15]) * cos1_16;
  308. p[13] = (pp[9] - pp[14]) * cos3_16;
  309. p[14] = (pp[10] - pp[13]) * cos5_16;
  310. p[15] = (pp[11] - pp[12]) * cos7_16;
  311. }
  313. {
  314. pp[0] = p[0] + p[3];
  315. pp[1] = p[1] + p[2];
  316. pp[2] = (p[0] - p[3]) * cos1_8;
  317. pp[3] = (p[1] - p[2]) * cos3_8;
  318. pp[4] = p[4] + p[7];
  319. pp[5] = p[5] + p[6];
  320. pp[6] = (p[4] - p[7]) * cos1_8;
  321. pp[7] = (p[5] - p[6]) * cos3_8;
  322. pp[8] = p[8] + p[11];
  323. pp[9] = p[9] + p[10];
  324. pp[10] = (p[8] - p[11]) * cos1_8;
  325. pp[11] = (p[9] - p[10]) * cos3_8;
  326. pp[12] = p[12] + p[15];
  327. pp[13] = p[13] + p[14];
  328. pp[14] = (p[12] - p[15]) * cos1_8;
  329. pp[15] = (p[13] - p[14]) * cos3_8;
  330. }
  332. {
  333. p[0] = pp[0] + pp[1];
  334. p[1] = (pp[0] - pp[1]) * cos1_4;
  335. p[2] = pp[2] + pp[3];
  336. p[3] = (pp[2] - pp[3]) * cos1_4;
  337. p[4] = pp[4] + pp[5];
  338. p[5] = (pp[4] - pp[5]) * cos1_4;
  339. p[6] = pp[6] + pp[7];
  340. p[7] = (pp[6] - pp[7]) * cos1_4;
  341. p[8] = pp[8] + pp[9];
  342. p[9] = (pp[8] - pp[9]) * cos1_4;
  343. p[10] = pp[10] + pp[11];
  344. p[11] = (pp[10] - pp[11]) * cos1_4;
  345. p[12] = pp[12] + pp[13];
  346. p[13] = (pp[12] - pp[13]) * cos1_4;
  347. p[14] = pp[14] + pp[15];
  348. p[15] = (pp[14] - pp[15]) * cos1_4;
  349. }
  351. {
  352. // this is pretty insane coding
  353. register real tmp1;
  354. new_v[36-17] = -(new_v[4] = (new_v[12] = p[7]) + p[5]) - p[6];
  355. new_v[44-17] = -p[6] - p[7] - p[4];
  356. new_v[6] = (new_v[10] = (new_v[14] = p[15]) + p[11]) + p[13];
  357. new_v[34-17] = -(new_v[2] = p[15] + p[13] + p[9]) - p[14];
  358. new_v[38-17] = (tmp1 = -p[14] - p[15] - p[10] - p[11]) - p[13];
  359. new_v[46-17] = -p[14] - p[15] - p[12] - p[8];
  360. new_v[42-17] = tmp1 - p[12];
  361. new_v[48-17] = -p[0];
  362. new_v[0] = p[1];
  363. new_v[40-17] = -(new_v[8] = p[3]) - p[2];
  364. }
  366. {
  367. register real *x1 = samples;
  369. p[0] = (x1[0] - x1[31]) * cos1_64;
  370. p[1] = (x1[1] - x1[30]) * cos3_64;
  371. p[2] = (x1[2] - x1[29]) * cos5_64;
  372. p[3] = (x1[3] - x1[28]) * cos7_64;
  373. p[4] = (x1[4] - x1[27]) * cos9_64;
  374. p[5] = (x1[5] - x1[26]) * cos11_64;
  375. p[6] = (x1[6] - x1[25]) * cos13_64;
  376. p[7] = (x1[7] - x1[24]) * cos15_64;
  377. p[8] = (x1[8] - x1[23]) * cos17_64;
  378. p[9] = (x1[9] - x1[22]) * cos19_64;
  379. p[10] = (x1[10] - x1[21]) * cos21_64;
  380. p[11] = (x1[11] - x1[20]) * cos23_64;
  381. p[12] = (x1[12] - x1[19]) * cos25_64;
  382. p[13] = (x1[13] - x1[18]) * cos27_64;
  383. p[14] = (x1[14] - x1[17]) * cos29_64;
  384. p[15] = (x1[15] - x1[16]) * cos31_64;
  385. }
  387. {
  388. pp[0] = p[0] + p[15];
  389. pp[1] = p[1] + p[14];
  390. pp[2] = p[2] + p[13];
  391. pp[3] = p[3] + p[12];
  392. pp[4] = p[4] + p[11];
  393. pp[5] = p[5] + p[10];
  394. pp[6] = p[6] + p[9];
  395. pp[7] = p[7] + p[8];
  396. pp[8] = (p[0] - p[15]) * cos1_32;
  397. pp[9] = (p[1] - p[14]) * cos3_32;
  398. pp[10] = (p[2] - p[13]) * cos5_32;
  399. pp[11] = (p[3] - p[12]) * cos7_32;
  400. pp[12] = (p[4] - p[11]) * cos9_32;
  401. pp[13] = (p[5] - p[10]) * cos11_32;
  402. pp[14] = (p[6] - p[9]) * cos13_32;
  403. pp[15] = (p[7] - p[8]) * cos15_32;
  404. }
  406. {
  407. p[0] = pp[0] + pp[7];
  408. p[1] = pp[1] + pp[6];
  409. p[2] = pp[2] + pp[5];
  410. p[3] = pp[3] + pp[4];
  411. p[4] = (pp[0] - pp[7]) * cos1_16;
  412. p[5] = (pp[1] - pp[6]) * cos3_16;
  413. p[6] = (pp[2] - pp[5]) * cos5_16;
  414. p[7] = (pp[3] - pp[4]) * cos7_16;
  415. p[8] = pp[8] + pp[15];
  416. p[9] = pp[9] + pp[14];
  417. p[10] = pp[10] + pp[13];
  418. p[11] = pp[11] + pp[12];
  419. p[12] = (pp[8] - pp[15]) * cos1_16;
  420. p[13] = (pp[9] - pp[14]) * cos3_16;
  421. p[14] = (pp[10] - pp[13]) * cos5_16;
  422. p[15] = (pp[11] - pp[12]) * cos7_16;
  423. }
  425. {
  426. pp[0] = p[0] + p[3];
  427. pp[1] = p[1] + p[2];
  428. pp[2] = (p[0] - p[3]) * cos1_8;
  429. pp[3] = (p[1] - p[2]) * cos3_8;
  430. pp[4] = p[4] + p[7];
  431. pp[5] = p[5] + p[6];
  432. pp[6] = (p[4] - p[7]) * cos1_8;
  433. pp[7] = (p[5] - p[6]) * cos3_8;
  434. pp[8] = p[8] + p[11];
  435. pp[9] = p[9] + p[10];
  436. pp[10] = (p[8] - p[11]) * cos1_8;
  437. pp[11] = (p[9] - p[10]) * cos3_8;
  438. pp[12] = p[12] + p[15];
  439. pp[13] = p[13] + p[14];
  440. pp[14] = (p[12] - p[15]) * cos1_8;
  441. pp[15] = (p[13] - p[14]) * cos3_8;
  442. }
  444. {
  445. p[0] = pp[0] + pp[1];
  446. p[1] = (pp[0] - pp[1]) * cos1_4;
  447. p[2] = pp[2] + pp[3];
  448. p[3] = (pp[2] - pp[3]) * cos1_4;
  449. p[4] = pp[4] + pp[5];
  450. p[5] = (pp[4] - pp[5]) * cos1_4;
  451. p[6] = pp[6] + pp[7];
  452. p[7] = (pp[6] - pp[7]) * cos1_4;
  453. p[8] = pp[8] + pp[9];
  454. p[9] = (pp[8] - pp[9]) * cos1_4;
  455. p[10] = pp[10] + pp[11];
  456. p[11] = (pp[10] - pp[11]) * cos1_4;
  457. p[12] = pp[12] + pp[13];
  458. p[13] = (pp[12] - pp[13]) * cos1_4;
  459. p[14] = pp[14] + pp[15];
  460. p[15] = (pp[14] - pp[15]) * cos1_4;
  461. }
  463. {
  464. // manually doing something that a compiler should handle sucks
  465. // coding like this is hard to read
  466. register real tmp1, tmp2;
  467. new_v[5] = (new_v[11] = (new_v[13] = (new_v[15] = p[15]) + p[7]) + p[11])
  468. + p[5] + p[13];
  469. new_v[7] = (new_v[9] = p[15] + p[11] + p[3]) + p[13];
  470. new_v[33-17] = -(new_v[1] = (tmp1 = p[13] + p[15] + p[9]) + p[1]) - p[14];
  471. new_v[35-17] = -(new_v[3] = tmp1 + p[5] + p[7]) - p[6] - p[14];
  473. new_v[39-17] = (tmp1 = -p[10] - p[11] - p[14] - p[15])
  474. - p[13] - p[2] - p[3];
  475. new_v[37-17] = tmp1 - p[13] - p[5] - p[6] - p[7];
  476. new_v[41-17] = tmp1 - p[12] - p[2] - p[3];
  477. new_v[43-17] = tmp1 - p[12] - (tmp2 = p[4] + p[6] + p[7]);
  478. new_v[47-17] = (tmp1 = -p[8] - p[12] - p[14] - p[15]) - p[0];
  479. new_v[45-17] = tmp1 - tmp2;
  480. }
  482. {
  483. // insert V[0-15] (== new_v[0-15]) into actual v:
  484. register real *x1 = new_v;
  485. register real *x2 = actual_v + actual_write_pos;
  486. x2[0] = x1[0];
  487. x2[16] = x1[1];
  488. x2[32] = x1[2];
  489. x2[48] = x1[3];
  490. x2[64] = x1[4];
  491. x2[80] = x1[5];
  492. x2[96] = x1[6];
  493. x2[112] = x1[7];
  494. x2[128] = x1[8];
  495. x2[144] = x1[9];
  496. x2[160] = x1[10];
  497. x2[176] = x1[11];
  498. x2[192] = x1[12];
  499. x2[208] = x1[13];
  500. x2[224] = x1[14];
  501. x2[240] = x1[15];
  503. // V[16] is always 0.0:
  504. x2[256] = 0.0f;
  506. // insert V[17-31] (== -new_v[15-1]) into actual v:
  507. x2[272] = -x1[15];
  508. x2[288] = -x1[14];
  509. x2[304] = -x1[13];
  510. x2[320] = -x1[12];
  511. x2[336] = -x1[11];
  512. x2[352] = -x1[10];
  513. x2[368] = -x1[9];
  514. x2[384] = -x1[8];
  515. x2[400] = -x1[7];
  516. x2[416] = -x1[6];
  517. x2[432] = -x1[5];
  518. x2[448] = -x1[4];
  519. x2[464] = -x1[3];
  520. x2[480] = -x1[2];
  521. x2[496] = -x1[1];
  523. // insert V[32] (== -new_v[0]) into other v:
  524. x2 = (actual_v == v1 ? v2 : v1) + actual_write_pos;
  525. x2[0] = -x1[0];
  527. // insert V[33-48] (== new_v[16-31]) into other v:
  528. x2[16] = x1[16];
  529. x2[32] = x1[17];
  530. x2[48] = x1[18];
  531. x2[64] = x1[19];
  532. x2[80] = x1[20];
  533. x2[96] = x1[21];
  534. x2[112] = x1[22];
  535. x2[128] = x1[23];
  536. x2[144] = x1[24];
  537. x2[160] = x1[25];
  538. x2[176] = x1[26];
  539. x2[192] = x1[27];
  540. x2[208] = x1[28];
  541. x2[224] = x1[29];
  542. x2[240] = x1[30];
  543. x2[256] = x1[31];
  545. // insert V[49-63] (== new_v[30-16]) into other v:
  546. x2[272] = x1[30];
  547. x2[288] = x1[29];
  548. x2[304] = x1[28];
  549. x2[320] = x1[27];
  550. x2[336] = x1[26];
  551. x2[352] = x1[25];
  552. x2[368] = x1[24];
  553. x2[384] = x1[23];
  554. x2[400] = x1[22];
  555. x2[416] = x1[21];
  556. x2[432] = x1[20];
  557. x2[448] = x1[19];
  558. x2[464] = x1[18];
  559. x2[480] = x1[17];
  560. x2[496] = x1[16];
  561. }
  562. }
  564. void SynthesisFilter::compute_pcm_samples(Obuffer *buffer)
  565. {
  566. // scoping variables makes it easier to optimize for the compiler
  567. register real *vp = actual_v;
  568. //real a[32];
  569. //real * pa=a;
  570. // switch depending on the value for actual_write_pos
  571. switch (actual_write_pos) {
  572. case 0: {
  573. // fat chance of having this loop unroll
  574. for( register const real *dp = d;
  575.  dp < d + (32 * 16);
  576.  dp += 16, vp += 16 ) {
  577. register real pcm_sample;
  579. pcm_sample = (real)(((vp[0] * dp[0]) +
  580. (vp[15] * dp[1]) +
  581. (vp[14] * dp[2]) +
  582. (vp[13] * dp[3]) +
  583. (vp[12] * dp[4]) +
  584. (vp[11] * dp[5]) +
  585. (vp[10] * dp[6]) +
  586. (vp[9] * dp[7]) +
  587. (vp[8] * dp[8]) +
  588. (vp[7] * dp[9]) +
  589. (vp[6] * dp[10]) +
  590. (vp[5] * dp[11]) +
  591. (vp[4] * dp[12]) +
  592. (vp[3] * dp[13]) +
  593. (vp[2] * dp[14]) +
  594. (vp[1] * dp[15])
  595. ) /* * scalefactor */);
  596.             buffer->append (channel, (int16)pcm_sample);
  597. } // for
  598. } break; // case 0:
  600. case 1: {
  601. // fat chance of having this loop unroll
  602. for( register const real *dp = d;
  603.  dp < d + (32 * 16);
  604.  dp += 16, vp += 16 ) {
  605. register real pcm_sample;
  607. pcm_sample = (real)(((vp[1] * dp[0]) +
  608. (vp[0] * dp[1]) +
  609. (vp[15] * dp[2]) +
  610. (vp[14] * dp[3]) +
  611. (vp[13] * dp[4]) +
  612. (vp[12] * dp[5]) +
  613. (vp[11] * dp[6]) +
  614. (vp[10] * dp[7]) +
  615. (vp[9] * dp[8]) +
  616. (vp[8] * dp[9]) +
  617. (vp[7] * dp[10]) +
  618. (vp[6] * dp[11]) +
  619. (vp[5] * dp[12]) +
  620. (vp[4] * dp[13]) +
  621. (vp[3] * dp[14]) +
  622. (vp[2] * dp[15])
  623. ) /* * scalefactor */);
  625.             buffer->append (channel, (int16)pcm_sample);
  626.          } // for
  627. } break; // case 1:
  629. case 2: {
  630. // fat chance of having this loop unroll
  631. for( register const real *dp = d;
  632.  dp < d + (32 * 16);
  633.  dp += 16, vp += 16 ) {
  634. register real pcm_sample;
  636. pcm_sample = (real)(((vp[2] * dp[0]) +
  637. (vp[1] * dp[1]) +
  638. (vp[0] * dp[2]) +
  639. (vp[15] * dp[3]) +
  640. (vp[14] * dp[4]) +
  641. (vp[13] * dp[5]) +
  642. (vp[12] * dp[6]) +
  643. (vp[11] * dp[7]) +
  644. (vp[10] * dp[8]) +
  645. (vp[9] * dp[9]) +
  646. (vp[8] * dp[10]) +
  647. (vp[7] * dp[11]) +
  648. (vp[6] * dp[12]) +
  649. (vp[5] * dp[13]) +
  650. (vp[4] * dp[14]) +
  651. (vp[3] * dp[15])
  652. ) /* * scalefactor */);
  654.             buffer->append (channel, (int16)pcm_sample);
  655. } // for
  656. } break; // case 2:
  658. case 3: {
  659. // fat chance of having this loop unroll
  660. for( register const real *dp = d;
  661.  dp < d + (32 * 16);
  662.  dp += 16, vp += 16 ) {
  663. register real pcm_sample;
  665. pcm_sample = (real)(((vp[3] * dp[0]) +
  666. (vp[2] * dp[1]) +
  667. (vp[1] * dp[2]) +
  668. (vp[0] * dp[3]) +
  669. (vp[15] * dp[4]) +
  670. (vp[14] * dp[5]) +
  671. (vp[13] * dp[6]) +
  672. (vp[12] * dp[7]) +
  673. (vp[11] * dp[8]) +
  674. (vp[10] * dp[9]) +
  675. (vp[9] * dp[10]) +
  676. (vp[8] * dp[11]) +
  677. (vp[7] * dp[12]) +
  678. (vp[6] * dp[13]) +
  679. (vp[5] * dp[14]) +
  680. (vp[4] * dp[15])
  681. ) /* * scalefactor */);
  683. buffer->append (channel, (int16)pcm_sample);
  684. } // for
  685. } break; // case 3:
  687. case 4: {
  688. // fat chance of having this loop unroll
  689. for( register const real *dp = d;
  690.  dp < d + (32 * 16);
  691.  dp += 16, vp += 16 ) {
  692. register real pcm_sample;
  694. pcm_sample = (real)(((vp[4] * dp[0]) +
  695. (vp[3] * dp[1]) +
  696. (vp[2] * dp[2]) +
  697. (vp[1] * dp[3]) +
  698. (vp[0] * dp[4]) +
  699. (vp[15] * dp[5]) +
  700. (vp[14] * dp[6]) +
  701. (vp[13] * dp[7]) +
  702. (vp[12] * dp[8]) +
  703. (vp[11] * dp[9]) +
  704. (vp[10] * dp[10]) +
  705. (vp[9] * dp[11]) +
  706. (vp[8] * dp[12]) +
  707. (vp[7] * dp[13]) +
  708. (vp[6] * dp[14]) +
  709. (vp[5] * dp[15])
  710. ) /* * scalefactor */);
  712. buffer->append (channel, (int16)pcm_sample);
  713. } // for
  714. } break; // case 4:
  716. case 5: {
  717. // fat chance of having this loop unroll
  718. for( register const real *dp = d;
  719.  dp < d + (32 * 16);
  720.  dp += 16, vp += 16 ) {
  721. register real pcm_sample;
  723. pcm_sample = (real)(((vp[5] * dp[0]) +
  724. (vp[4] * dp[1]) +
  725. (vp[3] * dp[2]) +
  726. (vp[2] * dp[3]) +
  727. (vp[1] * dp[4]) +
  728. (vp[0] * dp[5]) +
  729. (vp[15] * dp[6]) +
  730. (vp[14] * dp[7]) +
  731. (vp[13] * dp[8]) +
  732. (vp[12] * dp[9]) +
  733. (vp[11] * dp[10]) +
  734. (vp[10] * dp[11]) +
  735. (vp[9] * dp[12]) +
  736. (vp[8] * dp[13]) +
  737. (vp[7] * dp[14]) +
  738. (vp[6] * dp[15])
  739. ) /* * scalefactor */);
  741. buffer->append (channel, (int16)pcm_sample);
  742. } // for
  743. } break; // case 5:
  745. case 6: {
  746. // fat chance of having this loop unroll
  747. for( register const real *dp = d;
  748.  dp < d + (32 * 16);
  749.  dp += 16, vp += 16 ) {
  750. register real pcm_sample;
  752. pcm_sample = (real)(((vp[6] * dp[0]) +
  753. (vp[5] * dp[1]) +
  754. (vp[4] * dp[2]) +
  755. (vp[3] * dp[3]) +
  756. (vp[2] * dp[4]) +
  757. (vp[1] * dp[5]) +
  758. (vp[0] * dp[6]) +
  759. (vp[15] * dp[7]) +
  760. (vp[14] * dp[8]) +
  761. (vp[13] * dp[9]) +
  762. (vp[12] * dp[10]) +
  763. (vp[11] * dp[11]) +
  764. (vp[10] * dp[12]) +
  765. (vp[9] * dp[13]) +
  766. (vp[8] * dp[14]) +
  767. (vp[7] * dp[15])
  768. ) /* * scalefactor */);
  770.             buffer->append (channel, (int16)pcm_sample);
  771. } // for
  772. } break; // case 6:
  774. case 7: {
  775. // fat chance of having this loop unroll
  776. for( register const real *dp = d;
  777.  dp < d + (32 * 16);
  778.  dp += 16, vp += 16 ) {
  779. register real pcm_sample;
  781. pcm_sample = (real)(((vp[7] * dp[0]) +
  782. (vp[6] * dp[1]) +
  783. (vp[5] * dp[2]) +
  784. (vp[4] * dp[3]) +
  785. (vp[3] * dp[4]) +
  786. (vp[2] * dp[5]) +
  787. (vp[1] * dp[6]) +
  788. (vp[0] * dp[7]) +
  789. (vp[15] * dp[8]) +
  790. (vp[14] * dp[9]) +
  791. (vp[13] * dp[10]) +
  792. (vp[12] * dp[11]) +
  793. (vp[11] * dp[12]) +
  794. (vp[10] * dp[13]) +
  795. (vp[9] * dp[14]) +
  796. (vp[8] * dp[15])
  797. ) /* * scalefactor */);
  799.             buffer->append (channel, (int16)pcm_sample);
  800. } // for
  801. } break; // case 7:
  803. case 8: {
  804. // fat chance of having this loop unroll
  805. for( register const real *dp = d;
  806.  dp < d + (32 * 16);
  807.  dp += 16, vp += 16 ) {
  808. register real pcm_sample;
  810. pcm_sample = (real)(((vp[8] * dp[0]) +
  811. (vp[7] * dp[1]) +
  812. (vp[6] * dp[2]) +
  813. (vp[5] * dp[3]) +
  814. (vp[4] * dp[4]) +
  815. (vp[3] * dp[5]) +
  816. (vp[2] * dp[6]) +
  817. (vp[1] * dp[7]) +
  818. (vp[0] * dp[8]) +
  819. (vp[15] * dp[9]) +
  820. (vp[14] * dp[10]) +
  821. (vp[13] * dp[11]) +
  822. (vp[12] * dp[12]) +
  823. (vp[11] * dp[13]) +
  824. (vp[10] * dp[14]) +
  825. (vp[9] * dp[15])
  826. ) /* * scalefactor */);
  828.             buffer->append (channel, (int16)pcm_sample);
  829. } // for
  830. } break; // case 8:
  832. case 9: {
  833. // fat chance of having this loop unroll
  834. for( register const real *dp = d;
  835.  dp < d + (32 * 16);
  836.  dp += 16, vp += 16 ) {
  837. register real pcm_sample;
  839. pcm_sample = (real)(((vp[9] * dp[0]) +
  840. (vp[8] * dp[1]) +
  841. (vp[7] * dp[2]) +
  842. (vp[6] * dp[3]) +
  843. (vp[5] * dp[4]) +
  844. (vp[4] * dp[5]) +
  845. (vp[3] * dp[6]) +
  846. (vp[2] * dp[7]) +
  847. (vp[1] * dp[8]) +
  848. (vp[0] * dp[9]) +
  849. (vp[15] * dp[10]) +
  850. (vp[14] * dp[11]) +
  851. (vp[13] * dp[12]) +
  852. (vp[12] * dp[13]) +
  853. (vp[11] * dp[14]) +
  854. (vp[10] * dp[15])
  855. ) /* * scalefactor */);
  857.             buffer->append (channel, (int16)pcm_sample);
  858. } // for
  859. } break; // case 9:
  861. case 10: {
  862. // fat chance of having this loop unroll
  863. for( register const real *dp = d;
  864.  dp < d + (32 * 16);
  865.  dp += 16, vp += 16 ) {
  866. register real pcm_sample;
  868. pcm_sample = (real)(((vp[10] * dp[0]) +
  869. (vp[9] * dp[1]) +
  870. (vp[8] * dp[2]) +
  871. (vp[7] * dp[3]) +
  872. (vp[6] * dp[4]) +
  873. (vp[5] * dp[5]) +
  874. (vp[4] * dp[6]) +
  875. (vp[3] * dp[7]) +
  876. (vp[2] * dp[8]) +
  877. (vp[1] * dp[9]) +
  878. (vp[0] * dp[10]) +
  879. (vp[15] * dp[11]) +
  880. (vp[14] * dp[12]) +
  881. (vp[13] * dp[13]) +
  882. (vp[12] * dp[14]) +
  883. (vp[11] * dp[15])
  884. ) /* * scalefactor */);
  886.             buffer->append (channel, (int16)pcm_sample);
  887. } // for
  888. } break; // case 10:
  890. case 11: {
  891. // fat chance of having this loop unroll
  892. for( register const real *dp = d;
  893.  dp < d + (32 * 16);
  894.  dp += 16, vp += 16 ) {
  895. register real pcm_sample;
  897. pcm_sample = (real)(((vp[11] * dp[0]) +
  898. (vp[10] * dp[1]) +
  899. (vp[9] * dp[2]) +
  900. (vp[8] * dp[3]) +
  901. (vp[7] * dp[4]) +
  902. (vp[6] * dp[5]) +
  903. (vp[5] * dp[6]) +
  904. (vp[4] * dp[7]) +
  905. (vp[3] * dp[8]) +
  906. (vp[2] * dp[9]) +
  907. (vp[1] * dp[10]) +
  908. (vp[0] * dp[11]) +
  909. (vp[15] * dp[12]) +
  910. (vp[14] * dp[13]) +
  911. (vp[13] * dp[14]) +
  912. (vp[12] * dp[15])
  913. ) /* * scalefactor */);
  915.             buffer->append (channel, (int16)pcm_sample);
  916. } // for
  917. } break; // case 11:
  919. case 12: {
  920. // fat chance of having this loop unroll
  921. for( register const real *dp = d;
  922.  dp < d + (32 * 16);
  923.  dp += 16, vp += 16 ) {
  924. register real pcm_sample;
  926. pcm_sample = (real)(((vp[12] * dp[0]) +
  927. (vp[11] * dp[1]) +
  928. (vp[10] * dp[2]) +
  929. (vp[9] * dp[3]) +
  930. (vp[8] * dp[4]) +
  931. (vp[7] * dp[5]) +
  932. (vp[6] * dp[6]) +
  933. (vp[5] * dp[7]) +
  934. (vp[4] * dp[8]) +
  935. (vp[3] * dp[9]) +
  936. (vp[2] * dp[10]) +
  937. (vp[1] * dp[11]) +
  938. (vp[0] * dp[12]) +
  939. (vp[15] * dp[13]) +
  940. (vp[14] * dp[14]) +
  941. (vp[13] * dp[15])
  942. ) /* * scalefactor */);
  944.             buffer->append (channel, (int16)pcm_sample);
  945. } // for
  946. } break; // case 12:
  948. case 13: {
  949. // fat chance of having this loop unroll
  950. for( register const real *dp = d;
  951.  dp < d + (32 * 16);
  952.  dp += 16, vp += 16 ) {
  953. register real pcm_sample;
  955. pcm_sample = (real)(((vp[13] * dp[0]) +
  956. (vp[12] * dp[1]) +
  957. (vp[11] * dp[2]) +
  958. (vp[10] * dp[3]) +
  959. (vp[9] * dp[4]) +
  960. (vp[8] * dp[5]) +
  961. (vp[7] * dp[6]) +
  962. (vp[6] * dp[7]) +
  963. (vp[5] * dp[8]) +
  964. (vp[4] * dp[9]) +
  965. (vp[3] * dp[10]) +
  966. (vp[2] * dp[11]) +
  967. (vp[1] * dp[12]) +
  968. (vp[0] * dp[13]) +
  969. (vp[15] * dp[14]) +
  970. (vp[14] * dp[15])
  971. ) /* * scalefactor */);
  973.             buffer->append (channel, (int16)pcm_sample);
  974. } // for
  975. } break; // case 13:
  977. case 14: {
  978. // fat chance of having this loop unroll
  979. for( register const real *dp = d;
  980.  dp < d + (32 * 16);
  981.  dp += 16, vp += 16 ) {
  982. register real pcm_sample;
  984. pcm_sample = (real)(((vp[14] * dp[0]) +
  985. (vp[13] * dp[1]) +
  986. (vp[12] * dp[2]) +
  987. (vp[11] * dp[3]) +
  988. (vp[10] * dp[4]) +
  989. (vp[9] * dp[5]) +
  990. (vp[8] * dp[6]) +
  991. (vp[7] * dp[7]) +
  992. (vp[6] * dp[8]) +
  993. (vp[5] * dp[9]) +
  994. (vp[4] * dp[10]) +
  995. (vp[3] * dp[11]) +
  996. (vp[2] * dp[12]) +
  997. (vp[1] * dp[13]) +
  998. (vp[0] * dp[14]) +
  999. (vp[15] * dp[15])
  1000. ) /* * scalefactor */);
  1002.             buffer->append (channel, (int16)pcm_sample);
  1003. } // for
  1004. } break; // case 14:
  1006. case 15: {
  1007. // fat chance of having this loop unroll
  1008. for( register const real *dp = d;
  1009.  dp < d + (32 * 16);
  1010.  dp += 16, vp += 16 ) {
  1011. register real pcm_sample;
  1013. pcm_sample = (real)(((vp[15] * dp[0]) +
  1014. (vp[14] * dp[1]) +
  1015. (vp[13] * dp[2]) +
  1016. (vp[12] * dp[3]) +
  1017. (vp[11] * dp[4]) +
  1018. (vp[10] * dp[5]) +
  1019. (vp[9] * dp[6]) +
  1020. (vp[8] * dp[7]) +
  1021. (vp[7] * dp[8]) +
  1022. (vp[6] * dp[9]) +
  1023. (vp[5] * dp[10]) +
  1024. (vp[4] * dp[11]) +
  1025. (vp[3] * dp[12]) +
  1026. (vp[2] * dp[13]) +
  1027. (vp[1] * dp[14]) +
  1028. (vp[0] * dp[15])
  1029. ) /* * scalefactor */);
  1031.             buffer->append (channel, (int16)pcm_sample);
  1032. } // for
  1033. } break; // case 15:
  1034. }; // switch (actual_write_pos)
  1036. //buffer->appendblock (channel, a, 32);
  1038. }
  1040. void SynthesisFilter::calculate_pcm_samples(Obuffer *buffer)
  1041. {
  1042. compute_new_v();
  1043. compute_pcm_samples(buffer);
  1045.    actual_write_pos = (actual_write_pos + 1) & 0xf;
  1046. actual_v = (actual_v == v1) ? v2 : v1;
  1048. // initialize samples[]:
  1049.   for (register real *floatp = samples + 32; floatp > samples; )
  1050.  *--floatp = 0.0f;
  1051. }