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