80 #define PFILTER(name, type, sin, cos, cc) \ 81 static void pfilter_channel_## name(AVFilterContext *ctx, \ 83 AVFrame *in, AVFrame *out) \ 85 AFreqShift *s = ctx->priv; \ 86 const int nb_samples = in->nb_samples; \ 87 const type *src = (const type *)in->extended_data[ch]; \ 88 type *dst = (type *)out->extended_data[ch]; \ 89 type *i1 = (type *)s->i1->extended_data[ch]; \ 90 type *o1 = (type *)s->o1->extended_data[ch]; \ 91 type *i2 = (type *)s->i2->extended_data[ch]; \ 92 type *o2 = (type *)s->o2->extended_data[ch]; \ 93 const type *c = s->cc; \ 94 const type level = s->level; \ 95 type shift = s->shift * M_PI; \ 96 type cos_theta = cos(shift); \ 97 type sin_theta = sin(shift); \ 99 for (int n = 0; n < nb_samples; n++) { \ 100 type xn1 = src[n], xn2 = src[n]; \ 103 for (int j = 0; j < NB_COEFS / 2; j++) { \ 104 I = c[j] * (xn1 + o2[j]) - i2[j]; \ 112 for (int j = NB_COEFS / 2; j < NB_COEFS; j++) { \ 113 Q = c[j] * (xn2 + o2[j]) - i2[j]; \ 120 Q = o2[NB_COEFS - 1]; \ 122 dst[n] = (I * cos_theta - Q * sin_theta) * level; \ 129 #define FFILTER(name, type, sin, cos, fmod, cc) \ 130 static void ffilter_channel_## name(AVFilterContext *ctx, \ 132 AVFrame *in, AVFrame *out) \ 134 AFreqShift *s = ctx->priv; \ 135 const int nb_samples = in->nb_samples; \ 136 const type *src = (const type *)in->extended_data[ch]; \ 137 type *dst = (type *)out->extended_data[ch]; \ 138 type *i1 = (type *)s->i1->extended_data[ch]; \ 139 type *o1 = (type *)s->o1->extended_data[ch]; \ 140 type *i2 = (type *)s->i2->extended_data[ch]; \ 141 type *o2 = (type *)s->o2->extended_data[ch]; \ 142 const type *c = s->cc; \ 143 const type level = s->level; \ 144 type ts = 1. / in->sample_rate; \ 145 type shift = s->shift; \ 146 int64_t N = s->in_samples; \ 148 for (int n = 0; n < nb_samples; n++) { \ 149 type xn1 = src[n], xn2 = src[n]; \ 152 for (int j = 0; j < NB_COEFS / 2; j++) { \ 153 I = c[j] * (xn1 + o2[j]) - i2[j]; \ 161 for (int j = NB_COEFS / 2; j < NB_COEFS; j++) { \ 162 Q = c[j] * (xn2 + o2[j]) - i2[j]; \ 169 Q = o2[NB_COEFS - 1]; \ 171 theta = 2. * M_PI * fmod(shift * (N + n) * ts, 1.); \ 172 dst[n] = (I * cos(theta) - Q * sin(theta)) * level; \ 181 double kksqrt, e, e2, e4, k, q;
183 k = tan((1. - transition * 2.) *
M_PI / 4.);
185 kksqrt = pow(1 - k * k, 0.25);
186 e = 0.5 * (1. - kksqrt) / (1. + kksqrt);
189 q = e * (1. + e4 * (2. + e4 * (15. + 150. * e4)));
195 static double ipowp(
double x, int64_t n)
217 q_ii1 =
ipowp(q, i * (i + 1));
218 q_ii1 *= sin((i * 2 + 1) * c *
M_PI / order) * j;
223 }
while (
fabs(q_ii1) > 1e-100);
236 q_i2 =
ipowp(q, i * i);
237 q_i2 *= cos(i * 2 * c *
M_PI / order) * j;
242 }
while (
fabs(q_i2) > 1e-100);
249 const int c = index + 1;
252 const double ww = num / den;
253 const double wwsq = ww * ww;
255 const double x = sqrt((1 - wwsq * k) * (1 - wwsq / k)) / (1 + wwsq);
256 const double coef = (1 - x) / (1 + x);
261 static void compute_coefs(
double *coef_arrd,
float *coef_arrf,
int nbr_coefs,
double transition)
263 const int order = nbr_coefs * 2 + 1;
268 for (
int n = 0; n < nbr_coefs; n++) {
269 const int idx = (n / 2) + (n & 1) * nbr_coefs / 2;
272 coef_arrf[idx] = coef_arrd[idx];
287 if (!s->
i1 || !s->
o1 || !s->
i2 || !s->
o2)
315 const int start = (in->
channels * jobnr) / nb_jobs;
316 const int end = (in->
channels * (jobnr+1)) / nb_jobs;
318 for (
int ch = start; ch < end; ch++)
364 #define OFFSET(x) offsetof(AFreqShift, x) 365 #define FLAGS AV_OPT_FLAG_AUDIO_PARAM|AV_OPT_FLAG_FILTERING_PARAM|AV_OPT_FLAG_RUNTIME_PARAM 394 .
name =
"afreqshift",
398 .priv_class = &afreqshift_class,
416 .
name =
"aphaseshift",
420 .priv_class = &aphaseshift_class,
static int query_formats(AVFilterContext *ctx)
This structure describes decoded (raw) audio or video data.
AVFILTER_DEFINE_CLASS(afreqshift)
Main libavfilter public API header.
AVFilter ff_af_afreqshift
static double compute_acc_num(double q, int order, int c)
static const AVOption aphaseshift_options[]
static const AVOption afreqshift_options[]
#define AVFILTER_FLAG_SUPPORT_TIMELINE_GENERIC
Some filters support a generic "enable" expression option that can be used to enable or disable a fil...
const char * name
Pad name.
int ff_filter_frame(AVFilterLink *link, AVFrame *frame)
Send a frame of data to the next filter.
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
#define PFILTER(name, type, sin, cos, cc)
static const AVFilterPad outputs[]
static double ipowp(double x, int64_t n)
A filter pad used for either input or output.
A link between two filters.
int sample_rate
samples per second
AVFrame * ff_get_audio_buffer(AVFilterLink *link, int nb_samples)
Request an audio samples buffer with a specific set of permissions.
static double compute_acc_den(double q, int order, int c)
static __device__ float fabs(float a)
void av_frame_free(AVFrame **frame)
Free the frame and any dynamically allocated objects in it, e.g.
#define NULL_IF_CONFIG_SMALL(x)
Return NULL if CONFIG_SMALL is true, otherwise the argument without modification. ...
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...
void * priv
private data for use by the filter
#define AVFILTER_FLAG_SLICE_THREADS
The filter supports multithreading by splitting frames into multiple parts and processing them concur...
static double compute_coef(int index, double k, double q, int order)
#define FFILTER(name, type, sin, cos, fmod, cc)
static int config_input(AVFilterLink *inlink)
static void compute_transition_param(double *K, double *Q, double transition)
int channels
number of audio channels, only used for audio.
audio channel layout utility functions
int ff_filter_get_nb_threads(AVFilterContext *ctx)
Get number of threads for current filter instance.
static int process_command(AVFilterContext *ctx, const char *cmd, const char *args, char *res, int res_len, int flags)
typedef void(APIENTRY *FF_PFNGLACTIVETEXTUREPROC)(GLenum texture)
static const AVFilterPad inputs[]
static int filter_channels(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
static av_cold void uninit(AVFilterContext *ctx)
int format
agreed upon media format
A list of supported channel layouts.
void(* filter_channel)(AVFilterContext *ctx, int channel, AVFrame *in, AVFrame *out)
AVSampleFormat
Audio sample formats.
int av_frame_is_writable(AVFrame *frame)
Check if the frame data is writable.
Used for passing data between threads.
static void compute_coefs(double *coef_arrd, float *coef_arrf, int nbr_coefs, double transition)
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)\n", 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.
const char * name
Filter name.
AVFilterLink ** outputs
array of pointers to output links
enum MovChannelLayoutTag * layouts
#define flags(name, subs,...)
AVFilterInternal * internal
An opaque struct for libavfilter internal use.
static int filter_frame(AVFilterLink *inlink, AVFrame *in)
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
internal math functions header
channel
Use these values when setting the channel map with ebur128_set_channel().
int channels
Number of channels.
avfilter_execute_func * execute
AVFilterContext * dst
dest filter
static enum AVSampleFormat sample_fmts[]
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
int nb_samples
number of audio samples (per channel) described by this frame
const AVFilter * filter
the AVFilter of which this is an instance
int av_frame_copy_props(AVFrame *dst, const AVFrame *src)
Copy only "metadata" fields from src to dst.
AVFilter ff_af_aphaseshift