Update SILK code using the CELT range coder
[opus.git] / src_SigProc_FLP / SKP_Silk_decimate2_coarse_FLP.c
similarity index 56%
rename from src/SKP_Silk_resample_1_2_coarsest.c
rename to src_SigProc_FLP/SKP_Silk_decimate2_coarse_FLP.c
index 8e72965..4dffbc8 100644 (file)
@@ -26,46 +26,47 @@ OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
 ***********************************************************************/\r
 \r
 /*                                                                      *\r
- * SKP_Silk_resample_1_2_coarsest.c                                   *\r
+ * SKP_Silk_decimate2_coarse.c                                        *\r
  *                                                                      *\r
- * Downsample by a factor 2, coarsest                                   *\r
+ * downsample by a factor 2, coarser                                    *\r
  *                                                                      *\r
- * Copyright 2006 (c), Skype Limited                                    *\r
+ * Copyright 2006 (c), Camino Networks                                  *\r
  * Date: 060221                                                         *\r
  *                                                                      */\r
-#include "SKP_Silk_SigProc_FIX.h"\r
+#include "SKP_Silk_SigProc_FLP.h"\r
 \r
+/* coefficients for coarser 2-fold resampling */\r
+static SKP_float A20c_FLP[ 2 ] = {0.064666748046875f, 0.508514404296875f};\r
+static SKP_float A21c_FLP[ 2 ] = {0.245666503906250f, 0.819732666015625f};\r
 \r
-/* Coefficients for coarsest 2-fold resampling */\r
-static SKP_int16 A20cst[ 1 ] = {  3786 };\r
-static SKP_int16 A21cst[ 1 ] = { 17908 };\r
-\r
-/* Downsample by a factor 2, coarsest */\r
-void SKP_Silk_resample_1_2_coarsest(\r
-    const SKP_int16     *in,                /* I:   16 kHz signal [2*len]   */\r
-    SKP_int32           *S,                 /* I/O: State vector [2]        */\r
-    SKP_int16           *out,               /* O:   8 kHz signal [len]      */\r
-    SKP_int32           *scratch,           /* I:   Scratch memory [3*len]  */\r
-    const SKP_int32     len                 /* I:   Number of OUTPUT samples*/\r
+/* downsample by a factor 2, coarser */\r
+void SKP_Silk_decimate2_coarse_FLP(\r
+    const SKP_float     *in,        /* I:   16 kHz signal [2*len]       */\r
+    SKP_float           *S,         /* I/O: state vector [4]            */\r
+    SKP_float           *out,       /* O:   8 kHz signal [len]          */\r
+    SKP_float           *scratch,   /* I:   scratch memory [3*len]      */\r
+    const SKP_int32     len         /* I:   number of OUTPUT samples    */\r
 )\r
 {\r
-    SKP_int32 k, idx;\r
+    SKP_int32 k;\r
 \r
-    /* De-interleave allpass inputs, and convert Q15 -> Q25 */\r
-    for( k = 0; k < len; k++ ) {\r
-        idx = SKP_LSHIFT( k, 1 );\r
-        scratch[ k ]       = SKP_LSHIFT( (SKP_int32)in[ idx     ], 10 );\r
-        scratch[ k + len ] = SKP_LSHIFT( (SKP_int32)in[ idx + 1 ], 10 );\r
+    /* de-interleave allpass inputs */\r
+    for ( k = 0; k < len; k++) {\r
+        scratch[ k ]       = in[ 2 * k ];\r
+        scratch[ k + len ] = in[ 2 * k + 1 ];\r
     }\r
 \r
-    idx = SKP_LSHIFT( len, 1 );\r
-    /* Allpass filters */\r
-    SKP_Silk_allpass_int( scratch,       S,     A21cst[ 0 ], scratch + idx, len );\r
-    SKP_Silk_allpass_int( scratch + len, S + 1, A20cst[ 0 ], scratch,       len );\r
-\r
-    /* Add two allpass outputs */\r
-    for( k = 0; k < len; k++ ) {\r
-        out[ k ] = (SKP_int16)SKP_SAT16( SKP_RSHIFT_ROUND( scratch[ k ] + scratch[ k + idx ], 11 ) );\r
-    }\r
+    /* allpass filters */\r
+    SKP_Silk_allpass_int_FLP( scratch, S + 0, A21c_FLP[ 0 ], scratch + 2 * len, len );\r
+    SKP_Silk_allpass_int_FLP( scratch + 2 * len, S + 1, A21c_FLP[ 1 ], scratch, len );\r
+    \r
+    SKP_Silk_allpass_int_FLP( scratch + len, S + 2, A20c_FLP[ 0 ], scratch + 2 * len, len );\r
+    SKP_Silk_allpass_int_FLP( scratch + 2 * len, S + 3, A20c_FLP[ 1 ], scratch + len, len );\r
+    \r
+    /* add two allpass outputs */\r
+    for ( k = 0; k < len; k++ ) {\r
+        out[ k ] = 0.5f * ( scratch[ k ] + scratch[ k + len ] );\r
+    }       \r
 }\r
 \r
+\r