Update SILK code using the CELT range coder
[opus.git] / src_SigProc_FIX / SKP_Silk_resampler_up2.c
similarity index 54%
rename from src/SKP_Silk_resample_1_2_coarse.c
rename to src_SigProc_FIX/SKP_Silk_resampler_up2.c
index f90b059..835e113 100644 (file)
@@ -26,49 +26,50 @@ OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
 ***********************************************************************/\r
 \r
 /*                                                                      *\r
- * SKP_Silk_resample_1_2_coarse.c                                     *\r
+ * SKP_Silk_resampler_up2.c                                           *\r
  *                                                                      *\r
- * Downsample by a factor 2, coarser                                    *\r
+ * Upsample by a factor 2, low quality                                  *\r
  *                                                                      *\r
- * Copyright 2006 (c), Skype Limited                                    *\r
- * Date: 060221                                                         *\r
+ * Copyright 2010 (c), Skype Limited                                    *\r
  *                                                                      */\r
+\r
 #include "SKP_Silk_SigProc_FIX.h"\r
+#include "SKP_Silk_resampler_rom.h"\r
 \r
-/* downsample by a factor 2, coarser */\r
-void SKP_Silk_resample_1_2_coarse(\r
-    const SKP_int16      *in,            /* I:   16 kHz signal [2*len]   */\r
-    SKP_int32            *S,             /* I/O: State vector [4]        */\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
+/* Upsample by a factor 2, low quality */\r
+void SKP_Silk_resampler_up2(\r
+    SKP_int32                           *S,         /* I/O: State vector [ 2 ]                  */\r
+    SKP_int16                           *out,       /* O:   Output signal [ 2 * len ]           */\r
+    const SKP_int16                     *in,        /* I:   Input signal [ len ]                */\r
+    SKP_int32                           len         /* I:   Number of input samples             */\r
 )\r
 {\r
-    SKP_int32 k, idx;\r
-    \r
-    /* Coefficients for coarser 2-fold resampling */\r
-    const SKP_int16 A20c[ 2 ] = { 2119, 16663 };\r
-    const SKP_int16 A21c[ 2 ] = { 8050, 26861 };\r
+    SKP_int32 k;\r
+    SKP_int32 in32, out32, Y, X;\r
 \r
-    /* De-interleave allpass inputs, and convert Q15 -> Q25 */\r
+    SKP_assert( SKP_Silk_resampler_up2_lq_0 > 0 );\r
+    SKP_assert( SKP_Silk_resampler_up2_lq_1 < 0 );\r
+    /* Internal variables and state are in Q10 format */\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
-    }\r
+        /* Convert to Q10 */\r
+        in32 = SKP_LSHIFT( (SKP_int32)in[ k ], 10 );\r
 \r
-    idx = SKP_LSHIFT( len, 1 );\r
-    /* Allpass filters */\r
-    SKP_Silk_allpass_int( scratch,       S,     A21c[ 0 ], scratch + idx, len );\r
-    SKP_Silk_allpass_int( scratch + idx, S + 1, A21c[ 1 ], scratch,       len );\r
+        /* All-pass section for even output sample */\r
+        Y      = SKP_SUB32( in32, S[ 0 ] );\r
+        X      = SKP_SMULWB( Y, SKP_Silk_resampler_up2_lq_0 );\r
+        out32  = SKP_ADD32( S[ 0 ], X );\r
+        S[ 0 ] = SKP_ADD32( in32, X );\r
 \r
-    SKP_Silk_allpass_int( scratch + len, S + 2, A20c[ 0 ], scratch + idx, len );\r
-    SKP_Silk_allpass_int( scratch + idx, S + 3, A20c[ 1 ], scratch + len, len );\r
+        /* Convert back to int16 and store to output */\r
+        out[ 2 * k ] = (SKP_int16)SKP_SAT16( SKP_RSHIFT_ROUND( out32, 10 ) );\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 + len ], 11 ) );\r
+        /* All-pass section for odd output sample */\r
+        Y      = SKP_SUB32( in32, S[ 1 ] );\r
+        X      = SKP_SMLAWB( Y, Y, SKP_Silk_resampler_up2_lq_1 );\r
+        out32  = SKP_ADD32( S[ 1 ], X );\r
+        S[ 1 ] = SKP_ADD32( in32, X );\r
+\r
+        /* Convert back to int16 and store to output */\r
+        out[ 2 * k + 1 ] = (SKP_int16)SKP_SAT16( SKP_RSHIFT_ROUND( out32, 10 ) );\r
     }\r
 }\r
-\r
-\r