[deb_ffmpeg.git] / ffmpeg / libavfilter / vf_noise.c

/*
 * Copyright (c) 2002 Michael Niedermayer <michaelni@gmx.at>
 * Copyright (c) 2013 Paul B Mahol
 *
 * This file is part of FFmpeg.
 *
 * FFmpeg is free software; you can redistribute it and/or
 * modify it under the terms of the GNU Lesser General Public
 * License as published by the Free Software Foundation; either
 * version 2.1 of the License, or (at your option) any later version.
 *
 * FFmpeg is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 * Lesser General Public License for more details.
 *
 * You should have received a copy of the GNU Lesser General Public
 * License along with FFmpeg; if not, write to the Free Software
 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
 */

/**
 * @file
 * noise generator
 */

#include "libavutil/opt.h"
#include "libavutil/imgutils.h"
#include "libavutil/lfg.h"
#include "libavutil/parseutils.h"
#include "libavutil/pixdesc.h"
#include "avfilter.h"
#include "formats.h"
#include "internal.h"
#include "vf_noise.h"
#include "video.h"

typedef struct ThreadData {
    AVFrame *in, *out;
} ThreadData;

#define OFFSET(x) offsetof(NoiseContext, x)
#define FLAGS AV_OPT_FLAG_FILTERING_PARAM|AV_OPT_FLAG_VIDEO_PARAM

