niklase@google.com | 470e71d | 2011-07-07 08:21:25 +0000 | [diff] [blame] | 1 | /* |
kma@webrtc.org | beb1851 | 2012-03-01 20:03:26 +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 | |
| 11 | /* |
| 12 | * This file contains the implementation of functions |
kma@webrtc.org | 0221b78 | 2012-09-08 00:09:26 +0000 | [diff] [blame] | 13 | * WebRtcSpl_MaxAbsValueW16C() |
| 14 | * WebRtcSpl_MaxAbsValueW32C() |
| 15 | * WebRtcSpl_MaxValueW16C() |
| 16 | * WebRtcSpl_MaxValueW32C() |
| 17 | * WebRtcSpl_MinValueW16C() |
| 18 | * WebRtcSpl_MinValueW32C() |
kma@webrtc.org | 33ced9c | 2012-04-11 17:40:40 +0000 | [diff] [blame] | 19 | * WebRtcSpl_MaxAbsIndexW16() |
| 20 | * WebRtcSpl_MaxIndexW16() |
| 21 | * WebRtcSpl_MaxIndexW32() |
| 22 | * WebRtcSpl_MinIndexW16() |
niklase@google.com | 470e71d | 2011-07-07 08:21:25 +0000 | [diff] [blame] | 23 | * WebRtcSpl_MinIndexW32() |
| 24 | * |
niklase@google.com | 470e71d | 2011-07-07 08:21:25 +0000 | [diff] [blame] | 25 | */ |
| 26 | |
kma@webrtc.org | beb1851 | 2012-03-01 20:03:26 +0000 | [diff] [blame] | 27 | #include <stdlib.h> |
| 28 | |
Mirko Bonadei | 92ea95e | 2017-09-15 06:47:31 +0200 | [diff] [blame] | 29 | #include "rtc_base/checks.h" |
| 30 | #include "common_audio/signal_processing/include/signal_processing_library.h" |
Peter Kasting | 1380e26 | 2015-08-28 17:31:03 -0700 | [diff] [blame] | 31 | |
kma@webrtc.org | 33ced9c | 2012-04-11 17:40:40 +0000 | [diff] [blame] | 32 | // TODO(bjorn/kma): Consolidate function pairs (e.g. combine |
kma@webrtc.org | 0221b78 | 2012-09-08 00:09:26 +0000 | [diff] [blame] | 33 | // WebRtcSpl_MaxAbsValueW16C and WebRtcSpl_MaxAbsIndexW16 into a single one.) |
| 34 | // TODO(kma): Move the next six functions into min_max_operations_c.c. |
kma@webrtc.org | 33ced9c | 2012-04-11 17:40:40 +0000 | [diff] [blame] | 35 | |
kma@webrtc.org | 0221b78 | 2012-09-08 00:09:26 +0000 | [diff] [blame] | 36 | // Maximum absolute value of word16 vector. C version for generic platforms. |
Peter Kasting | dce40cf | 2015-08-24 14:52:23 -0700 | [diff] [blame] | 37 | int16_t WebRtcSpl_MaxAbsValueW16C(const int16_t* vector, size_t length) { |
| 38 | size_t i = 0; |
| 39 | int absolute = 0, maximum = 0; |
kma@webrtc.org | 33ced9c | 2012-04-11 17:40:40 +0000 | [diff] [blame] | 40 | |
kwiberg | 1e8ed4a | 2016-08-26 04:33:34 -0700 | [diff] [blame] | 41 | RTC_DCHECK_GT(length, 0); |
niklase@google.com | 470e71d | 2011-07-07 08:21:25 +0000 | [diff] [blame] | 42 | |
kma@webrtc.org | beb1851 | 2012-03-01 20:03:26 +0000 | [diff] [blame] | 43 | for (i = 0; i < length; i++) { |
| 44 | absolute = abs((int)vector[i]); |
| 45 | |
| 46 | if (absolute > maximum) { |
| 47 | maximum = absolute; |
niklase@google.com | 470e71d | 2011-07-07 08:21:25 +0000 | [diff] [blame] | 48 | } |
kma@webrtc.org | beb1851 | 2012-03-01 20:03:26 +0000 | [diff] [blame] | 49 | } |
| 50 | |
| 51 | // Guard the case for abs(-32768). |
| 52 | if (maximum > WEBRTC_SPL_WORD16_MAX) { |
| 53 | maximum = WEBRTC_SPL_WORD16_MAX; |
| 54 | } |
| 55 | |
| 56 | return (int16_t)maximum; |
niklase@google.com | 470e71d | 2011-07-07 08:21:25 +0000 | [diff] [blame] | 57 | } |
| 58 | |
kma@webrtc.org | 0221b78 | 2012-09-08 00:09:26 +0000 | [diff] [blame] | 59 | // Maximum absolute value of word32 vector. C version for generic platforms. |
Peter Kasting | dce40cf | 2015-08-24 14:52:23 -0700 | [diff] [blame] | 60 | int32_t WebRtcSpl_MaxAbsValueW32C(const int32_t* vector, size_t length) { |
kma@webrtc.org | 33ced9c | 2012-04-11 17:40:40 +0000 | [diff] [blame] | 61 | // Use uint32_t for the local variables, to accommodate the return value |
| 62 | // of abs(0x80000000), which is 0x80000000. |
kma@webrtc.org | 95c3d40 | 2012-04-02 03:55:20 +0000 | [diff] [blame] | 63 | |
kma@webrtc.org | 33ced9c | 2012-04-11 17:40:40 +0000 | [diff] [blame] | 64 | uint32_t absolute = 0, maximum = 0; |
Peter Kasting | dce40cf | 2015-08-24 14:52:23 -0700 | [diff] [blame] | 65 | size_t i = 0; |
kma@webrtc.org | 33ced9c | 2012-04-11 17:40:40 +0000 | [diff] [blame] | 66 | |
kwiberg | 1e8ed4a | 2016-08-26 04:33:34 -0700 | [diff] [blame] | 67 | RTC_DCHECK_GT(length, 0); |
kma@webrtc.org | 33ced9c | 2012-04-11 17:40:40 +0000 | [diff] [blame] | 68 | |
| 69 | for (i = 0; i < length; i++) { |
| 70 | absolute = abs((int)vector[i]); |
| 71 | if (absolute > maximum) { |
| 72 | maximum = absolute; |
henrika@webrtc.org | afedb63 | 2012-04-02 07:12:08 +0000 | [diff] [blame] | 73 | } |
kma@webrtc.org | 33ced9c | 2012-04-11 17:40:40 +0000 | [diff] [blame] | 74 | } |
| 75 | |
| 76 | maximum = WEBRTC_SPL_MIN(maximum, WEBRTC_SPL_WORD32_MAX); |
| 77 | |
| 78 | return (int32_t)maximum; |
niklase@google.com | 470e71d | 2011-07-07 08:21:25 +0000 | [diff] [blame] | 79 | } |
| 80 | |
kma@webrtc.org | 0221b78 | 2012-09-08 00:09:26 +0000 | [diff] [blame] | 81 | // Maximum value of word16 vector. C version for generic platforms. |
Peter Kasting | dce40cf | 2015-08-24 14:52:23 -0700 | [diff] [blame] | 82 | int16_t WebRtcSpl_MaxValueW16C(const int16_t* vector, size_t length) { |
kma@webrtc.org | 33ced9c | 2012-04-11 17:40:40 +0000 | [diff] [blame] | 83 | int16_t maximum = WEBRTC_SPL_WORD16_MIN; |
Peter Kasting | dce40cf | 2015-08-24 14:52:23 -0700 | [diff] [blame] | 84 | size_t i = 0; |
niklase@google.com | 470e71d | 2011-07-07 08:21:25 +0000 | [diff] [blame] | 85 | |
kwiberg | 1e8ed4a | 2016-08-26 04:33:34 -0700 | [diff] [blame] | 86 | RTC_DCHECK_GT(length, 0); |
henrika@webrtc.org | afedb63 | 2012-04-02 07:12:08 +0000 | [diff] [blame] | 87 | |
kma@webrtc.org | 33ced9c | 2012-04-11 17:40:40 +0000 | [diff] [blame] | 88 | for (i = 0; i < length; i++) { |
| 89 | if (vector[i] > maximum) |
| 90 | maximum = vector[i]; |
| 91 | } |
| 92 | return maximum; |
henrika@webrtc.org | afedb63 | 2012-04-02 07:12:08 +0000 | [diff] [blame] | 93 | } |
| 94 | |
kma@webrtc.org | 0221b78 | 2012-09-08 00:09:26 +0000 | [diff] [blame] | 95 | // Maximum value of word32 vector. C version for generic platforms. |
Peter Kasting | dce40cf | 2015-08-24 14:52:23 -0700 | [diff] [blame] | 96 | int32_t WebRtcSpl_MaxValueW32C(const int32_t* vector, size_t length) { |
kma@webrtc.org | 33ced9c | 2012-04-11 17:40:40 +0000 | [diff] [blame] | 97 | int32_t maximum = WEBRTC_SPL_WORD32_MIN; |
Peter Kasting | dce40cf | 2015-08-24 14:52:23 -0700 | [diff] [blame] | 98 | size_t i = 0; |
henrika@webrtc.org | afedb63 | 2012-04-02 07:12:08 +0000 | [diff] [blame] | 99 | |
kwiberg | 1e8ed4a | 2016-08-26 04:33:34 -0700 | [diff] [blame] | 100 | RTC_DCHECK_GT(length, 0); |
niklase@google.com | 470e71d | 2011-07-07 08:21:25 +0000 | [diff] [blame] | 101 | |
kma@webrtc.org | 33ced9c | 2012-04-11 17:40:40 +0000 | [diff] [blame] | 102 | for (i = 0; i < length; i++) { |
| 103 | if (vector[i] > maximum) |
| 104 | maximum = vector[i]; |
| 105 | } |
| 106 | return maximum; |
henrika@webrtc.org | afedb63 | 2012-04-02 07:12:08 +0000 | [diff] [blame] | 107 | } |
kma@webrtc.org | 95c3d40 | 2012-04-02 03:55:20 +0000 | [diff] [blame] | 108 | |
kma@webrtc.org | 0221b78 | 2012-09-08 00:09:26 +0000 | [diff] [blame] | 109 | // Minimum value of word16 vector. C version for generic platforms. |
Peter Kasting | dce40cf | 2015-08-24 14:52:23 -0700 | [diff] [blame] | 110 | int16_t WebRtcSpl_MinValueW16C(const int16_t* vector, size_t length) { |
kma@webrtc.org | 33ced9c | 2012-04-11 17:40:40 +0000 | [diff] [blame] | 111 | int16_t minimum = WEBRTC_SPL_WORD16_MAX; |
Peter Kasting | dce40cf | 2015-08-24 14:52:23 -0700 | [diff] [blame] | 112 | size_t i = 0; |
henrika@webrtc.org | afedb63 | 2012-04-02 07:12:08 +0000 | [diff] [blame] | 113 | |
kwiberg | 1e8ed4a | 2016-08-26 04:33:34 -0700 | [diff] [blame] | 114 | RTC_DCHECK_GT(length, 0); |
kma@webrtc.org | 33ced9c | 2012-04-11 17:40:40 +0000 | [diff] [blame] | 115 | |
| 116 | for (i = 0; i < length; i++) { |
| 117 | if (vector[i] < minimum) |
| 118 | minimum = vector[i]; |
| 119 | } |
| 120 | return minimum; |
| 121 | } |
| 122 | |
kma@webrtc.org | 0221b78 | 2012-09-08 00:09:26 +0000 | [diff] [blame] | 123 | // Minimum value of word32 vector. C version for generic platforms. |
Peter Kasting | dce40cf | 2015-08-24 14:52:23 -0700 | [diff] [blame] | 124 | int32_t WebRtcSpl_MinValueW32C(const int32_t* vector, size_t length) { |
kma@webrtc.org | 33ced9c | 2012-04-11 17:40:40 +0000 | [diff] [blame] | 125 | int32_t minimum = WEBRTC_SPL_WORD32_MAX; |
Peter Kasting | dce40cf | 2015-08-24 14:52:23 -0700 | [diff] [blame] | 126 | size_t i = 0; |
kma@webrtc.org | 33ced9c | 2012-04-11 17:40:40 +0000 | [diff] [blame] | 127 | |
kwiberg | 1e8ed4a | 2016-08-26 04:33:34 -0700 | [diff] [blame] | 128 | RTC_DCHECK_GT(length, 0); |
kma@webrtc.org | 33ced9c | 2012-04-11 17:40:40 +0000 | [diff] [blame] | 129 | |
| 130 | for (i = 0; i < length; i++) { |
| 131 | if (vector[i] < minimum) |
| 132 | minimum = vector[i]; |
| 133 | } |
| 134 | return minimum; |
| 135 | } |
kma@webrtc.org | 33ced9c | 2012-04-11 17:40:40 +0000 | [diff] [blame] | 136 | |
| 137 | // Index of maximum absolute value in a word16 vector. |
Peter Kasting | 1380e26 | 2015-08-28 17:31:03 -0700 | [diff] [blame] | 138 | size_t WebRtcSpl_MaxAbsIndexW16(const int16_t* vector, size_t length) { |
kma@webrtc.org | 33ced9c | 2012-04-11 17:40:40 +0000 | [diff] [blame] | 139 | // Use type int for local variables, to accomodate the value of abs(-32768). |
| 140 | |
Peter Kasting | dce40cf | 2015-08-24 14:52:23 -0700 | [diff] [blame] | 141 | size_t i = 0, index = 0; |
| 142 | int absolute = 0, maximum = 0; |
kma@webrtc.org | 33ced9c | 2012-04-11 17:40:40 +0000 | [diff] [blame] | 143 | |
kwiberg | 1e8ed4a | 2016-08-26 04:33:34 -0700 | [diff] [blame] | 144 | RTC_DCHECK_GT(length, 0); |
kma@webrtc.org | 33ced9c | 2012-04-11 17:40:40 +0000 | [diff] [blame] | 145 | |
| 146 | for (i = 0; i < length; i++) { |
| 147 | absolute = abs((int)vector[i]); |
| 148 | |
| 149 | if (absolute > maximum) { |
| 150 | maximum = absolute; |
| 151 | index = i; |
henrika@webrtc.org | afedb63 | 2012-04-02 07:12:08 +0000 | [diff] [blame] | 152 | } |
kma@webrtc.org | 33ced9c | 2012-04-11 17:40:40 +0000 | [diff] [blame] | 153 | } |
| 154 | |
Peter Kasting | 1380e26 | 2015-08-28 17:31:03 -0700 | [diff] [blame] | 155 | return index; |
kma@webrtc.org | 33ced9c | 2012-04-11 17:40:40 +0000 | [diff] [blame] | 156 | } |
| 157 | |
| 158 | // Index of maximum value in a word16 vector. |
Peter Kasting | 1380e26 | 2015-08-28 17:31:03 -0700 | [diff] [blame] | 159 | size_t WebRtcSpl_MaxIndexW16(const int16_t* vector, size_t length) { |
Peter Kasting | dce40cf | 2015-08-24 14:52:23 -0700 | [diff] [blame] | 160 | size_t i = 0, index = 0; |
kma@webrtc.org | 33ced9c | 2012-04-11 17:40:40 +0000 | [diff] [blame] | 161 | int16_t maximum = WEBRTC_SPL_WORD16_MIN; |
| 162 | |
kwiberg | 1e8ed4a | 2016-08-26 04:33:34 -0700 | [diff] [blame] | 163 | RTC_DCHECK_GT(length, 0); |
kma@webrtc.org | 33ced9c | 2012-04-11 17:40:40 +0000 | [diff] [blame] | 164 | |
| 165 | for (i = 0; i < length; i++) { |
| 166 | if (vector[i] > maximum) { |
| 167 | maximum = vector[i]; |
| 168 | index = i; |
| 169 | } |
| 170 | } |
| 171 | |
Peter Kasting | 1380e26 | 2015-08-28 17:31:03 -0700 | [diff] [blame] | 172 | return index; |
kma@webrtc.org | 33ced9c | 2012-04-11 17:40:40 +0000 | [diff] [blame] | 173 | } |
| 174 | |
| 175 | // Index of maximum value in a word32 vector. |
Peter Kasting | 1380e26 | 2015-08-28 17:31:03 -0700 | [diff] [blame] | 176 | size_t WebRtcSpl_MaxIndexW32(const int32_t* vector, size_t length) { |
Peter Kasting | dce40cf | 2015-08-24 14:52:23 -0700 | [diff] [blame] | 177 | size_t i = 0, index = 0; |
kma@webrtc.org | 33ced9c | 2012-04-11 17:40:40 +0000 | [diff] [blame] | 178 | int32_t maximum = WEBRTC_SPL_WORD32_MIN; |
| 179 | |
kwiberg | 1e8ed4a | 2016-08-26 04:33:34 -0700 | [diff] [blame] | 180 | RTC_DCHECK_GT(length, 0); |
kma@webrtc.org | 33ced9c | 2012-04-11 17:40:40 +0000 | [diff] [blame] | 181 | |
| 182 | for (i = 0; i < length; i++) { |
| 183 | if (vector[i] > maximum) { |
| 184 | maximum = vector[i]; |
| 185 | index = i; |
| 186 | } |
| 187 | } |
| 188 | |
Peter Kasting | 1380e26 | 2015-08-28 17:31:03 -0700 | [diff] [blame] | 189 | return index; |
niklase@google.com | 470e71d | 2011-07-07 08:21:25 +0000 | [diff] [blame] | 190 | } |
| 191 | |
| 192 | // Index of minimum value in a word16 vector. |
Peter Kasting | 1380e26 | 2015-08-28 17:31:03 -0700 | [diff] [blame] | 193 | size_t WebRtcSpl_MinIndexW16(const int16_t* vector, size_t length) { |
Peter Kasting | dce40cf | 2015-08-24 14:52:23 -0700 | [diff] [blame] | 194 | size_t i = 0, index = 0; |
kma@webrtc.org | 33ced9c | 2012-04-11 17:40:40 +0000 | [diff] [blame] | 195 | int16_t minimum = WEBRTC_SPL_WORD16_MAX; |
niklase@google.com | 470e71d | 2011-07-07 08:21:25 +0000 | [diff] [blame] | 196 | |
kwiberg | 1e8ed4a | 2016-08-26 04:33:34 -0700 | [diff] [blame] | 197 | RTC_DCHECK_GT(length, 0); |
kma@webrtc.org | 33ced9c | 2012-04-11 17:40:40 +0000 | [diff] [blame] | 198 | |
| 199 | for (i = 0; i < length; i++) { |
| 200 | if (vector[i] < minimum) { |
| 201 | minimum = vector[i]; |
| 202 | index = i; |
niklase@google.com | 470e71d | 2011-07-07 08:21:25 +0000 | [diff] [blame] | 203 | } |
kma@webrtc.org | 33ced9c | 2012-04-11 17:40:40 +0000 | [diff] [blame] | 204 | } |
niklase@google.com | 470e71d | 2011-07-07 08:21:25 +0000 | [diff] [blame] | 205 | |
Peter Kasting | 1380e26 | 2015-08-28 17:31:03 -0700 | [diff] [blame] | 206 | return index; |
niklase@google.com | 470e71d | 2011-07-07 08:21:25 +0000 | [diff] [blame] | 207 | } |
| 208 | |
| 209 | // Index of minimum value in a word32 vector. |
Peter Kasting | 1380e26 | 2015-08-28 17:31:03 -0700 | [diff] [blame] | 210 | size_t WebRtcSpl_MinIndexW32(const int32_t* vector, size_t length) { |
Peter Kasting | dce40cf | 2015-08-24 14:52:23 -0700 | [diff] [blame] | 211 | size_t i = 0, index = 0; |
kma@webrtc.org | 33ced9c | 2012-04-11 17:40:40 +0000 | [diff] [blame] | 212 | int32_t minimum = WEBRTC_SPL_WORD32_MAX; |
niklase@google.com | 470e71d | 2011-07-07 08:21:25 +0000 | [diff] [blame] | 213 | |
kwiberg | 1e8ed4a | 2016-08-26 04:33:34 -0700 | [diff] [blame] | 214 | RTC_DCHECK_GT(length, 0); |
kma@webrtc.org | 33ced9c | 2012-04-11 17:40:40 +0000 | [diff] [blame] | 215 | |
| 216 | for (i = 0; i < length; i++) { |
| 217 | if (vector[i] < minimum) { |
| 218 | minimum = vector[i]; |
| 219 | index = i; |
niklase@google.com | 470e71d | 2011-07-07 08:21:25 +0000 | [diff] [blame] | 220 | } |
kma@webrtc.org | 33ced9c | 2012-04-11 17:40:40 +0000 | [diff] [blame] | 221 | } |
| 222 | |
Peter Kasting | 1380e26 | 2015-08-28 17:31:03 -0700 | [diff] [blame] | 223 | return index; |
niklase@google.com | 470e71d | 2011-07-07 08:21:25 +0000 | [diff] [blame] | 224 | } |