ref: 25aeed74f17bd0c99f549e91e117e3dd1d908099
parent: 89c5fd161bcf529e5e51add61c908d4cc4135b3b
author: Jean-Marc Valin <[email protected]>
date: Wed Mar 26 09:45:01 EDT 2008
changed 1*rcp(x) to just rcp(x)
--- a/libcelt/bands.c
+++ b/libcelt/bands.c
@@ -126,7 +126,7 @@
shift = celt_zlog2(bank[i*C+c])-13;
E = VSHR32(bank[i*C+c], shift);
if (E>0)
- g = EXTRACT16(celt_div(QCONST32(1.f,28),MULT16_16_Q14(E,sqrtC_1[C-1])));
+ g = EXTRACT16(celt_rcp(SHR32(MULT16_16(E,sqrtC_1[C-1]),11)));
else
g = 0;
for (j=B*eBands[i];j<B*eBands[i+1];j++)