198913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project/* Copyright (C) 2005 Jean-Marc Valin, CSIRO, Christopher Montgomery 298913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project File: vorbis_psy.h 398913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 498913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project Redistribution and use in source and binary forms, with or without 598913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project modification, are permitted provided that the following conditions 698913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project are met: 798913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 898913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project - Redistributions of source code must retain the above copyright 998913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project notice, this list of conditions and the following disclaimer. 1098913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 1198913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project - Redistributions in binary form must reproduce the above copyright 1298913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project notice, this list of conditions and the following disclaimer in the 1398913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project documentation and/or other materials provided with the distribution. 1498913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 1598913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project - Neither the name of the Xiph.org Foundation nor the names of its 1698913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project contributors may be used to endorse or promote products derived from 1798913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project this software without specific prior written permission. 1898913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 1998913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS 2098913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT 2198913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR 2298913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE FOUNDATION OR 2398913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, 2498913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, 2598913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR 2698913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF 2798913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING 2898913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS 2998913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. 3098913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project*/ 3198913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 3298913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#ifndef VORBIS_PSY_H 3398913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#define VORBIS_PSY_H 3498913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 3598913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#ifdef VORBIS_PSYCHO 3698913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 3798913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#include "smallft.h" 3898913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#define P_BANDS 17 /* 62Hz to 16kHz */ 3998913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#define NOISE_COMPAND_LEVELS 40 4098913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 4198913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 4298913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#define todB(x) ((x)>1e-13?log((x)*(x))*4.34294480f:-30) 4398913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#define fromdB(x) (exp((x)*.11512925f)) 4498913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 4598913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project/* The bark scale equations are approximations, since the original 4698913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project table was somewhat hand rolled. The below are chosen to have the 4798913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project best possible fit to the rolled tables, thus their somewhat odd 4898913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project appearance (these are more accurate and over a longer range than 4998913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project the oft-quoted bark equations found in the texts I have). The 5098913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project approximations are valid from 0 - 30kHz (nyquist) or so. 5198913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 5298913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project all f in Hz, z in Bark */ 5398913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 5498913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#define toBARK(n) (13.1f*atan(.00074f*(n))+2.24f*atan((n)*(n)*1.85e-8f)+1e-4f*(n)) 5598913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#define fromBARK(z) (102.f*(z)-2.f*pow(z,2.f)+.4f*pow(z,3.f)+pow(1.46f,z)-1.f) 5698913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 5798913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project/* Frequency to octave. We arbitrarily declare 63.5 Hz to be octave 5898913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 0.0 */ 5998913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 6098913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#define toOC(n) (log(n)*1.442695f-5.965784f) 6198913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#define fromOC(o) (exp(((o)+5.965784f)*.693147f)) 6298913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 6398913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 6498913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Projecttypedef struct { 6598913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 6698913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project float noisewindowlo; 6798913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project float noisewindowhi; 6898913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project int noisewindowlomin; 6998913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project int noisewindowhimin; 7098913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project int noisewindowfixed; 7198913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project float noiseoff[P_BANDS]; 7298913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project float noisecompand[NOISE_COMPAND_LEVELS]; 7398913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 7498913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project} VorbisPsyInfo; 7598913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 7698913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 7798913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 7898913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Projecttypedef struct { 7998913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project int n; 8098913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project int rate; 8198913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project struct drft_lookup lookup; 8298913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project VorbisPsyInfo *vi; 8398913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 8498913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project float *window; 8598913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project float *noiseoffset; 8698913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project long *bark; 8798913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 8898913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project} VorbisPsy; 8998913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 9098913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 9198913fed6520d8849fb2e246be943e04474aefaThe Android Open Source ProjectVorbisPsy *vorbis_psy_init(int rate, int size); 9298913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Projectvoid vorbis_psy_destroy(VorbisPsy *psy); 9398913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Projectvoid compute_curve(VorbisPsy *psy, float *audio, float *curve); 9498913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Projectvoid curve_to_lpc(VorbisPsy *psy, float *curve, float *awk1, float *awk2, int ord); 9598913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project 9698913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#endif 9798913fed6520d8849fb2e246be943e04474aefaThe Android Open Source Project#endif 98