1 /****************************************************************************** 2 * * 3 * Copyright (C) 2018 The Android Open Source Project 4 * 5 * Licensed under the Apache License, Version 2.0 (the "License"); 6 * you may not use this file except in compliance with the License. 7 * You may obtain a copy of the License at: 8 * 9 * http://www.apache.org/licenses/LICENSE-2.0 10 * 11 * Unless required by applicable law or agreed to in writing, software 12 * distributed under the License is distributed on an "AS IS" BASIS, 13 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. 14 * See the License for the specific language governing permissions and 15 * limitations under the License. 16 * 17 ***************************************************************************** 18 * Originally developed and contributed by Ittiam Systems Pvt. Ltd, Bangalore 19 */ 20 #include "ixheaacd_sbr_common.h" 21 #include <ixheaacd_type_def.h> 22 23 #include "ixheaacd_constants.h" 24 #include <ixheaacd_basic_ops32.h> 25 #include <ixheaacd_basic_ops16.h> 26 #include <ixheaacd_basic_ops40.h> 27 #include "ixheaacd_basic_ops.h" 28 #include "ixheaacd_common_rom.h" 29 #include "ixheaacd_dsp_fft32x32s.h" 30 #include "ixheaacd_intrinsics.h" 31 32 #include <ixheaacd_basic_op.h> 33 34 VOID ixheaacd_inv_dit_fft_8pt_dec(WORD32 *y, WORD32 *real, WORD32 *imag) { 35 WORD32 a0, a1, a2, a3, a00, a10, a20, a30; 36 WORD32 vr, vi; 37 38 WORD32 x[16]; 39 40 a00 = ixheaacd_add32_sat(y[0], y[8]); 41 a0 = ixheaacd_sub32_sat(y[0], y[8]); 42 43 a20 = ixheaacd_add32_sat(y[1], y[9]); 44 a3 = ixheaacd_sub32_sat(y[1], y[9]); 45 46 a10 = ixheaacd_add32_sat(y[4], y[12]); 47 a2 = ixheaacd_sub32_sat(y[4], y[12]); 48 49 a30 = ixheaacd_add32_sat(y[5], y[13]); 50 a1 = ixheaacd_sub32_sat(y[5], y[13]); 51 52 x[0] = ixheaacd_add32_sat(a00, a10); 53 x[4] = ixheaacd_sub32_sat(a00, a10); 54 x[1] = ixheaacd_add32_sat(a20, a30); 55 x[5] = ixheaacd_sub32_sat(a20, a30); 56 57 x[2] = ixheaacd_sub32_sat(a0, a1); 58 x[6] = ixheaacd_add32_sat(a0, a1); 59 x[3] = ixheaacd_add32_sat(a3, a2); 60 x[7] = ixheaacd_sub32_sat(a3, a2); 61 62 a00 = ixheaacd_add32_sat(y[2], y[10]); 63 a0 = ixheaacd_sub32_sat(y[2], y[10]); 64 65 a20 = ixheaacd_add32_sat(y[3], y[11]); 66 a3 = ixheaacd_sub32_sat(y[3], y[11]); 67 68 a10 = ixheaacd_add32_sat(y[6], y[14]); 69 a2 = ixheaacd_sub32_sat(y[6], y[14]); 70 71 a30 = ixheaacd_add32_sat(y[7], y[15]); 72 a1 = ixheaacd_sub32_sat(y[7], y[15]); 73 74 x[8] = ixheaacd_add32_sat(a00, a10); 75 x[12] = ixheaacd_sub32_sat(a00, a10); 76 x[9] = ixheaacd_add32_sat(a20, a30); 77 x[13] = ixheaacd_sub32_sat(a20, a30); 78 79 x[10] = ixheaacd_sub32_sat(a0, a1); 80 x[14] = ixheaacd_add32_sat(a0, a1); 81 x[11] = ixheaacd_add32_sat(a3, a2); 82 x[15] = ixheaacd_sub32_sat(a3, a2); 83 84 real[0] = ixheaacd_add32_sat(x[0], x[8]); 85 imag[0] = ixheaacd_add32_sat(x[1], x[9]); 86 a00 = ixheaacd_sub32_sat(x[0], x[8]); 87 a10 = ixheaacd_sub32_sat(x[1], x[9]); 88 89 a0 = ixheaacd_sub32_sat(x[4], x[13]); 90 a1 = ixheaacd_add32_sat(x[5], x[12]); 91 92 real[4] = ixheaacd_add32_sat(x[4], x[13]); 93 imag[4] = ixheaacd_sub32_sat(x[5], x[12]); 94 95 vr = ixheaacd_mult32x16in32_shl_sat(ixheaacd_sub32_sat(x[10], x[11]), 0x5A82); 96 vi = ixheaacd_mult32x16in32_shl_sat(ixheaacd_add32_sat(x[10], x[11]), 0x5A82); 97 98 real[1] = ixheaacd_add32_sat(x[2], vr); 99 imag[1] = ixheaacd_add32_sat(x[3], vi); 100 101 a2 = ixheaacd_sub32_sat(x[2], vr); 102 a3 = ixheaacd_sub32_sat(x[3], vi); 103 104 real[2] = ixheaacd_add32_sat(a0, a2); 105 imag[2] = ixheaacd_add32_sat(a1, a3); 106 vr = ixheaacd_mult32x16in32_shl_sat(ixheaacd_add32_sat(x[14], x[15]), 0x5A82); 107 vi = ixheaacd_mult32x16in32_shl_sat(ixheaacd_sub32_sat(x[14], x[15]), 0x5A82); 108 109 a20 = ixheaacd_sub32_sat(x[6], vr); 110 a30 = ixheaacd_add32_sat(x[7], vi); 111 112 real[3] = ixheaacd_add32_sat(a00, a20); 113 imag[3] = ixheaacd_add32_sat(a10, a30); 114 115 real[5] = ixheaacd_add32_sat(x[6], vr); 116 imag[5] = ixheaacd_sub32_sat(x[7], vi); 117 } 118