#define NOISE_PARAMS(name, x, param)                                                                                             \
    {#name"_seed", "set component #"#x" noise seed", OFFSET(param.seed), AV_OPT_TYPE_INT, {.i64=-1}, -1, INT_MAX, FLAGS},        \
    {#name"_strength", "set component #"#x" strength", OFFSET(param.strength), AV_OPT_TYPE_INT, {.i64=0}, 0, 100, FLAGS},        \
    {#name"s",         "set component #"#x" strength", OFFSET(param.strength), AV_OPT_TYPE_INT, {.i64=0}, 0, 100, FLAGS},        \
    {#name"_flags", "set component #"#x" flags", OFFSET(param.flags), AV_OPT_TYPE_FLAGS, {.i64=0}, 0, 31, FLAGS, #name"_flags"}, \
    {#name"f", "set component #"#x" flags", OFFSET(param.flags), AV_OPT_TYPE_FLAGS, {.i64=0}, 0, 31, FLAGS, #name"_flags"},      \
    {"a", "averaged noise", 0, AV_OPT_TYPE_CONST, {.i64=NOISE_AVERAGED}, 0, 0, FLAGS, #name"_flags"},                            \
    {"p", "(semi)regular pattern", 0, AV_OPT_TYPE_CONST, {.i64=NOISE_PATTERN},  0, 0, FLAGS, #name"_flags"},                     \
    {"t", "temporal noise", 0, AV_OPT_TYPE_CONST, {.i64=NOISE_TEMPORAL}, 0, 0, FLAGS, #name"_flags"},                            \
    {"u", "uniform noise",  0, AV_OPT_TYPE_CONST, {.i64=NOISE_UNIFORM},  0, 0, FLAGS, #name"_flags"},

static const AVOption noise_options[] = {
    NOISE_PARAMS(all, 0, all)
    NOISE_PARAMS(c0,  0, param[0])
    NOISE_PARAMS(c1,  1, param[1])
    NOISE_PARAMS(c2,  2, param[2])
    NOISE_PARAMS(c3,  3, param[3])
    {NULL}
};

AVFILTER_DEFINE_CLASS(noise);

static const int8_t patt[4] = { -1, 0, 1, 0 };

#define RAND_N(range) ((int) ((double) range * av_lfg_get(lfg) / (UINT_MAX + 1.0)))
static av_cold int init_noise(NoiseContext *n, int comp)
{
    int8_t *noise = av_malloc(MAX_NOISE * sizeof(int8_t));
    FilterParams *fp = &n->param[comp];
    AVLFG *lfg = &n->param[comp].lfg;
    int strength = fp->strength;
    int flags = fp->flags;
    int i, j;

    if (!noise)
        return AVERROR(ENOMEM);

    av_lfg_init(&fp->lfg, fp->seed + comp*31415U);

    for (i = 0, j = 0; i < MAX_NOISE; i++, j++) {
        if (flags & NOISE_UNIFORM) {
            if (flags & NOISE_AVERAGED) {
                if (flags & NOISE_PATTERN) {
                    noise[i] = (RAND_N(strength) - strength / 2) / 6
                        + patt[j % 4] * strength * 0.25 / 3;
                } else {
                    noise[i] = (RAND_N(strength) - strength / 2) / 3;
                }
            } else {
                if (flags & NOISE_PATTERN) {
                    noise[i] = (RAND_N(strength) - strength / 2) / 2
                        + patt[j % 4] * strength * 0.25;
                } else {
                    noise[i] = RAND_N(strength) - strength / 2;
                }
            }
        } else {
            double x1, x2, w, y1;
            do {
                x1 = 2.0 * av_lfg_get(lfg) / (float)UINT_MAX - 1.0;
                x2 = 2.0 * av_lfg_get(lfg) / (float)UINT_MAX - 1.0;
                w = x1 * x1 + x2 * x2;
            } while (w >= 1.0);

            w   = sqrt((-2.0 * log(w)) / w);
            y1  = x1 * w;
            y1 *= strength / sqrt(3.0);
            if (flags & NOISE_PATTERN) {
                y1 /= 2;
                y1 += patt[j % 4] * strength * 0.35;
            }
            y1 = av_clipf(y1, -128, 127);
            if (flags & NOISE_AVERAGED)
                y1 /= 3.0;
            noise[i] = (int)y1;
        }
        if (RAND_N(6) == 0)
            j--;
    }

    for (i = 0; i < MAX_RES; i++)
        for (j = 0; j < 3; j++)
            fp->prev_shift[i][j] = noise + (av_lfg_get(lfg) & (MAX_SHIFT - 1));

    fp->noise = noise;
    return 0;
}

static int query_formats(AVFilterContext *ctx)
{
    AVFilterFormats *formats = NULL;
    int fmt;

    for (fmt = 0; av_pix_fmt_desc_get(fmt); fmt++) {
        const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(fmt);
        if (desc->flags & AV_PIX_FMT_FLAG_PLANAR && !((desc->comp[0].depth_minus1 + 1) & 7))
            ff_add_format(&formats, fmt);
    }

    ff_set_common_formats(ctx, formats);
    return 0;
}

static int config_input(AVFilterLink *inlink)
{
    NoiseContext *n = inlink->dst->priv;
    const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(inlink->format);
    int ret;

    n->nb_planes = av_pix_fmt_count_planes(inlink->format);

    if ((ret = av_image_fill_linesizes(n->bytewidth, inlink->format, inlink->w)) < 0)
        return ret;

    n->height[1] = n->height[2] = FF_CEIL_RSHIFT(inlink->h, desc->log2_chroma_h);
    n->height[0] = n->height[3] = inlink->h;

    return 0;
}

void ff_line_noise_c(uint8_t *dst, const uint8_t *src, const int8_t *noise,
                     int len, int shift)
{
    int i;

    noise += shift;
    for (i = 0; i < len; i++) {
        int v = src[i] + noise[i];

        dst[i] = av_clip_uint8(v);
    }
}

void ff_line_noise_avg_c(uint8_t *dst, const uint8_t *src,
                         int len, const int8_t * const *shift)
{
    int i;
    const int8_t *src2 = (const int8_t*)src;

    for (i = 0; i < len; i++) {
        const int n = shift[0][i] + shift[1][i] + shift[2][i];
        dst[i] = src2[i] + ((n * src2[i]) >> 7);
    }
}

static void noise(uint8_t *dst, const uint8_t *src,
                  int dst_linesize, int src_linesize,
                  int width, int start, int end, NoiseContext *n, int comp)
{
    FilterParams *p = &n->param[comp];
    int8_t *noise = p->noise;
    const int flags = p->flags;
    int y;

    if (!noise) {
        if (dst != src)
            av_image_copy_plane(dst, dst_linesize, src, src_linesize, width, end - start);
        return;
    }

    for (y = start; y < end; y++) {
        const int ix = y & (MAX_RES - 1);
        int x;
        for (x=0; x < width; x+= MAX_RES) {
            int w = FFMIN(width - x, MAX_RES);
            int shift = p->rand_shift[ix];

            if (flags & NOISE_AVERAGED) {
                n->line_noise_avg(dst + x, src + x, w, (const int8_t**)p->prev_shift[ix]);
                p->prev_shift[ix][shift & 3] = noise + shift;
            } else {
                n->line_noise(dst + x, src + x, noise, w, shift);
            }
        }
        dst += dst_linesize;
        src += src_linesize;
    }
}

static int filter_slice(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
{
    NoiseContext *s = ctx->priv;
    ThreadData *td = arg;
    int plane;

    for (plane = 0; plane < s->nb_planes; plane++) {
        const int height = s->height[plane];
        const int start  = (height *  jobnr   ) / nb_jobs;
        const int end    = (height * (jobnr+1)) / nb_jobs;
        noise(td->out->data[plane] + start * td->out->linesize[plane],
              td->in->data[plane]  + start * td->in->linesize[plane],
              td->out->linesize[plane], td->in->linesize[plane],
              s->bytewidth[plane], start, end, s, plane);
    }
    return 0;
}

static int filter_frame(AVFilterLink *inlink, AVFrame *inpicref)
{
    AVFilterContext *ctx = inlink->dst;
    AVFilterLink *outlink = ctx->outputs[0];
    NoiseContext *n = ctx->priv;
    ThreadData td;
    AVFrame *out;
    int comp, i;

    if (av_frame_is_writable(inpicref)) {
        out = inpicref;
    } else {
        out = ff_get_video_buffer(outlink, outlink->w, outlink->h);
        if (!out) {
            av_frame_free(&inpicref);
            return AVERROR(ENOMEM);
        }
        av_frame_copy_props(out, inpicref);
    }

    for (comp = 0; comp < 4; comp++) {
        FilterParams *fp = &n->param[comp];

        if ((!fp->rand_shift_init || (fp->flags & NOISE_TEMPORAL)) && fp->strength) {

            for (i = 0; i < MAX_RES; i++) {
                fp->rand_shift[i] = av_lfg_get(&fp->lfg) & (MAX_SHIFT - 1);
            }
            fp->rand_shift_init = 1;
        }
    }

    td.in = inpicref; td.out = out;
    ctx->internal->execute(ctx, filter_slice, &td, NULL, FFMIN(n->height[0], ctx->graph->nb_threads));
    emms_c();

    if (inpicref != out)
        av_frame_free(&inpicref);
    return ff_filter_frame(outlink, out);
}

static av_cold int init(AVFilterContext *ctx)
{
    NoiseContext *n = ctx->priv;
    int ret, i;

    for (i = 0; i < 4; i++) {
        if (n->all.seed >= 0)
            n->param[i].seed = n->all.seed;
        else
            n->param[i].seed = 123457;
        if (n->all.strength)
            n->param[i].strength = n->all.strength;
        if (n->all.flags)
            n->param[i].flags = n->all.flags;
    }

    for (i = 0; i < 4; i++) {
        if (n->param[i].strength && ((ret = init_noise(n, i)) < 0))
            return ret;
    }

    n->line_noise     = ff_line_noise_c;
    n->line_noise_avg = ff_line_noise_avg_c;

    if (ARCH_X86)
        ff_noise_init_x86(n);

    return 0;
}

static av_cold void uninit(AVFilterContext *ctx)
{
    NoiseContext *n = ctx->priv;
    int i;

    for (i = 0; i < 4; i++)
        av_freep(&n->param[i].noise);
}

static const AVFilterPad noise_inputs[] = {
    {
        .name         = "default",
        .type         = AVMEDIA_TYPE_VIDEO,
        .filter_frame = filter_frame,
        .config_props = config_input,
    },
    { NULL }
};

static const AVFilterPad noise_outputs[] = {
    {
        .name = "default",
        .type = AVMEDIA_TYPE_VIDEO,
    },
    { NULL }
};

AVFilter ff_vf_noise = {
    .name          = "noise",
    .description   = NULL_IF_CONFIG_SMALL("Add noise."),
    .priv_size     = sizeof(NoiseContext),
    .init          = init,
    .uninit        = uninit,
    .query_formats = query_formats,
    .inputs        = noise_inputs,
    .outputs       = noise_outputs,
    .priv_class    = &noise_class,
    .flags         = AVFILTER_FLAG_SUPPORT_TIMELINE_GENERIC | AVFILTER_FLAG_SLICE_THREADS,
};
Commit	Line	Data
2ba45a60 DM	1	/*
	2	* Copyright (c) 2002 Michael Niedermayer <michaelni@gmx.at>
	3	* Copyright (c) 2013 Paul B Mahol
	4	*
	5	* This file is part of FFmpeg.
	6	*
	7	* FFmpeg is free software; you can redistribute it and/or
	8	* modify it under the terms of the GNU Lesser General Public
	9	* License as published by the Free Software Foundation; either
	10	* version 2.1 of the License, or (at your option) any later version.
	11	*
	12	* FFmpeg is distributed in the hope that it will be useful,
	13	* but WITHOUT ANY WARRANTY; without even the implied warranty of
	14	* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
	15	* Lesser General Public License for more details.
	16	*
	17	* You should have received a copy of the GNU Lesser General Public
	18	* License along with FFmpeg; if not, write to the Free Software
	19	* Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
	20	*/
	21
	22	/**
	23	* @file
	24	* noise generator
	25	*/
	26
	27	#include "libavutil/opt.h"
	28	#include "libavutil/imgutils.h"
	29	#include "libavutil/lfg.h"
	30	#include "libavutil/parseutils.h"
	31	#include "libavutil/pixdesc.h"
2ba45a60 DM	32	#include "avfilter.h"
	33	#include "formats.h"
	34	#include "internal.h"
f6fa7814	35	#include "vf_noise.h"
2ba45a60 DM	36	#include "video.h"
2ba45a60 DM	37
2ba45a60 DM	38	typedef struct ThreadData {
	39	AVFrame in, out;
	40	} ThreadData;
	41
	42	#define OFFSET(x) offsetof(NoiseContext, x)
	43	#define FLAGS AV_OPT_FLAG_FILTERING_PARAM\|AV_OPT_FLAG_VIDEO_PARAM
	44
	45	#define NOISE_PARAMS(name, x, param) \
	46	{#name"_seed", "set component #"#x" noise seed", OFFSET(param.seed), AV_OPT_TYPE_INT, {.i64=-1}, -1, INT_MAX, FLAGS}, \
	47	{#name"_strength", "set component #"#x" strength", OFFSET(param.strength), AV_OPT_TYPE_INT, {.i64=0}, 0, 100, FLAGS}, \
	48	{#name"s", "set component #"#x" strength", OFFSET(param.strength), AV_OPT_TYPE_INT, {.i64=0}, 0, 100, FLAGS}, \
	49	{#name"_flags", "set component #"#x" flags", OFFSET(param.flags), AV_OPT_TYPE_FLAGS, {.i64=0}, 0, 31, FLAGS, #name"_flags"}, \
	50	{#name"f", "set component #"#x" flags", OFFSET(param.flags), AV_OPT_TYPE_FLAGS, {.i64=0}, 0, 31, FLAGS, #name"_flags"}, \
	51	{"a", "averaged noise", 0, AV_OPT_TYPE_CONST, {.i64=NOISE_AVERAGED}, 0, 0, FLAGS, #name"_flags"}, \
	52	{"p", "(semi)regular pattern", 0, AV_OPT_TYPE_CONST, {.i64=NOISE_PATTERN}, 0, 0, FLAGS, #name"_flags"}, \
	53	{"t", "temporal noise", 0, AV_OPT_TYPE_CONST, {.i64=NOISE_TEMPORAL}, 0, 0, FLAGS, #name"_flags"}, \
	54	{"u", "uniform noise", 0, AV_OPT_TYPE_CONST, {.i64=NOISE_UNIFORM}, 0, 0, FLAGS, #name"_flags"},
	55
	56	static const AVOption noise_options[] = {
	57	NOISE_PARAMS(all, 0, all)
	58	NOISE_PARAMS(c0, 0, param[0])
	59	NOISE_PARAMS(c1, 1, param[1])
	60	NOISE_PARAMS(c2, 2, param[2])
	61	NOISE_PARAMS(c3, 3, param[3])
	62	{NULL}
	63	};
	64
	65	AVFILTER_DEFINE_CLASS(noise);
	66
	67	static const int8_t patt[4] = { -1, 0, 1, 0 };
	68
	69	#define RAND_N(range) ((int) ((double) range * av_lfg_get(lfg) / (UINT_MAX + 1.0)))
	70	static av_cold int init_noise(NoiseContext *n, int comp)
	71	{
	72	int8_t noise = av_malloc(MAX_NOISE sizeof(int8_t));
	73	FilterParams *fp = &n->param[comp];
	74	AVLFG *lfg = &n->param[comp].lfg;
	75	int strength = fp->strength;
	76	int flags = fp->flags;
	77	int i, j;
	78
	79	if (!noise)
	80	return AVERROR(ENOMEM);
	81
f6fa7814	82	av_lfg_init(&fp->lfg, fp->seed + comp*31415U);
2ba45a60 DM	83
	84	for (i = 0, j = 0; i < MAX_NOISE; i++, j++) {
	85	if (flags & NOISE_UNIFORM) {
	86	if (flags & NOISE_AVERAGED) {
	87	if (flags & NOISE_PATTERN) {
	88	noise[i] = (RAND_N(strength) - strength / 2) / 6
	89	+ patt[j % 4] * strength * 0.25 / 3;
	90	} else {
	91	noise[i] = (RAND_N(strength) - strength / 2) / 3;
	92	}
	93	} else {
	94	if (flags & NOISE_PATTERN) {
	95	noise[i] = (RAND_N(strength) - strength / 2) / 2
	96	+ patt[j % 4] * strength * 0.25;
	97	} else {
	98	noise[i] = RAND_N(strength) - strength / 2;
	99	}
	100	}
	101	} else {
	102	double x1, x2, w, y1;
	103	do {
	104	x1 = 2.0 * av_lfg_get(lfg) / (float)UINT_MAX - 1.0;
	105	x2 = 2.0 * av_lfg_get(lfg) / (float)UINT_MAX - 1.0;
	106	w = x1 * x1 + x2 * x2;
	107	} while (w >= 1.0);
	108
	109	w = sqrt((-2.0 * log(w)) / w);
	110	y1 = x1 * w;
	111	y1 *= strength / sqrt(3.0);
	112	if (flags & NOISE_PATTERN) {
	113	y1 /= 2;
	114	y1 += patt[j % 4] * strength * 0.35;
	115	}
	116	y1 = av_clipf(y1, -128, 127);
	117	if (flags & NOISE_AVERAGED)
	118	y1 /= 3.0;
	119	noise[i] = (int)y1;
	120	}
	121	if (RAND_N(6) == 0)
	122	j--;
	123	}
	124
	125	for (i = 0; i < MAX_RES; i++)
	126	for (j = 0; j < 3; j++)
	127	fp->prev_shift[i][j] = noise + (av_lfg_get(lfg) & (MAX_SHIFT - 1));
	128
2ba45a60 DM	129	fp->noise = noise;
	130	return 0;
	131	}
	132
	133	static int query_formats(AVFilterContext *ctx)
	134	{
	135	AVFilterFormats *formats = NULL;
	136	int fmt;
	137
	138	for (fmt = 0; av_pix_fmt_desc_get(fmt); fmt++) {
	139	const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(fmt);
	140	if (desc->flags & AV_PIX_FMT_FLAG_PLANAR && !((desc->comp[0].depth_minus1 + 1) & 7))
	141	ff_add_format(&formats, fmt);
	142	}
	143
	144	ff_set_common_formats(ctx, formats);
	145	return 0;
	146	}
	147
	148	static int config_input(AVFilterLink *inlink)
	149	{
	150	NoiseContext *n = inlink->dst->priv;
	151	const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(inlink->format);
	152	int ret;
	153
	154	n->nb_planes = av_pix_fmt_count_planes(inlink->format);
	155
	156	if ((ret = av_image_fill_linesizes(n->bytewidth, inlink->format, inlink->w)) < 0)
	157	return ret;
	158
	159	n->height[1] = n->height[2] = FF_CEIL_RSHIFT(inlink->h, desc->log2_chroma_h);
	160	n->height[0] = n->height[3] = inlink->h;
	161
	162	return 0;
	163	}
	164
f6fa7814 DM	165	void ff_line_noise_c(uint8_t dst, const uint8_t src, const int8_t *noise,
f6fa7814 DM	166	int len, int shift)
2ba45a60 DM	167	{
	168	int i;
	169
	170	noise += shift;
	171	for (i = 0; i < len; i++) {
	172	int v = src[i] + noise[i];
	173
	174	dst[i] = av_clip_uint8(v);
	175	}
	176	}
	177
f6fa7814 DM	178	void ff_line_noise_avg_c(uint8_t dst, const uint8_t src,
f6fa7814 DM	179	int len, const int8_t * const *shift)
2ba45a60 DM	180	{
	181	int i;
	182	const int8_t src2 = (const int8_t)src;
	183
	184	for (i = 0; i < len; i++) {
	185	const int n = shift[0][i] + shift[1][i] + shift[2][i];
	186	dst[i] = src2[i] + ((n * src2[i]) >> 7);
	187	}
	188	}
	189
2ba45a60 DM	190	static void noise(uint8_t dst, const uint8_t src,
	191	int dst_linesize, int src_linesize,
	192	int width, int start, int end, NoiseContext *n, int comp)
	193	{
	194	FilterParams *p = &n->param[comp];
	195	int8_t *noise = p->noise;
	196	const int flags = p->flags;
f6fa7814	197	int y;
2ba45a60 DM	198
	199	if (!noise) {
	200	if (dst != src)
	201	av_image_copy_plane(dst, dst_linesize, src, src_linesize, width, end - start);
	202	return;
	203	}
	204
	205	for (y = start; y < end; y++) {
	206	const int ix = y & (MAX_RES - 1);
f6fa7814 DM	207	int x;
	208	for (x=0; x < width; x+= MAX_RES) {
	209	int w = FFMIN(width - x, MAX_RES);
	210	int shift = p->rand_shift[ix];
2ba45a60	211
f6fa7814 DM	212	if (flags & NOISE_AVERAGED) {
	213	n->line_noise_avg(dst + x, src + x, w, (const int8_t**)p->prev_shift[ix]);
	214	p->prev_shift[ix][shift & 3] = noise + shift;
	215	} else {
	216	n->line_noise(dst + x, src + x, noise, w, shift);
	217	}
2ba45a60 DM	218	}
	219	dst += dst_linesize;
	220	src += src_linesize;
	221	}
	222	}
	223
	224	static int filter_slice(AVFilterContext ctx, void arg, int jobnr, int nb_jobs)
	225	{
	226	NoiseContext *s = ctx->priv;
	227	ThreadData *td = arg;
	228	int plane;
	229
	230	for (plane = 0; plane < s->nb_planes; plane++) {
	231	const int height = s->height[plane];
	232	const int start = (height * jobnr ) / nb_jobs;
	233	const int end = (height * (jobnr+1)) / nb_jobs;
	234	noise(td->out->data[plane] + start * td->out->linesize[plane],
	235	td->in->data[plane] + start * td->in->linesize[plane],
	236	td->out->linesize[plane], td->in->linesize[plane],
	237	s->bytewidth[plane], start, end, s, plane);
	238	}
	239	return 0;
	240	}
	241
	242	static int filter_frame(AVFilterLink inlink, AVFrame inpicref)
	243	{
	244	AVFilterContext *ctx = inlink->dst;
	245	AVFilterLink *outlink = ctx->outputs[0];
	246	NoiseContext *n = ctx->priv;
	247	ThreadData td;
	248	AVFrame *out;
f6fa7814	249	int comp, i;
2ba45a60 DM	250
	251	if (av_frame_is_writable(inpicref)) {
	252	out = inpicref;
	253	} else {
	254	out = ff_get_video_buffer(outlink, outlink->w, outlink->h);
	255	if (!out) {
	256	av_frame_free(&inpicref);
	257	return AVERROR(ENOMEM);
	258	}
	259	av_frame_copy_props(out, inpicref);
	260	}
	261
f6fa7814 DM	262	for (comp = 0; comp < 4; comp++) {
	263	FilterParams *fp = &n->param[comp];
	264
	265	if ((!fp->rand_shift_init \|\| (fp->flags & NOISE_TEMPORAL)) && fp->strength) {
	266
	267	for (i = 0; i < MAX_RES; i++) {
	268	fp->rand_shift[i] = av_lfg_get(&fp->lfg) & (MAX_SHIFT - 1);
	269	}
	270	fp->rand_shift_init = 1;
	271	}
	272	}
	273
2ba45a60 DM	274	td.in = inpicref; td.out = out;
	275	ctx->internal->execute(ctx, filter_slice, &td, NULL, FFMIN(n->height[0], ctx->graph->nb_threads));
	276	emms_c();
	277
	278	if (inpicref != out)
	279	av_frame_free(&inpicref);
	280	return ff_filter_frame(outlink, out);
	281	}
	282
	283	static av_cold int init(AVFilterContext *ctx)
	284	{
	285	NoiseContext *n = ctx->priv;
	286	int ret, i;
2ba45a60 DM	287
	288	for (i = 0; i < 4; i++) {
	289	if (n->all.seed >= 0)
	290	n->param[i].seed = n->all.seed;
	291	else
	292	n->param[i].seed = 123457;
	293	if (n->all.strength)
	294	n->param[i].strength = n->all.strength;
	295	if (n->all.flags)
	296	n->param[i].flags = n->all.flags;
	297	}
	298
	299	for (i = 0; i < 4; i++) {
	300	if (n->param[i].strength && ((ret = init_noise(n, i)) < 0))
	301	return ret;
	302	}
	303
f6fa7814 DM	304	n->line_noise = ff_line_noise_c;
f6fa7814 DM	305	n->line_noise_avg = ff_line_noise_avg_c;
2ba45a60	306
f6fa7814 DM	307	if (ARCH_X86)
f6fa7814 DM	308	ff_noise_init_x86(n);
2ba45a60 DM	309
	310	return 0;
	311	}
	312
	313	static av_cold void uninit(AVFilterContext *ctx)
	314	{
	315	NoiseContext *n = ctx->priv;
	316	int i;
	317
	318	for (i = 0; i < 4; i++)
	319	av_freep(&n->param[i].noise);
	320	}
	321
	322	static const AVFilterPad noise_inputs[] = {
	323	{
	324	.name = "default",
	325	.type = AVMEDIA_TYPE_VIDEO,
	326	.filter_frame = filter_frame,
	327	.config_props = config_input,
	328	},
	329	{ NULL }
	330	};
	331
	332	static const AVFilterPad noise_outputs[] = {
	333	{
	334	.name = "default",
	335	.type = AVMEDIA_TYPE_VIDEO,
	336	},
	337	{ NULL }
	338	};
	339
	340	AVFilter ff_vf_noise = {
	341	.name = "noise",
	342	.description = NULL_IF_CONFIG_SMALL("Add noise."),
	343	.priv_size = sizeof(NoiseContext),
	344	.init = init,
	345	.uninit = uninit,
	346	.query_formats = query_formats,
	347	.inputs = noise_inputs,
	348	.outputs = noise_outputs,
	349	.priv_class = &noise_class,
	350	.flags = AVFILTER_FLAG_SUPPORT_TIMELINE_GENERIC \| AVFILTER_FLAG_SLICE_THREADS,
	351	};