Go to the documentation of this file.
40 #define MAX_BANDS MAX_SPLITS + 1
84 #define OFFSET(x) offsetof(AudioCrossoverContext, x)
85 #define AF AV_OPT_FLAG_AUDIO_PARAM | AV_OPT_FLAG_FILTERING_PARAM
110 char *p, *
arg, *saveptr =
NULL;
130 if (
c[0] ==
'd' &&
c[1] ==
'B')
145 char *p, *
arg, *saveptr =
NULL;
170 if (
i > 0 && freq <= s->splits[
i-1]) {
184 for (
i = 0;
i <=
s->nb_splits;
i++) {
203 double omega = 2. *
M_PI *
fc / sr;
204 double cosine = cos(omega);
205 double alpha = sin(omega) / (2. * q);
207 double b0 = (1. - cosine) / 2.;
208 double b1 = 1. - cosine;
209 double b2 = (1. - cosine) / 2.;
211 double a1 = -2. * cosine;
229 double omega = 2. *
M_PI *
fc / sr;
230 double cosine = cos(omega);
231 double alpha = sin(omega) / (2. * q);
233 double b0 = (1. + cosine) / 2.;
234 double b1 = -1. - cosine;
235 double b2 = (1. + cosine) / 2.;
237 double a1 = -2. * cosine;
255 double omega = 2. *
M_PI *
fc / sr;
256 double cosine = cos(omega);
257 double alpha = sin(omega) / (2. * q);
260 double a1 = -2. * cosine;
281 double omega = 2. *
M_PI *
fc / sr;
298 double n = order / 2.;
300 for (
int i = 0;
i < n / 2;
i++)
301 q[
i] = 1. / (-2. * cos(
M_PI * (2. * (
i + 1) + n - 1.) / (2. * n)));
304 #define BIQUAD_PROCESS(name, type) \
305 static void biquad_process_## name(const type *const c, \
307 type *dst, const type *src, \
310 const type b0 = c[B0]; \
311 const type b1 = c[B1]; \
312 const type b2 = c[B2]; \
313 const type a1 = c[A1]; \
314 const type a2 = c[A2]; \
318 for (int n = 0; n + 1 < nb_samples; n++) { \
322 out = in * b0 + z1; \
323 z1 = b1 * in + z2 + a1 * out; \
324 z2 = b2 * in + a2 * out; \
329 out = in * b0 + z1; \
330 z1 = b1 * in + z2 + a1 * out; \
331 z2 = b2 * in + a2 * out; \
335 if (nb_samples & 1) { \
336 const int n = nb_samples - 1; \
337 const type in = src[n]; \
340 out = in * b0 + z1; \
341 z1 = b1 * in + z2 + a1 * out; \
342 z2 = b2 * in + a2 * out; \
353 #define XOVER_PROCESS(name, type, one, ff) \
354 static int filter_channels_## name(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs) \
356 AudioCrossoverContext *s = ctx->priv; \
357 AVFrame *in = s->input_frame; \
358 AVFrame **frames = s->frames; \
359 const int start = (in->channels * jobnr) / nb_jobs; \
360 const int end = (in->channels * (jobnr+1)) / nb_jobs; \
361 const int nb_samples = in->nb_samples; \
362 const int nb_outs = ctx->nb_outputs; \
363 const int first_order = s->first_order; \
365 for (int ch = start; ch < end; ch++) { \
366 const type *src = (const type *)in->extended_data[ch]; \
367 type *xover = (type *)s->xover->extended_data[ch]; \
369 s->fdsp->vector_## ff ##mul_scalar((type *)frames[0]->extended_data[ch], src, \
370 s->level_in, FFALIGN(nb_samples, sizeof(type))); \
372 for (int band = 0; band < nb_outs; band++) { \
373 for (int f = 0; band + 1 < nb_outs && f < s->filter_count; f++) { \
374 const type *prv = (const type *)frames[band]->extended_data[ch]; \
375 type *dst = (type *)frames[band + 1]->extended_data[ch]; \
376 const type *hsrc = f == 0 ? prv : dst; \
377 type *hp = xover + nb_outs * 20 + band * 20 + f * 2; \
378 const type *const hpc = (type *)&s->hp[band][f].c ## ff; \
380 biquad_process_## name(hpc, hp, dst, hsrc, nb_samples); \
383 for (int f = 0; band + 1 < nb_outs && f < s->filter_count; f++) { \
384 type *dst = (type *)frames[band]->extended_data[ch]; \
385 const type *lsrc = dst; \
386 type *lp = xover + band * 20 + f * 2; \
387 const type *const lpc = (type *)&s->lp[band][f].c ## ff; \
389 biquad_process_## name(lpc, lp, dst, lsrc, nb_samples); \
392 for (int aband = band + 1; aband + 1 < nb_outs; aband++) { \
394 const type *asrc = (const type *)frames[band]->extended_data[ch]; \
395 type *dst = (type *)frames[band]->extended_data[ch]; \
396 type *ap = xover + nb_outs * 40 + (aband * nb_outs + band) * 20; \
397 const type *const apc = (type *)&s->ap[aband][0].c ## ff; \
399 biquad_process_## name(apc, ap, dst, asrc, nb_samples); \
402 for (int f = first_order; f < s->ap_filter_count; f++) { \
403 const type *asrc = (const type *)frames[band]->extended_data[ch]; \
404 type *dst = (type *)frames[band]->extended_data[ch]; \
405 type *ap = xover + nb_outs * 40 + (aband * nb_outs + band) * 20 + f * 2;\
406 const type *const apc = (type *)&s->ap[aband][f].c ## ff; \
408 biquad_process_## name(apc, ap, dst, asrc, nb_samples); \
413 for (int band = 0; band < nb_outs; band++) { \
414 const type gain = s->gains[band] * ((band & 1 && first_order) ? -one : one); \
415 type *dst = (type *)frames[band]->extended_data[ch]; \
417 s->fdsp->vector_## ff ##mul_scalar(dst, dst, gain, \
418 FFALIGN(nb_samples, sizeof(type))); \
435 s->order = (
s->order_opt + 1) * 2;
436 s->filter_count =
s->order / 2;
437 s->first_order =
s->filter_count & 1;
438 s->ap_filter_count =
s->filter_count / 2 +
s->first_order;
441 for (
int band = 0; band <=
s->nb_splits; band++) {
442 if (
s->first_order) {
447 for (
int n =
s->first_order; n < s->filter_count; n++) {
448 const int idx =
s->filter_count / 2 - ((n +
s->first_order) / 2 -
s->first_order) - 1;
457 for (
int n =
s->first_order; n < s->ap_filter_count; n++) {
458 const int idx = (
s->filter_count / 2 - ((n * 2 +
s->first_order) / 2 -
s->first_order) - 1);
471 ctx->nb_outputs *
ctx->nb_outputs * 10));
485 for (
i = 0;
i <
ctx->nb_outputs;
i++) {
503 for (
i = 0;
i <
ctx->nb_outputs;
i++) {
511 for (
i = 0;
i <
ctx->nb_outputs;
i++)
514 s->input_frame =
NULL;
537 .
name =
"acrossover",
540 .priv_class = &acrossover_class,
AVFrame * ff_get_audio_buffer(AVFilterLink *link, int nb_samples)
Request an audio samples buffer with a specific set of permissions.
@ AV_SAMPLE_FMT_FLTP
float, planar
it s the only field you need to keep assuming you have a context There is some magic you don t need to care about around this just let it vf default minimum maximum flags name is the option name
Filter the word “frame” indicates either a video frame or a group of audio as stored in an AVFrame structure Format for each input and each output the list of supported formats For video that means pixel format For audio that means channel sample they are references to shared objects When the negotiation mechanism computes the intersection of the formats supported at each end of a all references to both lists are replaced with a reference to the intersection And when a single format is eventually chosen for a link amongst the remaining all references to the list are updated That means that if a filter requires that its input and output have the same format amongst a supported all it has to do is use a reference to the same list of formats query_formats can leave some formats unset and return AVERROR(EAGAIN) to cause the negotiation mechanism toagain later. That can be used by filters with complex requirements to use the format negotiated on one link to set the formats supported on another. Frame references ownership and permissions
BiquadCoeffs lp[MAX_BANDS][20]
int ff_filter_frame(AVFilterLink *link, AVFrame *frame)
Send a frame of data to the next filter.
The exact code depends on how similar the blocks are and how related they are to the and needs to apply these operations to the correct inlink or outlink if there are several Macros are available to factor that when no extra processing is inlink
char * av_asprintf(const char *fmt,...)
BiquadCoeffs hp[MAX_BANDS][20]
void av_frame_free(AVFrame **frame)
Free the frame and any dynamically allocated objects in it, e.g.
This structure describes decoded (raw) audio or video data.
int64_t pts
Presentation timestamp in time_base units (time when frame should be shown to user).
AVFILTER_DEFINE_CLASS(acrossover)
#define fc(width, name, range_min, range_max)
#define BIQUAD_PROCESS(name, type)
const char * name
Filter name.
AVFrame * frames[MAX_BANDS]
A link between two filters.
int(* filter_channels)(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
static double b1(void *priv, double x, double y)
static const AVOption acrossover_options[]
if it could not because there are no more frames
static int parse_gains(AVFilterContext *ctx)
static void set_hp(BiquadCoeffs *b, double fc, double q, double sr)
A filter pad used for either input or output.
#define AV_LOG_ERROR
Something went wrong and cannot losslessly be recovered.
char * av_strtok(char *s, const char *delim, char **saveptr)
Split the string into several tokens which can be accessed by successive calls to av_strtok().
const AVFilter ff_af_acrossover
#define FILTER_INPUTS(array)
int av_sscanf(const char *string, const char *format,...)
See libc sscanf manual for more information.
Describe the class of an AVClass context structure.
static int filter_frame(AVFilterLink *inlink, AVFrame *in)
Undefined Behavior In the C some operations are like signed integer dereferencing freed accessing outside allocated Undefined Behavior must not occur in a C it is not safe even if the output of undefined operations is unused The unsafety may seem nit picking but Optimizing compilers have in fact optimized code on the assumption that no undefined Behavior occurs Optimizing code based on wrong assumptions can and has in some cases lead to effects beyond the output of computations The signed integer overflow problem in speed critical code Code which is highly optimized and works with signed integers sometimes has the problem that often the output of the computation does not c
BiquadCoeffs ap[MAX_BANDS][20]
#define AVFILTER_FLAG_DYNAMIC_OUTPUTS
The number of the filter outputs is not determined just by AVFilter.outputs.
#define NULL_IF_CONFIG_SMALL(x)
Return NULL if CONFIG_SMALL is true, otherwise the argument without modification.
static av_cold void uninit(AVFilterContext *ctx)
static double b2(void *priv, double x, double y)
static void calc_q_factors(int order, double *q)
int nb_samples
number of audio samples (per channel) described by this frame
#define i(width, name, range_min, range_max)
int ff_filter_get_nb_threads(AVFilterContext *ctx)
Get number of threads for current filter instance.
const char * name
Pad name.
static void set_lp(BiquadCoeffs *b, double fc, double q, double sr)
static const AVFilterPad inputs[]
static void set_ap(BiquadCoeffs *b, double fc, double q, double sr)
enum AVMediaType type
AVFilterPad type.
static void set_ap1(BiquadCoeffs *b, double fc, double sr)
static av_cold int init(AVFilterContext *ctx)
@ AV_SAMPLE_FMT_DBLP
double, planar
#define AVFILTER_FLAG_SLICE_THREADS
The filter supports multithreading by splitting frames into multiple parts and processing them concur...
static int config_input(AVFilterLink *inlink)
static const int16_t alpha[]
int ff_append_outpad_free_name(AVFilterContext *f, AVFilterPad *p)
av_cold AVFloatDSPContext * avpriv_float_dsp_alloc(int bit_exact)
Allocate a float DSP context.
#define XOVER_PROCESS(name, type, one, ff)
#define AVERROR_BUG
Internal bug, also see AVERROR_BUG2.
static double b0(void *priv, double x, double y)
static av_always_inline int ff_filter_execute(AVFilterContext *ctx, avfilter_action_func *func, void *arg, int *ret, int nb_jobs)
#define FILTER_SAMPLEFMTS(...)