FFmpeg
af_anlms.c
Go to the documentation of this file.
1 /*
2  * Copyright (c) 2019 Paul B Mahol
3  *
4  * This file is part of FFmpeg.
5  *
6  * FFmpeg is free software; you can redistribute it and/or
7  * modify it under the terms of the GNU Lesser General Public
8  * License as published by the Free Software Foundation; either
9  * version 2.1 of the License, or (at your option) any later version.
10  *
11  * FFmpeg is distributed in the hope that it will be useful,
12  * but WITHOUT ANY WARRANTY; without even the implied warranty of
13  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
14  * Lesser General Public License for more details.
15  *
16  * You should have received a copy of the GNU Lesser General Public
17  * License along with FFmpeg; if not, write to the Free Software
18  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
19  */
20 
22 #include "libavutil/common.h"
23 #include "libavutil/float_dsp.h"
24 #include "libavutil/mem.h"
25 #include "libavutil/opt.h"
26 
27 #include "audio.h"
28 #include "avfilter.h"
29 #include "filters.h"
30 #include "formats.h"
31 
32 enum OutModes {
39 };
40 
41 typedef struct AudioNLMSContext {
42  const AVClass *class;
43 
44  int order;
45  float mu;
46  float eps;
47  float leakage;
49  int precision;
50 
56 
58 
59  int anlmf;
60 
61  int (*filter_channels)(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs);
62 
65 
66 #define OFFSET(x) offsetof(AudioNLMSContext, x)
67 #define A AV_OPT_FLAG_AUDIO_PARAM|AV_OPT_FLAG_FILTERING_PARAM
68 #define AT AV_OPT_FLAG_AUDIO_PARAM|AV_OPT_FLAG_FILTERING_PARAM|AV_OPT_FLAG_RUNTIME_PARAM
69 
70 static const AVOption anlms_options[] = {
71  { "order", "set the filter order", OFFSET(order), AV_OPT_TYPE_INT, {.i64=256}, 1, INT16_MAX, A },
72  { "mu", "set the filter mu", OFFSET(mu), AV_OPT_TYPE_FLOAT, {.dbl=0.75}, 0, 2, AT },
73  { "eps", "set the filter eps", OFFSET(eps), AV_OPT_TYPE_FLOAT, {.dbl=1}, 0, 1, AT },
74  { "leakage", "set the filter leakage", OFFSET(leakage), AV_OPT_TYPE_FLOAT, {.dbl=0}, 0, 1, AT },
75  { "out_mode", "set output mode", OFFSET(output_mode), AV_OPT_TYPE_INT, {.i64=OUT_MODE}, 0, NB_OMODES-1, AT, .unit = "mode" },
76  { "i", "input", 0, AV_OPT_TYPE_CONST, {.i64=IN_MODE}, 0, 0, AT, .unit = "mode" },
77  { "d", "desired", 0, AV_OPT_TYPE_CONST, {.i64=DESIRED_MODE}, 0, 0, AT, .unit = "mode" },
78  { "o", "output", 0, AV_OPT_TYPE_CONST, {.i64=OUT_MODE}, 0, 0, AT, .unit = "mode" },
79  { "n", "noise", 0, AV_OPT_TYPE_CONST, {.i64=NOISE_MODE}, 0, 0, AT, .unit = "mode" },
80  { "e", "error", 0, AV_OPT_TYPE_CONST, {.i64=ERROR_MODE}, 0, 0, AT, .unit = "mode" },
81  { "precision", "set processing precision", OFFSET(precision), AV_OPT_TYPE_INT, {.i64=0}, 0, 2, A, .unit = "precision" },
82  { "auto", "set auto processing precision", 0, AV_OPT_TYPE_CONST, {.i64=0}, 0, 0, A, .unit = "precision" },
83  { "float", "set single-floating point processing precision", 0, AV_OPT_TYPE_CONST, {.i64=1}, 0, 0, A, .unit = "precision" },
84  { "double","set double-floating point processing precision", 0, AV_OPT_TYPE_CONST, {.i64=2}, 0, 0, A, .unit = "precision" },
85  { NULL }
86 };
87 
88 AVFILTER_DEFINE_CLASS_EXT(anlms, "anlm(f|s)", anlms_options);
89 
90 static int query_formats(const AVFilterContext *ctx,
91  AVFilterFormatsConfig **cfg_in,
92  AVFilterFormatsConfig **cfg_out)
93 {
94  const AudioNLMSContext *s = ctx->priv;
95  static const enum AVSampleFormat sample_fmts[3][3] = {
99  };
100  int ret;
101 
102  if ((ret = ff_set_common_formats_from_list2(ctx, cfg_in, cfg_out,
103  sample_fmts[s->precision])) < 0)
104  return ret;
105 
106  return 0;
107 }
108 
110 {
111  AudioNLMSContext *s = ctx->priv;
112  int i, ret, status;
113  int nb_samples;
114  int64_t pts;
115 
117 
118  nb_samples = FFMIN(ff_inlink_queued_samples(ctx->inputs[0]),
119  ff_inlink_queued_samples(ctx->inputs[1]));
120  for (i = 0; i < ctx->nb_inputs && nb_samples > 0; i++) {
121  if (s->frame[i])
122  continue;
123 
124  if (ff_inlink_check_available_samples(ctx->inputs[i], nb_samples) > 0) {
125  ret = ff_inlink_consume_samples(ctx->inputs[i], nb_samples, nb_samples, &s->frame[i]);
126  if (ret < 0)
127  return ret;
128  }
129  }
130 
131  if (s->frame[0] && s->frame[1]) {
132  AVFrame *out;
133 
134  out = ff_get_audio_buffer(ctx->outputs[0], s->frame[0]->nb_samples);
135  if (!out) {
136  av_frame_free(&s->frame[0]);
137  av_frame_free(&s->frame[1]);
138  return AVERROR(ENOMEM);
139  }
140 
141  ff_filter_execute(ctx, s->filter_channels, out, NULL,
142  FFMIN(ctx->outputs[0]->ch_layout.nb_channels, ff_filter_get_nb_threads(ctx)));
143 
144  out->pts = s->frame[0]->pts;
145  out->duration = s->frame[0]->duration;
146 
147  av_frame_free(&s->frame[0]);
148  av_frame_free(&s->frame[1]);
149 
150  ret = ff_filter_frame(ctx->outputs[0], out);
151  if (ret < 0)
152  return ret;
153  }
154 
155  if (!nb_samples) {
156  for (i = 0; i < 2; i++) {
157  if (ff_inlink_acknowledge_status(ctx->inputs[i], &status, &pts)) {
158  ff_outlink_set_status(ctx->outputs[0], status, pts);
159  return 0;
160  }
161  }
162  }
163 
164  if (ff_outlink_frame_wanted(ctx->outputs[0])) {
165  for (i = 0; i < 2; i++) {
166  if (s->frame[i] || ff_inlink_queued_samples(ctx->inputs[i]) > 0)
167  continue;
168  ff_inlink_request_frame(ctx->inputs[i]);
169  return 0;
170  }
171  }
172  return 0;
173 }
174 
175 #define DEPTH 32
176 #include "anlms_template.c"
177 
178 #undef DEPTH
179 #define DEPTH 64
180 #include "anlms_template.c"
181 
182 static int config_output(AVFilterLink *outlink)
183 {
184  AVFilterContext *ctx = outlink->src;
185  AudioNLMSContext *s = ctx->priv;
186 
187  s->anlmf = !strcmp(ctx->filter->name, "anlmf");
188  s->kernel_size = FFALIGN(s->order, 16);
189 
190  if (!s->offset)
191  s->offset = ff_get_audio_buffer(outlink, 1);
192  if (!s->delay)
193  s->delay = ff_get_audio_buffer(outlink, 2 * s->kernel_size);
194  if (!s->coeffs)
195  s->coeffs = ff_get_audio_buffer(outlink, 2 * s->kernel_size);
196  if (!s->tmp)
197  s->tmp = ff_get_audio_buffer(outlink, s->kernel_size);
198  if (!s->delay || !s->coeffs || !s->offset || !s->tmp)
199  return AVERROR(ENOMEM);
200 
201  switch (outlink->format) {
202  case AV_SAMPLE_FMT_DBLP:
203  s->filter_channels = filter_channels_double;
204  break;
205  case AV_SAMPLE_FMT_FLTP:
206  s->filter_channels = filter_channels_float;
207  break;
208  }
209 
210  return 0;
211 }
212 
214 {
215  AudioNLMSContext *s = ctx->priv;
216 
217  s->fdsp = avpriv_float_dsp_alloc(0);
218  if (!s->fdsp)
219  return AVERROR(ENOMEM);
220 
221  return 0;
222 }
223 
225 {
226  AudioNLMSContext *s = ctx->priv;
227 
228  av_freep(&s->fdsp);
229  av_frame_free(&s->delay);
230  av_frame_free(&s->coeffs);
231  av_frame_free(&s->offset);
232  av_frame_free(&s->tmp);
233 }
234 
235 static const AVFilterPad inputs[] = {
236  {
237  .name = "input",
238  .type = AVMEDIA_TYPE_AUDIO,
239  },
240  {
241  .name = "desired",
242  .type = AVMEDIA_TYPE_AUDIO,
243  },
244 };
245 
246 static const AVFilterPad outputs[] = {
247  {
248  .name = "default",
249  .type = AVMEDIA_TYPE_AUDIO,
250  .config_props = config_output,
251  },
252 };
253 
255  .name = "anlms",
256  .description = NULL_IF_CONFIG_SMALL("Apply Normalized Least-Mean-Squares algorithm to first audio stream."),
257  .priv_size = sizeof(AudioNLMSContext),
258  .priv_class = &anlms_class,
259  .init = init,
260  .uninit = uninit,
261  .activate = activate,
267  .process_command = ff_filter_process_command,
268 };
269 
271  .name = "anlmf",
272  .description = NULL_IF_CONFIG_SMALL("Apply Normalized Least-Mean-Fourth algorithm to first audio stream."),
273  .priv_size = sizeof(AudioNLMSContext),
274  .priv_class = &anlms_class,
275  .init = init,
276  .uninit = uninit,
277  .activate = activate,
283  .process_command = ff_filter_process_command,
284 };
ff_get_audio_buffer
AVFrame * ff_get_audio_buffer(AVFilterLink *link, int nb_samples)
Request an audio samples buffer with a specific set of permissions.
Definition: audio.c:98
AV_SAMPLE_FMT_FLTP
@ AV_SAMPLE_FMT_FLTP
float, planar
Definition: samplefmt.h:66
AVERROR
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
opt.h
out
FILE * out
Definition: movenc.c:55
ff_filter_frame
int ff_filter_frame(AVFilterLink *link, AVFrame *frame)
Send a frame of data to the next filter.
Definition: avfilter.c:1061
sample_fmts
static enum AVSampleFormat sample_fmts[]
Definition: adpcmenc.c:948
anlms_template.c
int64_t
long long int64_t
Definition: coverity.c:34
av_frame_free
void av_frame_free(AVFrame **frame)
Free the frame and any dynamically allocated objects in it, e.g.
Definition: frame.c:162
OutModes
OutModes
Definition: af_aap.c:32
FILTER_INPUTS
#define FILTER_INPUTS(array)
Definition: filters.h:262
AVFrame
This structure describes decoded (raw) audio or video data.
Definition: frame.h:389
AVOption
AVOption.
Definition: opt.h:429
OUT_MODE
@ OUT_MODE
Definition: af_anlms.c:35
AudioNLMSContext::mu
float mu
Definition: af_anlms.c:45
init
static av_cold int init(AVFilterContext *ctx)
Definition: af_anlms.c:213
AVFilter::name
const char * name
Filter name.
Definition: avfilter.h:205
IN_MODE
@ IN_MODE
Definition: af_anlms.c:33
ff_af_anlms
const AVFilter ff_af_anlms
Definition: af_anlms.c:254
formats.h
FF_FILTER_FORWARD_STATUS_BACK_ALL
#define FF_FILTER_FORWARD_STATUS_BACK_ALL(outlink, filter)
Forward the status on an output link to all input links.
Definition: filters.h:447
AudioNLMSContext::fdsp
AVFloatDSPContext * fdsp
Definition: af_anlms.c:63
anlms_options
static const AVOption anlms_options[]
Definition: af_anlms.c:70
pts
static int64_t pts
Definition: transcode_aac.c:644
AVFilterPad
A filter pad used for either input or output.
Definition: filters.h:38
OFFSET
#define OFFSET(x)
Definition: af_anlms.c:66
ff_inlink_check_available_samples
int ff_inlink_check_available_samples(AVFilterLink *link, unsigned min)
Test if enough samples are available on the link.
Definition: avfilter.c:1471
av_cold
#define av_cold
Definition: attributes.h:90
AudioNLMSContext::offset
AVFrame * offset
Definition: af_anlms.c:52
DESIRED_MODE
@ DESIRED_MODE
Definition: af_anlms.c:34
ff_outlink_set_status
static void ff_outlink_set_status(AVFilterLink *link, int status, int64_t pts)
Set the status field of a link from the source filter.
Definition: filters.h:424
ff_inlink_request_frame
void ff_inlink_request_frame(AVFilterLink *link)
Mark that a frame is wanted on the link.
Definition: avfilter.c:1593
AudioNLMSContext::tmp
AVFrame * tmp
Definition: af_anlms.c:55
s
#define s(width, name)
Definition: cbs_vp9.c:198
AudioNLMSContext::frame
AVFrame * frame[2]
Definition: af_anlms.c:57
AudioNLMSContext::kernel_size
int kernel_size
Definition: af_anlms.c:51
AVMEDIA_TYPE_AUDIO
@ AVMEDIA_TYPE_AUDIO
Definition: avutil.h:202
filters.h
ctx
AVFormatContext * ctx
Definition: movenc.c:49
AudioNLMSContext::leakage
float leakage
Definition: af_anlms.c:47
FILTER_OUTPUTS
#define FILTER_OUTPUTS(array)
Definition: filters.h:263
arg
const char * arg
Definition: jacosubdec.c:67
AVClass
Describe the class of an AVClass context structure.
Definition: log.h:75
ff_inlink_consume_samples
int ff_inlink_consume_samples(AVFilterLink *link, unsigned min, unsigned max, AVFrame **rframe)
Take samples from the link's FIFO and update the link's stats.
Definition: avfilter.c:1510
NULL
#define NULL
Definition: coverity.c:32
uninit
static av_cold void uninit(AVFilterContext *ctx)
Definition: af_anlms.c:224
AudioNLMSContext::anlmf
int anlmf
Definition: af_anlms.c:59
AudioNLMSContext::eps
float eps
Definition: af_anlms.c:46
ff_inlink_acknowledge_status
int ff_inlink_acknowledge_status(AVFilterLink *link, int *rstatus, int64_t *rpts)
Test and acknowledge the change of status on the link.
Definition: avfilter.c:1437
AVFilterFormatsConfig
Lists of formats / etc.
Definition: avfilter.h:111
float_dsp.h
ff_af_anlmf
const AVFilter ff_af_anlmf
Definition: af_anlms.c:270
FILTER_SINGLE_SAMPLEFMT
#define FILTER_SINGLE_SAMPLEFMT(sample_fmt_)
Definition: filters.h:255
NULL_IF_CONFIG_SMALL
#define NULL_IF_CONFIG_SMALL(x)
Return NULL if CONFIG_SMALL is true, otherwise the argument without modification.
Definition: internal.h:94
AV_SAMPLE_FMT_NONE
@ AV_SAMPLE_FMT_NONE
Definition: samplefmt.h:56
query_formats
static int query_formats(const AVFilterContext *ctx, AVFilterFormatsConfig **cfg_in, AVFilterFormatsConfig **cfg_out)
Definition: af_anlms.c:90
AVFloatDSPContext
Definition: float_dsp.h:24
ff_filter_process_command
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:900
AudioNLMSContext::coeffs
AVFrame * coeffs
Definition: af_anlms.c:54
NB_OMODES
@ NB_OMODES
Definition: af_anlms.c:38
AV_OPT_TYPE_FLOAT
@ AV_OPT_TYPE_FLOAT
Underlying C type is float.
Definition: opt.h:271
i
#define i(width, name, range_min, range_max)
Definition: cbs_h2645.c:256
AudioNLMSContext::order
int order
Definition: af_anlms.c:44
activate
static int activate(AVFilterContext *ctx)
Definition: af_anlms.c:109
common.h
ff_filter_get_nb_threads
int ff_filter_get_nb_threads(AVFilterContext *ctx)
Get number of threads for current filter instance.
Definition: avfilter.c:840
AVSampleFormat
AVSampleFormat
Audio sample formats.
Definition: samplefmt.h:55
FILTER_QUERY_FUNC2
#define FILTER_QUERY_FUNC2(func)
Definition: filters.h:239
FFMIN
#define FFMIN(a, b)
Definition: macros.h:49
AudioNLMSContext::filter_channels
int(* filter_channels)(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
Definition: af_anlms.c:61
AudioNLMSContext::output_mode
int output_mode
Definition: af_anlms.c:48
AVFilterPad::name
const char * name
Pad name.
Definition: filters.h:44
ff_inlink_queued_samples
int ff_inlink_queued_samples(AVFilterLink *link)
Definition: avfilter.c:1465
AVFilter
Filter definition.
Definition: avfilter.h:201
ret
ret
Definition: filter_design.txt:187
ERROR_MODE
@ ERROR_MODE
Definition: af_anlms.c:37
NOISE_MODE
@ NOISE_MODE
Definition: af_anlms.c:36
AudioNLMSContext
Definition: af_anlms.c:41
AudioNLMSContext::delay
AVFrame * delay
Definition: af_anlms.c:53
status
ov_status_e status
Definition: dnn_backend_openvino.c:100
ff_set_common_formats_from_list2
int ff_set_common_formats_from_list2(const AVFilterContext *ctx, AVFilterFormatsConfig **cfg_in, AVFilterFormatsConfig **cfg_out, const int *fmts)
Definition: formats.c:1016
channel_layout.h
ff_filter_execute
int ff_filter_execute(AVFilterContext *ctx, avfilter_action_func *func, void *arg, int *ret, int nb_jobs)
Definition: avfilter.c:1666
inputs
static const AVFilterPad inputs[]
Definition: af_anlms.c:235
AV_OPT_TYPE_INT
@ AV_OPT_TYPE_INT
Underlying C type is int.
Definition: opt.h:259
avfilter.h
A
#define A
Definition: af_anlms.c:67
AudioNLMSContext::precision
int precision
Definition: af_anlms.c:49
AV_SAMPLE_FMT_DBLP
@ AV_SAMPLE_FMT_DBLP
double, planar
Definition: samplefmt.h:67
AVFILTER_DEFINE_CLASS_EXT
AVFILTER_DEFINE_CLASS_EXT(anlms, "anlm(f|s)", anlms_options)
config_output
static int config_output(AVFilterLink *outlink)
Definition: af_anlms.c:182
AVFilterContext
An instance of a filter.
Definition: avfilter.h:457
AVFILTER_FLAG_SLICE_THREADS
#define AVFILTER_FLAG_SLICE_THREADS
The filter supports multithreading by splitting frames into multiple parts and processing them concur...
Definition: avfilter.h:152
mem.h
audio.h
FFALIGN
#define FFALIGN(x, a)
Definition: macros.h:78
av_freep
#define av_freep(p)
Definition: tableprint_vlc.h:34
avpriv_float_dsp_alloc
av_cold AVFloatDSPContext * avpriv_float_dsp_alloc(int bit_exact)
Allocate a float DSP context.
Definition: float_dsp.c:146
AVFILTER_FLAG_SUPPORT_TIMELINE_INTERNAL
#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:190
outputs
static const AVFilterPad outputs[]
Definition: af_anlms.c:246
ff_outlink_frame_wanted
the definition of that something depends on the semantic of the filter The callback must examine the status of the filter s links and proceed accordingly The status of output links is stored in the status_in and status_out fields and tested by the ff_outlink_frame_wanted() function. If this function returns true
AT
#define AT
Definition: af_anlms.c:68
AV_OPT_TYPE_CONST
@ AV_OPT_TYPE_CONST
Special option type for declaring named constants.
Definition: opt.h:299