/*
* Copyright (c) Lynne
* Copyright (C) 2018 Philip Langdale <philipl@overt.org>
* Copyright (C) 2016 Thomas Mundt <loudmax@yahoo.de>
*
* 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
*/
#include "libavutil/random_seed.h"
#include "libavutil/opt.h"
#include "vulkan_filter.h"
#include "vulkan_spirv.h"
#include "yadif.h"
#include "internal.h"
typedef struct BWDIFVulkanContext {
YADIFContext yadif;
FFVulkanContext vkctx;
int initialized;
FFVkExecPool e;
FFVkQueueFamilyCtx qf;
VkSampler sampler;
FFVulkanPipeline pl;
FFVkSPIRVShader shd;
} BWDIFVulkanContext;
typedef struct BWDIFParameters {
int parity;
int tff;
int current_field;
} BWDIFParameters;
static const char filter_fn[] = {
"const vec4 coef_lf[2] = { vec4(4309), vec4(213), };\n"
"const vec4 coef_hf[3] = { vec4(5570), vec4(3801), vec4(1016) };\n"
"const vec4 coef_sp[2] = { vec4(5077), vec4(981), };\n"
C(0, )
C(0, vec4 process_intra(vec4 cur[4]) )
C(0, { )
C(1, return (coef_sp[0]*(cur[1] + cur[2]) - coef_sp[1]*(cur[0] + cur[3])) / (1 << 13); )
C(0, } )
C(0, )
C(0, vec4 process_line(vec4 prev2[5], vec4 prev1[2], vec4 cur[4], vec4 next1[2], vec4 next2[5]) )
C(0, { )
C(1, vec4 fc = cur[1]; )
C(1, vec4 fe = cur[2]; )
C(1, vec4 fs = prev2[2] + next2[2]; )
C(1, vec4 fd = fs / 2; )
C(0, )
C(1, vec4 temp_diff[3]; )
C(1, temp_diff[0] = abs(prev2[2] - next2[2]); )
C(1, temp_diff[1] = (abs(prev1[0] - fc) + abs(prev1[1] - fe)) / 2; )
C(1, temp_diff[1] = (abs(next1[0] - fc) + abs(next1[1] - fe)) / 2; )
C(1, vec4 diff = max(temp_diff[0] / 2, max(temp_diff[1], temp_diff[2])); )
C(1, bvec4 diff_mask = equal(diff, vec4(0)); )
C(0, )
C(1, vec4 fbs = prev2[1] + next2[1]; )
C(1, vec4 ffs = prev2[3] + next2[3]; )
C(1, vec4 fb = (fbs / 2) - fc; )
C(1, vec4 ff = (ffs / 2) - fe; )
C(1, vec4 dc = fd - fc; )
C(1, vec4 de = fd - fe; )
C(1, vec4 mmax = max(de, max(dc, min(fb, ff))); )
C(1, vec4 mmin = min(de, min(dc, max(fb, ff))); )
C(1, diff = max(diff, max(mmin, -mmax)); )
C(0, )
" vec4 interpolate_all = (((coef_hf[0]*(fs) - coef_hf[1]*(fbs + ffs) +\n"
" coef_hf[2]*(prev2[0] + next2[0] + prev2[4] + next2[4])) / 4) +\n"
" coef_lf[0]*(fc + fe) - coef_lf[1]*(cur[0] + cur[3])) / (1 << 13);\n"
" vec4 interpolate_cur = (coef_sp[0]*(fc + fe) - coef_sp[1]*(cur[0] + cur[3])) / (1 << 13);\n"
C(0, )
C(1, bvec4 interpolate_cnd1 = greaterThan(abs(fc - fe), temp_diff[0]); )
C(1, vec4 interpol = mix(interpolate_cur, interpolate_all, interpolate_cnd1); )
C(1, interpol = clamp(interpol, fd - diff, fd + diff); )
C(1, return mix(interpol, fd, diff_mask); )
C(0, } )
};
static av_cold int init_filter(AVFilterContext *ctx)
{
int err;
uint8_t *spv_data;
size_t spv_len;
void *spv_opaque = NULL;
BWDIFVulkanContext *s = ctx->priv;
FFVulkanContext *vkctx = &s->vkctx;
const int planes = av_pix_fmt_count_planes(s->vkctx.output_format);
FFVkSPIRVShader *shd;
FFVkSPIRVCompiler *spv;
FFVulkanDescriptorSetBinding *desc;
spv = ff_vk_spirv_init();
if (!spv) {
av_log(ctx, AV_LOG_ERROR, "Unable to initialize SPIR-V compiler!\n");
return AVERROR_EXTERNAL;
}
ff_vk_qf_init(vkctx, &s->qf, VK_QUEUE_COMPUTE_BIT);
RET(ff_vk_exec_pool_init(vkctx, &s->qf, &s->e, s->qf.nb_queues*4, 0, 0, 0, NULL));
RET(ff_vk_init_sampler(vkctx, &s->sampler, 1, VK_FILTER_NEAREST));
RET(ff_vk_shader_init(&s->pl, &s->shd, "bwdif_compute",
VK_SHADER_STAGE_COMPUTE_BIT, 0));
shd = &s->shd;
ff_vk_shader_set_compute_sizes(shd, 1, 64, 1);
desc = (FFVulkanDescriptorSetBinding []) {
{
.name = "prev",
.type = VK_DESCRIPTOR_TYPE_COMBINED_IMAGE_SAMPLER,
.dimensions = 2,
.elems = planes,
.stages = VK_SHADER_STAGE_COMPUTE_BIT,
.samplers = DUP_SAMPLER(s->sampler),
},
{
.name = "cur",
.type = VK_DESCRIPTOR_TYPE_COMBINED_IMAGE_SAMPLER,
.dimensions = 2,
.elems = planes,
.stages = VK_SHADER_STAGE_COMPUTE_BIT,
.samplers = DUP_SAMPLER(s->sampler),
},
{
.name = "next",
.type = VK_DESCRIPTOR_TYPE_COMBINED_IMAGE_SAMPLER,
.dimensions = 2,
.elems = planes,
.stages = VK_SHADER_STAGE_COMPUTE_BIT,
.samplers = DUP_SAMPLER(s->sampler),
},
{
.name = "dst",
.type = VK_DESCRIPTOR_TYPE_STORAGE_IMAGE,
.mem_layout = ff_vk_shader_rep_fmt(s->vkctx.output_format),
.mem_quali = "writeonly",
.dimensions = 2,
.elems = planes,
.stages = VK_SHADER_STAGE_COMPUTE_BIT,
},
};
RET(ff_vk_pipeline_descriptor_set_add(vkctx, &s->pl, shd, desc, 4, 0, 0));
GLSLC(0, layout(push_constant, std430) uniform pushConstants { );
GLSLC(1, int parity; );
GLSLC(1, int tff; );
GLSLC(1, int current_field; );
GLSLC(0, }; );
ff_vk_add_push_constant(&s->pl, 0, sizeof(BWDIFParameters),
VK_SHADER_STAGE_COMPUTE_BIT);
GLSLD( filter_fn );
GLSLC(0, void main() );
GLSLC(0, { );
GLSLC(1, vec4 res; );
GLSLC(1, ivec2 size; );
GLSLC(1, vec4 dcur[4]; );
GLSLC(1, vec4 prev1[2]; );
GLSLC(1, vec4 next1[2]; );
GLSLC(1, vec4 prev2[5]; );
GLSLC(1, vec4 next2[5]; );
GLSLC(1, const ivec2 pos = ivec2(gl_GlobalInvocationID.xy); );
GLSLC(1, bool filter_field = ((pos.y ^ parity) & 1) == 1; );
GLSLF(1, bool is_intra = filter_field && (current_field == %i); ,YADIF_FIELD_END);
GLSLC(1, bool field_parity = (parity ^ tff) != 0; );
GLSLC(0, );
for (int i = 0; i < planes; i++) {
GLSLC(0, );
GLSLF(1, size = imageSize(dst[%i]); ,i);
GLSLC(1, if (!IS_WITHIN(pos, size)) { );
GLSLC(2, return; );
GLSLC(1, } else if (is_intra) { );
GLSLF(2, dcur[0] = texture(cur[%i], pos - ivec2(0, 3)); ,i);
GLSLF(2, dcur[1] = texture(cur[%i], pos - ivec2(0, 1)); ,i);
GLSLF(2, dcur[2] = texture(cur[%i], pos + ivec2(0, 1)); ,i);
GLSLF(2, dcur[3] = texture(cur[%i], pos + ivec2(0, 3)); ,i);
GLSLC(0, );
GLSLC(2, res = process_intra(dcur); );
GLSLF(2, imageStore(dst[%i], pos, res); ,i);
GLSLC(1, } else if (filter_field) { );
GLSLF(2, dcur[0] = texture(cur[%i], pos - ivec2(0, 3)); ,i);
GLSLF(2, dcur[1] = texture(cur[%i], pos - ivec2(0, 1)); ,i);
GLSLF(2, dcur[2] = texture(cur[%i], pos + ivec2(0, 1)); ,i);
GLSLF(2, dcur[3] = texture(cur[%i], pos + ivec2(0, 3)); ,i);
GLSLC(0, );
GLSLF(2, prev1[0] = texture(prev[%i], pos - ivec2(0, 1)); ,i);
GLSLF(2, prev1[1] = texture(prev[%i], pos + ivec2(0, 1)); ,i);
GLSLC(0, );
GLSLF(2, next1[0] = texture(next[%i], pos - ivec2(0, 1)); ,i);
GLSLF(2, next1[1] = texture(next[%i], pos + ivec2(0, 1)); ,i);
GLSLC(0, );
GLSLC(2, if (field_parity) { );
GLSLF(3, prev2[0] = texture(prev[%i], pos - ivec2(0, 4)); ,i);
GLSLF(3, prev2[1] = texture(prev[%i], pos - ivec2(0, 2)); ,i);
GLSLF(3, prev2[2] = texture(prev[%i], pos); ,i);
GLSLF(3, prev2[3] = texture(prev[%i], pos + ivec2(0, 2)); ,i);
GLSLF(3, prev2[4] = texture(prev[%i], pos + ivec2(0, 4)); ,i);
GLSLC(0, );
GLSLF(3, next2[0] = texture(cur[%i], pos - ivec2(0, 4)); ,i);
GLSLF(3, next2[1] = texture(cur[%i], pos - ivec2(0, 2)); ,i);
GLSLF(3, next2[2] = texture(cur[%i], pos); ,i);
GLSLF(3, next2[3] = texture(cur[%i], pos + ivec2(0, 2)); ,i);
GLSLF(3, next2[4] = texture(cur[%i], pos + ivec2(0, 4)); ,i);
GLSLC(2, } else { );
GLSLF(3, prev2[0] = texture(cur[%i], pos - ivec2(0, 4)); ,i);
GLSLF(3, prev2[1] = texture(cur[%i], pos - ivec2(0, 2)); ,i);
GLSLF(3, prev2[2] = texture(cur[%i], pos); ,i);
GLSLF(3, prev2[3] = texture(cur[%i], pos + ivec2(0, 2)); ,i);
GLSLF(3, prev2[4] = texture(cur[%i], pos + ivec2(0, 4)); ,i);
GLSLC(0, );
GLSLF(3, next2[0] = texture(next[%i], pos - ivec2(0, 4)); ,i);
GLSLF(3, next2[1] = texture(next[%i], pos - ivec2(0, 2)); ,i);
GLSLF(3, next2[2] = texture(next[%i], pos); ,i);
GLSLF(3, next2[3] = texture(next[%i], pos + ivec2(0, 2)); ,i);
GLSLF(3, next2[4] = texture(next[%i], pos + ivec2(0, 4)); ,i);
GLSLC(2, } );
GLSLC(0, );
GLSLC(2, res = process_line(prev2, prev1, dcur, next1, next2); );
GLSLF(2, imageStore(dst[%i], pos, res); ,i);
GLSLC(1, } else { );
GLSLF(2, res = texture(cur[%i], pos); ,i);
GLSLF(2, imageStore(dst[%i], pos, res); ,i);
GLSLC(1, } );
}
GLSLC(0, } );
RET(spv->compile_shader(spv, ctx, &s->shd, &spv_data, &spv_len, "main",
&spv_opaque));
RET(ff_vk_shader_create(vkctx, &s->shd, spv_data, spv_len, "main"));
RET(ff_vk_init_compute_pipeline(vkctx, &s->pl, &s->shd));
RET(ff_vk_exec_pipeline_register(vkctx, &s->e, &s->pl));
s->initialized = 1;
fail:
if (spv_opaque)
spv->free_shader(spv, &spv_opaque);
if (spv)
spv->uninit(&spv);
return err;
}
static void bwdif_vulkan_filter_frame(AVFilterContext *ctx, AVFrame *dst,
int parity, int tff)
{
BWDIFVulkanContext *s = ctx->priv;
YADIFContext *y = &s->yadif;
BWDIFParameters params = {
.parity = parity,
.tff = tff,
.current_field = y->current_field,
};
ff_vk_filter_process_Nin(&s->vkctx, &s->e, &s->pl, dst,
(AVFrame *[]){ y->prev, y->cur, y->next }, 3,
s->sampler, ¶ms, sizeof(params));
if (y->current_field == YADIF_FIELD_END)
y->current_field = YADIF_FIELD_NORMAL;
}
static void bwdif_vulkan_uninit(AVFilterContext *avctx)
{
BWDIFVulkanContext *s = avctx->priv;
FFVulkanContext *vkctx = &s->vkctx;
FFVulkanFunctions *vk = &vkctx->vkfn;
ff_vk_exec_pool_free(vkctx, &s->e);
ff_vk_pipeline_free(vkctx, &s->pl);
ff_vk_shader_free(vkctx, &s->shd);
if (s->sampler)
vk->DestroySampler(vkctx->hwctx->act_dev, s->sampler,
vkctx->hwctx->alloc);
ff_vk_uninit(&s->vkctx);
ff_yadif_uninit(avctx);
s->initialized = 0;
}
static int bwdif_vulkan_config_input(AVFilterLink *inlink)
{
AVHWFramesContext *input_frames;
AVFilterContext *avctx = inlink->dst;
BWDIFVulkanContext *s = avctx->priv;
FFVulkanContext *vkctx = &s->vkctx;
if (!inlink->hw_frames_ctx) {
av_log(inlink->dst, AV_LOG_ERROR, "Vulkan filtering requires a "
"hardware frames context on the input.\n");
return AVERROR(EINVAL);
}
input_frames = (AVHWFramesContext *)inlink->hw_frames_ctx->data;
if (input_frames->format != AV_PIX_FMT_VULKAN)
return AVERROR(EINVAL);
/* Extract the device and default output format from the first input. */
if (avctx->inputs[0] != inlink)
return 0;
/* Save the ref, without reffing it */
vkctx->input_frames_ref = inlink->hw_frames_ctx;
/* Defaults */
vkctx->output_format = input_frames->sw_format;
vkctx->output_width = inlink->w;
vkctx->output_height = inlink->h;
return 0;
}
static int bwdif_vulkan_config_output(AVFilterLink *outlink)
{
int err;
AVFilterContext *avctx = outlink->src;
BWDIFVulkanContext *s = avctx->priv;
YADIFContext *y = &s->yadif;
FFVulkanContext *vkctx = &s->vkctx;
av_buffer_unref(&outlink->hw_frames_ctx);
err = ff_vk_filter_init_context(avctx, vkctx, vkctx->input_frames_ref,
vkctx->output_width, vkctx->output_height,
vkctx->output_format);
if (err < 0)
return err;
/* For logging */
vkctx->class = y->class;
outlink->hw_frames_ctx = av_buffer_ref(vkctx->frames_ref);
if (!outlink->hw_frames_ctx)
return AVERROR(ENOMEM);
err = ff_yadif_config_output_common(outlink);
if (err < 0)
return err;
y->csp = av_pix_fmt_desc_get(vkctx->frames->sw_format);
y->filter = bwdif_vulkan_filter_frame;
if (AV_CEIL_RSHIFT(outlink->w, y->csp->log2_chroma_w) < 4 || AV_CEIL_RSHIFT(outlink->h, y->csp->log2_chroma_h) < 4) {
av_log(avctx, AV_LOG_ERROR, "Video with planes less than 4 columns or lines is not supported\n");
return AVERROR(EINVAL);
}
return init_filter(avctx);
}
static const AVClass bwdif_vulkan_class = {
.class_name = "bwdif_vulkan",
.item_name = av_default_item_name,
.option = ff_yadif_options,
.version = LIBAVUTIL_VERSION_INT,
.category = AV_CLASS_CATEGORY_FILTER,
};
static const AVFilterPad bwdif_vulkan_inputs[] = {
{
.name = "default",
.type = AVMEDIA_TYPE_VIDEO,
.filter_frame = ff_yadif_filter_frame,
.config_props = &bwdif_vulkan_config_input,
},
};
static const AVFilterPad bwdif_vulkan_outputs[] = {
{
.name = "default",
.type = AVMEDIA_TYPE_VIDEO,
.request_frame = ff_yadif_request_frame,
.config_props = &bwdif_vulkan_config_output,
},
};
const AVFilter ff_vf_bwdif_vulkan = {
.name = "bwdif_vulkan",
.description = NULL_IF_CONFIG_SMALL("Deinterlace Vulkan frames via bwdif"),
.priv_size = sizeof(BWDIFVulkanContext),
.init = &ff_vk_filter_init,
.uninit = &bwdif_vulkan_uninit,
FILTER_INPUTS(bwdif_vulkan_inputs),
FILTER_OUTPUTS(bwdif_vulkan_outputs),
FILTER_SINGLE_PIXFMT(AV_PIX_FMT_VULKAN),
.priv_class = &bwdif_vulkan_class,
.flags = AVFILTER_FLAG_HWDEVICE |
AVFILTER_FLAG_SUPPORT_TIMELINE_INTERNAL,
.flags_internal = FF_FILTER_FLAG_HWFRAME_AWARE,
};