FFmpeg  4.3.9
af_crossfeed.c
Go to the documentation of this file.
1 /*
2  * This file is part of FFmpeg.
3  *
4  * FFmpeg is free software; you can redistribute it and/or
5  * modify it under the terms of the GNU Lesser General Public
6  * License as published by the Free Software Foundation; either
7  * version 2.1 of the License, or (at your option) any later version.
8  *
9  * FFmpeg is distributed in the hope that it will be useful,
10  * but WITHOUT ANY WARRANTY; without even the implied warranty of
11  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
12  * Lesser General Public License for more details.
13  *
14  * You should have received a copy of the GNU Lesser General Public
15  * License along with FFmpeg; if not, write to the Free Software
16  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
17  */
18 
20 #include "libavutil/ffmath.h"
21 #include "libavutil/opt.h"
22 #include "avfilter.h"
23 #include "audio.h"
24 #include "formats.h"
25 
26 typedef struct CrossfeedContext {
27  const AVClass *class;
28 
29  double range;
30  double strength;
31  double slope;
32  double level_in;
33  double level_out;
34 
35  double a0, a1, a2;
36  double b0, b1, b2;
37 
38  double i1, i2;
39  double o1, o2;
41 
43 {
46  int ret;
47 
48  if ((ret = ff_add_format (&formats, AV_SAMPLE_FMT_DBL )) < 0 ||
49  (ret = ff_set_common_formats (ctx , formats )) < 0 ||
50  (ret = ff_add_channel_layout (&layout , AV_CH_LAYOUT_STEREO)) < 0 ||
51  (ret = ff_set_common_channel_layouts (ctx , layout )) < 0 ||
52  (ret = ff_set_common_samplerates (ctx , ff_all_samplerates())) < 0)
53  return ret;
54 
55  return 0;
56 }
57 
58 static int config_input(AVFilterLink *inlink)
59 {
60  AVFilterContext *ctx = inlink->dst;
61  CrossfeedContext *s = ctx->priv;
62  double A = ff_exp10(s->strength * -30 / 40);
63  double w0 = 2 * M_PI * (1. - s->range) * 2100 / inlink->sample_rate;
64  double alpha;
65 
66  alpha = sin(w0) / 2 * sqrt((A + 1 / A) * (1 / s->slope - 1) + 2);
67 
68  s->a0 = (A + 1) + (A - 1) * cos(w0) + 2 * sqrt(A) * alpha;
69  s->a1 = -2 * ((A - 1) + (A + 1) * cos(w0));
70  s->a2 = (A + 1) + (A - 1) * cos(w0) - 2 * sqrt(A) * alpha;
71  s->b0 = A * ((A + 1) - (A - 1) * cos(w0) + 2 * sqrt(A) * alpha);
72  s->b1 = 2 * A * ((A - 1) - (A + 1) * cos(w0));
73  s->b2 = A * ((A + 1) - (A - 1) * cos(w0) - 2 * sqrt(A) * alpha);
74 
75  s->a1 /= s->a0;
76  s->a2 /= s->a0;
77  s->b0 /= s->a0;
78  s->b1 /= s->a0;
79  s->b2 /= s->a0;
80 
81  return 0;
82 }
83 
84 static int filter_frame(AVFilterLink *inlink, AVFrame *in)
85 {
86  AVFilterContext *ctx = inlink->dst;
87  AVFilterLink *outlink = ctx->outputs[0];
88  CrossfeedContext *s = ctx->priv;
89  const double *src = (const double *)in->data[0];
90  const double level_in = s->level_in;
91  const double level_out = s->level_out;
92  const double b0 = s->b0;
93  const double b1 = s->b1;
94  const double b2 = s->b2;
95  const double a1 = s->a1;
96  const double a2 = s->a2;
97  AVFrame *out;
98  double *dst;
99  int n;
100 
101  if (av_frame_is_writable(in)) {
102  out = in;
103  } else {
104  out = ff_get_audio_buffer(outlink, in->nb_samples);
105  if (!out) {
106  av_frame_free(&in);
107  return AVERROR(ENOMEM);
108  }
110  }
111  dst = (double *)out->data[0];
112 
113  for (n = 0; n < out->nb_samples; n++, src += 2, dst += 2) {
114  double mid = (src[0] + src[1]) * level_in * .5;
115  double side = (src[0] - src[1]) * level_in * .5;
116  double oside = side * b0 + s->i1 * b1 + s->i2 * b2 - s->o1 * a1 - s->o2 * a2;
117 
118  s->i2 = s->i1;
119  s->i1 = side;
120  s->o2 = s->o1;
121  s->o1 = oside;
122 
123  if (ctx->is_disabled) {
124  dst[0] = src[0];
125  dst[1] = src[1];
126  } else {
127  dst[0] = (mid + oside) * level_out;
128  dst[1] = (mid - oside) * level_out;
129  }
130  }
131 
132  if (out != in)
133  av_frame_free(&in);
134  return ff_filter_frame(outlink, out);
135 }
136 
137 static int process_command(AVFilterContext *ctx, const char *cmd, const char *args,
138  char *res, int res_len, int flags)
139 {
140  int ret;
141 
142  ret = ff_filter_process_command(ctx, cmd, args, res, res_len, flags);
143  if (ret < 0)
144  return ret;
145 
146  return config_input(ctx->inputs[0]);
147 }
148 
149 #define OFFSET(x) offsetof(CrossfeedContext, x)
150 #define FLAGS AV_OPT_FLAG_AUDIO_PARAM|AV_OPT_FLAG_FILTERING_PARAM|AV_OPT_FLAG_RUNTIME_PARAM
151 
152 static const AVOption crossfeed_options[] = {
153  { "strength", "set crossfeed strength", OFFSET(strength), AV_OPT_TYPE_DOUBLE, {.dbl=.2}, 0, 1, FLAGS },
154  { "range", "set soundstage wideness", OFFSET(range), AV_OPT_TYPE_DOUBLE, {.dbl=.5}, 0, 1, FLAGS },
155  { "slope", "set curve slope", OFFSET(slope), AV_OPT_TYPE_DOUBLE, {.dbl=.5}, .01, 1, FLAGS },
156  { "level_in", "set level in", OFFSET(level_in), AV_OPT_TYPE_DOUBLE, {.dbl=.9}, 0, 1, FLAGS },
157  { "level_out", "set level out", OFFSET(level_out), AV_OPT_TYPE_DOUBLE, {.dbl=1.}, 0, 1, FLAGS },
158  { NULL }
159 };
160 
161 AVFILTER_DEFINE_CLASS(crossfeed);
162 
163 static const AVFilterPad inputs[] = {
164  {
165  .name = "default",
166  .type = AVMEDIA_TYPE_AUDIO,
167  .filter_frame = filter_frame,
168  .config_props = config_input,
169  },
170  { NULL }
171 };
172 
173 static const AVFilterPad outputs[] = {
174  {
175  .name = "default",
176  .type = AVMEDIA_TYPE_AUDIO,
177  },
178  { NULL }
179 };
180 
182  .name = "crossfeed",
183  .description = NULL_IF_CONFIG_SMALL("Apply headphone crossfeed filter."),
184  .query_formats = query_formats,
185  .priv_size = sizeof(CrossfeedContext),
186  .priv_class = &crossfeed_class,
187  .inputs = inputs,
188  .outputs = outputs,
191 };
#define NULL
Definition: coverity.c:32
int ff_set_common_channel_layouts(AVFilterContext *ctx, AVFilterChannelLayouts *layouts)
A helper for query_formats() which sets all links to the same list of channel layouts/sample rates...
Definition: formats.c:586
This structure describes decoded (raw) audio or video data.
Definition: frame.h:300
AVOption.
Definition: opt.h:246
static int config_input(AVFilterLink *inlink)
Definition: af_crossfeed.c:58
Main libavfilter public API header.
static const AVOption crossfeed_options[]
Definition: af_crossfeed.c:152
#define FLAGS
Definition: af_crossfeed.c:150
#define AV_CH_LAYOUT_STEREO
int is_disabled
the enabled state from the last expression evaluation
Definition: avfilter.h:385
const char * name
Pad name.
Definition: internal.h:60
AVFilterLink ** inputs
array of pointers to input links
Definition: avfilter.h:346
int ff_filter_frame(AVFilterLink *link, AVFrame *frame)
Send a frame of data to the next filter.
Definition: avfilter.c:1075
static int filter_frame(AVFilterLink *inlink, AVFrame *in)
Definition: af_crossfeed.c:84
AVOptions.
static const AVFilterPad inputs[]
Definition: af_crossfeed.c:163
#define A(x)
Definition: vp56_arith.h:28
static const AVFilterPad outputs[]
Definition: af_crossfeed.c:173
A filter pad used for either input or output.
Definition: internal.h:54
#define src
Definition: vp8dsp.c:254
static av_always_inline double ff_exp10(double x)
Compute 10^x for floating point values.
Definition: ffmath.h:42
int ff_set_common_formats(AVFilterContext *ctx, AVFilterFormats *formats)
A helper for query_formats() which sets all links to the same list of formats.
Definition: formats.c:605
int ff_add_channel_layout(AVFilterChannelLayouts **l, uint64_t channel_layout)
Definition: formats.c:356
AVFrame * ff_get_audio_buffer(AVFilterLink *link, int nb_samples)
Request an audio samples buffer with a specific set of permissions.
Definition: audio.c:86
#define AVERROR(e)
Definition: error.h:43
void av_frame_free(AVFrame **frame)
Free the frame and any dynamically allocated objects in it, e.g.
Definition: frame.c:203
#define NULL_IF_CONFIG_SMALL(x)
Return NULL if CONFIG_SMALL is true, otherwise the argument without modification. ...
Definition: internal.h:188
int ff_filter_process_command(AVFilterContext *ctx, const char *cmd, const char *arg, char *res, int res_len, int flags)
Generic processing of user supplied commands that are set in the same way as the filter options...
Definition: avfilter.c:869
void * priv
private data for use by the filter
Definition: avfilter.h:353
AVFILTER_DEFINE_CLASS(crossfeed)
int ff_add_format(AVFilterFormats **avff, int64_t fmt)
Add fmt to the list of media formats contained in *avff.
Definition: formats.c:350
audio channel layout utility functions
AVFormatContext * ctx
Definition: movenc.c:48
#define s(width, name)
Definition: cbs_vp9.c:257
A list of supported channel layouts.
Definition: formats.h:85
int av_frame_is_writable(AVFrame *frame)
Check if the frame data is writable.
Definition: frame.c:595
static const int16_t alpha[]
Definition: ilbcdata.h:55
static int query_formats(AVFilterContext *ctx)
Definition: af_crossfeed.c:42
uint8_t pi<< 24) CONV_FUNC_GROUP(AV_SAMPLE_FMT_FLT, float, AV_SAMPLE_FMT_U8, uint8_t,(*(const uint8_t *) pi - 0x80) *(1.0f/(1<< 7))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_DBL, double, AV_SAMPLE_FMT_U8, uint8_t,(*(const uint8_t *) pi - 0x80) *(1.0/(1<< 7))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_S16, int16_t,(*(const int16_t *) pi >> 8)+0x80) CONV_FUNC_GROUP(AV_SAMPLE_FMT_FLT, float, AV_SAMPLE_FMT_S16, int16_t, *(const int16_t *) pi *(1.0f/(1<< 15))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_DBL, double, AV_SAMPLE_FMT_S16, int16_t, *(const int16_t *) pi *(1.0/(1<< 15))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_S32, int32_t,(*(const int32_t *) pi >> 24)+0x80) CONV_FUNC_GROUP(AV_SAMPLE_FMT_FLT, float, AV_SAMPLE_FMT_S32, int32_t, *(const int32_t *) pi *(1.0f/(1U<< 31))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_DBL, double, AV_SAMPLE_FMT_S32, int32_t, *(const int32_t *) pi *(1.0/(1U<< 31))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_FLT, float, av_clip_uint8(lrintf(*(const float *) pi *(1<< 7))+0x80)) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S16, int16_t, AV_SAMPLE_FMT_FLT, float, av_clip_int16(lrintf(*(const float *) pi *(1<< 15)))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S32, int32_t, AV_SAMPLE_FMT_FLT, float, av_clipl_int32(llrintf(*(const float *) pi *(1U<< 31)))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_DBL, double, av_clip_uint8(lrint(*(const double *) pi *(1<< 7))+0x80)) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S16, int16_t, AV_SAMPLE_FMT_DBL, double, av_clip_int16(lrint(*(const double *) pi *(1<< 15)))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S32, int32_t, AV_SAMPLE_FMT_DBL, double, av_clipl_int32(llrint(*(const double *) pi *(1U<< 31)))) #define SET_CONV_FUNC_GROUP(ofmt, ifmt) static void set_generic_function(AudioConvert *ac) { } void ff_audio_convert_free(AudioConvert **ac) { if(! *ac) return;ff_dither_free(&(*ac) ->dc);av_freep(ac);} AudioConvert *ff_audio_convert_alloc(AVAudioResampleContext *avr, enum AVSampleFormat out_fmt, enum AVSampleFormat in_fmt, int channels, int sample_rate, int apply_map) { AudioConvert *ac;int in_planar, out_planar;ac=av_mallocz(sizeof(*ac));if(!ac) return NULL;ac->avr=avr;ac->out_fmt=out_fmt;ac->in_fmt=in_fmt;ac->channels=channels;ac->apply_map=apply_map;if(avr->dither_method !=AV_RESAMPLE_DITHER_NONE &&av_get_packed_sample_fmt(out_fmt)==AV_SAMPLE_FMT_S16 &&av_get_bytes_per_sample(in_fmt) > 2) { ac->dc=ff_dither_alloc(avr, out_fmt, in_fmt, channels, sample_rate, apply_map);if(!ac->dc) { av_free(ac);return NULL;} return ac;} in_planar=ff_sample_fmt_is_planar(in_fmt, channels);out_planar=ff_sample_fmt_is_planar(out_fmt, channels);if(in_planar==out_planar) { ac->func_type=CONV_FUNC_TYPE_FLAT;ac->planes=in_planar ? ac->channels :1;} else if(in_planar) ac->func_type=CONV_FUNC_TYPE_INTERLEAVE;else ac->func_type=CONV_FUNC_TYPE_DEINTERLEAVE;set_generic_function(ac);if(ARCH_AARCH64) ff_audio_convert_init_aarch64(ac);if(ARCH_ARM) ff_audio_convert_init_arm(ac);if(ARCH_X86) ff_audio_convert_init_x86(ac);return ac;} int ff_audio_convert(AudioConvert *ac, AudioData *out, AudioData *in) { int use_generic=1;int len=in->nb_samples;int p;if(ac->dc) { av_log(ac->avr, AV_LOG_TRACE, "%d samples - audio_convert: %s to %s (dithered)\", len, av_get_sample_fmt_name(ac->in_fmt), av_get_sample_fmt_name(ac->out_fmt));return ff_convert_dither(ac-> in
Describe the class of an AVClass context structure.
Definition: log.h:67
Filter definition.
Definition: avfilter.h:144
const char * name
Filter name.
Definition: avfilter.h:148
#define AVFILTER_FLAG_SUPPORT_TIMELINE_INTERNAL
Same as AVFILTER_FLAG_SUPPORT_TIMELINE_GENERIC, except that the filter will have its filter_frame() c...
Definition: avfilter.h:133
AVFilterLink ** outputs
array of pointers to output links
Definition: avfilter.h:350
AVFilterFormats * ff_all_samplerates(void)
Definition: formats.c:439
#define OFFSET(x)
Definition: af_crossfeed.c:149
#define flags(name, subs,...)
Definition: cbs_av1.c:576
uint8_t * data[AV_NUM_DATA_POINTERS]
pointer to the picture/channel planes.
Definition: frame.h:314
internal math functions header
if(ret< 0)
Definition: vf_mcdeint.c:279
AVFilter ff_af_crossfeed
Definition: af_crossfeed.c:181
static int process_command(AVFilterContext *ctx, const char *cmd, const char *args, char *res, int res_len, int flags)
Definition: af_crossfeed.c:137
A list of supported formats for one end of a filter link.
Definition: formats.h:64
uint64_t layout
An instance of a filter.
Definition: avfilter.h:338
FILE * out
Definition: movenc.c:54
#define M_PI
Definition: mathematics.h:52
formats
Definition: signature.h:48
int nb_samples
number of audio samples (per channel) described by this frame
Definition: frame.h:366
for(j=16;j >0;--j)
int ff_set_common_samplerates(AVFilterContext *ctx, AVFilterFormats *samplerates)
Definition: formats.c:593
int av_frame_copy_props(AVFrame *dst, const AVFrame *src)
Copy only "metadata" fields from src to dst.
Definition: frame.c:659