source: roaraudio/libroardsp/rms.c @ 5961:06e7fd9e4c25

Last change on this file since 5961:06e7fd9e4c25 was 5961:06e7fd9e4c25, checked in by phi, 10 years ago

Updates of copyright and license headers

File size: 4.1 KB
Line 
1//rms.c:
2
3/*
4 *      Copyright (C) Philipp 'ph3-der-loewe' Schafft - 2009-2014
5 *
6 *  This file is part of libroardsp a part of RoarAudio,
7 *  a cross-platform sound system for both, home and professional use.
8 *  See README for details.
9 *
10 *  This file is free software; you can redistribute it and/or modify
11 *  it under the terms of the GNU General Public License version 3
12 *  as published by the Free Software Foundation.
13 *
14 *  libroardsp is distributed in the hope that it will be useful,
15 *  but WITHOUT ANY WARRANTY; without even the implied warranty of
16 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
17 *  GNU General Public License for more details.
18 *
19 *  You should have received a copy of the GNU General Public License
20 *  along with this software; see the file COPYING.  If not, write to
21 *  the Free Software Foundation, 51 Franklin Street, Fifth Floor,
22 *  Boston, MA 02110-1301, USA.
23 *
24 */
25
26#include "libroardsp.h"
27
28int64_t roar_rms2_1_8  (int8_t  * data, size_t samples) {
29 register int64_t s = 0;
30 register size_t  i;
31
32 for (i = 0; i < samples; i++)
33  s += data[i] * data[i];
34
35 s /= samples;
36
37 return s;
38}
39
40int64_t roar_rms2_1_16 (int16_t * data, size_t samples) {
41 register int64_t s = 0;
42 register size_t  i;
43
44 for (i = 0; i < samples; i++)
45  s += data[i] * data[i];
46
47 s /= samples;
48
49 return s;
50}
51
52int64_t roar_rms2_1_32 (int32_t * data, size_t samples) {
53 register int64_t s = 0;
54 register size_t  i;
55
56 for (i = 0; i < samples; i++)
57  s += data[i] * data[i];
58
59 s /= samples;
60
61 return s;
62}
63
64
65int roar_rms2_1_8_2    (int8_t  * data, size_t samples, int64_t * rms) {
66 register int64_t a0 = 0, a1 = 0;
67 register size_t  i;
68
69 if ( samples == 0 )
70  return 0;
71
72 if ( data == NULL )
73  return -1;
74
75 if ( samples & 0x1 ) /* odd */
76  return -1;
77
78 for (i = 0; i < samples; i += 2) {
79  a0 += data[i+0] * data[i+0];
80  a1 += data[i+1] * data[i+1];
81 }
82
83 samples /= 2;
84
85 a0 /= samples;
86 a1 /= samples;
87
88 rms[0] = a0;
89 rms[1] = a1;
90
91 return 0;
92}
93
94int roar_rms2_1_16_2   (int16_t * data, size_t samples, int64_t * rms) {
95 register int64_t a0 = 0, a1 = 0;
96 register size_t  i;
97
98 if ( samples == 0 )
99  return 0;
100
101 if ( data == NULL )
102  return -1;
103
104 if ( samples & 0x1 ) /* odd */
105  return -1;
106
107 for (i = 0; i < samples; i += 2) {
108  a0 += data[i+0] * data[i+0];
109  a1 += data[i+1] * data[i+1];
110 }
111
112 samples /= 2;
113
114 a0 /= samples;
115 a1 /= samples;
116
117 rms[0] = a0;
118 rms[1] = a1;
119
120 return 0;
121}
122
123int roar_rms2_1_32_2   (int32_t * data, size_t samples, int64_t * rms) {
124 register int64_t a0 = 0, a1 = 0;
125 register size_t  i;
126
127 if ( samples == 0 )
128  return 0;
129
130 if ( data == NULL )
131  return -1;
132
133 if ( samples & 0x1 ) /* odd */
134  return -1;
135
136 for (i = 0; i < samples; i += 2) {
137  a0 += data[i+0] * data[i+0];
138  a1 += data[i+1] * data[i+1];
139 }
140
141 samples /= 2;
142
143 a0 /= samples;
144 a1 /= samples;
145
146 rms[0] = a0;
147 rms[1] = a1;
148
149 return 0;
150}
151
152
153int roar_rms2_1_8_n    (int8_t  * data, size_t samples, int64_t * rms, size_t n) {
154 if ( n == 0 )
155  return 0;
156
157 if ( rms == NULL )
158  return -1;
159
160 switch (n) {
161  case 1: return *rms = roar_rms2_1_8(data, samples); return *rms == -1 ? -1 : 0; break;
162  case 2: return roar_rms2_1_8_2(data, samples, rms); break;
163  default:
164   return -1;
165 }
166}
167
168int roar_rms2_1_16_n   (int16_t * data, size_t samples, int64_t * rms, size_t n) {
169 if ( n == 0 )
170  return 0;
171
172 if ( rms == NULL )
173  return -1;
174
175 switch (n) {
176  case 1: return *rms = roar_rms2_1_16(data, samples); return *rms == -1 ? -1 : 0; break;
177  case 2: return roar_rms2_1_16_2(data, samples, rms); break;
178  default:
179   return -1;
180 }
181}
182
183int roar_rms2_1_32_n   (int32_t * data, size_t samples, int64_t * rms, size_t n) {
184 if ( n == 0 )
185  return 0;
186
187 if ( rms == NULL )
188  return -1;
189
190 switch (n) {
191  case 1: return *rms = roar_rms2_1_32(data, samples); return *rms == -1 ? -1 : 0; break;
192  case 2: return roar_rms2_1_32_2(data, samples, rms); break;
193  default:
194   return -1;
195 }
196}
197
198
199int roar_rms2_1_b_n    (void    * data, size_t samples, int64_t * rms, size_t n, size_t bits) {
200 switch (bits) {
201  case  8: return roar_rms2_1_8_n(data, samples, rms, n);  break;
202  case 16: return roar_rms2_1_16_n(data, samples, rms, n); break;
203  case 32: return roar_rms2_1_32_n(data, samples, rms, n); break;
204  default:
205    return -1;
206 }
207}
208
209//ll
Note: See TracBrowser for help on using the repository browser.