niklase@google.com | 470e71d | 2011-07-07 08:21:25 +0000 | [diff] [blame] | 1 | /* |
kma@webrtc.org | 551fcc0 | 2012-02-07 18:03:11 +0000 | [diff] [blame] | 2 | * Copyright (c) 2012 The WebRTC project authors. All Rights Reserved. |
niklase@google.com | 470e71d | 2011-07-07 08:21:25 +0000 | [diff] [blame] | 3 | * |
| 4 | * Use of this source code is governed by a BSD-style license |
| 5 | * that can be found in the LICENSE file in the root of the source |
| 6 | * tree. An additional intellectual property rights grant can be found |
| 7 | * in the file PATENTS. All contributing project authors may |
| 8 | * be found in the AUTHORS file in the root of the source tree. |
| 9 | */ |
| 10 | |
pbos@webrtc.org | aa30bb7 | 2013-05-27 09:49:58 +0000 | [diff] [blame] | 11 | #include "webrtc/common_audio/signal_processing/include/signal_processing_library.h" |
niklase@google.com | 470e71d | 2011-07-07 08:21:25 +0000 | [diff] [blame] | 12 | |
kma@webrtc.org | 551fcc0 | 2012-02-07 18:03:11 +0000 | [diff] [blame] | 13 | // TODO(Bjornv): Change the function parameter order to WebRTC code style. |
kma@webrtc.org | 0221b78 | 2012-09-08 00:09:26 +0000 | [diff] [blame] | 14 | // C version of WebRtcSpl_DownsampleFast() for generic platforms. |
| 15 | int WebRtcSpl_DownsampleFastC(const int16_t* data_in, |
Peter Kasting | dce40cf | 2015-08-24 14:52:23 -0700 | [diff] [blame] | 16 | size_t data_in_length, |
kma@webrtc.org | 0221b78 | 2012-09-08 00:09:26 +0000 | [diff] [blame] | 17 | int16_t* data_out, |
Peter Kasting | dce40cf | 2015-08-24 14:52:23 -0700 | [diff] [blame] | 18 | size_t data_out_length, |
kma@webrtc.org | 0221b78 | 2012-09-08 00:09:26 +0000 | [diff] [blame] | 19 | const int16_t* __restrict coefficients, |
Peter Kasting | dce40cf | 2015-08-24 14:52:23 -0700 | [diff] [blame] | 20 | size_t coefficients_length, |
kma@webrtc.org | 0221b78 | 2012-09-08 00:09:26 +0000 | [diff] [blame] | 21 | int factor, |
Peter Kasting | dce40cf | 2015-08-24 14:52:23 -0700 | [diff] [blame] | 22 | size_t delay) { |
| 23 | size_t i = 0; |
| 24 | size_t j = 0; |
kma@webrtc.org | 551fcc0 | 2012-02-07 18:03:11 +0000 | [diff] [blame] | 25 | int32_t out_s32 = 0; |
Peter Kasting | dce40cf | 2015-08-24 14:52:23 -0700 | [diff] [blame] | 26 | size_t endpos = delay + factor * (data_out_length - 1) + 1; |
niklase@google.com | 470e71d | 2011-07-07 08:21:25 +0000 | [diff] [blame] | 27 | |
kma@webrtc.org | 551fcc0 | 2012-02-07 18:03:11 +0000 | [diff] [blame] | 28 | // Return error if any of the running conditions doesn't meet. |
Peter Kasting | dce40cf | 2015-08-24 14:52:23 -0700 | [diff] [blame] | 29 | if (data_out_length == 0 || coefficients_length == 0 |
kma@webrtc.org | 551fcc0 | 2012-02-07 18:03:11 +0000 | [diff] [blame] | 30 | || data_in_length < endpos) { |
| 31 | return -1; |
| 32 | } |
niklase@google.com | 470e71d | 2011-07-07 08:21:25 +0000 | [diff] [blame] | 33 | |
kma@webrtc.org | 551fcc0 | 2012-02-07 18:03:11 +0000 | [diff] [blame] | 34 | for (i = delay; i < endpos; i += factor) { |
| 35 | out_s32 = 2048; // Round value, 0.5 in Q12. |
| 36 | |
| 37 | for (j = 0; j < coefficients_length; j++) { |
| 38 | out_s32 += coefficients[j] * data_in[i - j]; // Q12. |
niklase@google.com | 470e71d | 2011-07-07 08:21:25 +0000 | [diff] [blame] | 39 | } |
| 40 | |
kma@webrtc.org | 551fcc0 | 2012-02-07 18:03:11 +0000 | [diff] [blame] | 41 | out_s32 >>= 12; // Q0. |
niklase@google.com | 470e71d | 2011-07-07 08:21:25 +0000 | [diff] [blame] | 42 | |
kma@webrtc.org | 551fcc0 | 2012-02-07 18:03:11 +0000 | [diff] [blame] | 43 | // Saturate and store the output. |
| 44 | *data_out++ = WebRtcSpl_SatW32ToW16(out_s32); |
| 45 | } |
niklase@google.com | 470e71d | 2011-07-07 08:21:25 +0000 | [diff] [blame] | 46 | |
kma@webrtc.org | 551fcc0 | 2012-02-07 18:03:11 +0000 | [diff] [blame] | 47 | return 0; |
niklase@google.com | 470e71d | 2011-07-07 08:21:25 +0000 | [diff] [blame] | 48 | } |