16-bit clean shift in lsp_to_lpc()
[speexdsp.git] / libspeex / misc.c
index 4b0223b..53bdd0b 100644 (file)
@@ -1,5 +1,5 @@
 /* Copyright (C) 2002-2005 Jean-Marc Valin 
-   File: mics.c
+   File: misc.c
    Various utility routines for Speex
 
    Redistribution and use in source and binary forms, with or without
 #include <stdio.h>
 #include "misc.h"
 
+#ifdef USER_MISC
+#include "user_misc.h"
+#endif
+
 #ifdef BFIN_ASM
 #include "misc_bfin.h"
 #endif
@@ -127,28 +131,40 @@ void speex_memset_bytes(char *dst, char src, int nbytes)
 }
 #endif
 
+#ifndef OVERRIDE_SPEEX_ALLOC
 void *speex_alloc (int size)
 {
    return calloc(size,1);
 }
+#endif
+
+#ifndef OVERRIDE_SPEEX_ALLOC_SCRATCH
 void *speex_alloc_scratch (int size)
 {
    return calloc(size,1);
 }
+#endif
 
+#ifndef OVERRIDE_SPEEX_REALLOC
 void *speex_realloc (void *ptr, int size)
 {
    return realloc(ptr, size);
 }
+#endif
 
+#ifndef OVERRIDE_SPEEX_FREE
 void speex_free (void *ptr)
 {
    free(ptr);
 }
+#endif
+
+#ifndef OVERRIDE_SPEEX_FREE_SCRATCH
 void speex_free_scratch (void *ptr)
 {
    free(ptr);
 }
+#endif
 
 #ifndef OVERRIDE_SPEEX_MOVE
 void *speex_move (void *dest, void *src, int n)
@@ -157,27 +173,35 @@ void *speex_move (void *dest, void *src, int n)
 }
 #endif
 
+#ifndef OVERRIDE_SPEEX_ERROR
 void speex_error(const char *str)
 {
    fprintf (stderr, "Fatal error: %s\n", str);
    exit(1);
 }
+#endif
 
+#ifndef OVERRIDE_SPEEX_WARNING
 void speex_warning(const char *str)
 {
    fprintf (stderr, "warning: %s\n", str);
 }
+#endif
 
+#ifndef OVERRIDE_SPEEX_WARNING_INT
 void speex_warning_int(const char *str, int val)
 {
    fprintf (stderr, "warning: %s %d\n", str, val);
 }
+#endif
 
 #ifdef FIXED_POINT
-spx_word32_t speex_rand(spx_word16_t std, spx_int32_t *seed)
+spx_word16_t speex_rand(spx_word16_t std, spx_int32_t *seed)
 {
+   spx_word32_t res;
    *seed = 1664525 * *seed + 1013904223;
-   return MULT16_16(EXTRACT16(SHR32(*seed,16)),std);
+   res = MULT16_16(EXTRACT16(SHR32(*seed,16)),std);
+   return PSHR32(SUB32(res, SHR(res, 3)),14);
 }
 #else
 spx_word16_t speex_rand(spx_word16_t std, spx_int32_t *seed)
@@ -187,26 +211,15 @@ spx_word16_t speex_rand(spx_word16_t std, spx_int32_t *seed)
    union {int i; float f;} ran;
    *seed = 1664525 * *seed + 1013904223;
    ran.i = jflone | (jflmsk & *seed);
-   ran.f -= 1;
-   return 1.7321*std*ran.f;
+   ran.f -= 1.5;
+   return 3.4642*std*ran.f;
 }
 #endif
 
-void speex_rand_vec(float std, spx_sig_t *data, int len)
-{
-   int i;
-   for (i=0;i<len;i++)
-      data[i]+=SIG_SCALING*3*std*((((float)rand())/RAND_MAX)-.5);
-}
-
-
-/*float speex_rand(float std)
-{
-   return 3*std*((((float)rand())/RAND_MAX)-.5);
-}*/
-
+#ifndef OVERRIDE_SPEEX_PUTC
 void _speex_putc(int ch, void *file)
 {
    FILE *f = (FILE *)file;
    fprintf(f, "%c", ch);
 }
+#endif