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