1 /* Copyright (c) 2013 The Chromium OS Authors. All rights reserved.
2 * Use of this source code is governed by a BSD-style license that can be
3 * found in the LICENSE file.
4 */
5
6 #include <stdio.h>
7 #include <stdlib.h>
8 #include <assert.h>
9 #include <time.h>
10 #include <memory.h>
11
12 #include "crossover.h"
13 #include "dsp_test_util.h"
14 #include "dsp_util.h"
15 #include "raw.h"
16
17 #ifndef min
18 #define min(a, b) ({ __typeof__(a) _a = (a); \
19 __typeof__(b) _b = (b); \
20 _a < _b ? _a : _b; })
21 #endif
22
tp_diff(struct timespec * tp2,struct timespec * tp1)23 static double tp_diff(struct timespec *tp2, struct timespec *tp1)
24 {
25 return (tp2->tv_sec - tp1->tv_sec)
26 + (tp2->tv_nsec - tp1->tv_nsec) * 1e-9;
27 }
28
process(struct crossover * xo,int count,float * data0,float * data1,float * data2)29 void process(struct crossover *xo, int count, float *data0, float *data1,
30 float *data2)
31 {
32 int start;
33 for (start = 0; start < count; start += 2048)
34 crossover_process(xo, min(2048, count - start),
35 data0 + start, data1 + start, data2 + start);
36 }
37
main(int argc,char ** argv)38 int main(int argc, char **argv)
39 {
40 size_t frames;
41 float *data0, *data1, *data2;
42 double NQ = 44100 / 2;
43 struct timespec tp1, tp2;
44 struct crossover xo;
45
46 if (argc != 3 && argc != 6) {
47 printf("Usage: crossover_test input.raw output.raw "
48 "[low.raw mid.raw high.raw]\n");
49 return 1;
50 }
51
52 dsp_enable_flush_denormal_to_zero();
53 dsp_util_clear_fp_exceptions();
54
55 data0 = read_raw(argv[1], &frames);
56 data1 = (float *)malloc(sizeof(float) * frames * 2);
57 data2 = (float *)malloc(sizeof(float) * frames * 2);
58
59 crossover_init(&xo, 400 / NQ, 4000 / NQ);
60 clock_gettime(CLOCK_THREAD_CPUTIME_ID, &tp1);
61 process(&xo, frames, data0, data1, data2);
62 clock_gettime(CLOCK_THREAD_CPUTIME_ID, &tp2);
63 printf("processing takes %g seconds for %zu samples\n",
64 tp_diff(&tp2, &tp1), frames);
65
66 crossover_init(&xo, 400 / NQ, 4000 / NQ);
67 clock_gettime(CLOCK_THREAD_CPUTIME_ID, &tp1);
68 process(&xo, frames, data0 + frames, data1 + frames,
69 data2 + frames);
70 clock_gettime(CLOCK_THREAD_CPUTIME_ID, &tp2);
71 printf("processing takes %g seconds for %zu samples\n",
72 tp_diff(&tp2, &tp1), frames);
73
74 if (argc == 6) {
75 write_raw(argv[3], data0, frames);
76 write_raw(argv[4], data1, frames);
77 write_raw(argv[5], data2, frames);
78 }
79
80 int i;
81 for (i = 0; i < frames * 2; i++)
82 data0[i] += data1[i] + data2[i];
83 write_raw(argv[2], data0, frames);
84
85 free(data0);
86 free(data1);
87 free(data2);
88
89 dsp_util_print_fp_exceptions();
90 return 0;
91 }
92