14f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/* ------------------------------------------------------------------
24f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * Copyright (C) 1998-2009 PacketVideo
34f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *
44f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * Licensed under the Apache License, Version 2.0 (the "License");
54f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * you may not use this file except in compliance with the License.
64f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * You may obtain a copy of the License at
74f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *
84f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *      http://www.apache.org/licenses/LICENSE-2.0
94f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *
104f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * Unless required by applicable law or agreed to in writing, software
114f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * distributed under the License is distributed on an "AS IS" BASIS,
124f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either
134f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * express or implied.
144f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * See the License for the specific language governing permissions
154f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * and limitations under the License.
164f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * -------------------------------------------------------------------
174f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber */
184f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/****************************************************************************************
194f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas HuberPortions of this file are derived from the following 3GPP standard:
204f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
214f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    3GPP TS 26.073
224f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    ANSI-C code for the Adaptive Multi-Rate (AMR) speech codec
234f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    Available from http://www.3gpp.org
244f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
254f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber(C) 2004, 3GPP Organizational Partners (ARIB, ATIS, CCSA, ETSI, TTA, TTC)
264f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas HuberPermission to distribute, modify and use this file under the standard license
274f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huberterms listed above has been obtained from the copyright holder.
284f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber****************************************************************************************/
294f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*
304f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------
314f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
324f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
334f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
344f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Pathname: ./audio/gsm-amr/c/src/int_lsf.c
354f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
364f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber     Date: 04/20/2000
374f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
384f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------
394f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber REVISION HISTORY
404f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
414f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Put file into template and first pass at optimization.
424f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
434f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Made changes based on comments from the review meeting. Used
444f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    pointers instead of index addressing in the arrays.
454f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
464f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Added type definition to the input/output section. Fixed tabs.
474f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber              Deleted pseudo-code.
484f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
494f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Synchronized file with UMTS versin 3.2.0. Updated coding
504f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber              template. Removed unnecessary include files.
514f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
524f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Made the following changes per comments from Phase 2/3 review:
534f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber              1. Modified FOR loops to count down.
544f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber              2. Made some cosmetic changes in the Pseudo-code section.
554f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
564f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Changed to pass in overflow flag pointer to the add() routine.
574f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
584f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description:  Replaced "int" and/or "char" with OSCL defined types.
594f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
604f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description:
614f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
624f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------
634f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber*/
644f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
654f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*----------------------------------------------------------------------------
664f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; INCLUDES
674f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/
684f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#include    "int_lsf.h"
694f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#include    "typedef.h"
704f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#include    "basic_op.h"
714f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#include    "cnst.h"
724f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
734f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*----------------------------------------------------------------------------
744f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; MACROS
754f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; Define module specific macros here
764f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/
774f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
784f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
794f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*----------------------------------------------------------------------------
804f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; DEFINES
814f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; Include all pre-processor statements here. Include conditional
824f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; compile variables also.
834f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/
844f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
854f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
864f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*----------------------------------------------------------------------------
874f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; LOCAL FUNCTION DEFINITIONS
884f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; Function Prototype declaration
894f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/
904f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
914f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
924f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*----------------------------------------------------------------------------
934f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; LOCAL STORE/BUFFER/POINTER DEFINITIONS
944f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; Variable declaration - defined here and used outside this module
954f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/
964f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
974f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
984f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*
994f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------
1004f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber FUNCTION NAME: Int_lsf
1014f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------
1024f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber INPUT AND OUTPUT DEFINITIONS
1034f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
1044f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Inputs:
1054f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    lsf_old = LSF vector at the 4th SF of past frame (Word16)
1064f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    lsf_new = LSF vector at the 4th SF of present frame (Word16)
1074f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    i_subfr = Current subframe (equal to 0,40,80 or 120) (Word16)
1084f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    lsf_out = interpolated LSF parameters for current subframe (Word16)
1094f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
1104f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Outputs:
1114f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    lsf_out   = new interpolated LSF parameters for current subframe
1124f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    pOverflow = pointer of type Flag * to overflow indicator.
1134f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
1144f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Returns:
1154f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    None.
1164f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
1174f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Global Variables Used:
1184f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    None.
1194f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
1204f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Local Variables Needed:
1214f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    None.
1224f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
1234f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------
1244f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber FUNCTION DESCRIPTION
1254f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
1264f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber This function interpolates the LSFs for selected subframe.
1274f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber The 20 ms speech frame is divided into 4 subframes. The LSFs are
1284f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber interpolated at the 1st, 2nd and 3rd subframe and only forwarded
1294f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber at the 4th subframe.
1304f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
1314f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber                      |------|------|------|------|
1324f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber                         sf1    sf2    sf3    sf4
1334f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber                   F0                          F1
1344f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
1354f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber                 sf1:   3/4 F0 + 1/4 F1         sf3:   1/4 F0 + 3/4 F1
1364f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber                 sf2:   1/2 F0 + 1/2 F1         sf4:       F1
1374f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
1384f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------
1394f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber REQUIREMENTS
1404f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
1414f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None.
1424f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
1434f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------
1444f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber REFERENCES
1454f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
1464f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber int_lsf.c, UMTS GSM AMR speech codec, R99 - Version 3.2.0, March 2, 2001
1474f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
1484f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------
1494f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber PSEUDO-CODE
1504f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
1514f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Hubervoid Int_lsf(
1524f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    Word16 lsf_old[], // i : LSF vector at the 4th SF of past frame
1534f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    Word16 lsf_new[], // i : LSF vector at the 4th SF of present frame
1544f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    Word16 i_subfr,   // i : Pointer to current sf (equal to 0,40,80 or 120)
1554f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    Word16 lsf_out[]  // o : interpolated LSF parameters for current sf
1564f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber)
1574f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber{
1584f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    Word16 i;
1594f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
1604f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    if ( i_subfr == 0 )
1614f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    {
1624f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber       for (i = 0; i < M; i++) {
1634f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber          lsf_out[i] = add(sub(lsf_old[i], shr(lsf_old[i], 2)),
1644f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber                           shr(lsf_new[i], 2));
1654f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber       }
1664f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    }
1674f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    else if ( sub(i_subfr, 40) == 0 )
1684f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    {
1694f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber       for (i = 0; i < M; i++) {
1704f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber          lsf_out[i] = add(shr(lsf_old[i],1), shr(lsf_new[i], 1) );
1714f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber       }
1724f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    }
1734f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    else if ( sub(i_subfr, 80) == 0 )
1744f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    {
1754f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber       for (i = 0; i < M; i++) {
1764f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber          lsf_out[i] = add(shr(lsf_old[i], 2),
1774f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber                           sub(lsf_new[i], shr(lsf_new[i], 2)));
1784f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber       }
1794f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    }
1804f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    else if ( sub(i_subfr, 120) == 0 )
1814f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    {
1824f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber       for (i = 0; i < M; i++) {
1834f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber          lsf_out[i] = lsf_new[i];
1844f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber       }
1854f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    }
1864f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
1874f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    return;
1884f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber}
1894f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
1904f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------
1914f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber RESOURCES USED [optional]
1924f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
1934f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber When the code is written for a specific target processor the
1944f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber the resources used should be documented below.
1954f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
1964f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber HEAP MEMORY USED: x bytes
1974f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
1984f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber STACK MEMORY USED: x bytes
1994f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
2004f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber CLOCK CYCLES: (cycle count equation for this function) + (variable
2014f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber                used to represent cycle count for each subroutine
2024f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber                called)
2034f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber     where: (cycle count variable) = cycle count for [subroutine
2044f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber                                     name]
2054f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
2064f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------
2074f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber CAUTION [optional]
2084f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber [State any special notes, constraints or cautions for users of this function]
2094f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
2104f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------
2114f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber*/
2124f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
2134f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Hubervoid Int_lsf(
2144f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    Word16 lsf_old[], /* i : LSF vector at the 4th SF of past frame         */
2154f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    Word16 lsf_new[], /* i : LSF vector at the 4th SF of present frame      */
2164f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    Word16 i_subfr,   /* i : Current sf (equal to 0,40,80 or 120)           */
2174f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    Word16 lsf_out[], /* o : interpolated LSF parameters for current sf     */
2184f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    Flag  *pOverflow  /* o : flag set if overflow occurs                    */
2194f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber)
2204f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber{
2212d0ac425564ff9882ebaac5267d1a04d4af67d00Bernhard Rosenkränzer    Word16 i;
2222d0ac425564ff9882ebaac5267d1a04d4af67d00Bernhard Rosenkränzer    Word16 temp1;
2232d0ac425564ff9882ebaac5267d1a04d4af67d00Bernhard Rosenkränzer    Word16 temp2;
2244f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
2254f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    if (i_subfr == 0)
2264f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    {
2274f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber        for (i = M - 1; i >= 0; i--)
2284f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber        {
2294f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            if (*(lsf_old + i) < 0)
2304f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            {
2314f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber                temp1 = ~(~(*(lsf_old + i)) >> 2);
2324f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            }
2334f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            else
2344f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            {
2354f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber                temp1 = *(lsf_old + i) >> 2;
2364f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            }
2374f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            if (*(lsf_new + i) < 0)
2384f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            {
2394f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber                temp2 = ~(~(*(lsf_new + i)) >> 2);
2404f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            }
2414f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            else
2424f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            {
2434f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber                temp2 = *(lsf_new + i) >> 2;
2444f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            }
2454f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            *(lsf_out + i) = add((Word16)(*(lsf_old + i) - temp1),
2464f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber                                 (Word16)temp2,
2474f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber                                 pOverflow);
2484f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber        }
2494f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    }
2504f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
2514f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    else if (i_subfr == 40)
2524f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    {
2534f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber        for (i = M - 1; i >= 0; i--)
2544f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber        {
2554f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            if (*(lsf_old + i) < 0)
2564f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            {
2574f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber                temp1 = ~(~(*(lsf_old + i)) >> 1);
2584f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            }
2594f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            else
2604f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            {
2614f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber                temp1 = *(lsf_old + i) >> 1;
2624f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            }
2634f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            if (*(lsf_new + i) < 0)
2644f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            {
2654f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber                temp2 = ~(~(*(lsf_new + i)) >> 1);
2664f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            }
2674f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            else
2684f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            {
2694f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber                temp2 = *(lsf_new + i) >> 1;
2704f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            }
2714f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            *(lsf_out + i) = add(
2724f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber                                 temp1,
2734f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber                                 temp2,
2744f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber                                 pOverflow);
2754f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber        }
2764f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    }
2774f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
2784f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    else if (i_subfr == 80)
2794f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    {
2804f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber        for (i = M - 1; i >= 0; i--)
2814f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber        {
2824f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            if (*(lsf_old + i) < 0)
2834f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            {
2844f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber                temp1 = ~(~(*(lsf_old + i)) >> 2);
2854f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            }
2864f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            else
2874f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            {
2884f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber                temp1 = *(lsf_old + i) >> 2;
2894f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            }
2904f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            if (*(lsf_new + i) < 0)
2914f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            {
2924f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber                temp2 = ~(~(*(lsf_new + i)) >> 2);
2934f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            }
2944f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            else
2954f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            {
2964f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber                temp2 = *(lsf_new + i) >> 2;
2974f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            }
2984f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            *(lsf_out + i) = add((Word16)temp1,
2994f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber                                 (Word16)(*(lsf_new + i) - temp2),
3004f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber                                 pOverflow);
3014f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
3024f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber        }
3034f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    }
3044f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
3054f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    else if (i_subfr == 120)
3064f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    {
3074f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber        for (i = M - 1; i >= 0; i--)
3084f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber        {
3094f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber            *(lsf_out + i) = *(lsf_new + i);
3104f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber        }
3114f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    }
3124f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
3134f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber    return;
3144f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber}
3154f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber
316