198913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project/* Copyright (C) 2006 Jean-Marc Valin */ 298913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project/** 398913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project @file filterbank.c 498913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project @brief Converting between psd and filterbank 598913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project */ 698913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project/* 798913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project Redistribution and use in source and binary forms, with or without 898913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project modification, are permitted provided that the following conditions are 998913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project met: 1098913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 1198913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 1. Redistributions of source code must retain the above copyright notice, 1298913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project this list of conditions and the following disclaimer. 1398913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 1498913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 2. Redistributions in binary form must reproduce the above copyright 1598913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project notice, this list of conditions and the following disclaimer in the 1698913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project documentation and/or other materials provided with the distribution. 1798913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 1898913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 3. The name of the author may not be used to endorse or promote products 1998913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project derived from this software without specific prior written permission. 2098913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 2198913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project THIS SOFTWARE IS PROVIDED BY THE AUTHOR ``AS IS'' AND ANY EXPRESS OR 2298913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES 2398913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE 2498913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project DISCLAIMED. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY DIRECT, 2598913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES 2698913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR 2798913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) 2898913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, 2998913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN 3098913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE 3198913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project POSSIBILITY OF SUCH DAMAGE. 3298913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project*/ 3398913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 3498913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#ifdef HAVE_CONFIG_H 3598913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#include "config.h" 3698913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#endif 3798913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 3898913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#include "filterbank.h" 3998913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#include "arch.h" 4098913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#include <math.h> 4198913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#include "math_approx.h" 4298913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#include "os_support.h" 4398913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 4498913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#ifdef FIXED_POINT 4598913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 4698913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#define toBARK(n) (MULT16_16(26829,spx_atan(SHR32(MULT16_16(97,n),2))) + MULT16_16(4588,spx_atan(MULT16_32_Q15(20,MULT16_16(n,n)))) + MULT16_16(3355,n)) 4798913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 4898913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#else 4998913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#define toBARK(n) (13.1f*atan(.00074f*(n))+2.24f*atan((n)*(n)*1.85e-8f)+1e-4f*(n)) 5098913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#endif 5198913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 5298913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#define toMEL(n) (2595.f*log10(1.f+(n)/700.f)) 5398913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 5498913fed6520d8849fb2e246be943e04474aefaThe Android Open Source ProjectFilterBank *filterbank_new(int banks, spx_word32_t sampling, int len, int type) 5598913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project{ 5698913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project FilterBank *bank; 5798913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project spx_word32_t df; 5898913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project spx_word32_t max_mel, mel_interval; 5998913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project int i; 6098913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project int id1; 6198913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project int id2; 6298913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project df = DIV32(SHL32(sampling,15),MULT16_16(2,len)); 6398913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project max_mel = toBARK(EXTRACT16(sampling/2)); 6498913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project mel_interval = PDIV32(max_mel,banks-1); 6598913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 6698913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project bank = (FilterBank*)speex_alloc(sizeof(FilterBank)); 6798913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project bank->nb_banks = banks; 6898913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project bank->len = len; 6998913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project bank->bank_left = (int*)speex_alloc(len*sizeof(int)); 7098913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project bank->bank_right = (int*)speex_alloc(len*sizeof(int)); 7198913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project bank->filter_left = (spx_word16_t*)speex_alloc(len*sizeof(spx_word16_t)); 7298913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project bank->filter_right = (spx_word16_t*)speex_alloc(len*sizeof(spx_word16_t)); 7398913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project /* Think I can safely disable normalisation that for fixed-point (and probably float as well) */ 7498913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#ifndef FIXED_POINT 7598913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project bank->scaling = (float*)speex_alloc(banks*sizeof(float)); 7698913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#endif 7798913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project for (i=0;i<len;i++) 7898913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project { 7998913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project spx_word16_t curr_freq; 8098913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project spx_word32_t mel; 8198913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project spx_word16_t val; 8298913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project curr_freq = EXTRACT16(MULT16_32_P15(i,df)); 8398913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project mel = toBARK(curr_freq); 8498913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project if (mel > max_mel) 8598913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project break; 8698913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#ifdef FIXED_POINT 8798913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project id1 = DIV32(mel,mel_interval); 8898913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#else 8998913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project id1 = (int)(floor(mel/mel_interval)); 9098913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#endif 9198913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project if (id1>banks-2) 9298913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project { 9398913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project id1 = banks-2; 9498913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project val = Q15_ONE; 9598913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project } else { 9698913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project val = DIV32_16(mel - id1*mel_interval,EXTRACT16(PSHR32(mel_interval,15))); 9798913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project } 9898913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project id2 = id1+1; 9998913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project bank->bank_left[i] = id1; 10098913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project bank->filter_left[i] = SUB16(Q15_ONE,val); 10198913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project bank->bank_right[i] = id2; 10298913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project bank->filter_right[i] = val; 10398913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project } 10498913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 10598913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project /* Think I can safely disable normalisation for fixed-point (and probably float as well) */ 10698913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#ifndef FIXED_POINT 10798913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project for (i=0;i<bank->nb_banks;i++) 10898913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project bank->scaling[i] = 0; 10998913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project for (i=0;i<bank->len;i++) 11098913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project { 11198913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project int id = bank->bank_left[i]; 11298913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project bank->scaling[id] += bank->filter_left[i]; 11398913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project id = bank->bank_right[i]; 11498913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project bank->scaling[id] += bank->filter_right[i]; 11598913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project } 11698913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project for (i=0;i<bank->nb_banks;i++) 11798913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project bank->scaling[i] = Q15_ONE/(bank->scaling[i]); 11898913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#endif 11998913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project return bank; 12098913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project} 12198913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 12298913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Projectvoid filterbank_destroy(FilterBank *bank) 12398913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project{ 12498913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project speex_free(bank->bank_left); 12598913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project speex_free(bank->bank_right); 12698913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project speex_free(bank->filter_left); 12798913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project speex_free(bank->filter_right); 12898913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#ifndef FIXED_POINT 12998913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project speex_free(bank->scaling); 13098913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#endif 13198913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project speex_free(bank); 13298913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project} 13398913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 13498913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Projectvoid filterbank_compute_bank32(FilterBank *bank, spx_word32_t *ps, spx_word32_t *mel) 13598913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project{ 13698913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project int i; 13798913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project for (i=0;i<bank->nb_banks;i++) 13898913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project mel[i] = 0; 13998913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 14098913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project for (i=0;i<bank->len;i++) 14198913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project { 14298913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project int id; 14398913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project id = bank->bank_left[i]; 14498913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project mel[id] += MULT16_32_P15(bank->filter_left[i],ps[i]); 14598913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project id = bank->bank_right[i]; 14698913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project mel[id] += MULT16_32_P15(bank->filter_right[i],ps[i]); 14798913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project } 14898913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project /* Think I can safely disable normalisation that for fixed-point (and probably float as well) */ 14998913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#ifndef FIXED_POINT 15098913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project /*for (i=0;i<bank->nb_banks;i++) 15198913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project mel[i] = MULT16_32_P15(Q15(bank->scaling[i]),mel[i]); 15298913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project */ 15398913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#endif 15498913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project} 15598913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 15698913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Projectvoid filterbank_compute_psd16(FilterBank *bank, spx_word16_t *mel, spx_word16_t *ps) 15798913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project{ 15898913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project int i; 15998913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project for (i=0;i<bank->len;i++) 16098913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project { 16198913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project spx_word32_t tmp; 16298913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project int id1, id2; 16398913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project id1 = bank->bank_left[i]; 16498913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project id2 = bank->bank_right[i]; 16598913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project tmp = MULT16_16(mel[id1],bank->filter_left[i]); 16698913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project tmp += MULT16_16(mel[id2],bank->filter_right[i]); 16798913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project ps[i] = EXTRACT16(PSHR32(tmp,15)); 16898913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project } 16998913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project} 17098913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 17198913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 17298913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#ifndef FIXED_POINT 17398913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Projectvoid filterbank_compute_bank(FilterBank *bank, float *ps, float *mel) 17498913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project{ 17598913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project int i; 17698913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project for (i=0;i<bank->nb_banks;i++) 17798913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project mel[i] = 0; 17898913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 17998913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project for (i=0;i<bank->len;i++) 18098913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project { 18198913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project int id = bank->bank_left[i]; 18298913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project mel[id] += bank->filter_left[i]*ps[i]; 18398913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project id = bank->bank_right[i]; 18498913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project mel[id] += bank->filter_right[i]*ps[i]; 18598913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project } 18698913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project for (i=0;i<bank->nb_banks;i++) 18798913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project mel[i] *= bank->scaling[i]; 18898913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project} 18998913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 19098913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Projectvoid filterbank_compute_psd(FilterBank *bank, float *mel, float *ps) 19198913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project{ 19298913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project int i; 19398913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project for (i=0;i<bank->len;i++) 19498913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project { 19598913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project int id = bank->bank_left[i]; 19698913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project ps[i] = mel[id]*bank->filter_left[i]; 19798913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project id = bank->bank_right[i]; 19898913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project ps[i] += mel[id]*bank->filter_right[i]; 19998913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project } 20098913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project} 20198913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 20298913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Projectvoid filterbank_psy_smooth(FilterBank *bank, float *ps, float *mask) 20398913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project{ 20498913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project /* Low freq slope: 14 dB/Bark*/ 20598913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project /* High freq slope: 9 dB/Bark*/ 20698913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project /* Noise vs tone: 5 dB difference */ 20798913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project /* FIXME: Temporary kludge */ 20898913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project float bark[100]; 20998913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project int i; 21098913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project /* Assumes 1/3 Bark resolution */ 21198913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project float decay_low = 0.34145f; 21298913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project float decay_high = 0.50119f; 21398913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project filterbank_compute_bank(bank, ps, bark); 21498913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project for (i=1;i<bank->nb_banks;i++) 21598913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project { 21698913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project /*float decay_high = 13-1.6*log10(bark[i-1]); 21798913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project decay_high = pow(10,(-decay_high/30.f));*/ 21898913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project bark[i] = bark[i] + decay_high*bark[i-1]; 21998913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project } 22098913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project for (i=bank->nb_banks-2;i>=0;i--) 22198913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project { 22298913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project bark[i] = bark[i] + decay_low*bark[i+1]; 22398913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project } 22498913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project filterbank_compute_psd(bank, bark, mask); 22598913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project} 22698913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 22798913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#endif 228