| /* ------------------------------------------------------------------ |
| * Copyright (C) 1998-2009 PacketVideo |
| * |
| * Licensed under the Apache License, Version 2.0 (the "License"); |
| * you may not use this file except in compliance with the License. |
| * You may obtain a copy of the License at |
| * |
| * http://www.apache.org/licenses/LICENSE-2.0 |
| * |
| * Unless required by applicable law or agreed to in writing, software |
| * distributed under the License is distributed on an "AS IS" BASIS, |
| * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either |
| * express or implied. |
| * See the License for the specific language governing permissions |
| * and limitations under the License. |
| * ------------------------------------------------------------------- |
| */ |
| /* |
| |
| Filename: pv_pow2.c |
| |
| |
| ------------------------------------------------------------------------------ |
| REVISION HISTORY |
| |
| |
| Who: Date: MM/DD/YYYY |
| Description: |
| |
| ------------------------------------------------------------------------------ |
| INPUT AND OUTPUT DEFINITIONS |
| |
| Input |
| Int32 x 32-bit integer input Q27 |
| |
| Output |
| Int32 32-bit integer in Q25 |
| |
| |
| ------------------------------------------------------------------------------ |
| FUNCTION DESCRIPTION |
| |
| Implement the power base 2 for positive numbers lesser than 5.999999 |
| ------------------------------------------------------------------------------ |
| REQUIREMENTS |
| |
| |
| ------------------------------------------------------------------------------ |
| REFERENCES |
| |
| ------------------------------------------------------------------------------ |
| PSEUDO-CODE |
| |
| ------------------------------------------------------------------------------ |
| */ |
| #ifdef AAC_PLUS |
| |
| /*---------------------------------------------------------------------------- |
| ; INCLUDES |
| ----------------------------------------------------------------------------*/ |
| #include "pv_pow2.h" |
| #include "fxp_mul32.h" |
| |
| /*---------------------------------------------------------------------------- |
| ; MACROS |
| ; Define module specific macros here |
| ----------------------------------------------------------------------------*/ |
| |
| |
| /*---------------------------------------------------------------------------- |
| ; DEFINES |
| ; Include all pre-processor statements here. Include conditional |
| ; compile variables also. |
| ----------------------------------------------------------------------------*/ |
| |
| /*---------------------------------------------------------------------------- |
| ; LOCAL FUNCTION DEFINITIONS |
| ; Function Prototype declaration |
| ----------------------------------------------------------------------------*/ |
| #define POW_2_TABLE_LENGTH 6 |
| #define POW_2_TABLE_LENGTH_m_2 (POW_2_TABLE_LENGTH - 2) |
| |
| /*---------------------------------------------------------------------------- |
| ; LOCAL STORE/BUFFER/POINTER DEFINITIONS |
| ; Variable declaration - defined here and used outside this module |
| ----------------------------------------------------------------------------*/ |
| |
| #define R_SHIFT 29 |
| #define Q_fmt(x) (Int32)(x*((Int32)1<<R_SHIFT) + (x>=0?0.5F:-0.5F)) |
| |
| #define Q27fmt(x) (Int32)(x*((Int32)1<<27) + (x>=0?0.5F:-0.5F)) |
| |
| const Int32 pow2_table[6] = |
| { |
| Q_fmt(0.00224510927441F), Q_fmt(0.00777943379416F), |
| Q_fmt(0.05737929218747F), Q_fmt(0.23918017179889F), |
| Q_fmt(0.69345251849351F), Q_fmt(0.99996347120248F) |
| }; |
| |
| |
| /*---------------------------------------------------------------------------- |
| ; EXTERNAL FUNCTION REFERENCES |
| ; Declare functions defined elsewhere and referenced in this module |
| ----------------------------------------------------------------------------*/ |
| |
| /*---------------------------------------------------------------------------- |
| ; EXTERNAL GLOBAL STORE/BUFFER/POINTER REFERENCES |
| ; Declare variables used in this module but defined elsewhere |
| ----------------------------------------------------------------------------*/ |
| |
| /*---------------------------------------------------------------------------- |
| ; FUNCTION CODE |
| ----------------------------------------------------------------------------*/ |
| |
| |
| /* |
| * z in Q27 format |
| */ |
| |
| Int32 pv_pow2(Int32 z) |
| { |
| const Int32 *pt_table = pow2_table; |
| Int32 multiplier = 0; |
| Int32 shift_factor; |
| Int32 i; |
| Int32 v_q; |
| Int32 y; |
| |
| |
| if (z > Q27fmt(1.0f)) |
| { |
| v_q = z - (z & 0xF8000000); |
| shift_factor = z >> 27; |
| } |
| else |
| { |
| v_q = z; |
| shift_factor = 0; |
| } |
| |
| if (v_q < Q27fmt(0.5f)) |
| { |
| v_q += Q27fmt(0.5f); |
| multiplier = Q_fmt(0.70710678118655F); |
| } |
| |
| v_q = v_q << 2; |
| |
| y = fxp_mul32_Q29(*(pt_table++), v_q); |
| |
| for (i = POW_2_TABLE_LENGTH_m_2; i != 0; i--) |
| { |
| y += *(pt_table++); |
| y = fxp_mul32_Q29(y, v_q); |
| } |
| y += *(pt_table++); |
| |
| if (multiplier) |
| { |
| y = fxp_mul32_Q29(y, multiplier); |
| } |
| |
| /* |
| * returns number on Q25 |
| */ |
| return (y >> (4 - shift_factor)); |
| |
| } |
| |
| #endif |