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

源码开发语言/平台
当前位置：首页 > 源码/资料 > 嵌入式/单片机编程 > Windows CE > 查看源码
synth_filter.c
资源名称：tcpmp.rar　[点击查看]
上传用户：wstnjxml
上传日期：2014-04-03
资源大小：7248k
文件大小：23k
源码类别：

Windows CE

开发平台：

C/C++

synth_filter.c：源码内容
  1. /*
  2.   Copyright (c) 2005, The Musepack Development Team
  3.   All rights reserved.
  5.   Redistribution and use in source and binary forms, with or without
  6.   modification, are permitted provided that the following conditions are
  7.   met:
  9.   * Redistributions of source code must retain the above copyright
  10.   notice, this list of conditions and the following disclaimer.
  12.   * Redistributions in binary form must reproduce the above
  13.   copyright notice, this list of conditions and the following
  14.   disclaimer in the documentation and/or other materials provided
  15.   with the distribution.
  17.   * Neither the name of the The Musepack Development Team nor the
  18.   names of its contributors may be used to endorse or promote
  19.   products derived from this software without specific prior
  20.   written permission.
  22.   THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
  23.   "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
  24.   LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
  25.   A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
  26.   OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
  27.   SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
  28.   LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
  29.   DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
  30.   THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
  31.   (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
  32.   OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
  33. */
  35. /// file synth_filter.c
  36. /// Synthesis functions.
  37. /// todo document me
  39. #include "musepack/musepack.h"
  40. #include "musepack/internal.h"
  42. //typedef mpc_int32_t ptrdiff_t;
  44. /* C O N S T A N T S */
  45. #undef _
  47. #define MPC_FIXED_POINT_SYNTH_FIX 2
  49. #ifdef MPC_FIXED_POINT
  50. #define _(value)  MPC_MAKE_FRACT_CONST((double)value/(double)(0x40000))
  51. #else
  52. #define _(value)  MAKE_MPC_SAMPLE((double)value/(double)(0x10000))
  53. #endif
  56. static const MPC_SAMPLE_FORMAT  Di_opt [32] [16] = {
  57.     { _(  0), _( -29), _( 213), _( -459), _( 2037), _(-5153), _(  6574), _(-37489), _(75038), _(37489), _(6574), _( 5153), _(2037), _( 459), _(213), _(29) },
  58.     { _( -1), _( -31), _( 218), _( -519), _( 2000), _(-5517), _(  5959), _(-39336), _(74992), _(35640), _(7134), _( 4788), _(2063), _( 401), _(208), _(26) },
  59.     { _( -1), _( -35), _( 222), _( -581), _( 1952), _(-5879), _(  5288), _(-41176), _(74856), _(33791), _(7640), _( 4425), _(2080), _( 347), _(202), _(24) },
  60.     { _( -1), _( -38), _( 225), _( -645), _( 1893), _(-6237), _(  4561), _(-43006), _(74630), _(31947), _(8092), _( 4063), _(2087), _( 294), _(196), _(21) },
  61.     { _( -1), _( -41), _( 227), _( -711), _( 1822), _(-6589), _(  3776), _(-44821), _(74313), _(30112), _(8492), _( 3705), _(2085), _( 244), _(190), _(19) },
  62.     { _( -1), _( -45), _( 228), _( -779), _( 1739), _(-6935), _(  2935), _(-46617), _(73908), _(28289), _(8840), _( 3351), _(2075), _( 197), _(183), _(17) },
  63.     { _( -1), _( -49), _( 228), _( -848), _( 1644), _(-7271), _(  2037), _(-48390), _(73415), _(26482), _(9139), _( 3004), _(2057), _( 153), _(176), _(16) },
  64.     { _( -2), _( -53), _( 227), _( -919), _( 1535), _(-7597), _(  1082), _(-50137), _(72835), _(24694), _(9389), _( 2663), _(2032), _( 111), _(169), _(14) },
  65.     { _( -2), _( -58), _( 224), _( -991), _( 1414), _(-7910), _(    70), _(-51853), _(72169), _(22929), _(9592), _( 2330), _(2001), _(  72), _(161), _(13) },
  66.     { _( -2), _( -63), _( 221), _(-1064), _( 1280), _(-8209), _(  -998), _(-53534), _(71420), _(21189), _(9750), _( 2006), _(1962), _(  36), _(154), _(11) },
  67.     { _( -2), _( -68), _( 215), _(-1137), _( 1131), _(-8491), _( -2122), _(-55178), _(70590), _(19478), _(9863), _( 1692), _(1919), _(   2), _(147), _(10) },
  68.     { _( -3), _( -73), _( 208), _(-1210), _(  970), _(-8755), _( -3300), _(-56778), _(69679), _(17799), _(9935), _( 1388), _(1870), _( -29), _(139), _( 9) },
  69.     { _( -3), _( -79), _( 200), _(-1283), _(  794), _(-8998), _( -4533), _(-58333), _(68692), _(16155), _(9966), _( 1095), _(1817), _( -57), _(132), _( 8) },
  70.     { _( -4), _( -85), _( 189), _(-1356), _(  605), _(-9219), _( -5818), _(-59838), _(67629), _(14548), _(9959), _(  814), _(1759), _( -83), _(125), _( 7) },
  71.     { _( -4), _( -91), _( 177), _(-1428), _(  402), _(-9416), _( -7154), _(-61289), _(66494), _(12980), _(9916), _(  545), _(1698), _(-106), _(117), _( 7) },
  72.     { _( -5), _( -97), _( 163), _(-1498), _(  185), _(-9585), _( -8540), _(-62684), _(65290), _(11455), _(9838), _(  288), _(1634), _(-127), _(111), _( 6) },
  73.     { _( -5), _(-104), _( 146), _(-1567), _(  -45), _(-9727), _( -9975), _(-64019), _(64019), _( 9975), _(9727), _(   45), _(1567), _(-146), _(104), _( 5) },
  74.     { _( -6), _(-111), _( 127), _(-1634), _( -288), _(-9838), _(-11455), _(-65290), _(62684), _( 8540), _(9585), _( -185), _(1498), _(-163), _( 97), _( 5) },
  75.     { _( -7), _(-117), _( 106), _(-1698), _( -545), _(-9916), _(-12980), _(-66494), _(61289), _( 7154), _(9416), _( -402), _(1428), _(-177), _( 91), _( 4) },
  76.     { _( -7), _(-125), _(  83), _(-1759), _( -814), _(-9959), _(-14548), _(-67629), _(59838), _( 5818), _(9219), _( -605), _(1356), _(-189), _( 85), _( 4) },
  77.     { _( -8), _(-132), _(  57), _(-1817), _(-1095), _(-9966), _(-16155), _(-68692), _(58333), _( 4533), _(8998), _( -794), _(1283), _(-200), _( 79), _( 3) },
  78.     { _( -9), _(-139), _(  29), _(-1870), _(-1388), _(-9935), _(-17799), _(-69679), _(56778), _( 3300), _(8755), _( -970), _(1210), _(-208), _( 73), _( 3) },
  79.     { _(-10), _(-147), _(  -2), _(-1919), _(-1692), _(-9863), _(-19478), _(-70590), _(55178), _( 2122), _(8491), _(-1131), _(1137), _(-215), _( 68), _( 2) },
  80.     { _(-11), _(-154), _( -36), _(-1962), _(-2006), _(-9750), _(-21189), _(-71420), _(53534), _(  998), _(8209), _(-1280), _(1064), _(-221), _( 63), _( 2) },
  81.     { _(-13), _(-161), _( -72), _(-2001), _(-2330), _(-9592), _(-22929), _(-72169), _(51853), _(  -70), _(7910), _(-1414), _( 991), _(-224), _( 58), _( 2) },
  82.     { _(-14), _(-169), _(-111), _(-2032), _(-2663), _(-9389), _(-24694), _(-72835), _(50137), _(-1082), _(7597), _(-1535), _( 919), _(-227), _( 53), _( 2) },
  83.     { _(-16), _(-176), _(-153), _(-2057), _(-3004), _(-9139), _(-26482), _(-73415), _(48390), _(-2037), _(7271), _(-1644), _( 848), _(-228), _( 49), _( 1) },
  84.     { _(-17), _(-183), _(-197), _(-2075), _(-3351), _(-8840), _(-28289), _(-73908), _(46617), _(-2935), _(6935), _(-1739), _( 779), _(-228), _( 45), _( 1) },
  85.     { _(-19), _(-190), _(-244), _(-2085), _(-3705), _(-8492), _(-30112), _(-74313), _(44821), _(-3776), _(6589), _(-1822), _( 711), _(-227), _( 41), _( 1) },
  86.     { _(-21), _(-196), _(-294), _(-2087), _(-4063), _(-8092), _(-31947), _(-74630), _(43006), _(-4561), _(6237), _(-1893), _( 645), _(-225), _( 38), _( 1) },
  87.     { _(-24), _(-202), _(-347), _(-2080), _(-4425), _(-7640), _(-33791), _(-74856), _(41176), _(-5288), _(5879), _(-1952), _( 581), _(-222), _( 35), _( 1) },
  88.     { _(-26), _(-208), _(-401), _(-2063), _(-4788), _(-7134), _(-35640), _(-74992), _(39336), _(-5959), _(5517), _(-2000), _( 519), _(-218), _( 31), _( 1) }
  89. };
  91. #undef  _
  93. static void Calculate_New_V ( const MPC_SAMPLE_FORMAT * Sample, MPC_SAMPLE_FORMAT * V )
  94. {
  95.     // Calculating new V-buffer values for left channel
  96.     // calculate new V-values (ISO-11172-3, p. 39)
  97.     // based upon fast-MDCT algorithm by Byeong Gi Lee
  98.     /*static*/ MPC_SAMPLE_FORMAT A00, A01, A02, A03, A04, A05, A06, A07, A08, A09, A10, A11, A12, A13, A14, A15;
  99.     /*static*/ MPC_SAMPLE_FORMAT B00, B01, B02, B03, B04, B05, B06, B07, B08, B09, B10, B11, B12, B13, B14, B15;
  100.     MPC_SAMPLE_FORMAT tmp;
  102.     A00 = Sample[ 0] + Sample[31];
  103.     A01 = Sample[ 1] + Sample[30];
  104.     A02 = Sample[ 2] + Sample[29];
  105.     A03 = Sample[ 3] + Sample[28];
  106.     A04 = Sample[ 4] + Sample[27];
  107.     A05 = Sample[ 5] + Sample[26];
  108.     A06 = Sample[ 6] + Sample[25];
  109.     A07 = Sample[ 7] + Sample[24];
  110.     A08 = Sample[ 8] + Sample[23];
  111.     A09 = Sample[ 9] + Sample[22];
  112.     A10 = Sample[10] + Sample[21];
  113.     A11 = Sample[11] + Sample[20];
  114.     A12 = Sample[12] + Sample[19];
  115.     A13 = Sample[13] + Sample[18];
  116.     A14 = Sample[14] + Sample[17];
  117.     A15 = Sample[15] + Sample[16];
  119.     B00 = A00 + A15;
  120.     B01 = A01 + A14;
  121.     B02 = A02 + A13;
  122.     B03 = A03 + A12;
  123.     B04 = A04 + A11;
  124.     B05 = A05 + A10;
  125.     B06 = A06 + A09;
  126.     B07 = A07 + A08;;
  127.     B08 = MPC_SCALE_CONST((A00 - A15) , 0.5024192929f , 31);
  128.     B09 = MPC_SCALE_CONST((A01 - A14) , 0.5224986076f , 31);
  129.     B10 = MPC_SCALE_CONST((A02 - A13) , 0.5669440627f , 31);
  130.     B11 = MPC_SCALE_CONST((A03 - A12) , 0.6468217969f , 31);
  131.     B12 = MPC_SCALE_CONST((A04 - A11) , 0.7881546021f , 31);
  132.     B13 = MPC_SCALE_CONST((A05 - A10) , 1.0606776476f , 30);
  133.     B14 = MPC_SCALE_CONST((A06 - A09) , 1.7224471569f , 30);
  134.     B15 = MPC_SCALE_CONST((A07 - A08) , 5.1011486053f , 28);
  136.     A00 =  B00 + B07;
  137.     A01 =  B01 + B06;
  138.     A02 =  B02 + B05;
  139.     A03 =  B03 + B04;
  140.     A04 = MPC_SCALE_CONST((B00 - B07) , 0.5097956061f , 31);
  141.     A05 = MPC_SCALE_CONST((B01 - B06) , 0.6013448834f , 31);
  142.     A06 = MPC_SCALE_CONST((B02 - B05) , 0.8999761939f , 31);
  143.     A07 = MPC_SCALE_CONST((B03 - B04) , 2.5629155636f , 29);
  144.     A08 =  B08 + B15;
  145.     A09 =  B09 + B14;
  146.     A10 =  B10 + B13;
  147.     A11 =  B11 + B12;
  148.     A12 = MPC_SCALE_CONST((B08 - B15) , 0.5097956061f , 31);
  149.     A13 = MPC_SCALE_CONST((B09 - B14) , 0.6013448834f , 31);
  150.     A14 = MPC_SCALE_CONST((B10 - B13) , 0.8999761939f , 31);
  151.     A15 = MPC_SCALE_CONST((B11 - B12) , 2.5629155636f , 29);
  153.     B00 =  A00 + A03;
  154.     B01 =  A01 + A02;
  155.     B02 = MPC_MULTIPLY_FRACT_CONST_FIX((A00 - A03) , 0.5411961079f , 1);
  156.     B03 = MPC_MULTIPLY_FRACT_CONST_FIX((A01 - A02) , 1.3065630198f , 2);
  157.     B04 =  A04 + A07;
  158.     B05 =  A05 + A06;
  159.     B06 = MPC_MULTIPLY_FRACT_CONST_FIX((A04 - A07) , 0.5411961079f , 1);
  160.     B07 = MPC_MULTIPLY_FRACT_CONST_FIX((A05 - A06) , 1.3065630198f , 2);
  161.     B08 =  A08 + A11;
  162.     B09 =  A09 + A10;
  163.     B10 = MPC_MULTIPLY_FRACT_CONST_FIX((A08 - A11) , 0.5411961079f , 1);
  164.     B11 = MPC_MULTIPLY_FRACT_CONST_FIX((A09 - A10) , 1.3065630198f , 2);
  165.     B12 =  A12 + A15;
  166.     B13 =  A13 + A14;
  167.     B14 = MPC_MULTIPLY_FRACT_CONST_FIX((A12 - A15) , 0.5411961079f , 1);
  168.     B15 = MPC_MULTIPLY_FRACT_CONST_FIX((A13 - A14) , 1.3065630198f , 2);
  170.     A00 =  B00 + B01;
  171.     A01 = MPC_MULTIPLY_FRACT_CONST_FIX((B00 - B01) , 0.7071067691f , 1);
  172.     A02 =  B02 + B03;
  173.     A03 = MPC_MULTIPLY_FRACT_CONST_FIX((B02 - B03) , 0.7071067691f , 1);
  174.     A04 =  B04 + B05;
  175.     A05 = MPC_MULTIPLY_FRACT_CONST_FIX((B04 - B05) , 0.7071067691f , 1);
  176.     A06 =  B06 + B07;
  177.     A07 = MPC_MULTIPLY_FRACT_CONST_FIX((B06 - B07) , 0.7071067691f , 1);
  178.     A08 =  B08 + B09;
  179.     A09 = MPC_MULTIPLY_FRACT_CONST_FIX((B08 - B09) , 0.7071067691f , 1);
  180.     A10 =  B10 + B11;
  181.     A11 = MPC_MULTIPLY_FRACT_CONST_FIX((B10 - B11) , 0.7071067691f , 1);
  182.     A12 =  B12 + B13;
  183.     A13 = MPC_MULTIPLY_FRACT_CONST_FIX((B12 - B13) , 0.7071067691f , 1);
  184.     A14 =  B14 + B15;
  185.     A15 = MPC_MULTIPLY_FRACT_CONST_FIX((B14 - B15) , 0.7071067691f , 1);
  187.     V[48] = -A00;
  188.     V[ 0] =  A01;
  189.     V[40] = -A02 - (V[ 8] = A03);
  190.     V[36] = -((V[ 4] = A05 + (V[12] = A07)) + A06);
  191.     V[44] = - A04 - A06 - A07;
  192.     V[ 6] = (V[10] = A11 + (V[14] = A15)) + A13;
  193.     V[38] = (V[34] = -(V[ 2] = A09 + A13 + A15) - A14) + A09 - A10 - A11;
  194.     V[46] = (tmp = -(A12 + A14 + A15)) - A08;
  195.     V[42] = tmp - A10 - A11;
  197.     A00 = MPC_MULTIPLY_FRACT_CONST_SHR((Sample[ 0] - Sample[31]) , 0.5006030202f ,     MPC_FIXED_POINT_SYNTH_FIX);
  198.     A01 = MPC_MULTIPLY_FRACT_CONST_SHR((Sample[ 1] - Sample[30]) , 0.5054709315f ,     MPC_FIXED_POINT_SYNTH_FIX);
  199.     A02 = MPC_MULTIPLY_FRACT_CONST_SHR((Sample[ 2] - Sample[29]) , 0.5154473186f ,     MPC_FIXED_POINT_SYNTH_FIX);
  200.     A03 = MPC_MULTIPLY_FRACT_CONST_SHR((Sample[ 3] - Sample[28]) , 0.5310425758f ,     MPC_FIXED_POINT_SYNTH_FIX);
  201.     A04 = MPC_MULTIPLY_FRACT_CONST_SHR((Sample[ 4] - Sample[27]) , 0.5531039238f ,     MPC_FIXED_POINT_SYNTH_FIX);
  202.     A05 = MPC_MULTIPLY_FRACT_CONST_SHR((Sample[ 5] - Sample[26]) , 0.5829349756f ,     MPC_FIXED_POINT_SYNTH_FIX);
  203.     A06 = MPC_MULTIPLY_FRACT_CONST_SHR((Sample[ 6] - Sample[25]) , 0.6225041151f ,     MPC_FIXED_POINT_SYNTH_FIX);
  204.     A07 = MPC_MULTIPLY_FRACT_CONST_SHR((Sample[ 7] - Sample[24]) , 0.6748083234f ,     MPC_FIXED_POINT_SYNTH_FIX);
  205.     A08 = MPC_MULTIPLY_FRACT_CONST_SHR((Sample[ 8] - Sample[23]) , 0.7445362806f ,     MPC_FIXED_POINT_SYNTH_FIX);
  206.     A09 = MPC_MULTIPLY_FRACT_CONST_SHR((Sample[ 9] - Sample[22]) , 0.8393496275f ,     MPC_FIXED_POINT_SYNTH_FIX);
  207.     A10 = MPC_MULTIPLY_FRACT_CONST_SHR((Sample[10] - Sample[21]) , 0.9725682139f ,     MPC_FIXED_POINT_SYNTH_FIX);
  208. #if MPC_FIXED_POINT_SYNTH_FIX>=2
  209.     A11 = MPC_MULTIPLY_FRACT_CONST_SHR((Sample[11] - Sample[20]) , 1.1694399118f ,     MPC_FIXED_POINT_SYNTH_FIX);
  210.     A12 = MPC_MULTIPLY_FRACT_CONST_SHR((Sample[12] - Sample[19]) , 1.4841645956f ,     MPC_FIXED_POINT_SYNTH_FIX);
  211. #else
  212.     A11 = MPC_SCALE_CONST_SHR         ((Sample[11] - Sample[20]) , 1.1694399118f , 30, MPC_FIXED_POINT_SYNTH_FIX);
  213.     A12 = MPC_SCALE_CONST_SHR         ((Sample[12] - Sample[19]) , 1.4841645956f , 30, MPC_FIXED_POINT_SYNTH_FIX);
  214. #endif
  215.     A13 = MPC_SCALE_CONST_SHR         ((Sample[13] - Sample[18]) , 2.0577809811f , 29, MPC_FIXED_POINT_SYNTH_FIX);
  216.     A14 = MPC_SCALE_CONST_SHR         ((Sample[14] - Sample[17]) , 3.4076085091f , 29, MPC_FIXED_POINT_SYNTH_FIX);
  217.     A15 = MPC_SCALE_CONST_SHR         ((Sample[15] - Sample[16]) , 10.1900081635f, 27 ,MPC_FIXED_POINT_SYNTH_FIX);
  219.     B00 =  A00 + A15;
  220.     B01 =  A01 + A14;
  221.     B02 =  A02 + A13;
  222.     B03 =  A03 + A12;
  223.     B04 =  A04 + A11;
  224.     B05 =  A05 + A10;
  225.     B06 =  A06 + A09;
  226.     B07 =  A07 + A08;
  227.     B08 = MPC_SCALE_CONST((A00 - A15) , 0.5024192929f , 31);
  228.     B09 = MPC_SCALE_CONST((A01 - A14) , 0.5224986076f , 31);
  229.     B10 = MPC_SCALE_CONST((A02 - A13) , 0.5669440627f , 31);
  230.     B11 = MPC_SCALE_CONST((A03 - A12) , 0.6468217969f , 31);
  231.     B12 = MPC_SCALE_CONST((A04 - A11) , 0.7881546021f , 31);
  232.     B13 = MPC_SCALE_CONST((A05 - A10) , 1.0606776476f , 30);
  233.     B14 = MPC_SCALE_CONST((A06 - A09) , 1.7224471569f , 30);
  234.     B15 = MPC_SCALE_CONST((A07 - A08) , 5.1011486053f , 28);
  236.     A00 =  B00 + B07;
  237.     A01 =  B01 + B06;
  238.     A02 =  B02 + B05;
  239.     A03 =  B03 + B04;
  240.     A04 = MPC_SCALE_CONST((B00 - B07) , 0.5097956061f , 31);
  241.     A05 = MPC_SCALE_CONST((B01 - B06) , 0.6013448834f , 31);
  242.     A06 = MPC_SCALE_CONST((B02 - B05) , 0.8999761939f , 31);
  243.     A07 = MPC_SCALE_CONST((B03 - B04) , 2.5629155636f , 29);
  244.     A08 =  B08 + B15;
  245.     A09 =  B09 + B14;
  246.     A10 =  B10 + B13;
  247.     A11 =  B11 + B12;
  248.     A12 = MPC_SCALE_CONST((B08 - B15) , 0.5097956061f , 31);
  249.     A13 = MPC_SCALE_CONST((B09 - B14) , 0.6013448834f , 31);
  250.     A14 = MPC_SCALE_CONST((B10 - B13) , 0.8999761939f , 31);
  251.     A15 = MPC_SCALE_CONST((B11 - B12) , 2.5629155636f , 29);
  253.     B00 =  A00 + A03;
  254.     B01 =  A01 + A02;
  255.     B02 = MPC_SCALE_CONST((A00 - A03) , 0.5411961079f , 31);
  256.     B03 = MPC_SCALE_CONST((A01 - A02) , 1.3065630198f , 30);
  257.     B04 =  A04 + A07;
  258.     B05 =  A05 + A06;
  259.     B06 = MPC_SCALE_CONST((A04 - A07) , 0.5411961079f , 31);
  260.     B07 = MPC_SCALE_CONST((A05 - A06) , 1.3065630198f , 30);
  261.     B08 =  A08 + A11;
  262.     B09 =  A09 + A10;
  263.     B10 = MPC_SCALE_CONST((A08 - A11) , 0.5411961079f , 31);
  264.     B11 = MPC_SCALE_CONST((A09 - A10) , 1.3065630198f , 30);
  265.     B12 =  A12 + A15;
  266.     B13 =  A13 + A14;
  267.     B14 = MPC_SCALE_CONST((A12 - A15) , 0.5411961079f , 31);
  268.     B15 = MPC_SCALE_CONST((A13 - A14) , 1.3065630198f , 30);
  270.     A00 = MPC_SHL(B00 + B01, MPC_FIXED_POINT_SYNTH_FIX);
  271.     A01 = MPC_SCALE_CONST_SHL((B00 - B01) , 0.7071067691f , 31, MPC_FIXED_POINT_SYNTH_FIX);
  272.     A02 = MPC_SHL(B02 + B03, MPC_FIXED_POINT_SYNTH_FIX);
  273.     A03 = MPC_SCALE_CONST_SHL((B02 - B03) , 0.7071067691f , 31, MPC_FIXED_POINT_SYNTH_FIX);
  274.     A04 = MPC_SHL(B04 + B05, MPC_FIXED_POINT_SYNTH_FIX);
  275.     A05 = MPC_SCALE_CONST_SHL((B04 - B05) , 0.7071067691f , 31, MPC_FIXED_POINT_SYNTH_FIX);
  276.     A06 = MPC_SHL(B06 + B07, MPC_FIXED_POINT_SYNTH_FIX);
  277.     A07 = MPC_SCALE_CONST_SHL((B06 - B07) , 0.7071067691f , 31, MPC_FIXED_POINT_SYNTH_FIX);
  278.     A08 = MPC_SHL(B08 + B09, MPC_FIXED_POINT_SYNTH_FIX);
  279.     A09 = MPC_SCALE_CONST_SHL((B08 - B09) , 0.7071067691f , 31, MPC_FIXED_POINT_SYNTH_FIX);
  280.     A10 = MPC_SHL(B10 + B11, MPC_FIXED_POINT_SYNTH_FIX);
  281.     A11 = MPC_SCALE_CONST_SHL((B10 - B11) , 0.7071067691f , 31, MPC_FIXED_POINT_SYNTH_FIX);
  282.     A12 = MPC_SHL(B12 + B13, MPC_FIXED_POINT_SYNTH_FIX);
  283.     A13 = MPC_SCALE_CONST_SHL((B12 - B13) , 0.7071067691f , 31, MPC_FIXED_POINT_SYNTH_FIX);
  284.     A14 = MPC_SHL(B14 + B15, MPC_FIXED_POINT_SYNTH_FIX);
  285.     A15 = MPC_SCALE_CONST_SHL((B14 - B15) , 0.7071067691f , 31, MPC_FIXED_POINT_SYNTH_FIX);
  287.     // mehrfach verwendete Ausdr點ke: A04+A06+A07, A09+A13+A15
  288.     V[ 5] = (V[11] = (V[13] = A07 + (V[15] = A15)) + A11) + A05 + A13;
  289.     V[ 7] = (V[ 9] = A03 + A11 + A15) + A13;
  290.     V[33] = -(V[ 1] = A01 + A09 + A13 + A15) - A14;
  291.     V[35] = -(V[ 3] = A05 + A07 + A09 + A13 + A15) - A06 - A14;
  292.     V[37] = (tmp = -(A10 + A11 + A13 + A14 + A15)) - A05 - A06 - A07;
  293.     V[39] = tmp - A02 - A03;                      // abh鋘gig vom Befehl dr黚er
  294.     V[41] = (tmp += A13 - A12) - A02 - A03;       // abh鋘gig vom Befehl 2 dr黚er
  295.     V[43] = tmp - A04 - A06 - A07;                // abh鋘gig von Befehlen 1 und 3 dr黚er
  296.     V[47] = (tmp = -(A08 + A12 + A14 + A15)) - A00;
  297.     V[45] = tmp - A04 - A06 - A07;                // abh鋘gig vom Befehl dr黚er
  299.     V[32] = -V[ 0];
  300.     V[31] = -V[ 1];
  301.     V[30] = -V[ 2];
  302.     V[29] = -V[ 3];
  303.     V[28] = -V[ 4];
  304.     V[27] = -V[ 5];
  305.     V[26] = -V[ 6];
  306.     V[25] = -V[ 7];
  307.     V[24] = -V[ 8];
  308.     V[23] = -V[ 9];
  309.     V[22] = -V[10];
  310.     V[21] = -V[11];
  311.     V[20] = -V[12];
  312.     V[19] = -V[13];
  313.     V[18] = -V[14];
  314.     V[17] = -V[15];
  316.     V[63] =  V[33];
  317.     V[62] =  V[34];
  318.     V[61] =  V[35];
  319.     V[60] =  V[36];
  320.     V[59] =  V[37];
  321.     V[58] =  V[38];
  322.     V[57] =  V[39];
  323.     V[56] =  V[40];
  324.     V[55] =  V[41];
  325.     V[54] =  V[42];
  326.     V[53] =  V[43];
  327.     V[52] =  V[44];
  328.     V[51] =  V[45];
  329.     V[50] =  V[46];
  330.     V[49] =  V[47];
  331. }
  333. static void Synthese_Filter_float_internal(MPC_SAMPLE_FORMAT * OutData,MPC_SAMPLE_FORMAT * V,const MPC_SAMPLE_FORMAT * Y)
  334. {
  335.     mpc_uint32_t n;
  336.     for ( n = 0; n < 36; n++, Y += 32 ) {
  337.         V -= 64;
  338.         Calculate_New_V ( Y, V );
  339.         {
  340.             MPC_SAMPLE_FORMAT * Data = OutData;
  341.             const MPC_SAMPLE_FORMAT *  D = (const MPC_SAMPLE_FORMAT *) &Di_opt;
  342.             mpc_int32_t           k;
  343.             //mpc_int32_t           tmp;
  345.             
  346.             
  347.             for ( k = 0; k < 32; k++, D += 16, V++ ) {
  348.                 *Data = MPC_SHL(
  349.                     MPC_MULTIPLY_FRACT(V[  0],D[ 0]) + MPC_MULTIPLY_FRACT(V[ 96],D[ 1]) + MPC_MULTIPLY_FRACT(V[128],D[ 2]) + MPC_MULTIPLY_FRACT(V[224],D[ 3])
  350.                     + MPC_MULTIPLY_FRACT(V[256],D[ 4]) + MPC_MULTIPLY_FRACT(V[352],D[ 5]) + MPC_MULTIPLY_FRACT(V[384],D[ 6]) + MPC_MULTIPLY_FRACT(V[480],D[ 7])
  351.                     + MPC_MULTIPLY_FRACT(V[512],D[ 8]) + MPC_MULTIPLY_FRACT(V[608],D[ 9]) + MPC_MULTIPLY_FRACT(V[640],D[10]) + MPC_MULTIPLY_FRACT(V[736],D[11])
  352.                     + MPC_MULTIPLY_FRACT(V[768],D[12]) + MPC_MULTIPLY_FRACT(V[864],D[13]) + MPC_MULTIPLY_FRACT(V[896],D[14]) + MPC_MULTIPLY_FRACT(V[992],D[15])
  353.                     , 2);
  354.                 
  355.                 Data += 2;
  356.             }
  357.             V -= 32;//bleh
  358.             OutData+=64;
  359.         }
  360.     }
  361. }
  363. void
  364. mpc_decoder_synthese_filter_float(mpc_decoder *d, MPC_SAMPLE_FORMAT* OutData) 
  365. {
  366.     /********* left channel ********/
  367.     memmove(d->V_L + MPC_V_MEM, d->V_L, 960 * sizeof(MPC_SAMPLE_FORMAT) );
  369.     Synthese_Filter_float_internal(
  370.         OutData,
  371.         (MPC_SAMPLE_FORMAT *)(d->V_L + MPC_V_MEM),
  372.         (MPC_SAMPLE_FORMAT *)(d->Y_L [0]));
  374.     /******** right channel ********/
  375.     memmove(d->V_R + MPC_V_MEM, d->V_R, 960 * sizeof(MPC_SAMPLE_FORMAT) );
  377.     Synthese_Filter_float_internal(
  378.         OutData + 1,
  379.         (MPC_SAMPLE_FORMAT *)(d->V_R + MPC_V_MEM),
  380.         (MPC_SAMPLE_FORMAT *)(d->Y_R [0]));
  381. }
  383. /*******************************************/
  384. /*                                         */
  385. /*            dithered synthesis           */
  386. /*                                         */
  387. /*******************************************/
  389. static const unsigned char    Parity [256] = {  // parity
  390.     0,1,1,0,1,0,0,1,1,0,0,1,0,1,1,0,1,0,0,1,0,1,1,0,0,1,1,0,1,0,0,1,
  391.     1,0,0,1,0,1,1,0,0,1,1,0,1,0,0,1,0,1,1,0,1,0,0,1,1,0,0,1,0,1,1,0,
  392.     1,0,0,1,0,1,1,0,0,1,1,0,1,0,0,1,0,1,1,0,1,0,0,1,1,0,0,1,0,1,1,0,
  393.     0,1,1,0,1,0,0,1,1,0,0,1,0,1,1,0,1,0,0,1,0,1,1,0,0,1,1,0,1,0,0,1,
  394.     1,0,0,1,0,1,1,0,0,1,1,0,1,0,0,1,0,1,1,0,1,0,0,1,1,0,0,1,0,1,1,0,
  395.     0,1,1,0,1,0,0,1,1,0,0,1,0,1,1,0,1,0,0,1,0,1,1,0,0,1,1,0,1,0,0,1,
  396.     0,1,1,0,1,0,0,1,1,0,0,1,0,1,1,0,1,0,0,1,0,1,1,0,0,1,1,0,1,0,0,1,
  397.     1,0,0,1,0,1,1,0,0,1,1,0,1,0,0,1,0,1,1,0,1,0,0,1,1,0,0,1,0,1,1,0
  398. };
  400. /*
  401.  *  This is a simple random number generator with good quality for audio purposes.
  402.  *  It consists of two polycounters with opposite rotation direction and different
  403.  *  periods. The periods are coprime, so the total period is the product of both.
  404.  *
  405.  *     -------------------------------------------------------------------------------------------------
  406.  * +-> |31:30:29:28:27:26:25:24:23:22:21:20:19:18:17:16:15:14:13:12:11:10: 9: 8: 7: 6: 5: 4: 3: 2: 1: 0|
  407.  * |   -------------------------------------------------------------------------------------------------
  408.  * |                                                                          |  |  |  |     |        |
  409.  * |                                                                          +--+--+--+-XOR-+--------+
  410.  * |                                                                                      |
  411.  * +--------------------------------------------------------------------------------------+
  412.  *
  413.  *     -------------------------------------------------------------------------------------------------
  414.  *     |31:30:29:28:27:26:25:24:23:22:21:20:19:18:17:16:15:14:13:12:11:10: 9: 8: 7: 6: 5: 4: 3: 2: 1: 0| <-+
  415.  *     -------------------------------------------------------------------------------------------------   |
  416.  *       |  |           |  |                                                                               |
  417.  *       +--+----XOR----+--+                                                                               |
  418.  *                |                                                                                        |
  419.  *                +----------------------------------------------------------------------------------------+
  420.  *
  421.  *
  422.  *  The first has an period of 3*5*17*257*65537, the second of 7*47*73*178481,
  423.  *  which gives a period of 18.410.713.077.675.721.215. The result is the
  424.  *  XORed values of both generators.
  425.  */
  426. mpc_uint32_t
  427. random_int(mpc_decoder *d) 
  428. {
  429. #if 1
  430.     mpc_uint32_t  t1, t2, t3, t4;
  432.     t3   = t1 = d->__r1;   t4   = t2 = d->__r2;  // Parity calculation is done via table lookup, this is also available
  433.     t1  &= 0xF5;        t2 >>= 25;               // on CPUs without parity, can be implemented in C and avoid unpredictable
  434.     t1   = Parity [t1]; t2  &= 0x63;             // jumps and slow rotate through the carry flag operations.
  435.     t1 <<= 31;          t2   = Parity [t2];
  437.     return (d->__r1 = (t3 >> 1) | t1 ) ^ (d->__r2 = (t4 + t4) | t2 );
  438. #else
  439.     return (d->__r1 = (d->__r1 >> 1) | ((mpc_uint32_t)Parity [d->__r1 & 0xF5] << 31) ) ^
  440.         (d->__r2 = (d->__r2 << 1) |  (mpc_uint32_t)Parity [(d->__r2 >> 25) & 0x63] );
  441. #endif
  442. }