123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147 |
- #ifdef HAVE_CONFIG_H
- #include "config.h"
- #endif
- #include "vq.h"
- #include "stack_alloc.h"
- #include "arch.h"
- #ifdef _USE_SSE
- #include <xmmintrin.h>
- #include "vq_sse.h"
- #elif defined(SHORTCUTS) && (defined(ARM4_ASM) || defined(ARM5E_ASM))
- #include "vq_arm4.h"
- #elif defined(BFIN_ASM)
- #include "vq_bfin.h"
- #endif
- int scal_quant(spx_word16_t in, const spx_word16_t *boundary, int entries)
- {
- int i=0;
- while (i<entries-1 && in>boundary[0])
- {
- boundary++;
- i++;
- }
- return i;
- }
- int scal_quant32(spx_word32_t in, const spx_word32_t *boundary, int entries)
- {
- int i=0;
- while (i<entries-1 && in>boundary[0])
- {
- boundary++;
- i++;
- }
- return i;
- }
- #ifndef OVERRIDE_VQ_NBEST
- void vq_nbest(spx_word16_t *in, const spx_word16_t *codebook, int len, int entries, spx_word32_t *E, int N, int *nbest, spx_word32_t *best_dist, char *stack)
- {
- int i,j,k,used;
- used = 0;
- for (i=0;i<entries;i++)
- {
- spx_word32_t dist=0;
- for (j=0;j<len;j++)
- dist = MAC16_16(dist,in[j],*codebook++);
- #ifdef FIXED_POINT
- dist=SUB32(SHR32(E[i],1),dist);
- #else
- dist=.5f*E[i]-dist;
- #endif
- if (i<N || dist<best_dist[N-1])
- {
- for (k=N-1; (k >= 1) && (k > used || dist < best_dist[k-1]); k--)
- {
- best_dist[k]=best_dist[k-1];
- nbest[k] = nbest[k-1];
- }
- best_dist[k]=dist;
- nbest[k]=i;
- used++;
- }
- }
- }
- #endif
- #ifndef OVERRIDE_VQ_NBEST_SIGN
- void vq_nbest_sign(spx_word16_t *in, const spx_word16_t *codebook, int len, int entries, spx_word32_t *E, int N, int *nbest, spx_word32_t *best_dist, char *stack)
- {
- int i,j,k, sign, used;
- used=0;
- for (i=0;i<entries;i++)
- {
- spx_word32_t dist=0;
- for (j=0;j<len;j++)
- dist = MAC16_16(dist,in[j],*codebook++);
- if (dist>0)
- {
- sign=0;
- dist=-dist;
- } else
- {
- sign=1;
- }
- #ifdef FIXED_POINT
- dist = ADD32(dist,SHR32(E[i],1));
- #else
- dist = ADD32(dist,.5f*E[i]);
- #endif
- if (i<N || dist<best_dist[N-1])
- {
- for (k=N-1; (k >= 1) && (k > used || dist < best_dist[k-1]); k--)
- {
- best_dist[k]=best_dist[k-1];
- nbest[k] = nbest[k-1];
- }
- best_dist[k]=dist;
- nbest[k]=i;
- used++;
- if (sign)
- nbest[k]+=entries;
- }
- }
- }
- #endif
|