1 | /* |
---|
2 | Copyright (C) 2003-2015 Paul Brossier <piem@aubio.org> |
---|
3 | |
---|
4 | This file is part of aubio. |
---|
5 | |
---|
6 | aubio is free software: you can redistribute it and/or modify |
---|
7 | it under the terms of the GNU General Public License as published by |
---|
8 | the Free Software Foundation, either version 3 of the License, or |
---|
9 | (at your option) any later version. |
---|
10 | |
---|
11 | aubio is distributed in the hope that it will be useful, |
---|
12 | but WITHOUT ANY WARRANTY; without even the implied warranty of |
---|
13 | MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the |
---|
14 | GNU General Public License for more details. |
---|
15 | |
---|
16 | You should have received a copy of the GNU General Public License |
---|
17 | along with aubio. If not, see <http://www.gnu.org/licenses/>. |
---|
18 | |
---|
19 | */ |
---|
20 | |
---|
21 | /** @file |
---|
22 | * Private include file |
---|
23 | * |
---|
24 | * This file is for inclusion from _within_ the library only. |
---|
25 | */ |
---|
26 | |
---|
27 | #ifndef AUBIO_PRIV_H |
---|
28 | #define AUBIO_PRIV_H |
---|
29 | |
---|
30 | /********************* |
---|
31 | * |
---|
32 | * External includes |
---|
33 | * |
---|
34 | */ |
---|
35 | |
---|
36 | #ifdef HAVE_CONFIG_H |
---|
37 | #include "config.h" |
---|
38 | #endif |
---|
39 | |
---|
40 | #ifdef HAVE_STDLIB_H |
---|
41 | #include <stdlib.h> |
---|
42 | #endif |
---|
43 | |
---|
44 | #ifdef HAVE_STDIO_H |
---|
45 | #include <stdio.h> |
---|
46 | #endif |
---|
47 | |
---|
48 | /* must be included before fftw3.h */ |
---|
49 | #ifdef HAVE_COMPLEX_H |
---|
50 | #include <complex.h> |
---|
51 | #endif |
---|
52 | |
---|
53 | #if defined(HAVE_FFTW3) || defined(HAVE_FFTW3F) |
---|
54 | #include <fftw3.h> |
---|
55 | #endif |
---|
56 | |
---|
57 | #ifdef HAVE_MATH_H |
---|
58 | #include <math.h> |
---|
59 | #endif |
---|
60 | |
---|
61 | #ifdef HAVE_STRING_H |
---|
62 | #include <string.h> |
---|
63 | #endif |
---|
64 | |
---|
65 | #ifdef HAVE_ERRNO_H |
---|
66 | #include <errno.h> |
---|
67 | #endif |
---|
68 | |
---|
69 | #ifdef HAVE_LIMITS_H |
---|
70 | #include <limits.h> // for CHAR_BIT, in C99 standard |
---|
71 | #endif |
---|
72 | |
---|
73 | #ifdef HAVE_STDARG_H |
---|
74 | #include <stdarg.h> |
---|
75 | #endif |
---|
76 | |
---|
77 | #if defined(HAVE_BLAS) // --enable-blas=true |
---|
78 | // check which cblas header we found |
---|
79 | #if defined(HAVE_ATLAS_CBLAS_H) |
---|
80 | #define HAVE_ATLAS 1 |
---|
81 | #include <atlas/cblas.h> |
---|
82 | #elif defined(HAVE_OPENBLAS_CBLAS_H) |
---|
83 | #include <openblas/cblas.h> |
---|
84 | #elif defined(HAVE_CBLAS_H) |
---|
85 | #include <cblas.h> |
---|
86 | #elif !defined(HAVE_ACCELERATE) |
---|
87 | #error "HAVE_BLAS was defined, but no blas header was found" |
---|
88 | #endif /* end of cblas includes */ |
---|
89 | #endif |
---|
90 | |
---|
91 | #if defined(HAVE_ACCELERATE) |
---|
92 | // include accelerate framework after blas |
---|
93 | #define HAVE_ATLAS 1 |
---|
94 | #define HAVE_BLAS 1 |
---|
95 | #include <Accelerate/Accelerate.h> |
---|
96 | |
---|
97 | #ifndef HAVE_AUBIO_DOUBLE |
---|
98 | #define aubio_vDSP_mmov vDSP_mmov |
---|
99 | #define aubio_vDSP_vmul vDSP_vmul |
---|
100 | #define aubio_vDSP_vsmul vDSP_vsmul |
---|
101 | #define aubio_vDSP_vsadd vDSP_vsadd |
---|
102 | #define aubio_vDSP_vfill vDSP_vfill |
---|
103 | #define aubio_vDSP_meanv vDSP_meanv |
---|
104 | #define aubio_vDSP_sve vDSP_sve |
---|
105 | #define aubio_vDSP_maxv vDSP_maxv |
---|
106 | #define aubio_vDSP_maxvi vDSP_maxvi |
---|
107 | #define aubio_vDSP_minv vDSP_minv |
---|
108 | #define aubio_vDSP_minvi vDSP_minvi |
---|
109 | #define aubio_vDSP_dotpr vDSP_dotpr |
---|
110 | #define aubio_vDSP_vclr vDSP_vclr |
---|
111 | #else /* HAVE_AUBIO_DOUBLE */ |
---|
112 | #define aubio_vDSP_mmov vDSP_mmovD |
---|
113 | #define aubio_vDSP_vmul vDSP_vmulD |
---|
114 | #define aubio_vDSP_vsmul vDSP_vsmulD |
---|
115 | #define aubio_vDSP_vsadd vDSP_vsaddD |
---|
116 | #define aubio_vDSP_vfill vDSP_vfillD |
---|
117 | #define aubio_vDSP_meanv vDSP_meanvD |
---|
118 | #define aubio_vDSP_sve vDSP_sveD |
---|
119 | #define aubio_vDSP_maxv vDSP_maxvD |
---|
120 | #define aubio_vDSP_maxvi vDSP_maxviD |
---|
121 | #define aubio_vDSP_minv vDSP_minvD |
---|
122 | #define aubio_vDSP_minvi vDSP_minviD |
---|
123 | #define aubio_vDSP_dotpr vDSP_dotprD |
---|
124 | #define aubio_vDSP_vclr vDSP_vclrD |
---|
125 | #endif /* HAVE_AUBIO_DOUBLE */ |
---|
126 | #endif /* HAVE_ACCELERATE */ |
---|
127 | |
---|
128 | #if defined(HAVE_BLAS) |
---|
129 | #ifndef HAVE_AUBIO_DOUBLE |
---|
130 | #ifdef HAVE_ATLAS |
---|
131 | #define aubio_catlas_set catlas_sset |
---|
132 | #endif /* HAVE_ATLAS */ |
---|
133 | #define aubio_cblas_copy cblas_scopy |
---|
134 | #define aubio_cblas_swap cblas_sswap |
---|
135 | #define aubio_cblas_dot cblas_sdot |
---|
136 | #define aubio_cblas__gemv cblas_sgemv |
---|
137 | #else /* HAVE_AUBIO_DOUBLE */ |
---|
138 | #ifdef HAVE_ATLAS |
---|
139 | #define aubio_catlas_set catlas_dset |
---|
140 | #endif /* HAVE_ATLAS */ |
---|
141 | #define aubio_cblas_copy cblas_dcopy |
---|
142 | #define aubio_cblas_swap cblas_dswap |
---|
143 | #define aubio_cblas_dot cblas_ddot |
---|
144 | #define aubio_cblas__gemv cblas_dgemv |
---|
145 | #endif /* HAVE_AUBIO_DOUBLE */ |
---|
146 | #endif /* HAVE_BLAS */ |
---|
147 | |
---|
148 | #if defined HAVE_INTEL_IPP |
---|
149 | #include <ippcore.h> |
---|
150 | #include <ippvm.h> |
---|
151 | #include <ipps.h> |
---|
152 | #ifndef HAVE_AUBIO_DOUBLE |
---|
153 | #define aubio_ippsSet ippsSet_32f |
---|
154 | #define aubio_ippsZero ippsZero_32f |
---|
155 | #define aubio_ippsCopy ippsCopy_32f |
---|
156 | #define aubio_ippsMul ippsMul_32f |
---|
157 | #define aubio_ippsMulC ippsMulC_32f |
---|
158 | #define aubio_ippsAddC ippsAddC_32f |
---|
159 | #define aubio_ippsLn ippsLn_32f_A21 |
---|
160 | #define aubio_ippsMean(a,b,c) ippsMean_32f(a, b, c, ippAlgHintFast) |
---|
161 | #define aubio_ippsSum(a,b,c) ippsSum_32f(a, b, c, ippAlgHintFast) |
---|
162 | #define aubio_ippsMax ippsMax_32f |
---|
163 | #define aubio_ippsMin ippsMin_32f |
---|
164 | #else /* HAVE_AUBIO_DOUBLE */ |
---|
165 | #define aubio_ippsSet ippsSet_64f |
---|
166 | #define aubio_ippsZero ippsZero_64f |
---|
167 | #define aubio_ippsCopy ippsCopy_64f |
---|
168 | #define aubio_ippsMul ippsMul_64f |
---|
169 | #define aubio_ippsMulC ippsMulC_64f |
---|
170 | #define aubio_ippsAddC ippsAddC_64f |
---|
171 | #define aubio_ippsLn ippsLn_64f_A26 |
---|
172 | #define aubio_ippsMean ippsMean_64f |
---|
173 | #define aubio_ippsSum ippsSum_64f |
---|
174 | #define aubio_ippsMax ippsMax_64f |
---|
175 | #define aubio_ippsMin ippsMin_64f |
---|
176 | #endif /* HAVE_AUBIO_DOUBLE */ |
---|
177 | #endif |
---|
178 | |
---|
179 | #if !defined(HAVE_MEMCPY_HACKS) && !defined(HAVE_ACCELERATE) && !defined(HAVE_ATLAS) && !defined(HAVE_INTEL_IPP) |
---|
180 | #define HAVE_NOOPT 1 |
---|
181 | #endif |
---|
182 | |
---|
183 | #include "types.h" |
---|
184 | |
---|
185 | #define AUBIO_UNSTABLE 1 |
---|
186 | |
---|
187 | #include "mathutils.h" |
---|
188 | |
---|
189 | /**** |
---|
190 | * |
---|
191 | * SYSTEM INTERFACE |
---|
192 | * |
---|
193 | */ |
---|
194 | |
---|
195 | /* Memory management */ |
---|
196 | #define AUBIO_MALLOC(_n) malloc(_n) |
---|
197 | #define AUBIO_REALLOC(_p,_n) realloc(_p,_n) |
---|
198 | #define AUBIO_NEW(_t) (_t*)calloc(sizeof(_t), 1) |
---|
199 | #define AUBIO_ARRAY(_t,_n) (_t*)calloc((_n)*sizeof(_t), 1) |
---|
200 | #define AUBIO_MEMCPY(_dst,_src,_n) memcpy(_dst,_src,_n) |
---|
201 | #define AUBIO_MEMSET(_dst,_src,_t) memset(_dst,_src,_t) |
---|
202 | #define AUBIO_FREE(_p) free(_p) |
---|
203 | |
---|
204 | |
---|
205 | /* file interface */ |
---|
206 | #define AUBIO_FOPEN(_f,_m) fopen(_f,_m) |
---|
207 | #define AUBIO_FCLOSE(_f) fclose(_f) |
---|
208 | #define AUBIO_FREAD(_p,_s,_n,_f) fread(_p,_s,_n,_f) |
---|
209 | #define AUBIO_FSEEK(_f,_n,_set) fseek(_f,_n,_set) |
---|
210 | |
---|
211 | /* strings */ |
---|
212 | #define AUBIO_STRLEN(_s) strlen(_s) |
---|
213 | #define AUBIO_STRCMP(_s,_t) strcmp(_s,_t) |
---|
214 | #define AUBIO_STRNCMP(_s,_t,_n) strncmp(_s,_t,_n) |
---|
215 | #define AUBIO_STRCPY(_dst,_src) strcpy(_dst,_src) |
---|
216 | #define AUBIO_STRCHR(_s,_c) strchr(_s,_c) |
---|
217 | #ifdef strdup |
---|
218 | #define AUBIO_STRDUP(s) strdup(s) |
---|
219 | #else |
---|
220 | #define AUBIO_STRDUP(s) AUBIO_STRCPY(AUBIO_MALLOC(AUBIO_STRLEN(s) + 1), s) |
---|
221 | #endif |
---|
222 | |
---|
223 | |
---|
224 | /* Error reporting */ |
---|
225 | typedef enum { |
---|
226 | AUBIO_OK = 0, |
---|
227 | AUBIO_FAIL = 1 |
---|
228 | } aubio_status; |
---|
229 | |
---|
230 | /* Logging */ |
---|
231 | |
---|
232 | #include "utils/log.h" |
---|
233 | |
---|
234 | /** internal logging function, defined in utils/log.c */ |
---|
235 | uint_t aubio_log(sint_t level, const char_t *fmt, ...); |
---|
236 | |
---|
237 | #ifdef HAVE_C99_VARARGS_MACROS |
---|
238 | #define AUBIO_ERR(...) aubio_log(AUBIO_LOG_ERR, "AUBIO ERROR: " __VA_ARGS__) |
---|
239 | #define AUBIO_INF(...) aubio_log(AUBIO_LOG_INF, "AUBIO INFO: " __VA_ARGS__) |
---|
240 | #define AUBIO_MSG(...) aubio_log(AUBIO_LOG_MSG, __VA_ARGS__) |
---|
241 | #define _AUBIO_DBG(...) aubio_log(AUBIO_LOG_DBG, __VA_ARGS__) |
---|
242 | #define AUBIO_WRN(...) aubio_log(AUBIO_LOG_WRN, "AUBIO WARNING: " __VA_ARGS__) |
---|
243 | #else |
---|
244 | #define AUBIO_ERR(format, args...) aubio_log(AUBIO_LOG_ERR, "AUBIO ERROR: " format , ##args) |
---|
245 | #define AUBIO_INF(format, args...) aubio_log(AUBIO_LOG_INF, "AUBIO INFO: " format , ##args) |
---|
246 | #define AUBIO_MSG(format, args...) aubio_log(AUBIO_LOG_MSG, format , ##args) |
---|
247 | #define _AUBIO_DBG(format, args...) aubio_log(AUBIO_LOG_DBG, format , ##args) |
---|
248 | #define AUBIO_WRN(format, args...) aubio_log(AUBIO_LOG_WRN, "AUBIO WARNING: " format, ##args) |
---|
249 | #endif |
---|
250 | |
---|
251 | #ifdef DEBUG |
---|
252 | #define AUBIO_DBG _AUBIO_DBG |
---|
253 | #else |
---|
254 | // disable debug output |
---|
255 | #ifdef HAVE_C99_VARARGS_MACROS |
---|
256 | #define AUBIO_DBG(...) {} |
---|
257 | #else |
---|
258 | #define AUBIO_DBG(format, args...) {} |
---|
259 | #endif |
---|
260 | #endif |
---|
261 | |
---|
262 | #define AUBIO_ERROR AUBIO_ERR |
---|
263 | |
---|
264 | #define AUBIO_QUIT(_s) exit(_s) |
---|
265 | #define AUBIO_SPRINTF sprintf |
---|
266 | |
---|
267 | #define AUBIO_MAX_SAMPLERATE (192000*8) |
---|
268 | #define AUBIO_MAX_CHANNELS 1024 |
---|
269 | |
---|
270 | /* pi and 2*pi */ |
---|
271 | #ifndef M_PI |
---|
272 | #define PI (3.14159265358979323846) |
---|
273 | #else |
---|
274 | #define PI (M_PI) |
---|
275 | #endif |
---|
276 | #define TWO_PI (PI*2.) |
---|
277 | |
---|
278 | #ifndef PATH_MAX |
---|
279 | #define PATH_MAX 1024 |
---|
280 | #endif |
---|
281 | |
---|
282 | /* aliases to math.h functions */ |
---|
283 | #if !HAVE_AUBIO_DOUBLE |
---|
284 | #define EXP expf |
---|
285 | #define COS cosf |
---|
286 | #define SIN sinf |
---|
287 | #define ABS fabsf |
---|
288 | #define POW powf |
---|
289 | #define SQRT sqrtf |
---|
290 | #define LOG10 log10f |
---|
291 | #define LOG logf |
---|
292 | #define FLOOR floorf |
---|
293 | #define CEIL ceilf |
---|
294 | #define ATAN atanf |
---|
295 | #define ATAN2 atan2f |
---|
296 | #else |
---|
297 | #define EXP exp |
---|
298 | #define COS cos |
---|
299 | #define SIN sin |
---|
300 | #define ABS fabs |
---|
301 | #define POW pow |
---|
302 | #define SQRT sqrt |
---|
303 | #define LOG10 log10 |
---|
304 | #define LOG log |
---|
305 | #define FLOOR floor |
---|
306 | #define CEIL ceil |
---|
307 | #define ATAN atan |
---|
308 | #define ATAN2 atan2 |
---|
309 | #endif |
---|
310 | #define ROUND(x) FLOOR(x+.5) |
---|
311 | |
---|
312 | /* aliases to complex.h functions */ |
---|
313 | #if HAVE_AUBIO_DOUBLE || !defined(HAVE_COMPLEX_H) || defined(WIN32) |
---|
314 | /* mingw32 does not know about c*f functions */ |
---|
315 | #define EXPC cexp |
---|
316 | /** complex = CEXPC(complex) */ |
---|
317 | #define CEXPC cexp |
---|
318 | /** sample = ARGC(complex) */ |
---|
319 | #define ARGC carg |
---|
320 | /** sample = ABSC(complex) norm */ |
---|
321 | #define ABSC cabs |
---|
322 | /** sample = REAL(complex) */ |
---|
323 | #define REAL creal |
---|
324 | /** sample = IMAG(complex) */ |
---|
325 | #define IMAG cimag |
---|
326 | #else |
---|
327 | /** sample = EXPC(complex) */ |
---|
328 | #define EXPC cexpf |
---|
329 | /** complex = CEXPC(complex) */ |
---|
330 | #define CEXPC cexp |
---|
331 | /** sample = ARGC(complex) */ |
---|
332 | #define ARGC cargf |
---|
333 | /** sample = ABSC(complex) norm */ |
---|
334 | #define ABSC cabsf |
---|
335 | /** sample = REAL(complex) */ |
---|
336 | #define REAL crealf |
---|
337 | /** sample = IMAG(complex) */ |
---|
338 | #define IMAG cimagf |
---|
339 | #endif |
---|
340 | |
---|
341 | /* avoid unresolved symbol with msvc 9 */ |
---|
342 | #if defined(_MSC_VER) && (_MSC_VER < 1900) |
---|
343 | #define isnan _isnan |
---|
344 | #endif |
---|
345 | |
---|
346 | #if !defined(_MSC_VER) |
---|
347 | #define AUBIO_STRERROR(errno,buf,len) strerror_r(errno, buf, len) |
---|
348 | #else |
---|
349 | #define AUBIO_STRERROR(errno,buf,len) strerror_s(buf, len, errno) |
---|
350 | #endif |
---|
351 | |
---|
352 | #ifdef HAVE_C99_VARARGS_MACROS |
---|
353 | #define AUBIO_STRERR(...) \ |
---|
354 | char errorstr[256]; \ |
---|
355 | AUBIO_STRERROR(errno, errorstr, sizeof(errorstr)); \ |
---|
356 | AUBIO_ERR(__VA_ARGS__) |
---|
357 | #else |
---|
358 | #define AUBIO_STRERR(format, args...) \ |
---|
359 | char errorstr[256]; \ |
---|
360 | AUBIO_STRERROR(errno, errorstr, sizeof(errorstr)); \ |
---|
361 | AUBIO_ERR(format, ##args) |
---|
362 | #endif |
---|
363 | |
---|
364 | /* handy shortcuts */ |
---|
365 | #define DB2LIN(g) (POW(10.0,(g)*0.05f)) |
---|
366 | #define LIN2DB(v) (20.0*LOG10(v)) |
---|
367 | #define SQR(_a) ((_a)*(_a)) |
---|
368 | |
---|
369 | #ifndef MAX |
---|
370 | #define MAX(a,b) (((a)>(b))?(a):(b)) |
---|
371 | #endif /* MAX */ |
---|
372 | #ifndef MIN |
---|
373 | #define MIN(a,b) (((a)<(b))?(a):(b)) |
---|
374 | #endif /* MIN */ |
---|
375 | |
---|
376 | #define ELEM_SWAP(a,b) { register smpl_t t=(a);(a)=(b);(b)=t; } |
---|
377 | |
---|
378 | #define VERY_SMALL_NUMBER 2.e-42 //1.e-37 |
---|
379 | |
---|
380 | /** if ABS(f) < VERY_SMALL_NUMBER, returns 1, else 0 */ |
---|
381 | #define IS_DENORMAL(f) ABS(f) < VERY_SMALL_NUMBER |
---|
382 | |
---|
383 | /** if ABS(f) < VERY_SMALL_NUMBER, returns 0., else f */ |
---|
384 | #define KILL_DENORMAL(f) IS_DENORMAL(f) ? 0. : f |
---|
385 | |
---|
386 | /** if f > VERY_SMALL_NUMBER, returns f, else returns VERY_SMALL_NUMBER */ |
---|
387 | #define CEIL_DENORMAL(f) f < VERY_SMALL_NUMBER ? VERY_SMALL_NUMBER : f |
---|
388 | |
---|
389 | #define SAFE_LOG10(f) LOG10(CEIL_DENORMAL(f)) |
---|
390 | #define SAFE_LOG(f) LOG(CEIL_DENORMAL(f)) |
---|
391 | |
---|
392 | /** silence unused parameter warning by adding an attribute */ |
---|
393 | #if defined(__GNUC__) |
---|
394 | #define UNUSED __attribute__((unused)) |
---|
395 | #else |
---|
396 | #define UNUSED |
---|
397 | #endif |
---|
398 | |
---|
399 | /* are we using gcc -std=c99 ? */ |
---|
400 | #if defined(__STRICT_ANSI__) |
---|
401 | #define strnlen(a,b) MIN(strlen(a),b) |
---|
402 | #if !HAVE_AUBIO_DOUBLE |
---|
403 | #define floorf floor |
---|
404 | #endif |
---|
405 | #endif /* __STRICT_ANSI__ */ |
---|
406 | |
---|
407 | #if defined(DEBUG) |
---|
408 | #include <assert.h> |
---|
409 | #define AUBIO_ASSERT(x) assert(x) |
---|
410 | #else |
---|
411 | #define AUBIO_ASSERT(x) |
---|
412 | #endif /* DEBUG */ |
---|
413 | |
---|
414 | // goto to failure if condition x is not true |
---|
415 | #define AUBIO_GOTO_FAILURE(x) if (!(x)) goto failure |
---|
416 | |
---|
417 | #define AUBIO_ASSERT_EQUAL_SHAPE(t1, t2) { \ |
---|
418 | AUBIO_ASSERT(t1 && t2); \ |
---|
419 | AUBIO_ASSERT(t1->ndim == t2->ndim); \ |
---|
420 | uint_t nn; \ |
---|
421 | for (nn = 0; nn < t1->ndim; nn++) \ |
---|
422 | AUBIO_ASSERT(t1->shape[nn] == t2->shape[nn]); \ |
---|
423 | } |
---|
424 | |
---|
425 | #endif /* AUBIO_PRIV_H */ |
---|