git.sesse.net Git - ffmpeg/blob - libavcodec/cfhdencdsp.c

   1 /*
   2  * Copyright (c) 2015-2016 Kieran Kunhya <kieran@kunhya.com>
   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
  21 #include "libavutil/attributes.h"
  22 #include "libavutil/common.h"
  23 #include "libavutil/avassert.h"
  24
  25 #include "cfhdencdsp.h"
  26
  27 static av_always_inline void filter(int16_t *input, ptrdiff_t in_stride,
  28                           int16_t *low, ptrdiff_t low_stride,
  29                           int16_t *high, ptrdiff_t high_stride,
  30                           int len)
  31 {
  32     low[(0>>1) * low_stride]   = av_clip_int16(input[0*in_stride] + input[1*in_stride]);
  33     high[(0>>1) * high_stride] = av_clip_int16((5 * input[0*in_stride] - 11 * input[1*in_stride] +
  34                                                 4 * input[2*in_stride] +  4 * input[3*in_stride] -
  35                                                 1 * input[4*in_stride] -  1 * input[5*in_stride] + 4) >> 3);
  36
  37     for (int i = 2; i < len - 2; i += 2) {
  38         low[(i>>1) * low_stride]   = av_clip_int16(input[i*in_stride] + input[(i+1)*in_stride]);
  39         high[(i>>1) * high_stride] = av_clip_int16(((-input[(i-2)*in_stride] - input[(i-1)*in_stride] +
  40                                                       input[(i+2)*in_stride] + input[(i+3)*in_stride] + 4) >> 3) +
  41                                                       input[(i+0)*in_stride] - input[(i+1)*in_stride]);
  42     }
  43
  44     low[((len-2)>>1) * low_stride]   = av_clip_int16(input[((len-2)+0)*in_stride] + input[((len-2)+1)*in_stride]);
  45     high[((len-2)>>1) * high_stride] = av_clip_int16((11* input[((len-2)+0)*in_stride] - 5 * input[((len-2)+1)*in_stride] -
  46                                                       4 * input[((len-2)-1)*in_stride] - 4 * input[((len-2)-2)*in_stride] +
  47                                                       1 * input[((len-2)-3)*in_stride] + 1 * input[((len-2)-4)*in_stride] + 4) >> 3);
  48 }
  49
  50 static void horiz_filter(int16_t *input, int16_t *low, int16_t *high,
  51                          ptrdiff_t in_stride, ptrdiff_t low_stride,
  52                          ptrdiff_t high_stride,
  53                          int width, int height)
  54 {
  55     for (int i = 0; i < height; i++) {
  56         filter(input, 1, low, 1, high, 1, width);
  57         input += in_stride;
  58         low += low_stride;
  59         high += high_stride;
  60     }
  61 }
  62
  63 static void vert_filter(int16_t *input, int16_t *low, int16_t *high,
  64                         ptrdiff_t in_stride, ptrdiff_t low_stride,
  65                         ptrdiff_t high_stride,
  66                         int width, int height)
  67 {
  68     for (int i = 0; i < width; i++)
  69         filter(&input[i], in_stride, &low[i], low_stride, &high[i], high_stride, height);
  70 }
  71
  72 av_cold void ff_cfhdencdsp_init(CFHDEncDSPContext *c)
  73 {
  74     c->horiz_filter = horiz_filter;
  75     c->vert_filter = vert_filter;
  76
  77     if (ARCH_X86)
  78         ff_cfhdencdsp_init_x86(c);
  79 }