aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorAndroid Build Coastguard Worker <android-build-coastguard-worker@google.com>2021-10-02 12:05:07 +0000
committerAndroid Build Coastguard Worker <android-build-coastguard-worker@google.com>2021-10-02 12:05:07 +0000
commitf9637b7a5f264e396b5f4532d55ee5e6b570f2eb (patch)
treea75cc0c4cbcda77b7acf7d351a0939f8d066afcb
parent53f6a8facd3466bf58f3210873f03d73a43dc122 (diff)
parentc636f3175eee0df7d427dbad010e17ba0de20526 (diff)
downloadtremolo-f9637b7a5f264e396b5f4532d55ee5e6b570f2eb.tar.gz
Snap for 7788333 from c636f3175eee0df7d427dbad010e17ba0de20526 to mainline-media-swcodec-release
Change-Id: Ic24b9e2c0ef36902cd576c91e17bf5b3cc83dccb
-rw-r--r--Tremolo/floor0.c32
1 files changed, 21 insertions, 11 deletions
diff --git a/Tremolo/floor0.c b/Tremolo/floor0.c
index 7e74d1e..17ba4cf 100644
--- a/Tremolo/floor0.c
+++ b/Tremolo/floor0.c
@@ -103,28 +103,34 @@ static inline ogg_int32_t vorbis_coslook2_i(long a){
(COS_LOOKUP_I_SHIFT-LSP_FRACBITS+14);
}
-static const ogg_uint16_t barklook[54]={
+/* Values in barklook are defined such that toBARK(x) is an approximation to
+ POW(2, 15) * ((13.1*ATAN(0.00074*(x)))+(2.24*ATAN((x)*(x)*0.0000000185))+(0.0001*(x))) */
+static const ogg_uint16_t barklook[]={
0,51,102,154, 206,258,311,365,
420,477,535,594, 656,719,785,854,
926,1002,1082,1166, 1256,1352,1454,1564,
1683,1812,1953,2107, 2276,2463,2670,2900,
3155,3440,3756,4106, 4493,4919,5387,5901,
6466,7094,7798,8599, 9528,10623,11935,13524,
- 15453,17775,20517,23667, 27183,31004
+ 15453,17775,20517,23667, 27183,31004,35069
};
/* used in init only; interpolate the long way */
-static inline ogg_int32_t toBARK(int n){
+static inline ogg_int32_t toBARK(ogg_uint16_t n){
int i;
- for(i=0;i<54;i++)
- if(n>=barklook[i] && n<barklook[i+1])break;
-
- if(i==54){
- return 54<<14;
- }else{
- return (i<<14)+(((n-barklook[i])*
- ((1UL<<31)/(barklook[i+1]-barklook[i])))>>17);
+ int barklook_size = (sizeof(barklook) / sizeof(barklook[0]));
+ for(i=1;i<barklook_size;i++){
+ if(n<barklook[i]){
+ i--;
+ return (i<<14)+(((n-barklook[i])*
+ ((1UL<<31)/(barklook[i+1]-barklook[i])))>>17);
+ }
}
+ /* for a valid input n, which is half of info->rate (i.e. max 32767
+ as info->rate is 16 bit unsigned value), loop above will return
+ an output. So the following return will be used only when toBARK()
+ is called with invalid value */
+ return (barklook_size-1)<<14;
}
static const unsigned char MLOOP_1[64]={
@@ -174,6 +180,8 @@ void vorbis_lsp_to_curve(ogg_int32_t *curve,int n,int ln,
#else
ogg_uint32_t nextbark=MULT31(imap>>1,tBnyq1);
#endif
+ /* nextbark is guaranteed to be less than 54 << 14 here and that ensures index
+ to barklook can at max be 53 and 54 here */
int nextf=barklook[nextbark>>14]+(((nextbark&0x3fff)*
(barklook[(nextbark>>14)+1]-barklook[nextbark>>14]))>>14);
@@ -340,6 +348,8 @@ void vorbis_lsp_to_curve(ogg_int32_t *curve,int n,int ln,
#else
nextbark=MULT31((map+1)*(imap>>1),tBnyq1);
#endif
+ /* nextbark is guaranteed to be less than 54 << 14 here and that ensures index
+ to barklook can at max be 53 and 54 here */
nextf=barklook[nextbark>>14]+
(((nextbark&0x3fff)*
(barklook[(nextbark>>14)+1]-barklook[nextbark>>14]))>>14);