1 /* ------------------------------------------------------------------ 2 * Copyright (C) 1998-2009 PacketVideo 3 * 4 * Licensed under the Apache License, Version 2.0 (the "License"); 5 * you may not use this file except in compliance with the License. 6 * You may obtain a copy of the License at 7 * 8 * http://www.apache.org/licenses/LICENSE-2.0 9 * 10 * Unless required by applicable law or agreed to in writing, software 11 * distributed under the License is distributed on an "AS IS" BASIS, 12 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either 13 * express or implied. 14 * See the License for the specific language governing permissions 15 * and limitations under the License. 16 * ------------------------------------------------------------------- 17 */ 18 /* 19 20 Filename: ps_all_pass_filter_coeff.c 21 22 ------------------------------------------------------------------------------ 23 REVISION HISTORY 24 25 26 Who: Date: MM/DD/YYYY 27 Description: 28 29 ------------------------------------------------------------------------------ 30 INPUT AND OUTPUT DEFINITIONS 31 32 33 34 ------------------------------------------------------------------------------ 35 FUNCTION DESCRIPTION 36 37 Decorrelation is achieved by means of all-pass filtering 38 39 40 _______ ________ 41 | | _______ | | 42 ->|Hybrid | LF ---- | |->| Hybrid |--> 43 | Anal. | | | | | Synth | QMF -> L 44 ------- o----------------------->| | -------- Synth 45 QMF | s_k(n) |Stereo |--------------> 46 Anal. -------------------------->| | 47 _______ | | | | ________ 48 | | HF --o | ----------- |Process| | | 49 ->| Delay | | ->| |-------->| |->| Hybrid |--> 50 ------- | |decorrelate| d_k(n) | | | Synth | QMF -> R 51 ---->| |-------->| | -------- Synth 52 ----------- |_______|--------------> 53 54 55 ------------------------------------------------------------------------------ 56 REQUIREMENTS 57 58 59 ------------------------------------------------------------------------------ 60 REFERENCES 61 62 SC 29 Software Copyright Licencing Disclaimer: 63 64 This software module was originally developed by 65 Coding Technologies 66 67 and edited by 68 - 69 70 in the course of development of the ISO/IEC 13818-7 and ISO/IEC 14496-3 71 standards for reference purposes and its performance may not have been 72 optimized. This software module is an implementation of one or more tools as 73 specified by the ISO/IEC 13818-7 and ISO/IEC 14496-3 standards. 74 ISO/IEC gives users free license to this software module or modifications 75 thereof for use in products claiming conformance to audiovisual and 76 image-coding related ITU Recommendations and/or ISO/IEC International 77 Standards. ISO/IEC gives users the same free license to this software module or 78 modifications thereof for research purposes and further ISO/IEC standardisation. 79 Those intending to use this software module in products are advised that its 80 use may infringe existing patents. ISO/IEC have no liability for use of this 81 software module or modifications thereof. Copyright is not released for 82 products that do not conform to audiovisual and image-coding related ITU 83 Recommendations and/or ISO/IEC International Standards. 84 The original developer retains full right to modify and use the code for its 85 own purpose, assign or donate the code to a third party and to inhibit third 86 parties from using the code for products that do not conform to audiovisual and 87 image-coding related ITU Recommendations and/or ISO/IEC International Standards. 88 This copyright notice must be included in all copies or derivative works. 89 Copyright (c) ISO/IEC 2003. 90 91 ------------------------------------------------------------------------------ 92 PSEUDO-CODE 93 94 ------------------------------------------------------------------------------ 95 */ 96 97 98 /*---------------------------------------------------------------------------- 99 ; INCLUDES 100 ----------------------------------------------------------------------------*/ 101 102 #ifdef AAC_PLUS 103 104 #ifdef PARAMETRICSTEREO 105 106 #include "pv_audio_type_defs.h" 107 #include "s_ps_dec.h" 108 #include "ps_all_pass_filter_coeff.h" 109 #include "ps_all_pass_fract_delay_filter.h" 110 111 /*---------------------------------------------------------------------------- 112 ; MACROS 113 ; Define module specific macros here 114 ----------------------------------------------------------------------------*/ 115 116 /*---------------------------------------------------------------------------- 117 ; DEFINES 118 ; Include all pre-processor statements here. Include conditional 119 ; compile variables also. 120 ----------------------------------------------------------------------------*/ 121 122 123 /*---------------------------------------------------------------------------- 124 ; LOCAL FUNCTION DEFINITIONS 125 ; Function Prototype declaration 126 ----------------------------------------------------------------------------*/ 127 128 /*---------------------------------------------------------------------------- 129 ; LOCAL STORE/BUFFER/POINTER DEFINITIONS 130 ; Variable declaration - defined here and used outside this module 131 ----------------------------------------------------------------------------*/ 132 133 134 /*---------------------------------------------------------------------------- 135 ; EXTERNAL FUNCTION REFERENCES 136 ; Declare functions defined elsewhere and referenced in this module 137 ----------------------------------------------------------------------------*/ 138 139 140 const Int16 aRevLinkDecaySerCoeff[NO_ALLPASS_CHANNELS][NO_SERIAL_ALLPASS_LINKS] = 141 { 142 143 144 { Qfmt15(0.74915491616071f), Qfmt15(0.64942584030892f), Qfmt15(0.56297290849050f) }, 145 { Qfmt15(0.71658296328416f), Qfmt15(0.62118993420853f), Qfmt15(0.53849582551265f) }, 146 { Qfmt15(0.68401101040761f), Qfmt15(0.59295402810815f), Qfmt15(0.51401874253480f) }, 147 { Qfmt15(0.65143905753106f), Qfmt15(0.56471812200776f), Qfmt15(0.97908331911390f) }, /* 3 */ 148 { Qfmt15(0.61886710465450f), Qfmt15(0.53648221590737f), Qfmt15(0.93012915315822f) }, 149 { Qfmt15(0.58629515177795f), Qfmt15(0.50824630980698f), Qfmt15(0.88117498720252f) }, 150 { Qfmt15(0.55372319890140f), Qfmt15(0.48001040370660f), Qfmt15(0.83222082124682f) }, 151 { Qfmt15(0.52115124602484f), Qfmt15(0.45177449760621f), Qfmt15(0.78326665529112f) }, 152 { Qfmt15(0.48857929314829f), Qfmt15(0.42353859150582f), Qfmt15(0.73431248933542f) }, 153 { Qfmt15(0.45600734027174f), Qfmt15(0.39530268540543f), Qfmt15(0.68535832337974f) }, 154 { Qfmt15(0.42343538739519f), Qfmt15(0.36706677930504f), Qfmt15(0.63640415742404f) }, 155 { Qfmt15(0.39086343451863f), Qfmt15(0.33883087320466f), Qfmt15(0.58744999146834f) }, 156 { Qfmt15(0.35829148164208f), Qfmt15(0.31059496710427f), Qfmt15(0.53849582551265f) }, 157 { Qfmt15(0.32571952876553f), Qfmt15(0.28235906100388f), Qfmt15(0.48954165955695f) }, 158 { Qfmt15(0.29314757588898f), Qfmt15(0.25412315490349f), Qfmt15(0.44058749360126f) }, 159 { Qfmt15(0.26057562301242f), Qfmt15(0.22588724880310f), Qfmt15(0.39163332764556f) }, 160 { Qfmt15(0.22800367013587f), Qfmt15(0.19765134270272f), Qfmt15(0.34267916168986f) }, 161 { Qfmt15(0.19543171725932f), Qfmt15(0.16941543660233f), Qfmt15(0.29372499573418f) }, 162 { Qfmt15(0.16285976438276f), Qfmt15(0.14117953050194f), Qfmt15(0.24477082977848f) }, 163 { Qfmt15(0.13028781150621f), Qfmt15(0.11294362440155f), Qfmt15(0.19581666382278f) }, 164 { Qfmt15(0.09771585862966f), Qfmt15(0.08470771830116f), Qfmt15(0.14686249786708f) }, 165 { Qfmt15(0.06514390575311f), Qfmt15(0.05647181220078f), Qfmt15(0.09790833191140f) }, 166 { Qfmt15(0.03257195287655f), Qfmt15(0.02823590610039f), Qfmt15(0.04895416595570f) } 167 168 }; 169 170 171 172 173 174 const Char groupBorders[NO_IID_GROUPS + 1] = 175 { 176 4, 5, 0, 1, 2, 3, 7, 6, 8, 9, 3, 4, 177 5, 6, 7, 8, 9, 11, 14, 18, 23, 35, 64 178 }; 179 180 const Char bins2groupMap[NO_IID_GROUPS] = 181 { 182 1, 0, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19 183 }; 184 185 186 187 188 /* 189 * q_phi = 0.39 190 * 191 * cos(pi*([3:22]+0.5)*q_phi) 192 */ 193 194 195 /* 196 * sin(-pi*([3:22]+0.5)*q_phi) 197 */ 198 199 200 const Int32 aFractDelayPhaseFactor[NO_QMF_ALLPASS_CHANNELS] = 201 { 202 0xCB5474A9, 0x5BEC5914, 0x72F3C7B0, 0xF1F480C6, 0x8389E21E, 203 0xB9BA6AFC, 0x4CDB665C, 0x7A57DA5D, 0x06088024, 0x89BECF04, 204 0xA9DB5EAC, 0x3BE5711F, 0x7EB9EDF7, 0x19F582A9, 0x92DDBD1F, 205 0x9C1B5008, 0x29767919, 0x7FFC0203, 0x2D3F8843, 0x9EABACDF 206 }; 207 208 /* 209 * q(m) = { 0.43, 0.75, 0.347 } 210 * 211 * cos(pi*([3:22]+0.5)*q(m)) cos(pi*([3:22]+0.5)'*q) 212 * 213 * sin(-pi*([3:22]+0.5)*q(m)) 214 * 215 */ 216 217 218 const Int32 aaFractDelayPhaseFactorSerQmf[NO_QMF_ALLPASS_CHANNELS][3] = 219 { 220 { 0x02037FFC, 0xCF0489BE, 0x9BFB4FE0 }, 221 { 0x7D5719F5, 0xCF047642, 0x18947D9E }, 222 { 0x34AD8B57, 0x7642CF04, 0x7ABF244A }, 223 { 0x99A4B325, 0x89BECF04, 0x58EFA3F1 }, 224 { 0x9EAB5321, 0x30FC7642, 0xD77E8694 }, 225 { 0x3BE5711F, 0x30FC89BE, 0x819CEBC7 }, 226 { 0x7B77DE39, 0x89BE30FC, 0xB3A166B8 }, 227 { 0xF9F88024, 0x764230FC, 0x37C57336 }, 228 { 0x81E8E9FE, 0xCF0489BE, 0x7FF103D2 }, 229 { 0xCF047642, 0xCF047642, 0x3E8B9052 }, 230 { 0x68B9499A, 0x7642CF04, 0xB9E594E8 }, 231 { 0x5EACA9DB, 0x89BECF04, 0x80A00CA5 }, 232 { 0xC09590D1, 0x30FC7642, 0xD05276CA }, 233 { 0x85A925A3, 0x30FC89BE, 0x53486134 }, 234 { 0x0A0B7F9B, 0x89BE30FC, 0x7CB2E319 }, 235 { 0x7EB91209, 0x764230FC, 0x20078412 }, 236 { 0x2D3F8843, 0xCF0489BE, 0xA0ECAA4D }, 237 { 0x9504B9BA, 0xCF047642, 0x880D2CAE }, 238 { 0xA4145914, 0x7642CF04, 0xF0287F04 }, 239 { 0x42E16D23, 0x89BECF04, 0x694C48C7 } 240 }; 241 242 /* 243 * Fractional delay vector 244 * 245 * phi_fract(k) = exp(-j*pi*q_phi*f_center(k)) 0<= k <= SUBQMF_GROUPS 246 * 247 * q_phi = 0.39 248 * f_center(k) frequency vector 249 * 250 * 251 * f_center(k) = {0.5/4, 1.5/4, 2.5/4, 3.5/4, 252 * -1.5/4, -0.5/4, 253 * 3.5/2, 2.5/2, 4.5/2, 5.5/2}; 254 */ 255 256 257 258 const Int32 aFractDelayPhaseFactorSubQmf[SUBQMF_GROUPS] = 259 { 260 0x7E80EC79, 0x72BAC73D, 0x5C45A749, 0x3D398F97, 0x72BA38C3, 261 0x7E801387, 0xBA919478, 0x05068019, 0x895DCFF2, 0x834E1CE7, 262 }; 263 264 265 266 267 268 /* 269 * Fractional delay length matrix 270 * 271 * Q_fract(k,m) = exp(-j*pi*q(m)*f_center(k)) 272 * 273 * q(m) = { 0.43, 0.75, 0.347 } 274 * f_center(k) frequency vector 275 * 276 * 277 * f_center(k) = { 0.5/4, 1.5/4, 2.5/4, 3.5/4, 278 * -1.5/4, -0.5/4, 279 * 3.5/2, 2.5/2, 4.5/2, 5.5/2}; 280 */ 281 282 const Int32 aaFractDelayPhaseFactorSerSubQmf[SUBQMF_GROUPS][3] = 283 { 284 285 { 0x7E2EEA7D, 0x7A7DDAD8, 0x7ED0EE9D }, 286 { 0x6FEDC1E5, 0x51349D0E, 0x7574CD1E }, 287 { 0x5506A052, 0x0C8C809E, 0x636CAF62 }, 288 { 0x3085898D, 0xC3A98F1D, 0x4A0D9799 }, 289 { 0x6FED3E1B, 0x513462F2, 0x757432E2 }, 290 { 0x7E2E1583, 0x7A7D2528, 0x7ED01163 }, 291 { 0xA4C8A634, 0xB8E36A6E, 0xD5AF8732 }, 292 { 0xF0F580E3, 0x8276E707, 0x1A7382C3 }, 293 { 0x80ABF2F4, 0x471D6A6E, 0x9D2FAEA4 }, 294 { 0x9478456F, 0x7D8AE707, 0x8152EDAB } 295 }; 296 297 298 /*---------------------------------------------------------------------------- 299 ; EXTERNAL GLOBAL STORE/BUFFER/POINTER REFERENCES 300 ; Declare variables used in this module but defined elsewhere 301 ----------------------------------------------------------------------------*/ 302 303 /*---------------------------------------------------------------------------- 304 ; FUNCTION CODE 305 ----------------------------------------------------------------------------*/ 306 307 #endif 308 309 310 #endif 311 312
