lavfi/vf_cropdetect: replace round by lrint
[ffmpeg.git] / libavfilter / vf_stereo3d.c
1 /*
2  * Copyright (c) 2010 Gordon Schmidt <gordon.schmidt <at> s2000.tu-chemnitz.de>
3  * Copyright (c) 2013-2015 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 General Public
9  * License as published by the Free Software Foundation; either
10  * version 2 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
15  * GNU General Public License for more details.
16  *
17  * You should have received a copy of the GNU General Public License along
18  * 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 #include "libavutil/avassert.h"
23 #include "libavutil/imgutils.h"
24 #include "libavutil/intreadwrite.h"
25 #include "libavutil/opt.h"
26 #include "libavutil/parseutils.h"
27 #include "libavutil/pixdesc.h"
28 #include "avfilter.h"
29 #include "drawutils.h"
30 #include "formats.h"
31 #include "internal.h"
32 #include "video.h"
33 #include "stereo3d.h"
34
35 enum StereoCode {
36     ANAGLYPH_RC_GRAY,   // anaglyph red/cyan gray
37     ANAGLYPH_RC_HALF,   // anaglyph red/cyan half colored
38     ANAGLYPH_RC_COLOR,  // anaglyph red/cyan colored
39     ANAGLYPH_RC_DUBOIS, // anaglyph red/cyan dubois
40     ANAGLYPH_GM_GRAY,   // anaglyph green/magenta gray
41     ANAGLYPH_GM_HALF,   // anaglyph green/magenta half colored
42     ANAGLYPH_GM_COLOR,  // anaglyph green/magenta colored
43     ANAGLYPH_GM_DUBOIS, // anaglyph green/magenta dubois
44     ANAGLYPH_YB_GRAY,   // anaglyph yellow/blue gray
45     ANAGLYPH_YB_HALF,   // anaglyph yellow/blue half colored
46     ANAGLYPH_YB_COLOR,  // anaglyph yellow/blue colored
47     ANAGLYPH_YB_DUBOIS, // anaglyph yellow/blue dubois
48     ANAGLYPH_RB_GRAY,   // anaglyph red/blue gray
49     ANAGLYPH_RG_GRAY,   // anaglyph red/green gray
50     MONO_L,             // mono output for debugging (left eye only)
51     MONO_R,             // mono output for debugging (right eye only)
52     INTERLEAVE_ROWS_LR, // row-interleave (left eye has top row)
53     INTERLEAVE_ROWS_RL, // row-interleave (right eye has top row)
54     SIDE_BY_SIDE_LR,    // side by side parallel (left eye left, right eye right)
55     SIDE_BY_SIDE_RL,    // side by side crosseye (right eye left, left eye right)
56     SIDE_BY_SIDE_2_LR,  // side by side parallel with half width resolution
57     SIDE_BY_SIDE_2_RL,  // side by side crosseye with half width resolution
58     ABOVE_BELOW_LR,     // above-below (left eye above, right eye below)
59     ABOVE_BELOW_RL,     // above-below (right eye above, left eye below)
60     ABOVE_BELOW_2_LR,   // above-below with half height resolution
61     ABOVE_BELOW_2_RL,   // above-below with half height resolution
62     ALTERNATING_LR,     // alternating frames (left eye first, right eye second)
63     ALTERNATING_RL,     // alternating frames (right eye first, left eye second)
64     CHECKERBOARD_LR,    // checkerboard pattern (left eye first, right eye second)
65     CHECKERBOARD_RL,    // checkerboard pattern (right eye first, left eye second)
66     INTERLEAVE_COLS_LR, // column-interleave (left eye first, right eye second)
67     INTERLEAVE_COLS_RL, // column-interleave (right eye first, left eye second)
68     HDMI,               // HDMI frame pack (left eye first, right eye second)
69     STEREO_CODE_COUNT   // TODO: needs autodetection
70 };
71
72 typedef struct StereoComponent {
73     int format;                 ///< StereoCode
74     int width, height;
75     int off_left, off_right;
76     int off_lstep, off_rstep;
77     int row_left, row_right;
78     int row_step;
79 } StereoComponent;
80
81 static const int ana_coeff[][3][6] = {
82   [ANAGLYPH_RB_GRAY]   =
83     {{19595, 38470,  7471,     0,     0,     0},
84      {    0,     0,     0,     0,     0,     0},
85      {    0,     0,     0, 19595, 38470,  7471}},
86   [ANAGLYPH_RG_GRAY]   =
87     {{19595, 38470,  7471,     0,     0,     0},
88      {    0,     0,     0, 19595, 38470,  7471},
89      {    0,     0,     0,     0,     0,     0}},
90   [ANAGLYPH_RC_GRAY]   =
91     {{19595, 38470,  7471,     0,     0,     0},
92      {    0,     0,     0, 19595, 38470,  7471},
93      {    0,     0,     0, 19595, 38470,  7471}},
94   [ANAGLYPH_RC_HALF]   =
95     {{19595, 38470,  7471,     0,     0,     0},
96      {    0,     0,     0,     0, 65536,     0},
97      {    0,     0,     0,     0,     0, 65536}},
98   [ANAGLYPH_RC_COLOR]  =
99     {{65536,     0,     0,     0,     0,     0},
100      {    0,     0,     0,     0, 65536,     0},
101      {    0,     0,     0,     0,     0, 65536}},
102   [ANAGLYPH_RC_DUBOIS] =
103     {{29891, 32800, 11559, -2849, -5763,  -102},
104      {-2627, -2479, -1033, 24804, 48080, -1209},
105      { -997, -1350,  -358, -4729, -7403, 80373}},
106   [ANAGLYPH_GM_GRAY]   =
107     {{    0,     0,     0, 19595, 38470,  7471},
108      {19595, 38470,  7471,     0,     0,     0},
109      {    0,     0,     0, 19595, 38470,  7471}},
110   [ANAGLYPH_GM_HALF]   =
111     {{    0,     0,     0, 65536,     0,     0},
112      {19595, 38470,  7471,     0,     0,     0},
113      {    0,     0,     0,     0,     0, 65536}},
114   [ANAGLYPH_GM_COLOR]  =
115     {{    0,     0,     0, 65536,     0,     0},
116      {    0, 65536,     0,     0,     0,     0},
117      {    0,     0,     0,     0,     0, 65536}},
118   [ANAGLYPH_GM_DUBOIS]  =
119     {{-4063,-10354, -2556, 34669, 46203,  1573},
120      {18612, 43778,  9372, -1049,  -983, -4260},
121      { -983, -1769,  1376,   590,  4915, 61407}},
122   [ANAGLYPH_YB_GRAY]   =
123     {{    0,     0,     0, 19595, 38470,  7471},
124      {    0,     0,     0, 19595, 38470,  7471},
125      {19595, 38470,  7471,     0,     0,     0}},
126   [ANAGLYPH_YB_HALF]   =
127     {{    0,     0,     0, 65536,     0,     0},
128      {    0,     0,     0,     0, 65536,     0},
129      {19595, 38470,  7471,     0,     0,     0}},
130   [ANAGLYPH_YB_COLOR]  =
131     {{    0,     0,     0, 65536,     0,     0},
132      {    0,     0,     0,     0, 65536,     0},
133      {    0,     0, 65536,     0,     0,     0}},
134   [ANAGLYPH_YB_DUBOIS] =
135     {{65535,-12650,18451,   -987, -7590, -1049},
136      {-1604, 56032, 4196,    370,  3826, -1049},
137      {-2345,-10676, 1358,   5801, 11416, 56217}},
138 };
139
140 typedef struct Stereo3DContext {
141     const AVClass *class;
142     StereoComponent in, out;
143     int width, height;
144     const int *ana_matrix[3];
145     int nb_planes;
146     int linesize[4];
147     int pheight[4];
148     int hsub, vsub;
149     int pixstep[4];
150     AVFrame *prev;
151     int blanks;
152     int in_off_left[4], in_off_right[4];
153     Stereo3DDSPContext dsp;
154 } Stereo3DContext;
155
156 #define OFFSET(x) offsetof(Stereo3DContext, x)
157 #define FLAGS AV_OPT_FLAG_FILTERING_PARAM|AV_OPT_FLAG_VIDEO_PARAM
158
159 static const AVOption stereo3d_options[] = {
160     { "in",    "set input format",  OFFSET(in.format),   AV_OPT_TYPE_INT,   {.i64=SIDE_BY_SIDE_LR}, INTERLEAVE_ROWS_LR, STEREO_CODE_COUNT-1, FLAGS, "in"},
161     { "ab2l",  "above below half height left first",  0, AV_OPT_TYPE_CONST, {.i64=ABOVE_BELOW_2_LR},   0, 0, FLAGS, "in" },
162     { "ab2r",  "above below half height right first", 0, AV_OPT_TYPE_CONST, {.i64=ABOVE_BELOW_2_RL},   0, 0, FLAGS, "in" },
163     { "abl",   "above below left first",              0, AV_OPT_TYPE_CONST, {.i64=ABOVE_BELOW_LR},     0, 0, FLAGS, "in" },
164     { "abr",   "above below right first",             0, AV_OPT_TYPE_CONST, {.i64=ABOVE_BELOW_RL},     0, 0, FLAGS, "in" },
165     { "al",    "alternating frames left first",       0, AV_OPT_TYPE_CONST, {.i64=ALTERNATING_LR},     0, 0, FLAGS, "in" },
166     { "ar",    "alternating frames right first",      0, AV_OPT_TYPE_CONST, {.i64=ALTERNATING_RL},     0, 0, FLAGS, "in" },
167     { "sbs2l", "side by side half width left first",  0, AV_OPT_TYPE_CONST, {.i64=SIDE_BY_SIDE_2_LR},  0, 0, FLAGS, "in" },
168     { "sbs2r", "side by side half width right first", 0, AV_OPT_TYPE_CONST, {.i64=SIDE_BY_SIDE_2_RL},  0, 0, FLAGS, "in" },
169     { "sbsl",  "side by side left first",             0, AV_OPT_TYPE_CONST, {.i64=SIDE_BY_SIDE_LR},    0, 0, FLAGS, "in" },
170     { "sbsr",  "side by side right first",            0, AV_OPT_TYPE_CONST, {.i64=SIDE_BY_SIDE_RL},    0, 0, FLAGS, "in" },
171     { "irl",   "interleave rows left first",          0, AV_OPT_TYPE_CONST, {.i64=INTERLEAVE_ROWS_LR}, 0, 0, FLAGS, "in" },
172     { "irr",   "interleave rows right first",         0, AV_OPT_TYPE_CONST, {.i64=INTERLEAVE_ROWS_RL}, 0, 0, FLAGS, "in" },
173     { "out",   "set output format", OFFSET(out.format),  AV_OPT_TYPE_INT,   {.i64=ANAGLYPH_RC_DUBOIS}, 0, STEREO_CODE_COUNT-1, FLAGS, "out"},
174     { "ab2l",  "above below half height left first",  0, AV_OPT_TYPE_CONST, {.i64=ABOVE_BELOW_2_LR},   0, 0, FLAGS, "out" },
175     { "ab2r",  "above below half height right first", 0, AV_OPT_TYPE_CONST, {.i64=ABOVE_BELOW_2_RL},   0, 0, FLAGS, "out" },
176     { "abl",   "above below left first",              0, AV_OPT_TYPE_CONST, {.i64=ABOVE_BELOW_LR},     0, 0, FLAGS, "out" },
177     { "abr",   "above below right first",             0, AV_OPT_TYPE_CONST, {.i64=ABOVE_BELOW_RL},     0, 0, FLAGS, "out" },
178     { "agmc",  "anaglyph green magenta color",        0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_GM_COLOR},  0, 0, FLAGS, "out" },
179     { "agmd",  "anaglyph green magenta dubois",       0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_GM_DUBOIS}, 0, 0, FLAGS, "out" },
180     { "agmg",  "anaglyph green magenta gray",         0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_GM_GRAY},   0, 0, FLAGS, "out" },
181     { "agmh",  "anaglyph green magenta half color",   0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_GM_HALF},   0, 0, FLAGS, "out" },
182     { "al",    "alternating frames left first",       0, AV_OPT_TYPE_CONST, {.i64=ALTERNATING_LR},     0, 0, FLAGS, "out" },
183     { "ar",    "alternating frames right first",      0, AV_OPT_TYPE_CONST, {.i64=ALTERNATING_RL},     0, 0, FLAGS, "out" },
184     { "arbg",  "anaglyph red blue gray",              0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_RB_GRAY},   0, 0, FLAGS, "out" },
185     { "arcc",  "anaglyph red cyan color",             0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_RC_COLOR},  0, 0, FLAGS, "out" },
186     { "arcd",  "anaglyph red cyan dubois",            0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_RC_DUBOIS}, 0, 0, FLAGS, "out" },
187     { "arcg",  "anaglyph red cyan gray",              0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_RC_GRAY},   0, 0, FLAGS, "out" },
188     { "arch",  "anaglyph red cyan half color",        0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_RC_HALF},   0, 0, FLAGS, "out" },
189     { "argg",  "anaglyph red green gray",             0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_RG_GRAY},   0, 0, FLAGS, "out" },
190     { "aybc",  "anaglyph yellow blue color",          0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_YB_COLOR},  0, 0, FLAGS, "out" },
191     { "aybd",  "anaglyph yellow blue dubois",         0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_YB_DUBOIS}, 0, 0, FLAGS, "out" },
192     { "aybg",  "anaglyph yellow blue gray",           0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_YB_GRAY},   0, 0, FLAGS, "out" },
193     { "aybh",  "anaglyph yellow blue half color",     0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_YB_HALF},   0, 0, FLAGS, "out" },
194     { "irl",   "interleave rows left first",          0, AV_OPT_TYPE_CONST, {.i64=INTERLEAVE_ROWS_LR}, 0, 0, FLAGS, "out" },
195     { "irr",   "interleave rows right first",         0, AV_OPT_TYPE_CONST, {.i64=INTERLEAVE_ROWS_RL}, 0, 0, FLAGS, "out" },
196     { "ml",    "mono left",                           0, AV_OPT_TYPE_CONST, {.i64=MONO_L},             0, 0, FLAGS, "out" },
197     { "mr",    "mono right",                          0, AV_OPT_TYPE_CONST, {.i64=MONO_R},             0, 0, FLAGS, "out" },
198     { "sbs2l", "side by side half width left first",  0, AV_OPT_TYPE_CONST, {.i64=SIDE_BY_SIDE_2_LR},  0, 0, FLAGS, "out" },
199     { "sbs2r", "side by side half width right first", 0, AV_OPT_TYPE_CONST, {.i64=SIDE_BY_SIDE_2_RL},  0, 0, FLAGS, "out" },
200     { "sbsl",  "side by side left first",             0, AV_OPT_TYPE_CONST, {.i64=SIDE_BY_SIDE_LR},    0, 0, FLAGS, "out" },
201     { "sbsr",  "side by side right first",            0, AV_OPT_TYPE_CONST, {.i64=SIDE_BY_SIDE_RL},    0, 0, FLAGS, "out" },
202     { "chl",   "checkerboard left first",             0, AV_OPT_TYPE_CONST, {.i64=CHECKERBOARD_LR},    0, 0, FLAGS, "out" },
203     { "chr",   "checkerboard right first",            0, AV_OPT_TYPE_CONST, {.i64=CHECKERBOARD_RL},    0, 0, FLAGS, "out" },
204     { "icl",   "interleave columns left first",       0, AV_OPT_TYPE_CONST, {.i64=INTERLEAVE_COLS_LR}, 0, 0, FLAGS, "out" },
205     { "icr",   "interleave columns right first",      0, AV_OPT_TYPE_CONST, {.i64=INTERLEAVE_COLS_RL}, 0, 0, FLAGS, "out" },
206     { "hdmi",  "HDMI frame pack",                     0, AV_OPT_TYPE_CONST, {.i64=HDMI},               0, 0, FLAGS, "out" },
207     { NULL }
208 };
209
210 AVFILTER_DEFINE_CLASS(stereo3d);
211
212 static const enum AVPixelFormat anaglyph_pix_fmts[] = {
213     AV_PIX_FMT_RGB24, AV_PIX_FMT_BGR24,
214     AV_PIX_FMT_NONE
215 };
216
217 static const enum AVPixelFormat other_pix_fmts[] = {
218     AV_PIX_FMT_RGB24, AV_PIX_FMT_BGR24,
219     AV_PIX_FMT_RGB48BE, AV_PIX_FMT_BGR48BE,
220     AV_PIX_FMT_RGB48LE, AV_PIX_FMT_BGR48LE,
221     AV_PIX_FMT_RGBA64BE, AV_PIX_FMT_BGRA64BE,
222     AV_PIX_FMT_RGBA64LE, AV_PIX_FMT_BGRA64LE,
223     AV_PIX_FMT_RGBA,  AV_PIX_FMT_BGRA,
224     AV_PIX_FMT_ARGB,  AV_PIX_FMT_ABGR,
225     AV_PIX_FMT_RGB0,  AV_PIX_FMT_BGR0,
226     AV_PIX_FMT_0RGB,  AV_PIX_FMT_0BGR,
227     AV_PIX_FMT_GBRP,
228     AV_PIX_FMT_GBRP9BE,  AV_PIX_FMT_GBRP9LE,
229     AV_PIX_FMT_GBRP10BE, AV_PIX_FMT_GBRP10LE,
230     AV_PIX_FMT_GBRP12BE, AV_PIX_FMT_GBRP12LE,
231     AV_PIX_FMT_GBRP14BE, AV_PIX_FMT_GBRP14LE,
232     AV_PIX_FMT_GBRP16BE, AV_PIX_FMT_GBRP16LE,
233     AV_PIX_FMT_YUV410P,
234     AV_PIX_FMT_YUV411P,
235     AV_PIX_FMT_YUV420P, AV_PIX_FMT_YUVA420P,
236     AV_PIX_FMT_YUV422P, AV_PIX_FMT_YUVA422P,
237     AV_PIX_FMT_YUV440P,
238     AV_PIX_FMT_YUV444P, AV_PIX_FMT_YUVA444P,
239     AV_PIX_FMT_YUVJ411P,
240     AV_PIX_FMT_YUVJ420P,
241     AV_PIX_FMT_YUVJ422P,
242     AV_PIX_FMT_YUVJ440P,
243     AV_PIX_FMT_YUVJ444P,
244     AV_PIX_FMT_YUV420P9LE,  AV_PIX_FMT_YUVA420P9LE,
245     AV_PIX_FMT_YUV420P9BE,  AV_PIX_FMT_YUVA420P9BE,
246     AV_PIX_FMT_YUV422P9LE,  AV_PIX_FMT_YUVA422P9LE,
247     AV_PIX_FMT_YUV422P9BE,  AV_PIX_FMT_YUVA422P9BE,
248     AV_PIX_FMT_YUV444P9LE,  AV_PIX_FMT_YUVA444P9LE,
249     AV_PIX_FMT_YUV444P9BE,  AV_PIX_FMT_YUVA444P9BE,
250     AV_PIX_FMT_YUV420P10LE, AV_PIX_FMT_YUVA420P10LE,
251     AV_PIX_FMT_YUV420P10BE, AV_PIX_FMT_YUVA420P10BE,
252     AV_PIX_FMT_YUV422P10LE, AV_PIX_FMT_YUVA422P10LE,
253     AV_PIX_FMT_YUV422P10BE, AV_PIX_FMT_YUVA422P10BE,
254     AV_PIX_FMT_YUV444P10LE, AV_PIX_FMT_YUVA444P10LE,
255     AV_PIX_FMT_YUV444P10BE, AV_PIX_FMT_YUVA444P10BE,
256     AV_PIX_FMT_YUV420P12BE,  AV_PIX_FMT_YUV420P12LE,
257     AV_PIX_FMT_YUV422P12BE,  AV_PIX_FMT_YUV422P12LE,
258     AV_PIX_FMT_YUV444P12BE,  AV_PIX_FMT_YUV444P12LE,
259     AV_PIX_FMT_YUV420P14BE,  AV_PIX_FMT_YUV420P14LE,
260     AV_PIX_FMT_YUV422P14BE,  AV_PIX_FMT_YUV422P14LE,
261     AV_PIX_FMT_YUV444P14BE,  AV_PIX_FMT_YUV444P14LE,
262     AV_PIX_FMT_YUV420P16LE, AV_PIX_FMT_YUVA420P16LE,
263     AV_PIX_FMT_YUV420P16BE, AV_PIX_FMT_YUVA420P16BE,
264     AV_PIX_FMT_YUV422P16LE, AV_PIX_FMT_YUVA422P16LE,
265     AV_PIX_FMT_YUV422P16BE, AV_PIX_FMT_YUVA422P16BE,
266     AV_PIX_FMT_YUV444P16LE, AV_PIX_FMT_YUVA444P16LE,
267     AV_PIX_FMT_YUV444P16BE, AV_PIX_FMT_YUVA444P16BE,
268     AV_PIX_FMT_NONE
269 };
270
271 static int query_formats(AVFilterContext *ctx)
272 {
273     Stereo3DContext *s = ctx->priv;
274     const enum AVPixelFormat *pix_fmts;
275     AVFilterFormats *fmts_list;
276
277     switch (s->out.format) {
278     case ANAGLYPH_GM_COLOR:
279     case ANAGLYPH_GM_DUBOIS:
280     case ANAGLYPH_GM_GRAY:
281     case ANAGLYPH_GM_HALF:
282     case ANAGLYPH_RB_GRAY:
283     case ANAGLYPH_RC_COLOR:
284     case ANAGLYPH_RC_DUBOIS:
285     case ANAGLYPH_RC_GRAY:
286     case ANAGLYPH_RC_HALF:
287     case ANAGLYPH_RG_GRAY:
288     case ANAGLYPH_YB_COLOR:
289     case ANAGLYPH_YB_DUBOIS:
290     case ANAGLYPH_YB_GRAY:
291     case ANAGLYPH_YB_HALF:
292         pix_fmts = anaglyph_pix_fmts;
293         break;
294     default:
295         pix_fmts = other_pix_fmts;
296     }
297
298     fmts_list = ff_make_format_list(pix_fmts);
299     if (!fmts_list)
300         return AVERROR(ENOMEM);
301     return ff_set_common_formats(ctx, fmts_list);
302 }
303
304 static inline uint8_t ana_convert(const int *coeff, const uint8_t *left, const uint8_t *right)
305 {
306     int sum;
307
308     sum  = coeff[0] * left[0] + coeff[3] * right[0]; //red in
309     sum += coeff[1] * left[1] + coeff[4] * right[1]; //green in
310     sum += coeff[2] * left[2] + coeff[5] * right[2]; //blue in
311
312     return av_clip_uint8(sum >> 16);
313 }
314
315 static void anaglyph(uint8_t *dst, uint8_t *lsrc, uint8_t *rsrc,
316                      ptrdiff_t dst_linesize, ptrdiff_t l_linesize, ptrdiff_t r_linesize,
317                      int width, int height,
318                      const int *ana_matrix_r, const int *ana_matrix_g, const int *ana_matrix_b)
319 {
320     int x, y, o;
321
322     for (y = 0; y < height; y++) {
323         for (o = 0, x = 0; x < width; x++, o+= 3) {
324             dst[o    ] = ana_convert(ana_matrix_r, lsrc + o, rsrc + o);
325             dst[o + 1] = ana_convert(ana_matrix_g, lsrc + o, rsrc + o);
326             dst[o + 2] = ana_convert(ana_matrix_b, lsrc + o, rsrc + o);
327         }
328
329         dst  += dst_linesize;
330         lsrc += l_linesize;
331         rsrc += r_linesize;
332     }
333 }
334
335 static int config_output(AVFilterLink *outlink)
336 {
337     AVFilterContext *ctx = outlink->src;
338     AVFilterLink *inlink = ctx->inputs[0];
339     Stereo3DContext *s = ctx->priv;
340     AVRational aspect = inlink->sample_aspect_ratio;
341     AVRational fps = inlink->frame_rate;
342     AVRational tb = inlink->time_base;
343     const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(outlink->format);
344     int ret;
345
346     switch (s->in.format) {
347     case SIDE_BY_SIDE_2_LR:
348     case SIDE_BY_SIDE_LR:
349     case SIDE_BY_SIDE_2_RL:
350     case SIDE_BY_SIDE_RL:
351         if (inlink->w & 1) {
352             av_log(ctx, AV_LOG_ERROR, "width must be even\n");
353             return AVERROR_INVALIDDATA;
354         }
355         break;
356     case INTERLEAVE_ROWS_LR:
357     case INTERLEAVE_ROWS_RL:
358     case ABOVE_BELOW_2_LR:
359     case ABOVE_BELOW_LR:
360     case ABOVE_BELOW_2_RL:
361     case ABOVE_BELOW_RL:
362         if (inlink->h & 1) {
363             av_log(ctx, AV_LOG_ERROR, "height must be even\n");
364             return AVERROR_INVALIDDATA;
365         }
366         break;
367     }
368
369     s->in.width     =
370     s->width        = inlink->w;
371     s->in.height    =
372     s->height       = inlink->h;
373     s->in.off_lstep =
374     s->in.off_rstep =
375     s->in.off_left  =
376     s->in.off_right =
377     s->in.row_left  =
378     s->in.row_right = 0;
379     s->in.row_step  = 1;
380
381     switch (s->in.format) {
382     case SIDE_BY_SIDE_2_LR:
383         aspect.num     *= 2;
384     case SIDE_BY_SIDE_LR:
385         s->width        = inlink->w / 2;
386         s->in.off_right = s->width;
387         break;
388     case SIDE_BY_SIDE_2_RL:
389         aspect.num     *= 2;
390     case SIDE_BY_SIDE_RL:
391         s->width        = inlink->w / 2;
392         s->in.off_left  = s->width;
393         break;
394     case ABOVE_BELOW_2_LR:
395         aspect.den     *= 2;
396     case ABOVE_BELOW_LR:
397         s->in.row_right =
398         s->height       = inlink->h / 2;
399         break;
400     case ABOVE_BELOW_2_RL:
401         aspect.den     *= 2;
402     case ABOVE_BELOW_RL:
403         s->in.row_left  =
404         s->height       = inlink->h / 2;
405         break;
406     case ALTERNATING_RL:
407     case ALTERNATING_LR:
408         fps.den        *= 2;
409         tb.num         *= 2;
410         break;
411     case INTERLEAVE_ROWS_LR:
412     case INTERLEAVE_ROWS_RL:
413         s->in.row_step  = 2;
414         if (s->in.format == INTERLEAVE_ROWS_RL)
415             s->in.off_lstep = 1;
416         else
417             s->in.off_rstep = 1;
418         if (s->out.format != CHECKERBOARD_LR &&
419             s->out.format != CHECKERBOARD_RL)
420             s->height   = inlink->h / 2;
421         break;
422     default:
423         av_log(ctx, AV_LOG_ERROR, "input format %d is not supported\n", s->in.format);
424         return AVERROR(EINVAL);
425     }
426
427     s->out.width     = s->width;
428     s->out.height    = s->height;
429     s->out.off_lstep =
430     s->out.off_rstep =
431     s->out.off_left  =
432     s->out.off_right =
433     s->out.row_left  =
434     s->out.row_right = 0;
435     s->out.row_step  = 1;
436
437     switch (s->out.format) {
438     case ANAGLYPH_RB_GRAY:
439     case ANAGLYPH_RG_GRAY:
440     case ANAGLYPH_RC_GRAY:
441     case ANAGLYPH_RC_HALF:
442     case ANAGLYPH_RC_COLOR:
443     case ANAGLYPH_RC_DUBOIS:
444     case ANAGLYPH_GM_GRAY:
445     case ANAGLYPH_GM_HALF:
446     case ANAGLYPH_GM_COLOR:
447     case ANAGLYPH_GM_DUBOIS:
448     case ANAGLYPH_YB_GRAY:
449     case ANAGLYPH_YB_HALF:
450     case ANAGLYPH_YB_COLOR:
451     case ANAGLYPH_YB_DUBOIS: {
452         uint8_t rgba_map[4];
453
454         ff_fill_rgba_map(rgba_map, outlink->format);
455         s->ana_matrix[rgba_map[0]] = &ana_coeff[s->out.format][0][0];
456         s->ana_matrix[rgba_map[1]] = &ana_coeff[s->out.format][1][0];
457         s->ana_matrix[rgba_map[2]] = &ana_coeff[s->out.format][2][0];
458         break;
459     }
460     case SIDE_BY_SIDE_2_LR:
461         aspect.den      *= 2;
462     case SIDE_BY_SIDE_LR:
463         s->out.width     = s->width * 2;
464         s->out.off_right = s->width;
465         break;
466     case SIDE_BY_SIDE_2_RL:
467         aspect.den      *= 2;
468     case SIDE_BY_SIDE_RL:
469         s->out.width     = s->width * 2;
470         s->out.off_left  = s->width;
471         break;
472     case ABOVE_BELOW_2_LR:
473         aspect.num      *= 2;
474     case ABOVE_BELOW_LR:
475         s->out.height    = s->height * 2;
476         s->out.row_right = s->height;
477         break;
478     case HDMI:
479         if (s->height != 720 && s->height != 1080) {
480             av_log(ctx, AV_LOG_ERROR, "Only 720 and 1080 height supported\n");
481             return AVERROR(EINVAL);
482         }
483
484         s->blanks = s->height / 24;
485         s->out.height    = s->height * 2 + s->blanks;
486         s->out.row_right = s->height + s->blanks;
487         break;
488     case ABOVE_BELOW_2_RL:
489         aspect.num      *= 2;
490     case ABOVE_BELOW_RL:
491         s->out.height    = s->height * 2;
492         s->out.row_left  = s->height;
493         break;
494     case INTERLEAVE_ROWS_LR:
495         s->in.row_step   = 1 + (s->in.format == INTERLEAVE_ROWS_RL);
496         s->out.row_step  = 2;
497         s->out.height    = s->height * 2;
498         s->out.off_rstep = 1;
499         break;
500     case INTERLEAVE_ROWS_RL:
501         s->in.row_step   = 1 + (s->in.format == INTERLEAVE_ROWS_LR);
502         s->out.row_step  = 2;
503         s->out.height    = s->height * 2;
504         s->out.off_lstep = 1;
505         break;
506     case MONO_R:
507         s->in.off_left   = s->in.off_right;
508         s->in.row_left   = s->in.row_right;
509         if (s->in.format == INTERLEAVE_ROWS_LR)
510             FFSWAP(int, s->in.off_lstep, s->in.off_rstep);
511         break;
512     case MONO_L:
513         if (s->in.format == INTERLEAVE_ROWS_RL)
514             FFSWAP(int, s->in.off_lstep, s->in.off_rstep);
515         break;
516     case ALTERNATING_RL:
517     case ALTERNATING_LR:
518         fps.num         *= 2;
519         tb.den          *= 2;
520         break;
521     case CHECKERBOARD_LR:
522     case CHECKERBOARD_RL:
523         s->out.width     = s->width * 2;
524         break;
525     case INTERLEAVE_COLS_LR:
526     case INTERLEAVE_COLS_RL:
527         s->out.width     = s->width * 2;
528         break;
529     default:
530         av_log(ctx, AV_LOG_ERROR, "output format %d is not supported\n", s->out.format);
531         return AVERROR(EINVAL);
532     }
533
534     outlink->w = s->out.width;
535     outlink->h = s->out.height;
536     outlink->frame_rate = fps;
537     outlink->time_base = tb;
538     outlink->sample_aspect_ratio = aspect;
539
540     if ((ret = av_image_fill_linesizes(s->linesize, outlink->format, s->width)) < 0)
541         return ret;
542     s->nb_planes = av_pix_fmt_count_planes(outlink->format);
543     av_image_fill_max_pixsteps(s->pixstep, NULL, desc);
544     s->pheight[1] = s->pheight[2] = FF_CEIL_RSHIFT(s->height, desc->log2_chroma_h);
545     s->pheight[0] = s->pheight[3] = s->height;
546     s->hsub = desc->log2_chroma_w;
547     s->vsub = desc->log2_chroma_h;
548
549     s->dsp.anaglyph = anaglyph;
550     if (ARCH_X86)
551         ff_stereo3d_init_x86(&s->dsp);
552
553     return 0;
554 }
555
556 typedef struct ThreadData {
557     AVFrame *ileft, *iright;
558     AVFrame *out;
559 } ThreadData;
560
561 static int filter_slice(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
562 {
563     Stereo3DContext *s = ctx->priv;
564     ThreadData *td = arg;
565     AVFrame *ileft = td->ileft;
566     AVFrame *iright = td->iright;
567     AVFrame *out = td->out;
568     int height = s->out.height;
569     int start = (height *  jobnr   ) / nb_jobs;
570     int end   = (height * (jobnr+1)) / nb_jobs;
571     const int **ana_matrix = s->ana_matrix;
572
573     s->dsp.anaglyph(out->data[0] + out->linesize[0] * start,
574              ileft ->data[0] + s->in_off_left [0]  + ileft->linesize[0] * start * s->in.row_step,
575              iright->data[0] + s->in_off_right[0] + iright->linesize[0] * start * s->in.row_step,
576              out->linesize[0],
577              ileft->linesize[0] * s->in.row_step,
578              iright->linesize[0] * s->in.row_step,
579              s->out.width, end - start,
580              ana_matrix[0], ana_matrix[1], ana_matrix[2]);
581
582     return 0;
583 }
584
585 static int filter_frame(AVFilterLink *inlink, AVFrame *inpicref)
586 {
587     AVFilterContext *ctx  = inlink->dst;
588     Stereo3DContext *s = ctx->priv;
589     AVFilterLink *outlink = ctx->outputs[0];
590     AVFrame *out, *oleft, *oright, *ileft, *iright;
591     int out_off_left[4], out_off_right[4];
592     int i, ret;
593
594     if (s->in.format == s->out.format)
595         return ff_filter_frame(outlink, inpicref);
596
597     switch (s->out.format) {
598     case ALTERNATING_LR:
599     case ALTERNATING_RL:
600         if (!s->prev) {
601             s->prev = inpicref;
602             return 0;
603         }
604         break;
605     };
606
607     switch (s->in.format) {
608     case ALTERNATING_LR:
609     case ALTERNATING_RL:
610         if (!s->prev) {
611             s->prev = inpicref;
612             return 0;
613         }
614         ileft  = s->prev;
615         iright = inpicref;
616         if (s->in.format == ALTERNATING_RL)
617             FFSWAP(AVFrame *, ileft, iright);
618         break;
619     default:
620         ileft = iright = inpicref;
621     };
622
623     if ((s->out.format == ALTERNATING_LR ||
624          s->out.format == ALTERNATING_RL) &&
625         (s->in.format == SIDE_BY_SIDE_LR ||
626          s->in.format == SIDE_BY_SIDE_RL ||
627          s->in.format == SIDE_BY_SIDE_2_LR ||
628          s->in.format == SIDE_BY_SIDE_2_RL ||
629          s->in.format == ABOVE_BELOW_LR ||
630          s->in.format == ABOVE_BELOW_RL ||
631          s->in.format == ABOVE_BELOW_2_LR ||
632          s->in.format == ABOVE_BELOW_2_RL ||
633          s->in.format == INTERLEAVE_ROWS_LR ||
634          s->in.format == INTERLEAVE_ROWS_RL)) {
635         oright = av_frame_clone(s->prev);
636         oleft  = av_frame_clone(s->prev);
637         if (!oright || !oleft) {
638             av_frame_free(&oright);
639             av_frame_free(&oleft);
640             av_frame_free(&s->prev);
641             av_frame_free(&inpicref);
642             return AVERROR(ENOMEM);
643         }
644     } else if ((s->out.format == MONO_L ||
645                 s->out.format == MONO_R) &&
646         (s->in.format == SIDE_BY_SIDE_LR ||
647          s->in.format == SIDE_BY_SIDE_RL ||
648          s->in.format == SIDE_BY_SIDE_2_LR ||
649          s->in.format == SIDE_BY_SIDE_2_RL ||
650          s->in.format == ABOVE_BELOW_LR ||
651          s->in.format == ABOVE_BELOW_RL ||
652          s->in.format == ABOVE_BELOW_2_LR ||
653          s->in.format == ABOVE_BELOW_2_RL ||
654          s->in.format == INTERLEAVE_ROWS_LR ||
655          s->in.format == INTERLEAVE_ROWS_RL)) {
656         out = oleft = oright = av_frame_clone(inpicref);
657         if (!out) {
658             av_frame_free(&s->prev);
659             av_frame_free(&inpicref);
660             return AVERROR(ENOMEM);
661         }
662     } else if ((s->out.format == MONO_L && s->in.format == ALTERNATING_LR) ||
663                (s->out.format == MONO_R && s->in.format == ALTERNATING_RL)) {
664         s->prev->pts /= 2;
665         ret = ff_filter_frame(outlink, s->prev);
666         av_frame_free(&inpicref);
667         s->prev = NULL;
668         return ret;
669     } else if ((s->out.format == MONO_L && s->in.format == ALTERNATING_RL) ||
670                (s->out.format == MONO_R && s->in.format == ALTERNATING_LR)) {
671         av_frame_free(&s->prev);
672         inpicref->pts /= 2;
673         return ff_filter_frame(outlink, inpicref);
674     } else if ((s->out.format == ALTERNATING_LR && s->in.format == ALTERNATING_RL) ||
675                (s->out.format == ALTERNATING_RL && s->in.format == ALTERNATING_LR)) {
676         FFSWAP(int64_t, s->prev->pts, inpicref->pts);
677         ff_filter_frame(outlink, inpicref);
678         ret = ff_filter_frame(outlink, s->prev);
679         s->prev = NULL;
680         return ret;
681     } else {
682         out = oleft = oright = ff_get_video_buffer(outlink, outlink->w, outlink->h);
683         if (!out) {
684             av_frame_free(&s->prev);
685             av_frame_free(&inpicref);
686             return AVERROR(ENOMEM);
687         }
688         av_frame_copy_props(out, inpicref);
689
690         if (s->out.format == ALTERNATING_LR ||
691             s->out.format == ALTERNATING_RL) {
692             oright = ff_get_video_buffer(outlink, outlink->w, outlink->h);
693             if (!oright) {
694                 av_frame_free(&oleft);
695                 av_frame_free(&s->prev);
696                 av_frame_free(&inpicref);
697                 return AVERROR(ENOMEM);
698             }
699             av_frame_copy_props(oright, s->prev);
700         }
701     }
702
703     for (i = 0; i < 4; i++) {
704         int hsub = i == 1 || i == 2 ? s->hsub : 0;
705         int vsub = i == 1 || i == 2 ? s->vsub : 0;
706         s->in_off_left[i]   = (FF_CEIL_RSHIFT(s->in.row_left,   vsub) + s->in.off_lstep)  * ileft->linesize[i]  + FF_CEIL_RSHIFT(s->in.off_left   * s->pixstep[i], hsub);
707         s->in_off_right[i]  = (FF_CEIL_RSHIFT(s->in.row_right,  vsub) + s->in.off_rstep)  * iright->linesize[i] + FF_CEIL_RSHIFT(s->in.off_right  * s->pixstep[i], hsub);
708         out_off_left[i]  = (FF_CEIL_RSHIFT(s->out.row_left,  vsub) + s->out.off_lstep) * oleft->linesize[i]  + FF_CEIL_RSHIFT(s->out.off_left  * s->pixstep[i], hsub);
709         out_off_right[i] = (FF_CEIL_RSHIFT(s->out.row_right, vsub) + s->out.off_rstep) * oright->linesize[i] + FF_CEIL_RSHIFT(s->out.off_right * s->pixstep[i], hsub);
710     }
711
712     switch (s->out.format) {
713     case ALTERNATING_LR:
714     case ALTERNATING_RL:
715         switch (s->in.format) {
716         case INTERLEAVE_ROWS_LR:
717         case INTERLEAVE_ROWS_RL:
718             for (i = 0; i < s->nb_planes; i++) {
719                 oleft->linesize[i]  *= 2;
720                 oright->linesize[i] *= 2;
721             }
722         case ABOVE_BELOW_LR:
723         case ABOVE_BELOW_RL:
724         case ABOVE_BELOW_2_LR:
725         case ABOVE_BELOW_2_RL:
726         case SIDE_BY_SIDE_LR:
727         case SIDE_BY_SIDE_RL:
728         case SIDE_BY_SIDE_2_LR:
729         case SIDE_BY_SIDE_2_RL:
730             oleft->width   = outlink->w;
731             oright->width  = outlink->w;
732             oleft->height  = outlink->h;
733             oright->height = outlink->h;
734
735             for (i = 0; i < s->nb_planes; i++) {
736                 oleft->data[i]  += s->in_off_left[i];
737                 oright->data[i] += s->in_off_right[i];
738             }
739             break;
740         default:
741             goto copy;
742             break;
743         }
744         break;
745     case HDMI:
746         for (i = 0; i < s->nb_planes; i++) {
747             int j, h = s->height >> ((i == 1 || i == 2) ? s->vsub : 0);
748             int b = (s->blanks) >> ((i == 1 || i == 2) ? s->vsub : 0);
749
750             for (j = h; j < h + b; j++)
751                 memset(oleft->data[i] + j * s->linesize[i], 0, s->linesize[i]);
752         }
753     case SIDE_BY_SIDE_LR:
754     case SIDE_BY_SIDE_RL:
755     case SIDE_BY_SIDE_2_LR:
756     case SIDE_BY_SIDE_2_RL:
757     case ABOVE_BELOW_LR:
758     case ABOVE_BELOW_RL:
759     case ABOVE_BELOW_2_LR:
760     case ABOVE_BELOW_2_RL:
761     case INTERLEAVE_ROWS_LR:
762     case INTERLEAVE_ROWS_RL:
763 copy:
764         for (i = 0; i < s->nb_planes; i++) {
765             av_image_copy_plane(oleft->data[i] + out_off_left[i],
766                                 oleft->linesize[i] * s->out.row_step,
767                                 ileft->data[i] + s->in_off_left[i],
768                                 ileft->linesize[i] * s->in.row_step,
769                                 s->linesize[i], s->pheight[i]);
770             av_image_copy_plane(oright->data[i] + out_off_right[i],
771                                 oright->linesize[i] * s->out.row_step,
772                                 iright->data[i] + s->in_off_right[i],
773                                 iright->linesize[i] * s->in.row_step,
774                                 s->linesize[i], s->pheight[i]);
775         }
776         break;
777     case MONO_L:
778         iright = ileft;
779     case MONO_R:
780         switch (s->in.format) {
781         case INTERLEAVE_ROWS_LR:
782         case INTERLEAVE_ROWS_RL:
783             for (i = 0; i < s->nb_planes; i++) {
784                 out->linesize[i] *= 2;
785             }
786         case ABOVE_BELOW_LR:
787         case ABOVE_BELOW_RL:
788         case ABOVE_BELOW_2_LR:
789         case ABOVE_BELOW_2_RL:
790         case SIDE_BY_SIDE_LR:
791         case SIDE_BY_SIDE_RL:
792         case SIDE_BY_SIDE_2_LR:
793         case SIDE_BY_SIDE_2_RL:
794             out->width  = outlink->w;
795             out->height = outlink->h;
796
797             for (i = 0; i < s->nb_planes; i++) {
798                 out->data[i] += s->in_off_left[i];
799             }
800             break;
801         default:
802             for (i = 0; i < s->nb_planes; i++) {
803                 av_image_copy_plane(out->data[i], out->linesize[i],
804                                     iright->data[i] + s->in_off_left[i],
805                                     iright->linesize[i] * s->in.row_step,
806                                     s->linesize[i], s->pheight[i]);
807             }
808             break;
809         }
810         break;
811     case ANAGLYPH_RB_GRAY:
812     case ANAGLYPH_RG_GRAY:
813     case ANAGLYPH_RC_GRAY:
814     case ANAGLYPH_RC_HALF:
815     case ANAGLYPH_RC_COLOR:
816     case ANAGLYPH_RC_DUBOIS:
817     case ANAGLYPH_GM_GRAY:
818     case ANAGLYPH_GM_HALF:
819     case ANAGLYPH_GM_COLOR:
820     case ANAGLYPH_GM_DUBOIS:
821     case ANAGLYPH_YB_GRAY:
822     case ANAGLYPH_YB_HALF:
823     case ANAGLYPH_YB_COLOR:
824     case ANAGLYPH_YB_DUBOIS: {
825         ThreadData td;
826
827         td.ileft = ileft; td.iright = iright; td.out = out;
828         ctx->internal->execute(ctx, filter_slice, &td, NULL,
829                                FFMIN(s->out.height, ctx->graph->nb_threads));
830         break;
831     }
832     case CHECKERBOARD_RL:
833     case CHECKERBOARD_LR:
834         for (i = 0; i < s->nb_planes; i++) {
835             int x, y;
836
837             for (y = 0; y < s->pheight[i]; y++) {
838                 uint8_t *dst = out->data[i] + out->linesize[i] * y;
839                 uint8_t *left = ileft->data[i] + ileft->linesize[i] * y + s->in_off_left[i];
840                 uint8_t *right = iright->data[i] + iright->linesize[i] * y + s->in_off_right[i];
841                 int p, b;
842
843                 if (s->out.format == CHECKERBOARD_RL)
844                     FFSWAP(uint8_t*, left, right);
845                 switch (s->pixstep[i]) {
846                 case 1:
847                     for (x = 0, b = 0, p = 0; x < s->linesize[i] * 2; x+=2, p++, b++) {
848                         dst[x  ] = (b&1) == (y&1) ? left[p] : right[p];
849                         dst[x+1] = (b&1) != (y&1) ? left[p] : right[p];
850                     }
851                     break;
852                 case 2:
853                     for (x = 0, b = 0, p = 0; x < s->linesize[i] * 2; x+=4, p+=2, b++) {
854                         AV_WN16(&dst[x  ], (b&1) == (y&1) ? AV_RN16(&left[p]) : AV_RN16(&right[p]));
855                         AV_WN16(&dst[x+2], (b&1) != (y&1) ? AV_RN16(&left[p]) : AV_RN16(&right[p]));
856                     }
857                     break;
858                 case 3:
859                     for (x = 0, b = 0, p = 0; x < s->linesize[i] * 2; x+=6, p+=3, b++) {
860                         AV_WB24(&dst[x  ], (b&1) == (y&1) ? AV_RB24(&left[p]) : AV_RB24(&right[p]));
861                         AV_WB24(&dst[x+3], (b&1) != (y&1) ? AV_RB24(&left[p]) : AV_RB24(&right[p]));
862                     }
863                     break;
864                 case 4:
865                     for (x = 0, b = 0, p = 0; x < s->linesize[i] * 2; x+=8, p+=4, b++) {
866                         AV_WN32(&dst[x  ], (b&1) == (y&1) ? AV_RN32(&left[p]) : AV_RN32(&right[p]));
867                         AV_WN32(&dst[x+4], (b&1) != (y&1) ? AV_RN32(&left[p]) : AV_RN32(&right[p]));
868                     }
869                     break;
870                 case 6:
871                     for (x = 0, b = 0, p = 0; x < s->linesize[i] * 2; x+=12, p+=6, b++) {
872                         AV_WB48(&dst[x  ], (b&1) == (y&1) ? AV_RB48(&left[p]) : AV_RB48(&right[p]));
873                         AV_WB48(&dst[x+6], (b&1) != (y&1) ? AV_RB48(&left[p]) : AV_RB48(&right[p]));
874                     }
875                     break;
876                 case 8:
877                     for (x = 0, b = 0, p = 0; x < s->linesize[i] * 2; x+=16, p+=8, b++) {
878                         AV_WN64(&dst[x  ], (b&1) == (y&1) ? AV_RN64(&left[p]) : AV_RN64(&right[p]));
879                         AV_WN64(&dst[x+8], (b&1) != (y&1) ? AV_RN64(&left[p]) : AV_RN64(&right[p]));
880                     }
881                     break;
882                 }
883             }
884         }
885         break;
886     case INTERLEAVE_COLS_LR:
887     case INTERLEAVE_COLS_RL:
888         for (i = 0; i < s->nb_planes; i++) {
889             int x, y;
890
891             for (y = 0; y < s->pheight[i]; y++) {
892                 uint8_t *dst = out->data[i] + out->linesize[i] * y;
893                 uint8_t *left = ileft->data[i] + ileft->linesize[i] * y * s->in.row_step + s->in_off_left[i];
894                 uint8_t *right = iright->data[i] + iright->linesize[i] * y * s->in.row_step + s->in_off_right[i];
895                 int p, b;
896
897                 if (s->out.format == INTERLEAVE_COLS_LR)
898                     FFSWAP(uint8_t*, left, right);
899
900                 switch (s->pixstep[i]) {
901                 case 1:
902                     for (x = 0, b = 0, p = 0; x < s->linesize[i] * 2; x+=2, p++, b+=2) {
903                         dst[x  ] =   b&1  ? left[p] : right[p];
904                         dst[x+1] = !(b&1) ? left[p] : right[p];
905                     }
906                     break;
907                 case 2:
908                     for (x = 0, b = 0, p = 0; x < s->linesize[i] * 2; x+=4, p+=2, b+=2) {
909                         AV_WN16(&dst[x  ],   b&1  ? AV_RN16(&left[p]) : AV_RN16(&right[p]));
910                         AV_WN16(&dst[x+2], !(b&1) ? AV_RN16(&left[p]) : AV_RN16(&right[p]));
911                     }
912                     break;
913                 case 3:
914                     for (x = 0, b = 0, p = 0; x < s->linesize[i] * 2; x+=6, p+=3, b+=2) {
915                         AV_WB24(&dst[x  ],   b&1  ? AV_RB24(&left[p]) : AV_RB24(&right[p]));
916                         AV_WB24(&dst[x+3], !(b&1) ? AV_RB24(&left[p]) : AV_RB24(&right[p]));
917                     }
918                     break;
919                 case 4:
920                     for (x = 0, b = 0, p = 0; x < s->linesize[i] * 2; x+=8, p+=4, b+=2) {
921                         AV_WN32(&dst[x  ],   b&1  ? AV_RN32(&left[p]) : AV_RN32(&right[p]));
922                         AV_WN32(&dst[x+4], !(b&1) ? AV_RN32(&left[p]) : AV_RN32(&right[p]));
923                     }
924                     break;
925                 case 6:
926                     for (x = 0, b = 0, p = 0; x < s->linesize[i] * 2; x+=12, p+=6, b+=2) {
927                         AV_WB48(&dst[x  ],   b&1  ? AV_RB48(&left[p]) : AV_RB48(&right[p]));
928                         AV_WB48(&dst[x+6], !(b&1) ? AV_RB48(&left[p]) : AV_RB48(&right[p]));
929                     }
930                     break;
931                 case 8:
932                     for (x = 0, b = 0, p = 0; x < s->linesize[i] * 2; x+=16, p+=8, b+=2) {
933                         AV_WN64(&dst[x  ],   b&1 ?  AV_RN64(&left[p]) : AV_RN64(&right[p]));
934                         AV_WN64(&dst[x+8], !(b&1) ? AV_RN64(&left[p]) : AV_RN64(&right[p]));
935                     }
936                     break;
937                 }
938             }
939         }
940         break;
941     default:
942         av_assert0(0);
943     }
944
945     if (oright != oleft) {
946         if (s->out.format == ALTERNATING_LR)
947             FFSWAP(AVFrame *, oleft, oright);
948         oright->pts = s->prev->pts * 2;
949         ff_filter_frame(outlink, oright);
950         out = oleft;
951         oleft->pts = s->prev->pts + inpicref->pts;
952         av_frame_free(&s->prev);
953         s->prev = inpicref;
954     } else if (s->in.format == ALTERNATING_LR ||
955                s->in.format == ALTERNATING_RL) {
956         out->pts = s->prev->pts / 2;
957         av_frame_free(&s->prev);
958         av_frame_free(&inpicref);
959     } else {
960         av_frame_free(&s->prev);
961         av_frame_free(&inpicref);
962     }
963     return ff_filter_frame(outlink, out);
964 }
965
966 static av_cold void uninit(AVFilterContext *ctx)
967 {
968     Stereo3DContext *s = ctx->priv;
969
970     av_frame_free(&s->prev);
971 }
972
973 static const AVFilterPad stereo3d_inputs[] = {
974     {
975         .name         = "default",
976         .type         = AVMEDIA_TYPE_VIDEO,
977         .filter_frame = filter_frame,
978     },
979     { NULL }
980 };
981
982 static const AVFilterPad stereo3d_outputs[] = {
983     {
984         .name         = "default",
985         .type         = AVMEDIA_TYPE_VIDEO,
986         .config_props = config_output,
987     },
988     { NULL }
989 };
990
991 AVFilter ff_vf_stereo3d = {
992     .name          = "stereo3d",
993     .description   = NULL_IF_CONFIG_SMALL("Convert video stereoscopic 3D view."),
994     .priv_size     = sizeof(Stereo3DContext),
995     .uninit        = uninit,
996     .query_formats = query_formats,
997     .inputs        = stereo3d_inputs,
998     .outputs       = stereo3d_outputs,
999     .priv_class    = &stereo3d_class,
1000     .flags         = AVFILTER_FLAG_SLICE_THREADS,
1001 };