| /* ------------------------------------------------------------------ |
| * 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. |
| * ------------------------------------------------------------------- |
| */ |
| /**************************************************************************************** |
| Portions of this file are derived from the following 3GPP standard: |
| |
| 3GPP TS 26.073 |
| ANSI-C code for the Adaptive Multi-Rate (AMR) speech codec |
| Available from http://www.3gpp.org |
| |
| (C) 2004, 3GPP Organizational Partners (ARIB, ATIS, CCSA, ETSI, TTA, TTC) |
| Permission to distribute, modify and use this file under the standard license |
| terms listed above has been obtained from the copyright holder. |
| ****************************************************************************************/ |
| /* |
| Pathname: ./audio/gsm-amr/c/src/pred_lt.c |
| |
| ------------------------------------------------------------------------------ |
| REVISION HISTORY |
| |
| Description: Updated template used to PV coding template. First attempt at |
| optimizing C code. |
| |
| Description: Deleted variables listed in the Local Stores Needed/Modified |
| sections. |
| |
| Description: Updated file per comments from Phase 2/3 review. |
| |
| Description: Synchronized file with UMTS version 3.2.0. Updated coding |
| template. Removed unnecessary include files. |
| |
| Description: Fixed typecasting issue with TI C compiler. Updated copyright |
| year. |
| |
| Description: |
| (1) Removed instance of static in the const table "inter_6" |
| (2) Changed Overflow from a global to a parameter passed via a pointer. |
| (3) Made numerous small changes to bring code more in line with PV standards. |
| |
| Description: For pred_ltp() |
| 1. Eliminated unused include files. |
| 2. Replaced array addressing by pointers |
| 3. Eliminated math operations that unnecessary checked for |
| saturation |
| 4. Unrolled loops to speed up processing, use decrement loops, |
| loaded into memory filter coefficient in linear order for |
| faster execution in main loop. |
| 5. Eliminated call to round by proper initialization |
| |
| Description: Replaced "int" and/or "char" with defined types. |
| Added proper casting (Word32) to some left shifting operations |
| |
| |
| Description: Changed round function name to pv_round to avoid conflict with |
| round function in C standard library. |
| |
| Who: Date: |
| Description: |
| |
| ------------------------------------------------------------------------------ |
| */ |
| |
| /*---------------------------------------------------------------------------- |
| ; INCLUDES |
| ----------------------------------------------------------------------------*/ |
| #include "pred_lt.h" |
| #include "cnst.h" |
| |
| /*---------------------------------------------------------------------------- |
| ; MACROS |
| ; Define module specific macros here |
| ----------------------------------------------------------------------------*/ |
| |
| /*---------------------------------------------------------------------------- |
| ; DEFINES |
| ; Include all pre-processor statements here. Include conditional |
| ; compile variables also. |
| ----------------------------------------------------------------------------*/ |
| #define UP_SAMP_MAX 6 |
| #define L_INTER10 (L_INTERPOL-1) |
| #define FIR_SIZE (UP_SAMP_MAX*L_INTER10+1) |
| |
| /*---------------------------------------------------------------------------- |
| ; LOCAL FUNCTION DEFINITIONS |
| ; Function Prototype declaration |
| ----------------------------------------------------------------------------*/ |
| |
| /*---------------------------------------------------------------------------- |
| ; LOCAL STORE/BUFFER/POINTER DEFINITIONS |
| ; Variable declaration - defined here and used outside this module |
| ----------------------------------------------------------------------------*/ |
| |
| /* 1/6 resolution interpolation filter (-3 dB at 3600 Hz) */ |
| /* Note: the 1/3 resolution filter is simply a subsampled |
| * version of the 1/6 resolution filter, i.e. it uses |
| * every second coefficient: |
| * |
| * inter_3l[k] = inter_6[2*k], 0 <= k <= 3*L_INTER10 |
| */ |
| |
| const Word16 inter_6_pred_lt[FIR_SIZE] = |
| { |
| 29443, |
| 28346, 25207, 20449, 14701, 8693, 3143, |
| -1352, -4402, -5865, -5850, -4673, -2783, |
| -672, 1211, 2536, 3130, 2991, 2259, |
| 1170, 0, -1001, -1652, -1868, -1666, |
| -1147, -464, 218, 756, 1060, 1099, |
| 904, 550, 135, -245, -514, -634, |
| -602, -451, -231, 0, 191, 308, |
| 340, 296, 198, 78, -36, -120, |
| -163, -165, -132, -79, -19, 34, |
| 73, 91, 89, 70, 38, 0 |
| }; |
| |
| |
| /* |
| ------------------------------------------------------------------------------ |
| FUNCTION NAME: Pred_lt_3or6 |
| ------------------------------------------------------------------------------ |
| INPUT AND OUTPUT DEFINITIONS |
| |
| Inputs: |
| exc = buffer containing the excitation (Word16) |
| T0 = integer pitch lag (Word16) |
| frac = fraction of lag (Word16) |
| L_subfr = number of samples per subframe (Word16) |
| flag3 = flag to indicate the upsampling rate; if set, upsampling |
| rate is 3, otherwise, upsampling rate is 6 (Word16) |
| |
| pOverflow = pointer to overflow (Flag) |
| |
| Returns: |
| None |
| |
| Outputs: |
| exc buffer contains the newly formed adaptive codebook excitation |
| pOverflow -> 1 if the add operation resulted in overflow |
| |
| Global Variables Used: |
| inter_6_pred_lt = (1/6) resolution interpolation filter table (Word16) |
| |
| Local Variables Needed: |
| None |
| |
| ------------------------------------------------------------------------------ |
| FUNCTION DESCRIPTION |
| |
| This function computes the result of long term prediction with fractional |
| interpolation of resolution 1/3 or 1/6. (Interpolated past excitation). |
| |
| The past excitation signal at the given delay is interpolated at |
| the given fraction to build the adaptive codebook excitation. |
| On return exc[0..L_subfr-1] contains the interpolated signal |
| (adaptive codebook excitation). |
| |
| ------------------------------------------------------------------------------ |
| REQUIREMENTS |
| |
| None |
| |
| ------------------------------------------------------------------------------ |
| REFERENCES |
| |
| pred_lt.c, UMTS GSM AMR speech codec, R99 - Version 3.2.0, March 2, 2001 |
| |
| ------------------------------------------------------------------------------ |
| PSEUDO-CODE |
| |
| void Pred_lt_3or6 ( |
| Word16 exc[], // in/out: excitation buffer |
| Word16 T0, // input : integer pitch lag |
| Word16 frac, // input : fraction of lag |
| Word16 L_subfr, // input : subframe size |
| Word16 flag3 // input : if set, upsampling rate = 3 (6 otherwise) |
| ) |
| { |
| Word16 i, j, k; |
| Word16 *pX0, *pX1, *pX2; |
| const Word16 *pC1, *pC2; |
| Word32 s; |
| |
| pX0 = &exc[-T0]; |
| |
| frac = negate (frac); |
| if (flag3 != 0) |
| { |
| frac = shl (frac, 1); // inter_3l[k] = inter_6[2*k] -> k' = 2*k |
| } |
| |
| if (frac < 0) |
| { |
| frac = add (frac, UP_SAMP_MAX); |
| pX0--; |
| } |
| |
| for (j = 0; j < L_subfr; j++) |
| { |
| pX1 = pX0++; |
| pX2 = pX0; |
| pC1 = &inter_6[frac]; |
| pC2 = &inter_6[sub (UP_SAMP_MAX, frac)]; |
| |
| s = 0; |
| for (i = 0, k = 0; i < L_INTER10; i++, k += UP_SAMP_MAX) |
| { |
| s = L_mac (s, pX1[-i], pC1[k]); |
| s = L_mac (s, pX2[i], pC2[k]); |
| } |
| |
| exc[j] = pv_round (s); |
| } |
| |
| return; |
| } |
| |
| ------------------------------------------------------------------------------ |
| RESOURCES USED [optional] |
| |
| When the code is written for a specific target processor the |
| the resources used should be documented below. |
| |
| HEAP MEMORY USED: x bytes |
| |
| STACK MEMORY USED: x bytes |
| |
| CLOCK CYCLES: (cycle count equation for this function) + (variable |
| used to represent cycle count for each subroutine |
| called) |
| where: (cycle count variable) = cycle count for [subroutine |
| name] |
| |
| ------------------------------------------------------------------------------ |
| CAUTION [optional] |
| [State any special notes, constraints or cautions for users of this function] |
| |
| ------------------------------------------------------------------------------ |
| */ |
| |
| void Pred_lt_3or6( |
| Word16 exc[], /* in/out: excitation buffer */ |
| Word16 T0, /* input : integer pitch lag */ |
| Word16 frac, /* input : fraction of lag */ |
| Word16 L_subfr, /* input : subframe size */ |
| Word16 flag3, /* input : if set, upsampling rate = 3 (6 otherwise) */ |
| Flag *pOverflow /* output: if set, overflow occurred in this function */ |
| ) |
| { |
| Word16 i; |
| Word16 j; |
| Word16 k; |
| |
| Word16 *pX0; |
| Word16 *pX2; |
| Word16 *pX3; |
| Word16 *p_exc; |
| Word16 *pC1; |
| const Word16 *pC1_ref; |
| const Word16 *pC2_ref; |
| |
| Word16 Coeff_1[(L_INTER10<<1)]; |
| |
| Word32 s1; |
| Word32 s2; |
| OSCL_UNUSED_ARG(pOverflow); |
| |
| pX0 = &(exc[-T0]); |
| |
| /* frac goes between -3 and 3 */ |
| |
| frac = -frac; |
| |
| if (flag3 != 0) |
| { |
| frac <<= 1; /* inter_3l[k] = inter_6[2*k] -> k' = 2*k */ |
| } |
| |
| if (frac < 0) |
| { |
| frac += UP_SAMP_MAX; |
| pX0--; |
| } |
| |
| pC1_ref = &inter_6_pred_lt[frac]; |
| pC2_ref = &inter_6_pred_lt[UP_SAMP_MAX-frac]; |
| |
| |
| pC1 = Coeff_1; |
| |
| k = 0; |
| |
| for (i = L_INTER10 >> 1; i > 0; i--) |
| { |
| *(pC1++) = pC1_ref[k]; |
| *(pC1++) = pC2_ref[k]; |
| k += UP_SAMP_MAX; |
| *(pC1++) = pC1_ref[k]; |
| *(pC1++) = pC2_ref[k]; |
| k += UP_SAMP_MAX; |
| |
| } |
| |
| p_exc = exc; |
| |
| for (j = (L_subfr >> 1); j != 0 ; j--) |
| { |
| pX0++; |
| pX2 = pX0; |
| pX3 = pX0++; |
| |
| pC1 = Coeff_1; |
| |
| s1 = 0x00004000L; |
| s2 = 0x00004000L; |
| |
| for (i = L_INTER10 >> 1; i > 0; i--) |
| { |
| s2 += ((Word32) * (pX3--)) * *(pC1); |
| s1 += ((Word32) * (pX3)) * *(pC1++); |
| s1 += ((Word32) * (pX2++)) * *(pC1); |
| s2 += ((Word32) * (pX2)) * *(pC1++); |
| s2 += ((Word32) * (pX3--)) * *(pC1); |
| s1 += ((Word32) * (pX3)) * *(pC1++); |
| s1 += ((Word32) * (pX2++)) * *(pC1); |
| s2 += ((Word32) * (pX2)) * *(pC1++); |
| |
| } /* for (i = L_INTER10>>1; i > 0; i--) */ |
| |
| *(p_exc++) = (Word16)(s1 >> 15); |
| *(p_exc++) = (Word16)(s2 >> 15); |
| |
| } /* for (j = (L_subfr>>1); j != 0 ; j--) */ |
| |
| return; |
| } |