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.173 224f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber ANSI-C code for the Adaptive Multi-Rate - Wideband (AMR-WB) speech codec 234f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Available from http://www.3gpp.org 244f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 254f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber(C) 2007, 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 Filename: oversamp_12k8_to_16k.cpp 354f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 364f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Date: 05/08/2004 374f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 384f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 394f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber REVISION HISTORY 404f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 414f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 424f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: 434f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 444f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 454f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber INPUT AND OUTPUT DEFINITIONS 464f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 474f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber int16 signal[], input signal / output is divided by 16 484f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber int16 lg, lenght of signal 494f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber int16 mem[] in/out: memory (size=30) 504f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber int16 x[] scratch mem ( size= 60) 514f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 524f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 534f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber FUNCTION DESCRIPTION 544f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 554f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Oversamp_16k : oversampling from 12.8kHz to 16kHz. 564f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 574f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 584f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 594f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber REQUIREMENTS 604f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 614f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 624f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 634f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber REFERENCES 644f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 654f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 664f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber PSEUDO-CODE 674f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 684f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 694f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber*/ 704f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 714f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 724f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*---------------------------------------------------------------------------- 734f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; INCLUDES 744f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/ 754f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 764f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#include "pv_amr_wb_type_defs.h" 774f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#include "pvamrwbdecoder_basic_op.h" 784f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#include "pvamrwbdecoder_acelp.h" 794f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#include "pvamrwbdecoder_cnst.h" 804f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 814f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*---------------------------------------------------------------------------- 824f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; MACROS 834f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; Define module specific macros here 844f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/ 854f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 864f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 874f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*---------------------------------------------------------------------------- 884f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; DEFINES 894f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; Include all pre-processor statements here. Include conditional 904f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; compile variables also. 914f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/ 924f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 934f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#define FAC4 4 944f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#define FAC5 5 954f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#define INV_FAC5 6554 /* 1/5 in Q15 */ 964f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#define DOWN_FAC 26215 /* 4/5 in Q15 */ 974f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#define UP_FAC 20480 /* 5/4 in Q14 */ 984f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#define NB_COEF_DOWN 15 994f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#define NB_COEF_UP 12 1004f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#define N_LOOP_COEF_UP 4 1014f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1024f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*---------------------------------------------------------------------------- 1034f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; LOCAL FUNCTION DEFINITIONS 1044f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; Function Prototype declaration 1054f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/ 1064f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1074f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#ifdef __cplusplus 1084f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huberextern "C" 1094f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber{ 1104f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#endif 1114f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1124f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1134f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* Local functions */ 1144f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1154f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber void AmrWbUp_samp( 1164f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber int16 * sig_d, /* input: signal to oversampling */ 1174f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber int16 * sig_u, /* output: oversampled signal */ 1184f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber int16 L_frame /* input: length of output */ 1194f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber ); 1204f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1214f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1224f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber int16 AmrWbInterpol( /* return result of interpolation */ 1234f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber int16 * x, /* input vector */ 1244f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber const int16 * fir, /* filter coefficient */ 1254f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber int16 nb_coef /* number of coefficients */ 1264f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber ); 1274f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1284f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1294f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#ifdef __cplusplus 1304f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber} 1314f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#endif 1324f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1334f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*---------------------------------------------------------------------------- 1344f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; LOCAL STORE/BUFFER/POINTER DEFINITIONS 1354f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; Variable declaration - defined here and used outside this module 1364f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/ 1374f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1384f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1394f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/* 1/5 resolution interpolation filter (in Q14) */ 1404f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/* -1.5dB @ 6kHz, -6dB @ 6.4kHz, -10dB @ 6.6kHz, 1414f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber -20dB @ 6.9kHz, -25dB @ 7kHz, -55dB @ 8kHz */ 1424f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1434f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1444f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huberconst int16 fir_up[4][24] = 1454f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber{ 1464f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1474f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 1484f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber -1, 12, -33, 68, -119, 191, 1494f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber -291, 430, -634, 963, -1616, 3792, 1504f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 15317, -2496, 1288, -809, 542, -369, 1514f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 247, -160, 96, -52, 23, -6, 1524f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber }, 1534f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 1544f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber -4, 24, -62, 124, -213, 338, 1554f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber -510, 752, -1111, 1708, -2974, 8219, 1564f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 12368, -3432, 1881, -1204, 812, -552, 1574f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 368, -235, 139, -73, 30, -7, 1584f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber }, 1594f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 1604f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber -7, 30, -73, 139, -235, 368, 1614f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber -552, 812, -1204, 1881, -3432, 12368, 1624f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 8219, -2974, 1708, -1111, 752, -510, 1634f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 338, -213, 124, -62, 24, -4, 1644f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber }, 1654f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 1664f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber -6, 23, -52, 96, -160, 247, 1674f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber -369, 542, -809, 1288, -2496, 15317, 1684f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3792, -1616, 963, -634, 430, -291, 1694f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 191, -119, 68, -33, 12, -1, 1704f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 1714f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber}; 1724f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1734f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*---------------------------------------------------------------------------- 1744f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; EXTERNAL FUNCTION REFERENCES 1754f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; Declare functions defined elsewhere and referenced in this module 1764f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/ 1774f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1784f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*---------------------------------------------------------------------------- 1794f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; EXTERNAL GLOBAL STORE/BUFFER/POINTER REFERENCES 1804f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; Declare variables used in this module but defined elsewhere 1814f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/ 1824f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1834f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*---------------------------------------------------------------------------- 1844f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; FUNCTION CODE 1854f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/ 1864f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1874f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1884f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/* output: memory (2*NB_COEF_UP) set to zeros */ 1894f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Hubervoid oversamp_12k8_to_16k_init(int16 mem[]) 1904f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber{ 1914f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber pv_memset((void *)mem, 0, (2*NB_COEF_UP)*sizeof(*mem)); 1924f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1934f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber} 1944f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1954f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*---------------------------------------------------------------------------- 1964f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; FUNCTION CODE 1974f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/ 1984f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1994f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Hubervoid oversamp_12k8_to_16k( 2004f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber int16 sig12k8[], /* input: signal to oversampling */ 2014f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber int16 lg, /* input: length of input */ 2024f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber int16 sig16k[], /* output: oversampled signal */ 2034f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber int16 mem[], /* in/out: memory (2*NB_COEF_UP) */ 2044f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber int16 signal[] 2054f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber) 2064f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber{ 2074f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber int16 lg_up; 2084f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2094f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber pv_memcpy((void *)signal, 2104f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber (void *)mem, 2114f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber (2*NB_COEF_UP)*sizeof(*mem)); 2124f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2134f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber pv_memcpy((void *)(signal + (2*NB_COEF_UP)), 2144f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber (void *)sig12k8, 2154f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber lg*sizeof(*sig12k8)); 2164f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2174f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber lg_up = lg + (lg >> 2); /* 5/4 of lg */ 2184f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2194f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber AmrWbUp_samp(signal + NB_COEF_UP, sig16k, lg_up); 2204f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2214f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber pv_memcpy((void *)mem, 2224f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber (void *)(signal + lg), 2234f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber (2*NB_COEF_UP)*sizeof(*signal)); 2244f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2254f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber return; 2264f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber} 2274f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2284f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2294f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2304f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*---------------------------------------------------------------------------- 2314f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; FUNCTION CODE 2324f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/ 2334f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2344f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2354f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Hubervoid AmrWbUp_samp( 2364f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber int16 * sig_d, /* input: signal to oversampling */ 2374f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber int16 * sig_u, /* output: oversampled signal */ 2384f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber int16 L_frame /* input: length of output */ 2394f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber) 2404f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber{ 2414f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2424f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber int32 i; 243f27422a71112e77af77de0da2e60d4deb94a9375Martin Storsjo int16 frac, j; 2444f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber int16 * pt_sig_u = sig_u; 2454f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2464f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber frac = 1; 247f27422a71112e77af77de0da2e60d4deb94a9375Martin Storsjo for (j = 0; j < L_frame; j++) 2484f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 2494f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber i = ((int32)j * INV_FAC5) >> 13; /* integer part = pos * 1/5 */ 2504f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2514f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber frac--; 2524f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if (frac) 2534f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 2544f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *(pt_sig_u++) = AmrWbInterpol(&sig_d[i], 2554f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber fir_up[(FAC5-1) - frac], 2564f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber N_LOOP_COEF_UP); 2574f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 2584f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber else 2594f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 2604f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *(pt_sig_u++) = sig_d[i+12 - NB_COEF_UP ]; 2614f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber frac = FAC5; 2624f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 2634f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 2644f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2654f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber} 2664f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2674f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*---------------------------------------------------------------------------- 2684f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; FUNCTION CODE 2694f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/ 2704f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2714f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2724f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/* Fractional interpolation of signal at position (frac/resol) */ 2734f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2744f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2754f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huberint16 AmrWbInterpol( /* return result of interpolation */ 2764f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber int16 * x, /* input vector */ 2774f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber const int16 *fir, /* filter coefficient */ 2784f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber int16 nb_coef /* number of coefficients */ 2794f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber) 2804f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber{ 2814f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber int32 L_sum; 2824f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber const int16 *pt_fir = fir; 2834f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2844f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber int16 tmp1, tmp2, tmp3, tmp4; 2854f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber int16 *pt_x = x - nb_coef - (nb_coef << 1) + 1; 2864f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2874f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2884f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber tmp1 = *(pt_x++); 2894f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber tmp2 = *(pt_x++); 2904f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber tmp3 = *(pt_x++); 2914f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber tmp4 = *(pt_x++); 2924f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_sum = fxp_mac_16by16(tmp1, *(pt_fir++), 0x00002000L); 2934f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_sum = fxp_mac_16by16(tmp2, *(pt_fir++), L_sum); 2944f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_sum = fxp_mac_16by16(tmp3, *(pt_fir++), L_sum); 2954f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_sum = fxp_mac_16by16(tmp4, *(pt_fir++), L_sum); 2964f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber tmp1 = *(pt_x++); 2974f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber tmp2 = *(pt_x++); 2984f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber tmp3 = *(pt_x++); 2994f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber tmp4 = *(pt_x++); 3004f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_sum = fxp_mac_16by16(tmp1, *(pt_fir++), L_sum); 3014f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_sum = fxp_mac_16by16(tmp2, *(pt_fir++), L_sum); 3024f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_sum = fxp_mac_16by16(tmp3, *(pt_fir++), L_sum); 3034f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_sum = fxp_mac_16by16(tmp4, *(pt_fir++), L_sum); 3044f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber tmp1 = *(pt_x++); 3054f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber tmp2 = *(pt_x++); 3064f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber tmp3 = *(pt_x++); 3074f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber tmp4 = *(pt_x++); 3084f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_sum = fxp_mac_16by16(tmp1, *(pt_fir++), L_sum); 3094f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_sum = fxp_mac_16by16(tmp2, *(pt_fir++), L_sum); 3104f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_sum = fxp_mac_16by16(tmp3, *(pt_fir++), L_sum); 3114f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_sum = fxp_mac_16by16(tmp4, *(pt_fir++), L_sum); 3124f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber tmp1 = *(pt_x++); 3134f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber tmp2 = *(pt_x++); 3144f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber tmp3 = *(pt_x++); 3154f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber tmp4 = *(pt_x++); 3164f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_sum = fxp_mac_16by16(tmp1, *(pt_fir++), L_sum); 3174f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_sum = fxp_mac_16by16(tmp2, *(pt_fir++), L_sum); 3184f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_sum = fxp_mac_16by16(tmp3, *(pt_fir++), L_sum); 3194f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_sum = fxp_mac_16by16(tmp4, *(pt_fir++), L_sum); 3204f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber tmp1 = *(pt_x++); 3214f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber tmp2 = *(pt_x++); 3224f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber tmp3 = *(pt_x++); 3234f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber tmp4 = *(pt_x++); 3244f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_sum = fxp_mac_16by16(tmp1, *(pt_fir++), L_sum); 3254f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_sum = fxp_mac_16by16(tmp2, *(pt_fir++), L_sum); 3264f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_sum = fxp_mac_16by16(tmp3, *(pt_fir++), L_sum); 3274f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_sum = fxp_mac_16by16(tmp4, *(pt_fir++), L_sum); 3284f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber tmp1 = *(pt_x++); 3294f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber tmp2 = *(pt_x++); 3304f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber tmp3 = *(pt_x++); 3314f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber tmp4 = *(pt_x++); 3324f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_sum = fxp_mac_16by16(tmp1, *(pt_fir++), L_sum); 3334f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_sum = fxp_mac_16by16(tmp2, *(pt_fir++), L_sum); 3344f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_sum = fxp_mac_16by16(tmp3, *(pt_fir++), L_sum); 3354f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_sum = fxp_mac_16by16(tmp4, *(pt_fir++), L_sum); 3364f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3374f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3384f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_sum = shl_int32(L_sum, 2); /* saturation can occur here */ 3394f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 340ffa34edd0fd50d558f9bc7a650d959b4d3d29123Martin Storsjo return ((int16)(L_sum >> 16)); 3414f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber} 3424f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 343