cannam@154
|
1 /* Copyright (c) 2007-2008 CSIRO
|
cannam@154
|
2 Copyright (c) 2007-2008 Xiph.Org Foundation
|
cannam@154
|
3 Written by Jean-Marc Valin */
|
cannam@154
|
4 /*
|
cannam@154
|
5 Redistribution and use in source and binary forms, with or without
|
cannam@154
|
6 modification, are permitted provided that the following conditions
|
cannam@154
|
7 are met:
|
cannam@154
|
8
|
cannam@154
|
9 - Redistributions of source code must retain the above copyright
|
cannam@154
|
10 notice, this list of conditions and the following disclaimer.
|
cannam@154
|
11
|
cannam@154
|
12 - Redistributions in binary form must reproduce the above copyright
|
cannam@154
|
13 notice, this list of conditions and the following disclaimer in the
|
cannam@154
|
14 documentation and/or other materials provided with the distribution.
|
cannam@154
|
15
|
cannam@154
|
16 THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
|
cannam@154
|
17 ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
|
cannam@154
|
18 LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
|
cannam@154
|
19 A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER
|
cannam@154
|
20 OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
|
cannam@154
|
21 EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
|
cannam@154
|
22 PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
|
cannam@154
|
23 PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
|
cannam@154
|
24 LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
|
cannam@154
|
25 NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
|
cannam@154
|
26 SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
|
cannam@154
|
27 */
|
cannam@154
|
28
|
cannam@154
|
29 /* This is a simple MDCT implementation that uses a N/4 complex FFT
|
cannam@154
|
30 to do most of the work. It should be relatively straightforward to
|
cannam@154
|
31 plug in pretty much and FFT here.
|
cannam@154
|
32
|
cannam@154
|
33 This replaces the Vorbis FFT (and uses the exact same API), which
|
cannam@154
|
34 was a bit too messy and that was ending up duplicating code
|
cannam@154
|
35 (might as well use the same FFT everywhere).
|
cannam@154
|
36
|
cannam@154
|
37 The algorithm is similar to (and inspired from) Fabrice Bellard's
|
cannam@154
|
38 MDCT implementation in FFMPEG, but has differences in signs, ordering
|
cannam@154
|
39 and scaling in many places.
|
cannam@154
|
40 */
|
cannam@154
|
41
|
cannam@154
|
42 #ifndef MDCT_H
|
cannam@154
|
43 #define MDCT_H
|
cannam@154
|
44
|
cannam@154
|
45 #include "opus_defines.h"
|
cannam@154
|
46 #include "kiss_fft.h"
|
cannam@154
|
47 #include "arch.h"
|
cannam@154
|
48
|
cannam@154
|
49 typedef struct {
|
cannam@154
|
50 int n;
|
cannam@154
|
51 int maxshift;
|
cannam@154
|
52 const kiss_fft_state *kfft[4];
|
cannam@154
|
53 const kiss_twiddle_scalar * OPUS_RESTRICT trig;
|
cannam@154
|
54 } mdct_lookup;
|
cannam@154
|
55
|
cannam@154
|
56 #if defined(HAVE_ARM_NE10)
|
cannam@154
|
57 #include "arm/mdct_arm.h"
|
cannam@154
|
58 #endif
|
cannam@154
|
59
|
cannam@154
|
60
|
cannam@154
|
61 int clt_mdct_init(mdct_lookup *l,int N, int maxshift, int arch);
|
cannam@154
|
62 void clt_mdct_clear(mdct_lookup *l, int arch);
|
cannam@154
|
63
|
cannam@154
|
64 /** Compute a forward MDCT and scale by 4/N, trashes the input array */
|
cannam@154
|
65 void clt_mdct_forward_c(const mdct_lookup *l, kiss_fft_scalar *in,
|
cannam@154
|
66 kiss_fft_scalar * OPUS_RESTRICT out,
|
cannam@154
|
67 const opus_val16 *window, int overlap,
|
cannam@154
|
68 int shift, int stride, int arch);
|
cannam@154
|
69
|
cannam@154
|
70 /** Compute a backward MDCT (no scaling) and performs weighted overlap-add
|
cannam@154
|
71 (scales implicitly by 1/2) */
|
cannam@154
|
72 void clt_mdct_backward_c(const mdct_lookup *l, kiss_fft_scalar *in,
|
cannam@154
|
73 kiss_fft_scalar * OPUS_RESTRICT out,
|
cannam@154
|
74 const opus_val16 * OPUS_RESTRICT window,
|
cannam@154
|
75 int overlap, int shift, int stride, int arch);
|
cannam@154
|
76
|
cannam@154
|
77 #if !defined(OVERRIDE_OPUS_MDCT)
|
cannam@154
|
78 /* Is run-time CPU detection enabled on this platform? */
|
cannam@154
|
79 #if defined(OPUS_HAVE_RTCD) && defined(HAVE_ARM_NE10)
|
cannam@154
|
80
|
cannam@154
|
81 extern void (*const CLT_MDCT_FORWARD_IMPL[OPUS_ARCHMASK+1])(
|
cannam@154
|
82 const mdct_lookup *l, kiss_fft_scalar *in,
|
cannam@154
|
83 kiss_fft_scalar * OPUS_RESTRICT out, const opus_val16 *window,
|
cannam@154
|
84 int overlap, int shift, int stride, int arch);
|
cannam@154
|
85
|
cannam@154
|
86 #define clt_mdct_forward(_l, _in, _out, _window, _overlap, _shift, _stride, _arch) \
|
cannam@154
|
87 ((*CLT_MDCT_FORWARD_IMPL[(arch)&OPUS_ARCHMASK])(_l, _in, _out, \
|
cannam@154
|
88 _window, _overlap, _shift, \
|
cannam@154
|
89 _stride, _arch))
|
cannam@154
|
90
|
cannam@154
|
91 extern void (*const CLT_MDCT_BACKWARD_IMPL[OPUS_ARCHMASK+1])(
|
cannam@154
|
92 const mdct_lookup *l, kiss_fft_scalar *in,
|
cannam@154
|
93 kiss_fft_scalar * OPUS_RESTRICT out, const opus_val16 *window,
|
cannam@154
|
94 int overlap, int shift, int stride, int arch);
|
cannam@154
|
95
|
cannam@154
|
96 #define clt_mdct_backward(_l, _in, _out, _window, _overlap, _shift, _stride, _arch) \
|
cannam@154
|
97 (*CLT_MDCT_BACKWARD_IMPL[(arch)&OPUS_ARCHMASK])(_l, _in, _out, \
|
cannam@154
|
98 _window, _overlap, _shift, \
|
cannam@154
|
99 _stride, _arch)
|
cannam@154
|
100
|
cannam@154
|
101 #else /* if defined(OPUS_HAVE_RTCD) && defined(HAVE_ARM_NE10) */
|
cannam@154
|
102
|
cannam@154
|
103 #define clt_mdct_forward(_l, _in, _out, _window, _overlap, _shift, _stride, _arch) \
|
cannam@154
|
104 clt_mdct_forward_c(_l, _in, _out, _window, _overlap, _shift, _stride, _arch)
|
cannam@154
|
105
|
cannam@154
|
106 #define clt_mdct_backward(_l, _in, _out, _window, _overlap, _shift, _stride, _arch) \
|
cannam@154
|
107 clt_mdct_backward_c(_l, _in, _out, _window, _overlap, _shift, _stride, _arch)
|
cannam@154
|
108
|
cannam@154
|
109 #endif /* end if defined(OPUS_HAVE_RTCD) && defined(HAVE_ARM_NE10) && !defined(FIXED_POINT) */
|
cannam@154
|
110 #endif /* end if !defined(OVERRIDE_OPUS_MDCT) */
|
cannam@154
|
111
|
cannam@154
|
112 #endif
|