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 #include <stdio.h>
6 #include <stdlib.h>
7 #include <time.h>
8
9 #include "dsp_test_util.h"
10 #include "dsp_util.h"
11 #include "drc.h"
12 #include "raw.h"
13
14 #ifndef min
15 #define min(a, b) \
16 ({ \
17 __typeof__(a) _a = (a); \
18 __typeof__(b) _b = (b); \
19 _a < _b ? _a : _b; \
20 })
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 drc * drc,float * buf,size_t frames)29 static void process(struct drc *drc, float *buf, size_t frames)
30 {
31 struct timespec tp1, tp2;
32 int start;
33 float *data[2];
34 int chunk;
35 clock_gettime(CLOCK_THREAD_CPUTIME_ID, &tp1);
36 for (start = 0; start < frames;) {
37 data[0] = buf + start;
38 data[1] = buf + start + frames;
39 chunk = min(DRC_PROCESS_MAX_FRAMES, frames - start);
40 drc_process(drc, data, chunk);
41 start += chunk;
42 }
43 clock_gettime(CLOCK_THREAD_CPUTIME_ID, &tp2);
44 printf("drc processing takes %g seconds for %zu samples\n",
45 tp_diff(&tp2, &tp1), frames * 2);
46 }
47
main(int argc,char ** argv)48 int main(int argc, char **argv)
49 {
50 double NQ = 44100 / 2; /* nyquist frequency */
51 struct drc *drc;
52 size_t frames;
53 float *buf;
54
55 if (argc != 3) {
56 printf("Usage: drc_test input.raw output.raw\n");
57 return 1;
58 }
59
60 dsp_enable_flush_denormal_to_zero();
61 dsp_util_clear_fp_exceptions();
62 drc = drc_new(44100);
63
64 drc->emphasis_disabled = 0;
65 drc_set_param(drc, 0, PARAM_CROSSOVER_LOWER_FREQ, 0);
66 drc_set_param(drc, 0, PARAM_ENABLED, 1);
67 drc_set_param(drc, 0, PARAM_THRESHOLD, -29);
68 drc_set_param(drc, 0, PARAM_KNEE, 3);
69 drc_set_param(drc, 0, PARAM_RATIO, 6.677);
70 drc_set_param(drc, 0, PARAM_ATTACK, 0.02);
71 drc_set_param(drc, 0, PARAM_RELEASE, 0.2);
72 drc_set_param(drc, 0, PARAM_POST_GAIN, -7);
73
74 drc_set_param(drc, 1, PARAM_CROSSOVER_LOWER_FREQ, 200 / NQ);
75 drc_set_param(drc, 1, PARAM_ENABLED, 1);
76 drc_set_param(drc, 1, PARAM_THRESHOLD, -32);
77 drc_set_param(drc, 1, PARAM_KNEE, 23);
78 drc_set_param(drc, 1, PARAM_RATIO, 12);
79 drc_set_param(drc, 1, PARAM_ATTACK, 0.02);
80 drc_set_param(drc, 1, PARAM_RELEASE, 0.2);
81 drc_set_param(drc, 1, PARAM_POST_GAIN, 0.7);
82
83 drc_set_param(drc, 2, PARAM_CROSSOVER_LOWER_FREQ, 1200 / NQ);
84 drc_set_param(drc, 2, PARAM_ENABLED, 1);
85 drc_set_param(drc, 2, PARAM_THRESHOLD, -24);
86 drc_set_param(drc, 2, PARAM_KNEE, 30);
87 drc_set_param(drc, 2, PARAM_RATIO, 1);
88 drc_set_param(drc, 2, PARAM_ATTACK, 0.001);
89 drc_set_param(drc, 2, PARAM_RELEASE, 1);
90 drc_set_param(drc, 2, PARAM_POST_GAIN, 0);
91
92 drc_init(drc);
93 buf = read_raw(argv[1], &frames);
94 process(drc, buf, frames);
95 write_raw(argv[2], buf, frames);
96 drc_free(drc);
97 free(buf);
98 dsp_util_print_fp_exceptions();
99 return 0;
100 }
101