frameworks/base/media/libstagefright/codecs/aacdec/dct16.cpp

   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  Filename: dct16.c
  20
  21 ------------------------------------------------------------------------------
  22  REVISION HISTORY
  23
  24
  25  Who:                                   Date: MM/DD/YYYY
  26  Description:
  27
  28 ------------------------------------------------------------------------------
  29  INPUT AND OUTPUT DEFINITIONS
  30
  31     Int32 x             32-bit integer input length 16
  32
  33     Int32 flag           1  forward dct16, 0 modified dct-16
  34
  35 ------------------------------------------------------------------------------
  36  FUNCTION DESCRIPTION
  37
  38     Implement dct of lenght 16
  39
  40 ------------------------------------------------------------------------------
  41  REQUIREMENTS
  42
  43
  44 ------------------------------------------------------------------------------
  45  REFERENCES
  46
  47 ------------------------------------------------------------------------------
  48  PSEUDO-CODE
  49
  50 ------------------------------------------------------------------------------
  51 */
  52
  53
  54 /*----------------------------------------------------------------------------
  55 ; INCLUDES
  56 ----------------------------------------------------------------------------*/
  57
  58 #ifdef AAC_PLUS
  59
  60 #include "dct16.h"
  61 #include "fxp_mul32.h"
  62
  63 /*----------------------------------------------------------------------------
  64 ; MACROS
  65 ; Define module specific macros here
  66 ----------------------------------------------------------------------------*/
  67
  68
  69 /*----------------------------------------------------------------------------
  70 ; DEFINES
  71 ; Include all pre-processor statements here. Include conditional
  72 ; compile variables also.
  73 ----------------------------------------------------------------------------*/
  74
  75 #define Qfmt_31(a)   (Int32)(a*0x7FFFFFFF + (a>=0?0.5F:-0.5F))
  76
  77 #define Qfmt15(x)   (Int16)(x*((Int32)1<<15) + (x>=0?0.5F:-0.5F))
  78
  79 /*----------------------------------------------------------------------------
  80 ; LOCAL FUNCTION DEFINITIONS
  81 ; Function Prototype declaration
  82 ----------------------------------------------------------------------------*/
  83
  84 /*----------------------------------------------------------------------------
  85 ; LOCAL STORE/BUFFER/POINTER DEFINITIONS
  86 ; Variable declaration - defined here and used outside this module
  87 ----------------------------------------------------------------------------*/
  88
  89 /*----------------------------------------------------------------------------
  90 ; EXTERNAL FUNCTION REFERENCES
  91 ; Declare functions defined elsewhere and referenced in this module
  92 ----------------------------------------------------------------------------*/
  93
  94 /*----------------------------------------------------------------------------
  95 ; EXTERNAL GLOBAL STORE/BUFFER/POINTER REFERENCES
  96 ; Declare variables used in this module but defined elsewhere
  97 ----------------------------------------------------------------------------*/
  98
  99 /*----------------------------------------------------------------------------
 100 ; FUNCTION CODE
 101 ----------------------------------------------------------------------------*/
 102
 103
 104 void dct_16(Int32 vec[], Int flag)
 105 {
 106     Int32 tmp0;
 107     Int32 tmp1;
 108     Int32 tmp2;
 109     Int32 tmp3;
 110     Int32 tmp4;
 111     Int32 tmp5;
 112     Int32 tmp6;
 113     Int32 tmp7;
 114     Int32 tmp_o0;
 115     Int32 tmp_o1;
 116     Int32 tmp_o2;
 117     Int32 tmp_o3;
 118     Int32 tmp_o4;
 119     Int32 tmp_o5;
 120     Int32 tmp_o6;
 121     Int32 tmp_o7;
 122     Int32 itmp_e0;
 123     Int32 itmp_e1;
 124     Int32 itmp_e2;
 125
 126     /*  split input vector */
 127
 128
 129     tmp_o0 = fxp_mul32_by_16((vec[ 0] - vec[15]), Qfmt15(0.50241928618816F));
 130     tmp0   =  vec[ 0] + vec[15];
 131
 132     tmp_o7 = fxp_mul32_Q31((vec[ 7] - vec[ 8]) << 3, Qfmt_31(0.63764357733614F));
 133     tmp7   =  vec[ 7] + vec[ 8];
 134
 135     itmp_e0 = (tmp0 + tmp7);
 136     tmp7    = fxp_mul32_by_16((tmp0 - tmp7), Qfmt15(0.50979557910416F));
 137
 138     tmp_o1 = fxp_mul32_by_16((vec[ 1] - vec[14]), Qfmt15(0.52249861493969F));
 139     tmp1   =  vec[ 1] + vec[14];
 140     tmp_o6 = fxp_mul32_by_16((vec[ 6] - vec[ 9]) << 1, Qfmt15(0.86122354911916F));
 141     tmp6   =  vec[ 6] + vec[ 9];
 142
 143     itmp_e1 = (tmp1 + tmp6);
 144     tmp6    = fxp_mul32_by_16((tmp1 - tmp6), Qfmt15(0.60134488693505F));
 145
 146     tmp_o2 = fxp_mul32_by_16((vec[ 2] - vec[13]), Qfmt15(0.56694403481636F));
 147     tmp2   =  vec[ 2] + vec[13];
 148     tmp_o5 = fxp_mul32_by_16((vec[ 5] - vec[10]) << 1, Qfmt15(0.53033884299517F));
 149     tmp5   =  vec[ 5] + vec[10];
 150
 151     itmp_e2 = (tmp2 + tmp5);
 152     tmp5    = fxp_mul32_by_16((tmp2 - tmp5), Qfmt15(0.89997622313642F));
 153
 154     tmp_o3 = fxp_mul32_by_16((vec[ 3] - vec[12]), Qfmt15(0.64682178335999F));
 155     tmp3   =  vec[ 3] + vec[12];
 156     tmp_o4 = fxp_mul32_by_16((vec[ 4] - vec[11]), Qfmt15(0.78815462345125F));
 157     tmp4   =  vec[ 4] + vec[11];
 158
 159     tmp1   = (tmp3 + tmp4);
 160     tmp4   =  fxp_mul32_Q31((tmp3 - tmp4) << 2, Qfmt_31(0.64072886193538F));
 161
 162     /*  split even part of tmp_e */
 163
 164     tmp0 = (itmp_e0 + tmp1);
 165     tmp1 = fxp_mul32_by_16((itmp_e0 - tmp1), Qfmt15(0.54119610014620F));
 166
 167
 168     tmp3 = fxp_mul32_by_16((itmp_e1 - itmp_e2) << 1, Qfmt15(0.65328148243819F));
 169     tmp2 = (itmp_e1 + itmp_e2);
 170
 171     vec[ 0]  = (tmp0 + tmp2) >> 1;
 172     vec[ 8]  = fxp_mul32_by_16((tmp0 - tmp2), Qfmt15(0.70710678118655F));
 173     vec[12]  = fxp_mul32_by_16((tmp1 - tmp3) << 1, Qfmt15(0.70710678118655F));
 174     vec[ 4]  =  tmp1 + tmp3;
 175     vec[ 4] +=  vec[12];
 176
 177     /*  split odd part of tmp_e */
 178
 179     tmp1 = fxp_mul32_by_16((tmp7 - tmp4) << 1, Qfmt15(0.54119610014620F));
 180     tmp7 += tmp4;
 181     tmp3 = fxp_mul32_Q31((tmp6 - tmp5) << 2, Qfmt_31(0.65328148243819F));
 182
 183     tmp6 += tmp5;
 184
 185     vec[10]  = fxp_mul32_by_16((tmp7 - tmp6) << 1, Qfmt15(0.70710678118655F));
 186     vec[ 2]  =  tmp7 + tmp6;
 187     vec[14]  = fxp_mul32_by_16((tmp1 - tmp3) << 1, Qfmt15(0.70710678118655F));
 188
 189     tmp1    +=  tmp3 + vec[14];
 190     vec[ 2] +=  tmp1;
 191     vec[ 6]  =  tmp1 + vec[10];
 192
 193     vec[10] += vec[14];
 194
 195
 196     // dct8;
 197
 198     tmp7 = tmp_o0 + tmp_o7;
 199     tmp_o7 = fxp_mul32_by_16((tmp_o0 - tmp_o7) << 1, Qfmt15(0.50979557910416F));
 200
 201     tmp6 = tmp_o1 + tmp_o6;
 202     tmp_o1 = fxp_mul32_by_16((tmp_o1 - tmp_o6) << 1, Qfmt15(0.60134488693505F));
 203
 204     tmp5 = tmp_o2 + tmp_o5;
 205     tmp_o5 = fxp_mul32_by_16((tmp_o2 - tmp_o5) << 1, Qfmt15(0.89997622313642F));
 206
 207     tmp4 = tmp_o3 + tmp_o4;
 208
 209     tmp_o3 = fxp_mul32_Q31((tmp_o3 - tmp_o4) << 3, Qfmt_31(0.6407288619354F));
 210
 211     if (!flag)
 212     {
 213         tmp7   = -tmp7;
 214         tmp_o7 = -tmp_o7;
 215         tmp6   = -tmp6;
 216         tmp_o1 = -tmp_o1;
 217         tmp5   = -tmp5;
 218         tmp_o5 = -tmp_o5;
 219         tmp4   = -tmp4;
 220         tmp_o3 = -tmp_o3;
 221     }
 222
 223     // even part
 224
 225     tmp1 = fxp_mul32_by_16((tmp7 - tmp4) << 1, Qfmt15(0.54119610014620F));
 226     tmp0 =  tmp7 + tmp4;
 227     tmp3 = fxp_mul32_Q31((tmp6 - tmp5) << 2, Qfmt_31(0.65328148243819F));
 228     tmp2 =  tmp6 + tmp5;
 229
 230     vec[ 9]  = fxp_mul32_Q31((tmp0 - tmp2) << 1, Qfmt_31(0.70710678118655F));
 231     vec[ 1]  =  tmp0 + tmp2;
 232     vec[13]  = fxp_mul32_Q31((tmp1 - tmp3) << 1, Qfmt_31(0.70710678118655F));
 233
 234     vec[ 5]  =  tmp1 + tmp3 + vec[13];
 235
 236     // odd part
 237
 238     tmp0 =  tmp_o7 + tmp_o3;
 239     tmp1 = fxp_mul32_by_16((tmp_o7 - tmp_o3) << 1, Qfmt15(0.54119610014620F));
 240     tmp2 =  tmp_o1 + tmp_o5;
 241     tmp3 = fxp_mul32_Q31((tmp_o1 - tmp_o5) << 2, Qfmt_31(0.65328148243819F));
 242
 243     vec[11]  = fxp_mul32_Q31((tmp0 - tmp2) << 1, Qfmt_31(0.70710678118655F));
 244     vec[ 3]  =  tmp0 + tmp2;
 245     vec[15]  = fxp_mul32_Q31((tmp1 - tmp3) << 1, Qfmt_31(0.70710678118655F));
 246     vec[ 7]  =  tmp1 + tmp3 + vec[15];
 247
 248
 249     vec[ 3] += vec[ 7];
 250     vec[ 7] += vec[11];
 251     vec[11] += vec[15];
 252
 253     vec[ 1] += vec[ 3];
 254     vec[ 3] += vec[ 5];
 255     vec[ 5] += vec[ 7];
 256     vec[ 7] += vec[ 9];
 257     vec[ 9] += vec[11];
 258     vec[11] += vec[13];
 259     vec[13] += vec[15];
 260
 261
 262 }
 263
 264
 265 #endif
 266