git.sesse.net Git - ffmpeg/blob - libavcodec/cfhddsp.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 "cfhddsp.h"
  26
  27 static av_always_inline void filter(int16_t *output, ptrdiff_t out_stride,
  28                           const int16_t *low, ptrdiff_t low_stride,
  29                           const int16_t *high, ptrdiff_t high_stride,
  30                           int len, int clip)
  31 {
  32     int16_t tmp;
  33     int i;
  34
  35     tmp = (11*low[0*low_stride] - 4*low[1*low_stride] + low[2*low_stride] + 4) >> 3;
  36     output[(2*0+0)*out_stride] = (tmp + high[0*high_stride]) >> 1;
  37     if (clip)
  38         output[(2*0+0)*out_stride] = av_clip_uintp2_c(output[(2*0+0)*out_stride], clip);
  39
  40     tmp = ( 5*low[0*low_stride] + 4*low[1*low_stride] - low[2*low_stride] + 4) >> 3;
  41     output[(2*0+1)*out_stride] = (tmp - high[0*high_stride]) >> 1;
  42     if (clip)
  43         output[(2*0+1)*out_stride] = av_clip_uintp2_c(output[(2*0+1)*out_stride], clip);
  44
  45     for (i = 1; i < len - 1; i++) {
  46         tmp = (low[(i-1)*low_stride] - low[(i+1)*low_stride] + 4) >> 3;
  47         output[(2*i+0)*out_stride] = (tmp + low[i*low_stride] + high[i*high_stride]) >> 1;
  48         if (clip)
  49             output[(2*i+0)*out_stride] = av_clip_uintp2_c(output[(2*i+0)*out_stride], clip);
  50
  51         tmp = (low[(i+1)*low_stride] - low[(i-1)*low_stride] + 4) >> 3;
  52         output[(2*i+1)*out_stride] = (tmp + low[i*low_stride] - high[i*high_stride]) >> 1;
  53         if (clip)
  54             output[(2*i+1)*out_stride] = av_clip_uintp2_c(output[(2*i+1)*out_stride], clip);
  55     }
  56
  57     tmp = ( 5*low[i*low_stride] + 4*low[(i-1)*low_stride] - low[(i-2)*low_stride] + 4) >> 3;
  58     output[(2*i+0)*out_stride] = (tmp + high[i*high_stride]) >> 1;
  59     if (clip)
  60         output[(2*i+0)*out_stride] = av_clip_uintp2_c(output[(2*i+0)*out_stride], clip);
  61
  62     tmp = (11*low[i*low_stride] - 4*low[(i-1)*low_stride] + low[(i-2)*low_stride] + 4) >> 3;
  63     output[(2*i+1)*out_stride] = (tmp - high[i*high_stride]) >> 1;
  64     if (clip)
  65         output[(2*i+1)*out_stride] = av_clip_uintp2_c(output[(2*i+1)*out_stride], clip);
  66 }
  67
  68 static void vert_filter(int16_t *output, ptrdiff_t out_stride,
  69                         const int16_t *low, ptrdiff_t low_stride,
  70                         const int16_t *high, ptrdiff_t high_stride,
  71                         int width, int height)
  72 {
  73     for (int i = 0; i < width; i++) {
  74         filter(output, out_stride, low, low_stride, high, high_stride, height, 0);
  75         low++;
  76         high++;
  77         output++;
  78     }
  79 }
  80
  81 static void horiz_filter(int16_t *output, ptrdiff_t ostride,
  82                          const int16_t *low, ptrdiff_t lstride,
  83                          const int16_t *high, ptrdiff_t hstride,
  84                          int width, int height)
  85 {
  86     for (int i = 0; i < height; i++) {
  87         filter(output, 1, low, 1, high, 1, width, 0);
  88         low    += lstride;
  89         high   += hstride;
  90         output += ostride * 2;
  91     }
  92 }
  93
  94 static void horiz_filter_clip(int16_t *output, const int16_t *low, const int16_t *high,
  95                               int width, int clip)
  96 {
  97     filter(output, 1, low, 1, high, 1, width, clip);
  98 }
  99
 100 static void horiz_filter_clip_bayer(int16_t *output, const int16_t *low, const int16_t *high,
 101                                     int width, int clip)
 102 {
 103     filter(output, 2, low, 1, high, 1, width, clip);
 104 }
 105
 106 av_cold void ff_cfhddsp_init(CFHDDSPContext *c, int depth, int bayer)
 107 {
 108     c->horiz_filter = horiz_filter;
 109     c->vert_filter = vert_filter;
 110
 111     if (bayer)
 112         c->horiz_filter_clip = horiz_filter_clip_bayer;
 113     else
 114         c->horiz_filter_clip = horiz_filter_clip;
 115
 116     if (ARCH_X86)
 117         ff_cfhddsp_init_x86(c, depth, bayer);
 118 }