1dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber/* ------------------------------------------------------------------ 2dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber * Copyright (C) 1998-2009 PacketVideo 3dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber * 4dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber * Licensed under the Apache License, Version 2.0 (the "License"); 5dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber * you may not use this file except in compliance with the License. 6dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber * You may obtain a copy of the License at 7dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber * 8dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber * http://www.apache.org/licenses/LICENSE-2.0 9dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber * 10dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber * Unless required by applicable law or agreed to in writing, software 11dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber * distributed under the License is distributed on an "AS IS" BASIS, 12dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either 13dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber * express or implied. 14dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber * See the License for the specific language governing permissions 15dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber * and limitations under the License. 16dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber * ------------------------------------------------------------------- 17dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber */ 18dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber/* 19dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 20dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber Filename: idct8.c 21dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 22dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber------------------------------------------------------------------------------ 23dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber REVISION HISTORY 24dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 25dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 26dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber Who: Date: MM/DD/YYYY 27dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber Description: 28dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 29dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber------------------------------------------------------------------------------ 30dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber INPUT AND OUTPUT DEFINITIONS 31dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 32dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber Int32 x 32-bit integer input length 8 33dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 34dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 35dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber------------------------------------------------------------------------------ 36dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber FUNCTION DESCRIPTION 37dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 38dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber Implement inverse discrete cosine transform of lenght 8 39dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 40dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber------------------------------------------------------------------------------ 41dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber REQUIREMENTS 42dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 43dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 44dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber------------------------------------------------------------------------------ 45dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber REFERENCES 46dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 47dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber------------------------------------------------------------------------------ 48dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber PSEUDO-CODE 49dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 50dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber------------------------------------------------------------------------------ 51dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber*/ 52dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 53dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber#ifdef AAC_PLUS 54dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 55dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber/*---------------------------------------------------------------------------- 56dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber; INCLUDES 57dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber----------------------------------------------------------------------------*/ 58dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber#include "idct8.h" 59dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 60dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber#include "fxp_mul32.h" 61dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 62dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber/*---------------------------------------------------------------------------- 63dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber; MACROS 64dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber; Define module specific macros here 65dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber----------------------------------------------------------------------------*/ 66dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 67dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 68dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber/*---------------------------------------------------------------------------- 69dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber; DEFINES 70dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber; Include all pre-processor statements here. Include conditional 71dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber; compile variables also. 72dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber----------------------------------------------------------------------------*/ 73dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 74dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 75dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber#define R_SHIFT 29 76dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber#define Qfmt(x) (Int32)(x*((Int32)1<<R_SHIFT) + (x>=0?0.5F:-0.5F)) 77dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 78dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber#define Qfmt15(x) (Int16)(x*((Int32)1<<15) + (x>=0?0.5F:-0.5F)) 79dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 80dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 81dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber/*---------------------------------------------------------------------------- 82dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber; LOCAL FUNCTION DEFINITIONS 83dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber; Function Prototype declaration 84dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber----------------------------------------------------------------------------*/ 85dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 86dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber/*---------------------------------------------------------------------------- 87dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber; LOCAL STORE/BUFFER/POINTER DEFINITIONS 88dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber; Variable declaration - defined here and used outside this module 89dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber----------------------------------------------------------------------------*/ 90dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 91dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber/*---------------------------------------------------------------------------- 92dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber; EXTERNAL FUNCTION REFERENCES 93dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber; Declare functions defined elsewhere and referenced in this module 94dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber----------------------------------------------------------------------------*/ 95dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 96dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber/*---------------------------------------------------------------------------- 97dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber; EXTERNAL GLOBAL STORE/BUFFER/POINTER REFERENCES 98dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber; Declare variables used in this module but defined elsewhere 99dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber----------------------------------------------------------------------------*/ 100dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 101dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber/*---------------------------------------------------------------------------- 102dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber; FUNCTION CODE 103dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber----------------------------------------------------------------------------*/ 104dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 105dacaa73ae5010b66f4224d70a520945e5b653544Andreas Hubervoid idct_8(Int32 vec[]) 106dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber{ 107dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 108dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber Int32 tmp0; 109dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber Int32 tmp1; 110dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber Int32 tmp2; 111dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber Int32 tmp3; 112dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber Int32 tmp4; 113dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber Int32 tmp5; 114dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber Int32 tmp6; 115dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber Int32 tmp7; 116dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber Int32 tmp8; 117dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 118dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 119dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber tmp5 = fxp_mul32_by_16(vec[4] << 1, Qfmt15(0.70710678118655F)); 120dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 121dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber tmp1 = vec[0] + tmp5; 122dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber tmp5 = vec[0] - tmp5; 123dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 124dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber tmp3 = fxp_mul32_by_16(vec[2] << 1, Qfmt15(0.54119610014620F)); /* (1/(2*cos(2*phi)));*/ 125dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber tmp7 = fxp_mul32_Q29(vec[6], Qfmt(1.30656296487638F)); /* (1/(2*cos(6*phi)));*/ 126dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 127dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber tmp0 = fxp_mul32_by_16((tmp3 - tmp7) << 1, Qfmt15(0.70710678118655F)); /* (1/(2*cos(2*phi))); */ 128dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber tmp7 = (tmp3 + tmp7) + tmp0; 129dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 130dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber vec[0] = tmp1 + tmp7; 131dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber tmp2 = fxp_mul32_by_16(vec[1] << 1, Qfmt15(0.50979557910416F)); /* (1/(2*cos( phi)));*/ 132dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber vec[1] = tmp5 + tmp0; 133dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber vec[2] = tmp5 - tmp0; 134dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber tmp4 = fxp_mul32_by_16(vec[3] << 1, Qfmt15(0.60134488693505F)); /* (1/(2*cos(3*phi)));*/ 135dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber vec[3] = tmp1 - tmp7; 136dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 137dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber tmp6 = fxp_mul32_by_16(vec[5] << 1, Qfmt15(0.89997622313642F)); /* (1/(2*cos(5*phi)));*/ 138dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber tmp8 = fxp_mul32_Q29(vec[7], Qfmt(2.56291544774151F)); /* (1/(2*cos(7*phi)));*/ 139dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 140dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber tmp7 = tmp2 + tmp8; 141dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber tmp5 = fxp_mul32_by_16((tmp2 - tmp8) << 1, Qfmt15(0.54119610014620F)); 142dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber tmp8 = tmp4 + tmp6; 143dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber tmp6 = fxp_mul32_Q29((tmp4 - tmp6), Qfmt(1.30656296487638F)); 144dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 145dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber tmp0 = tmp7 + tmp8; 146dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber tmp2 = fxp_mul32_by_16((tmp7 - tmp8) << 1, Qfmt15(0.70710678118655F)); 147dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 148dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber tmp3 = fxp_mul32_by_16((tmp5 - tmp6) << 1, Qfmt15(0.70710678118655F)); 149dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber tmp1 = (tmp5 + tmp6) + tmp3; 150dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 151dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber tmp5 = tmp0 + tmp1; 152dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber tmp6 = tmp1 + tmp2; 153dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber tmp7 = tmp2 + tmp3; 154dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 155dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber vec[7] = vec[0] - tmp5; 156dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber vec[0] += tmp5; 157dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber vec[6] = vec[1] - tmp6; 158dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber vec[1] += tmp6; 159dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber vec[5] = vec[2] - tmp7; 160dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber vec[2] += tmp7; 161dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber vec[4] = vec[3] - tmp3; 162dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber vec[3] += tmp3; 163dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 164dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber} 165dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 166dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 167dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber#endif 168dacaa73ae5010b66f4224d70a520945e5b653544Andreas Huber 169