00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00038
00039 #include "musepack/musepack.h"
00040 #include "musepack/internal.h"
00041
00042 typedef mpc_int32_t ptrdiff_t;
00043
00044
00045 #undef _
00046
00047 #define MPC_FIXED_POINT_SYNTH_FIX 2
00048
00049 #ifdef MPC_FIXED_POINT
00050 #define _(value) MPC_MAKE_FRACT_CONST((double)value/(double)(0x40000))
00051 #else
00052 #define _(value) MAKE_MPC_SAMPLE((double)value/(double)(0x10000))
00053 #endif
00054
00055
00056 static const MPC_SAMPLE_FORMAT Di_opt [32] [16] = {
00057 { _( 0), _( -29), _( 213), _( -459), _( 2037), _(-5153), _( 6574), _(-37489), _(75038), _(37489), _(6574), _( 5153), _(2037), _( 459), _(213), _(29) },
00058 { _( -1), _( -31), _( 218), _( -519), _( 2000), _(-5517), _( 5959), _(-39336), _(74992), _(35640), _(7134), _( 4788), _(2063), _( 401), _(208), _(26) },
00059 { _( -1), _( -35), _( 222), _( -581), _( 1952), _(-5879), _( 5288), _(-41176), _(74856), _(33791), _(7640), _( 4425), _(2080), _( 347), _(202), _(24) },
00060 { _( -1), _( -38), _( 225), _( -645), _( 1893), _(-6237), _( 4561), _(-43006), _(74630), _(31947), _(8092), _( 4063), _(2087), _( 294), _(196), _(21) },
00061 { _( -1), _( -41), _( 227), _( -711), _( 1822), _(-6589), _( 3776), _(-44821), _(74313), _(30112), _(8492), _( 3705), _(2085), _( 244), _(190), _(19) },
00062 { _( -1), _( -45), _( 228), _( -779), _( 1739), _(-6935), _( 2935), _(-46617), _(73908), _(28289), _(8840), _( 3351), _(2075), _( 197), _(183), _(17) },
00063 { _( -1), _( -49), _( 228), _( -848), _( 1644), _(-7271), _( 2037), _(-48390), _(73415), _(26482), _(9139), _( 3004), _(2057), _( 153), _(176), _(16) },
00064 { _( -2), _( -53), _( 227), _( -919), _( 1535), _(-7597), _( 1082), _(-50137), _(72835), _(24694), _(9389), _( 2663), _(2032), _( 111), _(169), _(14) },
00065 { _( -2), _( -58), _( 224), _( -991), _( 1414), _(-7910), _( 70), _(-51853), _(72169), _(22929), _(9592), _( 2330), _(2001), _( 72), _(161), _(13) },
00066 { _( -2), _( -63), _( 221), _(-1064), _( 1280), _(-8209), _( -998), _(-53534), _(71420), _(21189), _(9750), _( 2006), _(1962), _( 36), _(154), _(11) },
00067 { _( -2), _( -68), _( 215), _(-1137), _( 1131), _(-8491), _( -2122), _(-55178), _(70590), _(19478), _(9863), _( 1692), _(1919), _( 2), _(147), _(10) },
00068 { _( -3), _( -73), _( 208), _(-1210), _( 970), _(-8755), _( -3300), _(-56778), _(69679), _(17799), _(9935), _( 1388), _(1870), _( -29), _(139), _( 9) },
00069 { _( -3), _( -79), _( 200), _(-1283), _( 794), _(-8998), _( -4533), _(-58333), _(68692), _(16155), _(9966), _( 1095), _(1817), _( -57), _(132), _( 8) },
00070 { _( -4), _( -85), _( 189), _(-1356), _( 605), _(-9219), _( -5818), _(-59838), _(67629), _(14548), _(9959), _( 814), _(1759), _( -83), _(125), _( 7) },
00071 { _( -4), _( -91), _( 177), _(-1428), _( 402), _(-9416), _( -7154), _(-61289), _(66494), _(12980), _(9916), _( 545), _(1698), _(-106), _(117), _( 7) },
00072 { _( -5), _( -97), _( 163), _(-1498), _( 185), _(-9585), _( -8540), _(-62684), _(65290), _(11455), _(9838), _( 288), _(1634), _(-127), _(111), _( 6) },
00073 { _( -5), _(-104), _( 146), _(-1567), _( -45), _(-9727), _( -9975), _(-64019), _(64019), _( 9975), _(9727), _( 45), _(1567), _(-146), _(104), _( 5) },
00074 { _( -6), _(-111), _( 127), _(-1634), _( -288), _(-9838), _(-11455), _(-65290), _(62684), _( 8540), _(9585), _( -185), _(1498), _(-163), _( 97), _( 5) },
00075 { _( -7), _(-117), _( 106), _(-1698), _( -545), _(-9916), _(-12980), _(-66494), _(61289), _( 7154), _(9416), _( -402), _(1428), _(-177), _( 91), _( 4) },
00076 { _( -7), _(-125), _( 83), _(-1759), _( -814), _(-9959), _(-14548), _(-67629), _(59838), _( 5818), _(9219), _( -605), _(1356), _(-189), _( 85), _( 4) },
00077 { _( -8), _(-132), _( 57), _(-1817), _(-1095), _(-9966), _(-16155), _(-68692), _(58333), _( 4533), _(8998), _( -794), _(1283), _(-200), _( 79), _( 3) },
00078 { _( -9), _(-139), _( 29), _(-1870), _(-1388), _(-9935), _(-17799), _(-69679), _(56778), _( 3300), _(8755), _( -970), _(1210), _(-208), _( 73), _( 3) },
00079 { _(-10), _(-147), _( -2), _(-1919), _(-1692), _(-9863), _(-19478), _(-70590), _(55178), _( 2122), _(8491), _(-1131), _(1137), _(-215), _( 68), _( 2) },
00080 { _(-11), _(-154), _( -36), _(-1962), _(-2006), _(-9750), _(-21189), _(-71420), _(53534), _( 998), _(8209), _(-1280), _(1064), _(-221), _( 63), _( 2) },
00081 { _(-13), _(-161), _( -72), _(-2001), _(-2330), _(-9592), _(-22929), _(-72169), _(51853), _( -70), _(7910), _(-1414), _( 991), _(-224), _( 58), _( 2) },
00082 { _(-14), _(-169), _(-111), _(-2032), _(-2663), _(-9389), _(-24694), _(-72835), _(50137), _(-1082), _(7597), _(-1535), _( 919), _(-227), _( 53), _( 2) },
00083 { _(-16), _(-176), _(-153), _(-2057), _(-3004), _(-9139), _(-26482), _(-73415), _(48390), _(-2037), _(7271), _(-1644), _( 848), _(-228), _( 49), _( 1) },
00084 { _(-17), _(-183), _(-197), _(-2075), _(-3351), _(-8840), _(-28289), _(-73908), _(46617), _(-2935), _(6935), _(-1739), _( 779), _(-228), _( 45), _( 1) },
00085 { _(-19), _(-190), _(-244), _(-2085), _(-3705), _(-8492), _(-30112), _(-74313), _(44821), _(-3776), _(6589), _(-1822), _( 711), _(-227), _( 41), _( 1) },
00086 { _(-21), _(-196), _(-294), _(-2087), _(-4063), _(-8092), _(-31947), _(-74630), _(43006), _(-4561), _(6237), _(-1893), _( 645), _(-225), _( 38), _( 1) },
00087 { _(-24), _(-202), _(-347), _(-2080), _(-4425), _(-7640), _(-33791), _(-74856), _(41176), _(-5288), _(5879), _(-1952), _( 581), _(-222), _( 35), _( 1) },
00088 { _(-26), _(-208), _(-401), _(-2063), _(-4788), _(-7134), _(-35640), _(-74992), _(39336), _(-5959), _(5517), _(-2000), _( 519), _(-218), _( 31), _( 1) }
00089 };
00090
00091 #undef _
00092
00093 static void Calculate_New_V ( const MPC_SAMPLE_FORMAT * Sample, MPC_SAMPLE_FORMAT * V )
00094 {
00095
00096
00097
00098 MPC_SAMPLE_FORMAT A00, A01, A02, A03, A04, A05, A06, A07, A08, A09, A10, A11, A12, A13, A14, A15;
00099 MPC_SAMPLE_FORMAT B00, B01, B02, B03, B04, B05, B06, B07, B08, B09, B10, B11, B12, B13, B14, B15;
00100 MPC_SAMPLE_FORMAT tmp;
00101
00102 A00 = Sample[ 0] + Sample[31];
00103 A01 = Sample[ 1] + Sample[30];
00104 A02 = Sample[ 2] + Sample[29];
00105 A03 = Sample[ 3] + Sample[28];
00106 A04 = Sample[ 4] + Sample[27];
00107 A05 = Sample[ 5] + Sample[26];
00108 A06 = Sample[ 6] + Sample[25];
00109 A07 = Sample[ 7] + Sample[24];
00110 A08 = Sample[ 8] + Sample[23];
00111 A09 = Sample[ 9] + Sample[22];
00112 A10 = Sample[10] + Sample[21];
00113 A11 = Sample[11] + Sample[20];
00114 A12 = Sample[12] + Sample[19];
00115 A13 = Sample[13] + Sample[18];
00116 A14 = Sample[14] + Sample[17];
00117 A15 = Sample[15] + Sample[16];
00118
00119 B00 = A00 + A15;
00120 B01 = A01 + A14;
00121 B02 = A02 + A13;
00122 B03 = A03 + A12;
00123 B04 = A04 + A11;
00124 B05 = A05 + A10;
00125 B06 = A06 + A09;
00126 B07 = A07 + A08;;
00127 B08 = MPC_SCALE_CONST((A00 - A15) , 0.5024192929f , 31);
00128 B09 = MPC_SCALE_CONST((A01 - A14) , 0.5224986076f , 31);
00129 B10 = MPC_SCALE_CONST((A02 - A13) , 0.5669440627f , 31);
00130 B11 = MPC_SCALE_CONST((A03 - A12) , 0.6468217969f , 31);
00131 B12 = MPC_SCALE_CONST((A04 - A11) , 0.7881546021f , 31);
00132 B13 = MPC_SCALE_CONST((A05 - A10) , 1.0606776476f , 30);
00133 B14 = MPC_SCALE_CONST((A06 - A09) , 1.7224471569f , 30);
00134 B15 = MPC_SCALE_CONST((A07 - A08) , 5.1011486053f , 28);
00135
00136 A00 = B00 + B07;
00137 A01 = B01 + B06;
00138 A02 = B02 + B05;
00139 A03 = B03 + B04;
00140 A04 = MPC_SCALE_CONST((B00 - B07) , 0.5097956061f , 31);
00141 A05 = MPC_SCALE_CONST((B01 - B06) , 0.6013448834f , 31);
00142 A06 = MPC_SCALE_CONST((B02 - B05) , 0.8999761939f , 31);
00143 A07 = MPC_SCALE_CONST((B03 - B04) , 2.5629155636f , 29);
00144 A08 = B08 + B15;
00145 A09 = B09 + B14;
00146 A10 = B10 + B13;
00147 A11 = B11 + B12;
00148 A12 = MPC_SCALE_CONST((B08 - B15) , 0.5097956061f , 31);
00149 A13 = MPC_SCALE_CONST((B09 - B14) , 0.6013448834f , 31);
00150 A14 = MPC_SCALE_CONST((B10 - B13) , 0.8999761939f , 31);
00151 A15 = MPC_SCALE_CONST((B11 - B12) , 2.5629155636f , 29);
00152
00153 B00 = A00 + A03;
00154 B01 = A01 + A02;
00155 B02 = MPC_MULTIPLY_FRACT_CONST_FIX((A00 - A03) , 0.5411961079f , 1);
00156 B03 = MPC_MULTIPLY_FRACT_CONST_FIX((A01 - A02) , 1.3065630198f , 2);
00157 B04 = A04 + A07;
00158 B05 = A05 + A06;
00159 B06 = MPC_MULTIPLY_FRACT_CONST_FIX((A04 - A07) , 0.5411961079f , 1);
00160 B07 = MPC_MULTIPLY_FRACT_CONST_FIX((A05 - A06) , 1.3065630198f , 2);
00161 B08 = A08 + A11;
00162 B09 = A09 + A10;
00163 B10 = MPC_MULTIPLY_FRACT_CONST_FIX((A08 - A11) , 0.5411961079f , 1);
00164 B11 = MPC_MULTIPLY_FRACT_CONST_FIX((A09 - A10) , 1.3065630198f , 2);
00165 B12 = A12 + A15;
00166 B13 = A13 + A14;
00167 B14 = MPC_MULTIPLY_FRACT_CONST_FIX((A12 - A15) , 0.5411961079f , 1);
00168 B15 = MPC_MULTIPLY_FRACT_CONST_FIX((A13 - A14) , 1.3065630198f , 2);
00169
00170 A00 = B00 + B01;
00171 A01 = MPC_MULTIPLY_FRACT_CONST_FIX((B00 - B01) , 0.7071067691f , 1);
00172 A02 = B02 + B03;
00173 A03 = MPC_MULTIPLY_FRACT_CONST_FIX((B02 - B03) , 0.7071067691f , 1);
00174 A04 = B04 + B05;
00175 A05 = MPC_MULTIPLY_FRACT_CONST_FIX((B04 - B05) , 0.7071067691f , 1);
00176 A06 = B06 + B07;
00177 A07 = MPC_MULTIPLY_FRACT_CONST_FIX((B06 - B07) , 0.7071067691f , 1);
00178 A08 = B08 + B09;
00179 A09 = MPC_MULTIPLY_FRACT_CONST_FIX((B08 - B09) , 0.7071067691f , 1);
00180 A10 = B10 + B11;
00181 A11 = MPC_MULTIPLY_FRACT_CONST_FIX((B10 - B11) , 0.7071067691f , 1);
00182 A12 = B12 + B13;
00183 A13 = MPC_MULTIPLY_FRACT_CONST_FIX((B12 - B13) , 0.7071067691f , 1);
00184 A14 = B14 + B15;
00185 A15 = MPC_MULTIPLY_FRACT_CONST_FIX((B14 - B15) , 0.7071067691f , 1);
00186
00187 V[48] = -A00;
00188 V[ 0] = A01;
00189 V[40] = -A02 - (V[ 8] = A03);
00190 V[36] = -((V[ 4] = A05 + (V[12] = A07)) + A06);
00191 V[44] = - A04 - A06 - A07;
00192 V[ 6] = (V[10] = A11 + (V[14] = A15)) + A13;
00193 V[38] = (V[34] = -(V[ 2] = A09 + A13 + A15) - A14) + A09 - A10 - A11;
00194 V[46] = (tmp = -(A12 + A14 + A15)) - A08;
00195 V[42] = tmp - A10 - A11;
00196
00197 A00 = MPC_MULTIPLY_FRACT_CONST_SHR((Sample[ 0] - Sample[31]) , 0.5006030202f , MPC_FIXED_POINT_SYNTH_FIX);
00198 A01 = MPC_MULTIPLY_FRACT_CONST_SHR((Sample[ 1] - Sample[30]) , 0.5054709315f , MPC_FIXED_POINT_SYNTH_FIX);
00199 A02 = MPC_MULTIPLY_FRACT_CONST_SHR((Sample[ 2] - Sample[29]) , 0.5154473186f , MPC_FIXED_POINT_SYNTH_FIX);
00200 A03 = MPC_MULTIPLY_FRACT_CONST_SHR((Sample[ 3] - Sample[28]) , 0.5310425758f , MPC_FIXED_POINT_SYNTH_FIX);
00201 A04 = MPC_MULTIPLY_FRACT_CONST_SHR((Sample[ 4] - Sample[27]) , 0.5531039238f , MPC_FIXED_POINT_SYNTH_FIX);
00202 A05 = MPC_MULTIPLY_FRACT_CONST_SHR((Sample[ 5] - Sample[26]) , 0.5829349756f , MPC_FIXED_POINT_SYNTH_FIX);
00203 A06 = MPC_MULTIPLY_FRACT_CONST_SHR((Sample[ 6] - Sample[25]) , 0.6225041151f , MPC_FIXED_POINT_SYNTH_FIX);
00204 A07 = MPC_MULTIPLY_FRACT_CONST_SHR((Sample[ 7] - Sample[24]) , 0.6748083234f , MPC_FIXED_POINT_SYNTH_FIX);
00205 A08 = MPC_MULTIPLY_FRACT_CONST_SHR((Sample[ 8] - Sample[23]) , 0.7445362806f , MPC_FIXED_POINT_SYNTH_FIX);
00206 A09 = MPC_MULTIPLY_FRACT_CONST_SHR((Sample[ 9] - Sample[22]) , 0.8393496275f , MPC_FIXED_POINT_SYNTH_FIX);
00207 A10 = MPC_MULTIPLY_FRACT_CONST_SHR((Sample[10] - Sample[21]) , 0.9725682139f , MPC_FIXED_POINT_SYNTH_FIX);
00208 #if MPC_FIXED_POINT_SYNTH_FIX>=2
00209 A11 = MPC_MULTIPLY_FRACT_CONST_SHR((Sample[11] - Sample[20]) , 1.1694399118f , MPC_FIXED_POINT_SYNTH_FIX);
00210 A12 = MPC_MULTIPLY_FRACT_CONST_SHR((Sample[12] - Sample[19]) , 1.4841645956f , MPC_FIXED_POINT_SYNTH_FIX);
00211 #else
00212 A11 = MPC_SCALE_CONST_SHR ((Sample[11] - Sample[20]) , 1.1694399118f , 30, MPC_FIXED_POINT_SYNTH_FIX);
00213 A12 = MPC_SCALE_CONST_SHR ((Sample[12] - Sample[19]) , 1.4841645956f , 30, MPC_FIXED_POINT_SYNTH_FIX);
00214 #endif
00215 A13 = MPC_SCALE_CONST_SHR ((Sample[13] - Sample[18]) , 2.0577809811f , 29, MPC_FIXED_POINT_SYNTH_FIX);
00216 A14 = MPC_SCALE_CONST_SHR ((Sample[14] - Sample[17]) , 3.4076085091f , 29, MPC_FIXED_POINT_SYNTH_FIX);
00217 A15 = MPC_SCALE_CONST_SHR ((Sample[15] - Sample[16]) , 10.1900081635f, 27 ,MPC_FIXED_POINT_SYNTH_FIX);
00218
00219 B00 = A00 + A15;
00220 B01 = A01 + A14;
00221 B02 = A02 + A13;
00222 B03 = A03 + A12;
00223 B04 = A04 + A11;
00224 B05 = A05 + A10;
00225 B06 = A06 + A09;
00226 B07 = A07 + A08;
00227 B08 = MPC_SCALE_CONST((A00 - A15) , 0.5024192929f , 31);
00228 B09 = MPC_SCALE_CONST((A01 - A14) , 0.5224986076f , 31);
00229 B10 = MPC_SCALE_CONST((A02 - A13) , 0.5669440627f , 31);
00230 B11 = MPC_SCALE_CONST((A03 - A12) , 0.6468217969f , 31);
00231 B12 = MPC_SCALE_CONST((A04 - A11) , 0.7881546021f , 31);
00232 B13 = MPC_SCALE_CONST((A05 - A10) , 1.0606776476f , 30);
00233 B14 = MPC_SCALE_CONST((A06 - A09) , 1.7224471569f , 30);
00234 B15 = MPC_SCALE_CONST((A07 - A08) , 5.1011486053f , 28);
00235
00236 A00 = B00 + B07;
00237 A01 = B01 + B06;
00238 A02 = B02 + B05;
00239 A03 = B03 + B04;
00240 A04 = MPC_SCALE_CONST((B00 - B07) , 0.5097956061f , 31);
00241 A05 = MPC_SCALE_CONST((B01 - B06) , 0.6013448834f , 31);
00242 A06 = MPC_SCALE_CONST((B02 - B05) , 0.8999761939f , 31);
00243 A07 = MPC_SCALE_CONST((B03 - B04) , 2.5629155636f , 29);
00244 A08 = B08 + B15;
00245 A09 = B09 + B14;
00246 A10 = B10 + B13;
00247 A11 = B11 + B12;
00248 A12 = MPC_SCALE_CONST((B08 - B15) , 0.5097956061f , 31);
00249 A13 = MPC_SCALE_CONST((B09 - B14) , 0.6013448834f , 31);
00250 A14 = MPC_SCALE_CONST((B10 - B13) , 0.8999761939f , 31);
00251 A15 = MPC_SCALE_CONST((B11 - B12) , 2.5629155636f , 29);
00252
00253 B00 = A00 + A03;
00254 B01 = A01 + A02;
00255 B02 = MPC_SCALE_CONST((A00 - A03) , 0.5411961079f , 31);
00256 B03 = MPC_SCALE_CONST((A01 - A02) , 1.3065630198f , 30);
00257 B04 = A04 + A07;
00258 B05 = A05 + A06;
00259 B06 = MPC_SCALE_CONST((A04 - A07) , 0.5411961079f , 31);
00260 B07 = MPC_SCALE_CONST((A05 - A06) , 1.3065630198f , 30);
00261 B08 = A08 + A11;
00262 B09 = A09 + A10;
00263 B10 = MPC_SCALE_CONST((A08 - A11) , 0.5411961079f , 31);
00264 B11 = MPC_SCALE_CONST((A09 - A10) , 1.3065630198f , 30);
00265 B12 = A12 + A15;
00266 B13 = A13 + A14;
00267 B14 = MPC_SCALE_CONST((A12 - A15) , 0.5411961079f , 31);
00268 B15 = MPC_SCALE_CONST((A13 - A14) , 1.3065630198f , 30);
00269
00270 A00 = MPC_SHL(B00 + B01, MPC_FIXED_POINT_SYNTH_FIX);
00271 A01 = MPC_SCALE_CONST_SHL((B00 - B01) , 0.7071067691f , 31, MPC_FIXED_POINT_SYNTH_FIX);
00272 A02 = MPC_SHL(B02 + B03, MPC_FIXED_POINT_SYNTH_FIX);
00273 A03 = MPC_SCALE_CONST_SHL((B02 - B03) , 0.7071067691f , 31, MPC_FIXED_POINT_SYNTH_FIX);
00274 A04 = MPC_SHL(B04 + B05, MPC_FIXED_POINT_SYNTH_FIX);
00275 A05 = MPC_SCALE_CONST_SHL((B04 - B05) , 0.7071067691f , 31, MPC_FIXED_POINT_SYNTH_FIX);
00276 A06 = MPC_SHL(B06 + B07, MPC_FIXED_POINT_SYNTH_FIX);
00277 A07 = MPC_SCALE_CONST_SHL((B06 - B07) , 0.7071067691f , 31, MPC_FIXED_POINT_SYNTH_FIX);
00278 A08 = MPC_SHL(B08 + B09, MPC_FIXED_POINT_SYNTH_FIX);
00279 A09 = MPC_SCALE_CONST_SHL((B08 - B09) , 0.7071067691f , 31, MPC_FIXED_POINT_SYNTH_FIX);
00280 A10 = MPC_SHL(B10 + B11, MPC_FIXED_POINT_SYNTH_FIX);
00281 A11 = MPC_SCALE_CONST_SHL((B10 - B11) , 0.7071067691f , 31, MPC_FIXED_POINT_SYNTH_FIX);
00282 A12 = MPC_SHL(B12 + B13, MPC_FIXED_POINT_SYNTH_FIX);
00283 A13 = MPC_SCALE_CONST_SHL((B12 - B13) , 0.7071067691f , 31, MPC_FIXED_POINT_SYNTH_FIX);
00284 A14 = MPC_SHL(B14 + B15, MPC_FIXED_POINT_SYNTH_FIX);
00285 A15 = MPC_SCALE_CONST_SHL((B14 - B15) , 0.7071067691f , 31, MPC_FIXED_POINT_SYNTH_FIX);
00286
00287
00288 V[ 5] = (V[11] = (V[13] = A07 + (V[15] = A15)) + A11) + A05 + A13;
00289 V[ 7] = (V[ 9] = A03 + A11 + A15) + A13;
00290 V[33] = -(V[ 1] = A01 + A09 + A13 + A15) - A14;
00291 V[35] = -(V[ 3] = A05 + A07 + A09 + A13 + A15) - A06 - A14;
00292 V[37] = (tmp = -(A10 + A11 + A13 + A14 + A15)) - A05 - A06 - A07;
00293 V[39] = tmp - A02 - A03;
00294 V[41] = (tmp += A13 - A12) - A02 - A03;
00295 V[43] = tmp - A04 - A06 - A07;
00296 V[47] = (tmp = -(A08 + A12 + A14 + A15)) - A00;
00297 V[45] = tmp - A04 - A06 - A07;
00298
00299 V[32] = -V[ 0];
00300 V[31] = -V[ 1];
00301 V[30] = -V[ 2];
00302 V[29] = -V[ 3];
00303 V[28] = -V[ 4];
00304 V[27] = -V[ 5];
00305 V[26] = -V[ 6];
00306 V[25] = -V[ 7];
00307 V[24] = -V[ 8];
00308 V[23] = -V[ 9];
00309 V[22] = -V[10];
00310 V[21] = -V[11];
00311 V[20] = -V[12];
00312 V[19] = -V[13];
00313 V[18] = -V[14];
00314 V[17] = -V[15];
00315
00316 V[63] = V[33];
00317 V[62] = V[34];
00318 V[61] = V[35];
00319 V[60] = V[36];
00320 V[59] = V[37];
00321 V[58] = V[38];
00322 V[57] = V[39];
00323 V[56] = V[40];
00324 V[55] = V[41];
00325 V[54] = V[42];
00326 V[53] = V[43];
00327 V[52] = V[44];
00328 V[51] = V[45];
00329 V[50] = V[46];
00330 V[49] = V[47];
00331 }
00332
00333 static void Synthese_Filter_float_internal(MPC_SAMPLE_FORMAT * OutData,MPC_SAMPLE_FORMAT * V,const MPC_SAMPLE_FORMAT * Y)
00334 {
00335 mpc_uint32_t n;
00336 for ( n = 0; n < 36; n++, Y += 32 ) {
00337 V -= 64;
00338 Calculate_New_V ( Y, V );
00339 {
00340 MPC_SAMPLE_FORMAT * Data = OutData;
00341 const MPC_SAMPLE_FORMAT * D = (const MPC_SAMPLE_FORMAT *) &Di_opt;
00342 mpc_int32_t k;
00343
00344
00345
00346
00347 for ( k = 0; k < 32; k++, D += 16, V++ ) {
00348 *Data = MPC_SHL(
00349 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])
00350 + 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])
00351 + 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])
00352 + 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])
00353 , 2);
00354
00355 Data += 2;
00356 }
00357 V -= 32;
00358 OutData+=64;
00359 }
00360 }
00361 }
00362
00363 void
00364 mpc_decoder_synthese_filter_float(mpc_decoder *d, MPC_SAMPLE_FORMAT* OutData)
00365 {
00366
00367 memmove(d->V_L + MPC_V_MEM, d->V_L, 960 * sizeof(MPC_SAMPLE_FORMAT) );
00368
00369 Synthese_Filter_float_internal(
00370 OutData,
00371 (MPC_SAMPLE_FORMAT *)(d->V_L + MPC_V_MEM),
00372 (MPC_SAMPLE_FORMAT *)(d->Y_L [0]));
00373
00374
00375 memmove(d->V_R + MPC_V_MEM, d->V_R, 960 * sizeof(MPC_SAMPLE_FORMAT) );
00376
00377 Synthese_Filter_float_internal(
00378 OutData + 1,
00379 (MPC_SAMPLE_FORMAT *)(d->V_R + MPC_V_MEM),
00380 (MPC_SAMPLE_FORMAT *)(d->Y_R [0]));
00381 }
00382
00383
00384
00385
00386
00387
00388
00389 static const unsigned char Parity [256] = {
00390 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,
00391 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,
00392 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,
00393 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,
00394 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,
00395 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,
00396 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,
00397 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
00398 };
00399
00400
00401
00402
00403
00404
00405
00406
00407
00408
00409
00410
00411
00412
00413
00414
00415
00416
00417
00418
00419
00420
00421
00422
00423
00424
00425
00426 mpc_uint32_t
00427 random_int(mpc_decoder *d)
00428 {
00429 #if 1
00430 mpc_uint32_t t1, t2, t3, t4;
00431
00432 t3 = t1 = d->__r1; t4 = t2 = d->__r2;
00433 t1 &= 0xF5; t2 >>= 25;
00434 t1 = Parity [t1]; t2 &= 0x63;
00435 t1 <<= 31; t2 = Parity [t2];
00436
00437 return (d->__r1 = (t3 >> 1) | t1 ) ^ (d->__r2 = (t4 + t4) | t2 );
00438 #else
00439 return (d->__r1 = (d->__r1 >> 1) | ((mpc_uint32_t)Parity [d->__r1 & 0xF5] << 31) ) ^
00440 (d->__r2 = (d->__r2 << 1) | (mpc_uint32_t)Parity [(d->__r2 >> 25) & 0x63] );
00441 #endif
00442 }