eden-emu/opus
15
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

New VQ search nearly fixed-point ready

Browse source
This commit is contained in:
Jean-Marc Valin 2008-12-05 00:28:28 -05:00
parent c7635b45f4
commit 09dc5a18c5
1 changed files with 7 additions and 3 deletions
Download patch file Download diff file Expand all files Collapse all files

10
libcelt/vq.c
View file

@ -141,18 +141,22 @@ void alg_quant(celt_norm_t *X, celt_mask_t *W, int N, int K, celt_norm_t *P, ec_
#if 0 #if 0
if (K > (N>>1)) if (K > (N>>1))
{ {
float sum=0; celt_word32_t sum=0;
j=0; do { j=0; do {
sum += X[j]; sum += X[j];
} while (++j<N); } while (++j<N);
sum = K/(EPSILON+sum); sum = DIV32(SHL32(EXTEND32(K),15),EPSILON+sum);
j=0; do { j=0; do {
#ifdef FIXED_POINT
iy[j] = MULT16_32_Q15(X[j],sum);
#else
iy[j] = floor(sum*X[j]); iy[j] = floor(sum*X[j]);
#endif
y[j] = SHL16(iy[j],yshift); y[j] = SHL16(iy[j],yshift);
yy = MAC16_16(yy, y[j],y[j]); yy = MAC16_16(yy, y[j],y[j]);
xy = MAC16_16(xy, X[j],y[j]); xy = MAC16_16(xy, X[j],y[j]);
yp += P[j]*y[j]; yp += P[j]*y[j];
y[j] *= 2;; y[j] *= 2;
pulsesLeft -= iy[j]; pulsesLeft -= iy[j];
} while (++j<N); } while (++j<N);
} }