Vignesh Venkatasubramanian | 2bd8b54 | 2014-02-20 10:50:35 -0800 | [diff] [blame] | 1 | /* Copyright (c) 2007-2008 CSIRO |
| 2 | Copyright (c) 2007-2009 Xiph.Org Foundation |
| 3 | Written by Jean-Marc Valin */ |
| 4 | /** |
| 5 | @file pitch.h |
| 6 | @brief Pitch analysis |
| 7 | */ |
| 8 | |
| 9 | /* |
| 10 | Redistribution and use in source and binary forms, with or without |
| 11 | modification, are permitted provided that the following conditions |
| 12 | are met: |
| 13 | |
| 14 | - Redistributions of source code must retain the above copyright |
| 15 | notice, this list of conditions and the following disclaimer. |
| 16 | |
| 17 | - Redistributions in binary form must reproduce the above copyright |
| 18 | notice, this list of conditions and the following disclaimer in the |
| 19 | documentation and/or other materials provided with the distribution. |
| 20 | |
| 21 | THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS |
| 22 | ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT |
| 23 | LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR |
| 24 | A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER |
| 25 | OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, |
| 26 | EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, |
| 27 | PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR |
| 28 | PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF |
| 29 | LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING |
| 30 | NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS |
| 31 | SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. |
| 32 | */ |
| 33 | |
| 34 | #ifndef PITCH_H |
| 35 | #define PITCH_H |
| 36 | |
| 37 | #include "modes.h" |
| 38 | #include "cpu_support.h" |
| 39 | |
| 40 | #if defined(__SSE__) && !defined(FIXED_POINT) |
| 41 | #include "x86/pitch_sse.h" |
| 42 | #endif |
| 43 | |
| 44 | #if defined(OPUS_ARM_ASM) && defined(FIXED_POINT) |
| 45 | # include "arm/pitch_arm.h" |
| 46 | #endif |
| 47 | |
| 48 | void pitch_downsample(celt_sig * OPUS_RESTRICT x[], opus_val16 * OPUS_RESTRICT x_lp, |
| 49 | int len, int C, int arch); |
| 50 | |
| 51 | void pitch_search(const opus_val16 * OPUS_RESTRICT x_lp, opus_val16 * OPUS_RESTRICT y, |
| 52 | int len, int max_pitch, int *pitch, int arch); |
| 53 | |
| 54 | opus_val16 remove_doubling(opus_val16 *x, int maxperiod, int minperiod, |
| 55 | int N, int *T0, int prev_period, opus_val16 prev_gain); |
| 56 | |
| 57 | /* OPT: This is the kernel you really want to optimize. It gets used a lot |
| 58 | by the prefilter and by the PLC. */ |
| 59 | #ifndef OVERRIDE_XCORR_KERNEL |
| 60 | static OPUS_INLINE void xcorr_kernel(const opus_val16 * x, const opus_val16 * y, opus_val32 sum[4], int len) |
| 61 | { |
| 62 | int j; |
| 63 | opus_val16 y_0, y_1, y_2, y_3; |
| 64 | celt_assert(len>=3); |
| 65 | y_3=0; /* gcc doesn't realize that y_3 can't be used uninitialized */ |
| 66 | y_0=*y++; |
| 67 | y_1=*y++; |
| 68 | y_2=*y++; |
| 69 | for (j=0;j<len-3;j+=4) |
| 70 | { |
| 71 | opus_val16 tmp; |
| 72 | tmp = *x++; |
| 73 | y_3=*y++; |
| 74 | sum[0] = MAC16_16(sum[0],tmp,y_0); |
| 75 | sum[1] = MAC16_16(sum[1],tmp,y_1); |
| 76 | sum[2] = MAC16_16(sum[2],tmp,y_2); |
| 77 | sum[3] = MAC16_16(sum[3],tmp,y_3); |
| 78 | tmp=*x++; |
| 79 | y_0=*y++; |
| 80 | sum[0] = MAC16_16(sum[0],tmp,y_1); |
| 81 | sum[1] = MAC16_16(sum[1],tmp,y_2); |
| 82 | sum[2] = MAC16_16(sum[2],tmp,y_3); |
| 83 | sum[3] = MAC16_16(sum[3],tmp,y_0); |
| 84 | tmp=*x++; |
| 85 | y_1=*y++; |
| 86 | sum[0] = MAC16_16(sum[0],tmp,y_2); |
| 87 | sum[1] = MAC16_16(sum[1],tmp,y_3); |
| 88 | sum[2] = MAC16_16(sum[2],tmp,y_0); |
| 89 | sum[3] = MAC16_16(sum[3],tmp,y_1); |
| 90 | tmp=*x++; |
| 91 | y_2=*y++; |
| 92 | sum[0] = MAC16_16(sum[0],tmp,y_3); |
| 93 | sum[1] = MAC16_16(sum[1],tmp,y_0); |
| 94 | sum[2] = MAC16_16(sum[2],tmp,y_1); |
| 95 | sum[3] = MAC16_16(sum[3],tmp,y_2); |
| 96 | } |
| 97 | if (j++<len) |
| 98 | { |
| 99 | opus_val16 tmp = *x++; |
| 100 | y_3=*y++; |
| 101 | sum[0] = MAC16_16(sum[0],tmp,y_0); |
| 102 | sum[1] = MAC16_16(sum[1],tmp,y_1); |
| 103 | sum[2] = MAC16_16(sum[2],tmp,y_2); |
| 104 | sum[3] = MAC16_16(sum[3],tmp,y_3); |
| 105 | } |
| 106 | if (j++<len) |
| 107 | { |
| 108 | opus_val16 tmp=*x++; |
| 109 | y_0=*y++; |
| 110 | sum[0] = MAC16_16(sum[0],tmp,y_1); |
| 111 | sum[1] = MAC16_16(sum[1],tmp,y_2); |
| 112 | sum[2] = MAC16_16(sum[2],tmp,y_3); |
| 113 | sum[3] = MAC16_16(sum[3],tmp,y_0); |
| 114 | } |
| 115 | if (j<len) |
| 116 | { |
| 117 | opus_val16 tmp=*x++; |
| 118 | y_1=*y++; |
| 119 | sum[0] = MAC16_16(sum[0],tmp,y_2); |
| 120 | sum[1] = MAC16_16(sum[1],tmp,y_3); |
| 121 | sum[2] = MAC16_16(sum[2],tmp,y_0); |
| 122 | sum[3] = MAC16_16(sum[3],tmp,y_1); |
| 123 | } |
| 124 | } |
| 125 | #endif /* OVERRIDE_XCORR_KERNEL */ |
| 126 | |
| 127 | #ifndef OVERRIDE_DUAL_INNER_PROD |
| 128 | static OPUS_INLINE void dual_inner_prod(const opus_val16 *x, const opus_val16 *y01, const opus_val16 *y02, |
| 129 | int N, opus_val32 *xy1, opus_val32 *xy2) |
| 130 | { |
| 131 | int i; |
| 132 | opus_val32 xy01=0; |
| 133 | opus_val32 xy02=0; |
| 134 | for (i=0;i<N;i++) |
| 135 | { |
| 136 | xy01 = MAC16_16(xy01, x[i], y01[i]); |
| 137 | xy02 = MAC16_16(xy02, x[i], y02[i]); |
| 138 | } |
| 139 | *xy1 = xy01; |
| 140 | *xy2 = xy02; |
| 141 | } |
| 142 | #endif |
| 143 | |
| 144 | #ifdef FIXED_POINT |
| 145 | opus_val32 |
| 146 | #else |
| 147 | void |
| 148 | #endif |
| 149 | celt_pitch_xcorr_c(const opus_val16 *_x, const opus_val16 *_y, |
| 150 | opus_val32 *xcorr, int len, int max_pitch); |
| 151 | |
| 152 | #if !defined(OVERRIDE_PITCH_XCORR) |
| 153 | /*Is run-time CPU detection enabled on this platform?*/ |
| 154 | # if defined(OPUS_HAVE_RTCD) |
| 155 | extern |
| 156 | # if defined(FIXED_POINT) |
| 157 | opus_val32 |
| 158 | # else |
| 159 | void |
| 160 | # endif |
| 161 | (*const CELT_PITCH_XCORR_IMPL[OPUS_ARCHMASK+1])(const opus_val16 *, |
| 162 | const opus_val16 *, opus_val32 *, int, int); |
| 163 | |
| 164 | # define celt_pitch_xcorr(_x, _y, xcorr, len, max_pitch, arch) \ |
| 165 | ((*CELT_PITCH_XCORR_IMPL[(arch)&OPUS_ARCHMASK])(_x, _y, \ |
| 166 | xcorr, len, max_pitch)) |
| 167 | # else |
| 168 | # define celt_pitch_xcorr(_x, _y, xcorr, len, max_pitch, arch) \ |
| 169 | ((void)(arch),celt_pitch_xcorr_c(_x, _y, xcorr, len, max_pitch)) |
| 170 | # endif |
| 171 | #endif |
| 172 | |
| 173 | #endif |