fixed-point: Wrapper for the 32-bit complex FFT used in the MDCT so we can use
authorJean-Marc Valin <jean-marc.valin@usherbrooke.ca>
Fri, 21 Mar 2008 01:09:27 +0000 (12:09 +1100)
committerJean-Marc Valin <jean-marc.valin@usherbrooke.ca>
Fri, 21 Mar 2008 01:09:27 +0000 (12:09 +1100)
the TI dsplib FFT.

libcelt/Makefile.am
libcelt/kfft_double.h [new file with mode: 0644]
libcelt/mdct.c

index b7f43f4..29e828a 100644 (file)
@@ -25,7 +25,7 @@ libcelt_la_LDFLAGS = -version-info @CELT_LT_CURRENT@:@CELT_LT_REVISION@:@CELT_LT
 
 noinst_HEADERS = _kiss_fft_guts.h arch.h bands.h fixed_c5x.h \
        cwrs.h ecintrin.h entcode.h entdec.h entenc.h fixed_generic.h \
-       kfft_single.h kiss_fft.h kiss_fftr.h laplace.h mdct.h mfrngcod.h \
+       kfft_double.h kfft_single.h kiss_fft.h kiss_fftr.h laplace.h mdct.h mfrngcod.h \
        mathops.h modes.h os_support.h pgain_table.h pitch.h psy.h \
        quant_bands.h quant_pitch.h rate.h stack_alloc.h vq.h
 
diff --git a/libcelt/kfft_double.h b/libcelt/kfft_double.h
new file mode 100644 (file)
index 0000000..341d818
--- /dev/null
@@ -0,0 +1,67 @@
+/* (C) 2008 Jean-Marc Valin, CSIRO
+*/
+/*
+   Redistribution and use in source and binary forms, with or without
+   modification, are permitted provided that the following conditions
+   are met:
+   
+   - Redistributions of source code must retain the above copyright
+   notice, this list of conditions and the following disclaimer.
+   
+   - Redistributions in binary form must reproduce the above copyright
+   notice, this list of conditions and the following disclaimer in the
+   documentation and/or other materials provided with the distribution.
+   
+   - Neither the name of the Xiph.org Foundation nor the names of its
+   contributors may be used to endorse or promote products derived from
+   this software without specific prior written permission.
+   
+   THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
+   ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
+   LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
+   A PARTICULAR PURPOSE ARE DISCLAIMED.  IN NO EVENT SHALL THE FOUNDATION OR
+   CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
+   EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
+   PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
+   PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
+   LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
+   NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
+   SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+*/
+
+#ifndef KFFT_DOUBLE_H
+#define KFFT_DOUBLE_H
+
+#ifdef ENABLE_TI_DSPLIB
+
+#include "dsplib.h"
+
+#define cpx32_fft_alloc(length) NULL
+#define cpx32_fft_free(state)
+
+#define cpx32_fft(state, X, Y, nx)\
+    (\
+      cfft_SCALE(X,nx),\
+      cbrev(X,X,nx),\
+    )
+
+#define cpx32_ifft(state, X, Y, nx) \
+    (\
+      cifft_NOSCALE(X,nx),\
+      cbrev(X,Y,nx)\
+    )
+
+
+#else /* ENABLE_TI_DSPLIB */
+
+#include "kiss_fft.h"
+#include "_kiss_fft_guts.h"
+
+#define cpx32_fft_alloc(length) kiss_fft_alloc(length, 0, 0);
+#define cpx32_fft_free(state) kiss_fft_free(state)
+#define cpx32_fft(state, X, Y, nx) kiss_fft(state,X, Y)
+#define cpx32_ifft(state, X, Y, nx) kiss_ifft(state,X, Y)
+
+#endif /* !ENABLE_TI_DSPLIB */
+
+#endif /* KFFT_DOUBLE_H */
index 85355d1..5a60822 100644 (file)
 #endif
 
 #include "mdct.h"
-#include "kiss_fft.h"
+#include "kfft_double.h"
 #include <math.h>
 #include "os_support.h"
-#include "_kiss_fft_guts.h"
 #include "mathops.h"
 #include "stack_alloc.h"
 
@@ -65,7 +64,7 @@ void mdct_init(mdct_lookup *l,int N)
    l->n = N;
    N2 = N/2;
    N4 = N/4;
-   l->kfft = kiss_fft_alloc(N4, NULL, NULL);
+   l->kfft = cpx32_fft_alloc(N4);
    l->trig = (kiss_twiddle_scalar*)celt_alloc(N2*sizeof(kiss_twiddle_scalar));
    /* We have enough points that sine isn't necessary */
 #if defined(FIXED_POINT)
@@ -84,7 +83,7 @@ void mdct_init(mdct_lookup *l,int N)
 
 void mdct_clear(mdct_lookup *l)
 {
-   kiss_fft_free(l->kfft);
+   cpx32_fft_free(l->kfft);
    celt_free(l->trig);
 }
 
@@ -121,7 +120,7 @@ void mdct_forward(const mdct_lookup *l, kiss_fft_scalar *in, kiss_fft_scalar *ou
    }
 
    /* N/4 complex FFT, which should normally down-scale by 4/N (but doesn't now) */
-   kiss_fft(l->kfft, (const kiss_fft_cpx *)out, (kiss_fft_cpx *)f);
+   cpx32_fft(l->kfft, (const kiss_fft_cpx *)out, (kiss_fft_cpx *)f, N4);
 
    /* Post-rotate and apply the scaling if the FFT doesn't to it itself */
    for(i=0;i<N4;i++)
@@ -152,7 +151,7 @@ void mdct_backward(const mdct_lookup *l, kiss_fft_scalar *in, kiss_fft_scalar *o
    }
 
    /* Inverse N/4 complex FFT. This one should *not* downscale even in fixed-point */
-   kiss_ifft(l->kfft, (const kiss_fft_cpx *)out, (kiss_fft_cpx *)f);
+   cpx32_ifft(l->kfft, (const kiss_fft_cpx *)out, (kiss_fft_cpx *)f, N4);
    
    /* Post-rotate */
    for(i=0;i<N4;i++)