Mercurial > libavcodec.hg
annotate svq3.c @ 3716:e583349b9720 libavcodec
Proper armv5te instructions support detection.
patch by Siarhei Siamashka, siarhei dot siamashka gmail com
author | diego |
---|---|
date | Wed, 13 Sep 2006 21:39:41 +0000 |
parents | 0b546eab515d |
children | c8c591fe26f8 |
rev | line source |
---|---|
1234 | 1 /* |
2 * Copyright (c) 2003 The FFmpeg Project. | |
3 * | |
4 * This library is free software; you can redistribute it and/or | |
5 * modify it under the terms of the GNU Lesser General Public | |
6 * License as published by the Free Software Foundation; either | |
7 * version 2 of the License, or (at your option) any later version. | |
8 * | |
9 * This library is distributed in the hope that it will be useful, | |
10 * but WITHOUT ANY WARRANTY; without even the implied warranty of | |
11 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU | |
12 * Lesser General Public License for more details. | |
13 * | |
14 * You should have received a copy of the GNU Lesser General Public | |
15 * License along with this library; if not, write to the Free Software | |
3036
0b546eab515d
Update licensing information: The FSF changed postal address.
diego
parents:
2979
diff
changeset
|
16 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA |
1234 | 17 * |
18 * | |
19 * How to use this decoder: | |
20 * SVQ3 data is transported within Apple Quicktime files. Quicktime files | |
1235
8bb75c3c2f21
change the way the ImageDescription is passed to the decoder
tmmm
parents:
1234
diff
changeset
|
21 * have stsd atoms to describe media trak properties. A stsd atom for a |
8bb75c3c2f21
change the way the ImageDescription is passed to the decoder
tmmm
parents:
1234
diff
changeset
|
22 * video trak contains 1 or more ImageDescription atoms. These atoms begin |
8bb75c3c2f21
change the way the ImageDescription is passed to the decoder
tmmm
parents:
1234
diff
changeset
|
23 * with the 4-byte length of the atom followed by the codec fourcc. Some |
8bb75c3c2f21
change the way the ImageDescription is passed to the decoder
tmmm
parents:
1234
diff
changeset
|
24 * decoders need information in this atom to operate correctly. Such |
8bb75c3c2f21
change the way the ImageDescription is passed to the decoder
tmmm
parents:
1234
diff
changeset
|
25 * is the case with SVQ3. In order to get the best use out of this decoder, |
8bb75c3c2f21
change the way the ImageDescription is passed to the decoder
tmmm
parents:
1234
diff
changeset
|
26 * the calling app must make the SVQ3 ImageDescription atom available |
1234 | 27 * via the AVCodecContext's extradata[_size] field: |
28 * | |
2967 | 29 * AVCodecContext.extradata = pointer to ImageDescription, first characters |
1235
8bb75c3c2f21
change the way the ImageDescription is passed to the decoder
tmmm
parents:
1234
diff
changeset
|
30 * are expected to be 'S', 'V', 'Q', and '3', NOT the 4-byte atom length |
2967 | 31 * AVCodecContext.extradata_size = size of ImageDescription atom memory |
32 * buffer (which will be the same as the ImageDescription atom size field | |
1235
8bb75c3c2f21
change the way the ImageDescription is passed to the decoder
tmmm
parents:
1234
diff
changeset
|
33 * from the QT file, minus 4 bytes since the length is missing) |
8bb75c3c2f21
change the way the ImageDescription is passed to the decoder
tmmm
parents:
1234
diff
changeset
|
34 * |
8bb75c3c2f21
change the way the ImageDescription is passed to the decoder
tmmm
parents:
1234
diff
changeset
|
35 * You will know you have these parameters passed correctly when the decoder |
8bb75c3c2f21
change the way the ImageDescription is passed to the decoder
tmmm
parents:
1234
diff
changeset
|
36 * correctly decodes this file: |
8bb75c3c2f21
change the way the ImageDescription is passed to the decoder
tmmm
parents:
1234
diff
changeset
|
37 * ftp://ftp.mplayerhq.hu/MPlayer/samples/V-codecs/SVQ3/Vertical400kbit.sorenson3.mov |
1234 | 38 * |
39 */ | |
2967 | 40 |
1234 | 41 /** |
42 * @file svq3.c | |
43 * svq3 decoder. | |
44 */ | |
45 | |
2967 | 46 #define FULLPEL_MODE 1 |
47 #define HALFPEL_MODE 2 | |
1265 | 48 #define THIRDPEL_MODE 3 |
1330
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
49 #define PREDICT_MODE 4 |
2967 | 50 |
1252 | 51 /* dual scan (from some older h264 draft) |
52 o-->o-->o o | |
53 | /| | |
54 o o o / o | |
55 | / | |/ | | |
56 o o o o | |
2967 | 57 / |
1252 | 58 o-->o-->o-->o |
59 */ | |
1234 | 60 static const uint8_t svq3_scan[16]={ |
61 0+0*4, 1+0*4, 2+0*4, 2+1*4, | |
62 2+2*4, 3+0*4, 3+1*4, 3+2*4, | |
63 0+1*4, 0+2*4, 1+1*4, 1+2*4, | |
64 0+3*4, 1+3*4, 2+3*4, 3+3*4, | |
65 }; | |
66 | |
67 static const uint8_t svq3_pred_0[25][2] = { | |
68 { 0, 0 }, | |
69 { 1, 0 }, { 0, 1 }, | |
70 { 0, 2 }, { 1, 1 }, { 2, 0 }, | |
71 { 3, 0 }, { 2, 1 }, { 1, 2 }, { 0, 3 }, | |
72 { 0, 4 }, { 1, 3 }, { 2, 2 }, { 3, 1 }, { 4, 0 }, | |
73 { 4, 1 }, { 3, 2 }, { 2, 3 }, { 1, 4 }, | |
74 { 2, 4 }, { 3, 3 }, { 4, 2 }, | |
75 { 4, 3 }, { 3, 4 }, | |
76 { 4, 4 } | |
77 }; | |
78 | |
79 static const int8_t svq3_pred_1[6][6][5] = { | |
80 { { 2,-1,-1,-1,-1 }, { 2, 1,-1,-1,-1 }, { 1, 2,-1,-1,-1 }, | |
81 { 2, 1,-1,-1,-1 }, { 1, 2,-1,-1,-1 }, { 1, 2,-1,-1,-1 } }, | |
82 { { 0, 2,-1,-1,-1 }, { 0, 2, 1, 4, 3 }, { 0, 1, 2, 4, 3 }, | |
83 { 0, 2, 1, 4, 3 }, { 2, 0, 1, 3, 4 }, { 0, 4, 2, 1, 3 } }, | |
84 { { 2, 0,-1,-1,-1 }, { 2, 1, 0, 4, 3 }, { 1, 2, 4, 0, 3 }, | |
85 { 2, 1, 0, 4, 3 }, { 2, 1, 4, 3, 0 }, { 1, 2, 4, 0, 3 } }, | |
86 { { 2, 0,-1,-1,-1 }, { 2, 0, 1, 4, 3 }, { 1, 2, 0, 4, 3 }, | |
87 { 2, 1, 0, 4, 3 }, { 2, 1, 3, 4, 0 }, { 2, 4, 1, 0, 3 } }, | |
88 { { 0, 2,-1,-1,-1 }, { 0, 2, 1, 3, 4 }, { 1, 2, 3, 0, 4 }, | |
89 { 2, 0, 1, 3, 4 }, { 2, 1, 3, 0, 4 }, { 2, 0, 4, 3, 1 } }, | |
90 { { 0, 2,-1,-1,-1 }, { 0, 2, 4, 1, 3 }, { 1, 4, 2, 0, 3 }, | |
91 { 4, 2, 0, 1, 3 }, { 2, 0, 1, 4, 3 }, { 4, 2, 1, 0, 3 } }, | |
92 }; | |
93 | |
94 static const struct { uint8_t run; uint8_t level; } svq3_dct_tables[2][16] = { | |
95 { { 0, 0 }, { 0, 1 }, { 1, 1 }, { 2, 1 }, { 0, 2 }, { 3, 1 }, { 4, 1 }, { 5, 1 }, | |
96 { 0, 3 }, { 1, 2 }, { 2, 2 }, { 6, 1 }, { 7, 1 }, { 8, 1 }, { 9, 1 }, { 0, 4 } }, | |
97 { { 0, 0 }, { 0, 1 }, { 1, 1 }, { 0, 2 }, { 2, 1 }, { 0, 3 }, { 0, 4 }, { 0, 5 }, | |
98 { 3, 1 }, { 4, 1 }, { 1, 2 }, { 1, 3 }, { 0, 6 }, { 0, 7 }, { 0, 8 }, { 0, 9 } } | |
99 }; | |
100 | |
101 static const uint32_t svq3_dequant_coeff[32] = { | |
102 3881, 4351, 4890, 5481, 6154, 6914, 7761, 8718, | |
103 9781, 10987, 12339, 13828, 15523, 17435, 19561, 21873, | |
104 24552, 27656, 30847, 34870, 38807, 43747, 49103, 54683, | |
105 61694, 68745, 77615, 89113,100253,109366,126635,141533 | |
106 }; | |
107 | |
108 | |
109 static void svq3_luma_dc_dequant_idct_c(DCTELEM *block, int qp){ | |
110 const int qmul= svq3_dequant_coeff[qp]; | |
111 #define stride 16 | |
112 int i; | |
113 int temp[16]; | |
114 static const int x_offset[4]={0, 1*stride, 4* stride, 5*stride}; | |
115 static const int y_offset[4]={0, 2*stride, 8* stride, 10*stride}; | |
116 | |
117 for(i=0; i<4; i++){ | |
118 const int offset= y_offset[i]; | |
119 const int z0= 13*(block[offset+stride*0] + block[offset+stride*4]); | |
120 const int z1= 13*(block[offset+stride*0] - block[offset+stride*4]); | |
121 const int z2= 7* block[offset+stride*1] - 17*block[offset+stride*5]; | |
122 const int z3= 17* block[offset+stride*1] + 7*block[offset+stride*5]; | |
123 | |
124 temp[4*i+0]= z0+z3; | |
125 temp[4*i+1]= z1+z2; | |
126 temp[4*i+2]= z1-z2; | |
127 temp[4*i+3]= z0-z3; | |
128 } | |
129 | |
130 for(i=0; i<4; i++){ | |
131 const int offset= x_offset[i]; | |
132 const int z0= 13*(temp[4*0+i] + temp[4*2+i]); | |
133 const int z1= 13*(temp[4*0+i] - temp[4*2+i]); | |
134 const int z2= 7* temp[4*1+i] - 17*temp[4*3+i]; | |
135 const int z3= 17* temp[4*1+i] + 7*temp[4*3+i]; | |
136 | |
137 block[stride*0 +offset]= ((z0 + z3)*qmul + 0x80000)>>20; | |
138 block[stride*2 +offset]= ((z1 + z2)*qmul + 0x80000)>>20; | |
139 block[stride*8 +offset]= ((z1 - z2)*qmul + 0x80000)>>20; | |
140 block[stride*10+offset]= ((z0 - z3)*qmul + 0x80000)>>20; | |
141 } | |
142 } | |
143 #undef stride | |
144 | |
145 static void svq3_add_idct_c (uint8_t *dst, DCTELEM *block, int stride, int qp, int dc){ | |
146 const int qmul= svq3_dequant_coeff[qp]; | |
147 int i; | |
148 uint8_t *cm = cropTbl + MAX_NEG_CROP; | |
149 | |
150 if (dc) { | |
151 dc = 13*13*((dc == 1) ? 1538*block[0] : ((qmul*(block[0] >> 3)) / 2)); | |
152 block[0] = 0; | |
153 } | |
154 | |
155 for (i=0; i < 4; i++) { | |
156 const int z0= 13*(block[0 + 4*i] + block[2 + 4*i]); | |
157 const int z1= 13*(block[0 + 4*i] - block[2 + 4*i]); | |
158 const int z2= 7* block[1 + 4*i] - 17*block[3 + 4*i]; | |
159 const int z3= 17* block[1 + 4*i] + 7*block[3 + 4*i]; | |
160 | |
161 block[0 + 4*i]= z0 + z3; | |
162 block[1 + 4*i]= z1 + z2; | |
163 block[2 + 4*i]= z1 - z2; | |
164 block[3 + 4*i]= z0 - z3; | |
165 } | |
166 | |
167 for (i=0; i < 4; i++) { | |
168 const int z0= 13*(block[i + 4*0] + block[i + 4*2]); | |
169 const int z1= 13*(block[i + 4*0] - block[i + 4*2]); | |
170 const int z2= 7* block[i + 4*1] - 17*block[i + 4*3]; | |
171 const int z3= 17* block[i + 4*1] + 7*block[i + 4*3]; | |
172 const int rr= (dc + 0x80000); | |
173 | |
174 dst[i + stride*0]= cm[ dst[i + stride*0] + (((z0 + z3)*qmul + rr) >> 20) ]; | |
175 dst[i + stride*1]= cm[ dst[i + stride*1] + (((z1 + z2)*qmul + rr) >> 20) ]; | |
176 dst[i + stride*2]= cm[ dst[i + stride*2] + (((z1 - z2)*qmul + rr) >> 20) ]; | |
177 dst[i + stride*3]= cm[ dst[i + stride*3] + (((z0 - z3)*qmul + rr) >> 20) ]; | |
178 } | |
179 } | |
180 | |
181 static void pred4x4_down_left_svq3_c(uint8_t *src, uint8_t *topright, int stride){ | |
2967 | 182 LOAD_TOP_EDGE |
183 LOAD_LEFT_EDGE | |
1234 | 184 const __attribute__((unused)) int unu0= t0; |
185 const __attribute__((unused)) int unu1= l0; | |
186 | |
187 src[0+0*stride]=(l1 + t1)>>1; | |
188 src[1+0*stride]= | |
189 src[0+1*stride]=(l2 + t2)>>1; | |
190 src[2+0*stride]= | |
191 src[1+1*stride]= | |
192 src[0+2*stride]= | |
193 src[3+0*stride]= | |
194 src[2+1*stride]= | |
195 src[1+2*stride]= | |
196 src[0+3*stride]= | |
197 src[3+1*stride]= | |
198 src[2+2*stride]= | |
199 src[1+3*stride]= | |
200 src[3+2*stride]= | |
201 src[2+3*stride]= | |
202 src[3+3*stride]=(l3 + t3)>>1; | |
1282 | 203 } |
1234 | 204 |
205 static void pred16x16_plane_svq3_c(uint8_t *src, int stride){ | |
206 pred16x16_plane_compat_c(src, stride, 1); | |
207 } | |
208 | |
209 static inline int svq3_decode_block (GetBitContext *gb, DCTELEM *block, | |
2979 | 210 int index, const int type) { |
1234 | 211 |
212 static const uint8_t *const scan_patterns[4] = | |
213 { luma_dc_zigzag_scan, zigzag_scan, svq3_scan, chroma_dc_scan }; | |
214 | |
215 int run, level, sign, vlc, limit; | |
216 const int intra = (3 * type) >> 2; | |
217 const uint8_t *const scan = scan_patterns[type]; | |
218 | |
219 for (limit=(16 >> intra); index < 16; index=limit, limit+=8) { | |
220 for (; (vlc = svq3_get_ue_golomb (gb)) != 0; index++) { | |
221 | |
222 if (vlc == INVALID_VLC) | |
2979 | 223 return -1; |
1234 | 224 |
225 sign = (vlc & 0x1) - 1; | |
226 vlc = (vlc + 1) >> 1; | |
227 | |
228 if (type == 3) { | |
2979 | 229 if (vlc < 3) { |
230 run = 0; | |
231 level = vlc; | |
232 } else if (vlc < 4) { | |
233 run = 1; | |
234 level = 1; | |
235 } else { | |
236 run = (vlc & 0x3); | |
237 level = ((vlc + 9) >> 2) - run; | |
238 } | |
1234 | 239 } else { |
2979 | 240 if (vlc < 16) { |
241 run = svq3_dct_tables[intra][vlc].run; | |
242 level = svq3_dct_tables[intra][vlc].level; | |
243 } else if (intra) { | |
244 run = (vlc & 0x7); | |
245 level = (vlc >> 3) + ((run == 0) ? 8 : ((run < 2) ? 2 : ((run < 5) ? 0 : -1))); | |
246 } else { | |
247 run = (vlc & 0xF); | |
248 level = (vlc >> 4) + ((run == 0) ? 4 : ((run < 3) ? 2 : ((run < 10) ? 1 : 0))); | |
249 } | |
1234 | 250 } |
251 | |
252 if ((index += run) >= limit) | |
2979 | 253 return -1; |
1234 | 254 |
255 block[scan[index]] = (level ^ sign) - sign; | |
256 } | |
257 | |
258 if (type != 2) { | |
259 break; | |
260 } | |
261 } | |
262 | |
263 return 0; | |
264 } | |
265 | |
1319 | 266 static inline void svq3_mc_dir_part (MpegEncContext *s, |
2979 | 267 int x, int y, int width, int height, |
268 int mx, int my, int dxy, | |
269 int thirdpel, int dir, int avg) { | |
1319 | 270 |
271 const Picture *pic = (dir == 0) ? &s->last_picture : &s->next_picture; | |
1234 | 272 uint8_t *src, *dest; |
273 int i, emu = 0; | |
1267
85b71f9f7450
moving the svq3 motion compensation stuff to dsputil (this also means that existing optimized halfpel code is used now ...)
michaelni
parents:
1265
diff
changeset
|
274 int blocksize= 2 - (width>>3); //16->0, 8->1, 4->2 |
1234 | 275 |
1265 | 276 mx += x; |
277 my += y; | |
2967 | 278 |
1268 | 279 if (mx < 0 || mx >= (s->h_edge_pos - width - 1) || |
280 my < 0 || my >= (s->v_edge_pos - height - 1)) { | |
1234 | 281 |
282 if ((s->flags & CODEC_FLAG_EMU_EDGE)) { | |
283 emu = 1; | |
284 } | |
285 | |
1268 | 286 mx = clip (mx, -16, (s->h_edge_pos - width + 15)); |
287 my = clip (my, -16, (s->v_edge_pos - height + 15)); | |
1234 | 288 } |
289 | |
290 /* form component predictions */ | |
291 dest = s->current_picture.data[0] + x + y*s->linesize; | |
1319 | 292 src = pic->data[0] + mx + my*s->linesize; |
1234 | 293 |
294 if (emu) { | |
1321 | 295 ff_emulated_edge_mc (s->edge_emu_buffer, src, s->linesize, (width + 1), (height + 1), |
2979 | 296 mx, my, s->h_edge_pos, s->v_edge_pos); |
1234 | 297 src = s->edge_emu_buffer; |
298 } | |
1267
85b71f9f7450
moving the svq3 motion compensation stuff to dsputil (this also means that existing optimized halfpel code is used now ...)
michaelni
parents:
1265
diff
changeset
|
299 if(thirdpel) |
1319 | 300 (avg ? s->dsp.avg_tpel_pixels_tab : s->dsp.put_tpel_pixels_tab)[dxy](dest, src, s->linesize, width, height); |
1267
85b71f9f7450
moving the svq3 motion compensation stuff to dsputil (this also means that existing optimized halfpel code is used now ...)
michaelni
parents:
1265
diff
changeset
|
301 else |
1319 | 302 (avg ? s->dsp.avg_pixels_tab : s->dsp.put_pixels_tab)[blocksize][dxy](dest, src, s->linesize, height); |
1234 | 303 |
304 if (!(s->flags & CODEC_FLAG_GRAY)) { | |
2979 | 305 mx = (mx + (mx < (int) x)) >> 1; |
306 my = (my + (my < (int) y)) >> 1; | |
1234 | 307 width = (width >> 1); |
308 height = (height >> 1); | |
1267
85b71f9f7450
moving the svq3 motion compensation stuff to dsputil (this also means that existing optimized halfpel code is used now ...)
michaelni
parents:
1265
diff
changeset
|
309 blocksize++; |
1234 | 310 |
311 for (i=1; i < 3; i++) { | |
312 dest = s->current_picture.data[i] + (x >> 1) + (y >> 1)*s->uvlinesize; | |
1319 | 313 src = pic->data[i] + mx + my*s->uvlinesize; |
1234 | 314 |
315 if (emu) { | |
1321 | 316 ff_emulated_edge_mc (s->edge_emu_buffer, src, s->uvlinesize, (width + 1), (height + 1), |
2979 | 317 mx, my, (s->h_edge_pos >> 1), (s->v_edge_pos >> 1)); |
1321 | 318 src = s->edge_emu_buffer; |
1234 | 319 } |
1267
85b71f9f7450
moving the svq3 motion compensation stuff to dsputil (this also means that existing optimized halfpel code is used now ...)
michaelni
parents:
1265
diff
changeset
|
320 if(thirdpel) |
1319 | 321 (avg ? s->dsp.avg_tpel_pixels_tab : s->dsp.put_tpel_pixels_tab)[dxy](dest, src, s->uvlinesize, width, height); |
1267
85b71f9f7450
moving the svq3 motion compensation stuff to dsputil (this also means that existing optimized halfpel code is used now ...)
michaelni
parents:
1265
diff
changeset
|
322 else |
1319 | 323 (avg ? s->dsp.avg_pixels_tab : s->dsp.put_pixels_tab)[blocksize][dxy](dest, src, s->uvlinesize, height); |
1234 | 324 } |
325 } | |
326 } | |
327 | |
1330
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
328 static inline int svq3_mc_dir (H264Context *h, int size, int mode, int dir, int avg) { |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
329 |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
330 int i, j, k, mx, my, dx, dy, x, y; |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
331 MpegEncContext *const s = (MpegEncContext *) h; |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
332 const int part_width = ((size & 5) == 4) ? 4 : 16 >> (size & 1); |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
333 const int part_height = 16 >> ((unsigned) (size + 1) / 3); |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
334 const int extra_width = (mode == PREDICT_MODE) ? -16*6 : 0; |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
335 const int h_edge_pos = 6*(s->h_edge_pos - part_width ) - extra_width; |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
336 const int v_edge_pos = 6*(s->v_edge_pos - part_height) - extra_width; |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
337 |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
338 for (i=0; i < 16; i+=part_height) { |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
339 for (j=0; j < 16; j+=part_width) { |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
340 const int b_xy = (4*s->mb_x+(j>>2)) + (4*s->mb_y+(i>>2))*h->b_stride; |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
341 int dxy; |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
342 x = 16*s->mb_x + j; |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
343 y = 16*s->mb_y + i; |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
344 k = ((j>>2)&1) + ((i>>1)&2) + ((j>>1)&4) + (i&8); |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
345 |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
346 if (mode != PREDICT_MODE) { |
2979 | 347 pred_motion (h, k, (part_width >> 2), dir, 1, &mx, &my); |
1330
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
348 } else { |
2979 | 349 mx = s->next_picture.motion_val[0][b_xy][0]<<1; |
350 my = s->next_picture.motion_val[0][b_xy][1]<<1; | |
1330
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
351 |
2979 | 352 if (dir == 0) { |
353 mx = ((mx * h->frame_num_offset) / h->prev_frame_num_offset + 1)>>1; | |
354 my = ((my * h->frame_num_offset) / h->prev_frame_num_offset + 1)>>1; | |
355 } else { | |
356 mx = ((mx * (h->frame_num_offset - h->prev_frame_num_offset)) / h->prev_frame_num_offset + 1)>>1; | |
357 my = ((my * (h->frame_num_offset - h->prev_frame_num_offset)) / h->prev_frame_num_offset + 1)>>1; | |
358 } | |
1330
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
359 } |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
360 |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
361 /* clip motion vector prediction to frame border */ |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
362 mx = clip (mx, extra_width - 6*x, h_edge_pos - 6*x); |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
363 my = clip (my, extra_width - 6*y, v_edge_pos - 6*y); |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
364 |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
365 /* get (optional) motion vector differential */ |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
366 if (mode == PREDICT_MODE) { |
2979 | 367 dx = dy = 0; |
1330
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
368 } else { |
2979 | 369 dy = svq3_get_se_golomb (&s->gb); |
370 dx = svq3_get_se_golomb (&s->gb); | |
1330
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
371 |
2979 | 372 if (dx == INVALID_VLC || dy == INVALID_VLC) { |
2452 | 373 av_log(h->s.avctx, AV_LOG_ERROR, "invalid MV vlc\n"); |
2979 | 374 return -1; |
375 } | |
1330
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
376 } |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
377 |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
378 /* compute motion vector */ |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
379 if (mode == THIRDPEL_MODE) { |
2979 | 380 int fx, fy; |
381 mx = ((mx + 1)>>1) + dx; | |
382 my = ((my + 1)>>1) + dy; | |
383 fx= ((unsigned)(mx + 0x3000))/3 - 0x1000; | |
384 fy= ((unsigned)(my + 0x3000))/3 - 0x1000; | |
385 dxy= (mx - 3*fx) + 4*(my - 3*fy); | |
1330
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
386 |
2979 | 387 svq3_mc_dir_part (s, x, y, part_width, part_height, fx, fy, dxy, 1, dir, avg); |
388 mx += mx; | |
389 my += my; | |
1330
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
390 } else if (mode == HALFPEL_MODE || mode == PREDICT_MODE) { |
2979 | 391 mx = ((unsigned)(mx + 1 + 0x3000))/3 + dx - 0x1000; |
392 my = ((unsigned)(my + 1 + 0x3000))/3 + dy - 0x1000; | |
393 dxy= (mx&1) + 2*(my&1); | |
1330
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
394 |
2979 | 395 svq3_mc_dir_part (s, x, y, part_width, part_height, mx>>1, my>>1, dxy, 0, dir, avg); |
396 mx *= 3; | |
397 my *= 3; | |
1330
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
398 } else { |
2979 | 399 mx = ((unsigned)(mx + 3 + 0x6000))/6 + dx - 0x1000; |
400 my = ((unsigned)(my + 3 + 0x6000))/6 + dy - 0x1000; | |
1330
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
401 |
2979 | 402 svq3_mc_dir_part (s, x, y, part_width, part_height, mx, my, 0, 0, dir, avg); |
403 mx *= 6; | |
404 my *= 6; | |
1330
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
405 } |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
406 |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
407 /* update mv_cache */ |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
408 if (mode != PREDICT_MODE) { |
2979 | 409 int32_t mv = pack16to32(mx,my); |
1330
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
410 |
2979 | 411 if (part_height == 8 && i < 8) { |
412 *(int32_t *) h->mv_cache[dir][scan8[k] + 1*8] = mv; | |
1330
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
413 |
2979 | 414 if (part_width == 8 && j < 8) { |
415 *(int32_t *) h->mv_cache[dir][scan8[k] + 1 + 1*8] = mv; | |
416 } | |
417 } | |
418 if (part_width == 8 && j < 8) { | |
419 *(int32_t *) h->mv_cache[dir][scan8[k] + 1] = mv; | |
420 } | |
421 if (part_width == 4 || part_height == 4) { | |
422 *(int32_t *) h->mv_cache[dir][scan8[k]] = mv; | |
423 } | |
1330
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
424 } |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
425 |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
426 /* write back motion vectors */ |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
427 fill_rectangle(s->current_picture.motion_val[dir][b_xy], part_width>>2, part_height>>2, h->b_stride, pack16to32(mx,my), 4); |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
428 } |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
429 } |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
430 |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
431 return 0; |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
432 } |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
433 |
1234 | 434 static int svq3_decode_mb (H264Context *h, unsigned int mb_type) { |
1330
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
435 int i, j, k, m, dir, mode; |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
436 int cbp = 0; |
1234 | 437 uint32_t vlc; |
438 int8_t *top, *left; | |
439 MpegEncContext *const s = (MpegEncContext *) h; | |
440 const int mb_xy = s->mb_x + s->mb_y*s->mb_stride; | |
441 const int b_xy = 4*s->mb_x + 4*s->mb_y*h->b_stride; | |
442 | |
2979 | 443 h->top_samples_available = (s->mb_y == 0) ? 0x33FF : 0xFFFF; |
444 h->left_samples_available = (s->mb_x == 0) ? 0x5F5F : 0xFFFF; | |
445 h->topright_samples_available = 0xFFFF; | |
1234 | 446 |
2979 | 447 if (mb_type == 0) { /* SKIP */ |
1330
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
448 if (s->pict_type == P_TYPE || s->next_picture.mb_type[mb_xy] == -1) { |
1319 | 449 svq3_mc_dir_part (s, 16*s->mb_x, 16*s->mb_y, 16, 16, 0, 0, 0, 0, 0, 0); |
450 | |
1330
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
451 if (s->pict_type == B_TYPE) { |
2979 | 452 svq3_mc_dir_part (s, 16*s->mb_x, 16*s->mb_y, 16, 16, 0, 0, 0, 0, 1, 1); |
1330
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
453 } |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
454 |
1319 | 455 mb_type = MB_TYPE_SKIP; |
456 } else { | |
2088
19d4d5c6a8f4
1000l to myself for FFMIN(...,0) instead of ...,6)
michael
parents:
2086
diff
changeset
|
457 mb_type= FFMIN(s->next_picture.mb_type[mb_xy], 6); |
19d4d5c6a8f4
1000l to myself for FFMIN(...,0) instead of ...,6)
michael
parents:
2086
diff
changeset
|
458 if(svq3_mc_dir (h, mb_type, PREDICT_MODE, 0, 0) < 0) |
19d4d5c6a8f4
1000l to myself for FFMIN(...,0) instead of ...,6)
michael
parents:
2086
diff
changeset
|
459 return -1; |
19d4d5c6a8f4
1000l to myself for FFMIN(...,0) instead of ...,6)
michael
parents:
2086
diff
changeset
|
460 if(svq3_mc_dir (h, mb_type, PREDICT_MODE, 1, 1) < 0) |
19d4d5c6a8f4
1000l to myself for FFMIN(...,0) instead of ...,6)
michael
parents:
2086
diff
changeset
|
461 return -1; |
1319 | 462 |
463 mb_type = MB_TYPE_16x16; | |
464 } | |
2979 | 465 } else if (mb_type < 8) { /* INTER */ |
1234 | 466 if (h->thirdpel_flag && h->halfpel_flag == !get_bits (&s->gb, 1)) { |
1265 | 467 mode = THIRDPEL_MODE; |
1234 | 468 } else if (h->halfpel_flag && h->thirdpel_flag == !get_bits (&s->gb, 1)) { |
1265 | 469 mode = HALFPEL_MODE; |
1234 | 470 } else { |
1265 | 471 mode = FULLPEL_MODE; |
1234 | 472 } |
473 | |
474 /* fill caches */ | |
1319 | 475 /* note ref_cache should contain here: |
1252 | 476 ???????? |
477 ???11111 | |
478 N??11111 | |
479 N??11111 | |
480 N??11111 | |
481 N | |
482 */ | |
2967 | 483 |
1319 | 484 for (m=0; m < 2; m++) { |
485 if (s->mb_x > 0 && h->intra4x4_pred_mode[mb_xy - 1][0] != -1) { | |
2979 | 486 for (i=0; i < 4; i++) { |
487 *(uint32_t *) h->mv_cache[m][scan8[0] - 1 + i*8] = *(uint32_t *) s->current_picture.motion_val[m][b_xy - 1 + i*h->b_stride]; | |
488 } | |
1319 | 489 } else { |
2979 | 490 for (i=0; i < 4; i++) { |
491 *(uint32_t *) h->mv_cache[m][scan8[0] - 1 + i*8] = 0; | |
492 } | |
1234 | 493 } |
1319 | 494 if (s->mb_y > 0) { |
2979 | 495 memcpy (h->mv_cache[m][scan8[0] - 1*8], s->current_picture.motion_val[m][b_xy - h->b_stride], 4*2*sizeof(int16_t)); |
496 memset (&h->ref_cache[m][scan8[0] - 1*8], (h->intra4x4_pred_mode[mb_xy - s->mb_stride][4] == -1) ? PART_NOT_AVAILABLE : 1, 4); | |
1234 | 497 |
2979 | 498 if (s->mb_x < (s->mb_width - 1)) { |
499 *(uint32_t *) h->mv_cache[m][scan8[0] + 4 - 1*8] = *(uint32_t *) s->current_picture.motion_val[m][b_xy - h->b_stride + 4]; | |
500 h->ref_cache[m][scan8[0] + 4 - 1*8] = | |
501 (h->intra4x4_pred_mode[mb_xy - s->mb_stride + 1][0] == -1 || | |
502 h->intra4x4_pred_mode[mb_xy - s->mb_stride][4] == -1) ? PART_NOT_AVAILABLE : 1; | |
503 }else | |
504 h->ref_cache[m][scan8[0] + 4 - 1*8] = PART_NOT_AVAILABLE; | |
505 if (s->mb_x > 0) { | |
506 *(uint32_t *) h->mv_cache[m][scan8[0] - 1 - 1*8] = *(uint32_t *) s->current_picture.motion_val[m][b_xy - h->b_stride - 1]; | |
507 h->ref_cache[m][scan8[0] - 1 - 1*8] = (h->intra4x4_pred_mode[mb_xy - s->mb_stride - 1][3] == -1) ? PART_NOT_AVAILABLE : 1; | |
508 }else | |
509 h->ref_cache[m][scan8[0] - 1 - 1*8] = PART_NOT_AVAILABLE; | |
1252 | 510 }else |
2979 | 511 memset (&h->ref_cache[m][scan8[0] - 1*8 - 1], PART_NOT_AVAILABLE, 8); |
1319 | 512 |
513 if (s->pict_type != B_TYPE) | |
2979 | 514 break; |
1319 | 515 } |
1234 | 516 |
517 /* decode motion vector(s) and form prediction(s) */ | |
1319 | 518 if (s->pict_type == P_TYPE) { |
2088
19d4d5c6a8f4
1000l to myself for FFMIN(...,0) instead of ...,6)
michael
parents:
2086
diff
changeset
|
519 if(svq3_mc_dir (h, (mb_type - 1), mode, 0, 0) < 0) |
19d4d5c6a8f4
1000l to myself for FFMIN(...,0) instead of ...,6)
michael
parents:
2086
diff
changeset
|
520 return -1; |
2979 | 521 } else { /* B_TYPE */ |
1330
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
522 if (mb_type != 2) { |
2979 | 523 if(svq3_mc_dir (h, 0, mode, 0, 0) < 0) |
2088
19d4d5c6a8f4
1000l to myself for FFMIN(...,0) instead of ...,6)
michael
parents:
2086
diff
changeset
|
524 return -1; |
1330
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
525 } else { |
2979 | 526 for (i=0; i < 4; i++) { |
527 memset (s->current_picture.motion_val[0][b_xy + i*h->b_stride], 0, 4*2*sizeof(int16_t)); | |
528 } | |
1234 | 529 } |
1319 | 530 if (mb_type != 1) { |
2979 | 531 if(svq3_mc_dir (h, 0, mode, 1, (mb_type == 3)) < 0) |
2088
19d4d5c6a8f4
1000l to myself for FFMIN(...,0) instead of ...,6)
michael
parents:
2086
diff
changeset
|
532 return -1; |
1319 | 533 } else { |
2979 | 534 for (i=0; i < 4; i++) { |
535 memset (s->current_picture.motion_val[1][b_xy + i*h->b_stride], 0, 4*2*sizeof(int16_t)); | |
536 } | |
1319 | 537 } |
1234 | 538 } |
539 | |
540 mb_type = MB_TYPE_16x16; | |
2979 | 541 } else if (mb_type == 8 || mb_type == 33) { /* INTRA4x4 */ |
1234 | 542 memset (h->intra4x4_pred_mode_cache, -1, 8*5*sizeof(int8_t)); |
543 | |
1319 | 544 if (mb_type == 8) { |
545 if (s->mb_x > 0) { | |
2979 | 546 for (i=0; i < 4; i++) { |
547 h->intra4x4_pred_mode_cache[scan8[0] - 1 + i*8] = h->intra4x4_pred_mode[mb_xy - 1][i]; | |
548 } | |
549 if (h->intra4x4_pred_mode_cache[scan8[0] - 1] == -1) { | |
550 h->left_samples_available = 0x5F5F; | |
551 } | |
1234 | 552 } |
1319 | 553 if (s->mb_y > 0) { |
2979 | 554 h->intra4x4_pred_mode_cache[4+8*0] = h->intra4x4_pred_mode[mb_xy - s->mb_stride][4]; |
555 h->intra4x4_pred_mode_cache[5+8*0] = h->intra4x4_pred_mode[mb_xy - s->mb_stride][5]; | |
556 h->intra4x4_pred_mode_cache[6+8*0] = h->intra4x4_pred_mode[mb_xy - s->mb_stride][6]; | |
557 h->intra4x4_pred_mode_cache[7+8*0] = h->intra4x4_pred_mode[mb_xy - s->mb_stride][3]; | |
1319 | 558 |
2979 | 559 if (h->intra4x4_pred_mode_cache[4+8*0] == -1) { |
560 h->top_samples_available = 0x33FF; | |
561 } | |
1319 | 562 } |
1234 | 563 |
1319 | 564 /* decode prediction codes for luma blocks */ |
565 for (i=0; i < 16; i+=2) { | |
2979 | 566 vlc = svq3_get_ue_golomb (&s->gb); |
1234 | 567 |
2979 | 568 if (vlc >= 25){ |
2452 | 569 av_log(h->s.avctx, AV_LOG_ERROR, "luma prediction:%d\n", vlc); |
2979 | 570 return -1; |
2452 | 571 } |
1319 | 572 |
2979 | 573 left = &h->intra4x4_pred_mode_cache[scan8[i] - 1]; |
574 top = &h->intra4x4_pred_mode_cache[scan8[i] - 8]; | |
1234 | 575 |
2979 | 576 left[1] = svq3_pred_1[top[0] + 1][left[0] + 1][svq3_pred_0[vlc][0]]; |
577 left[2] = svq3_pred_1[top[1] + 1][left[1] + 1][svq3_pred_0[vlc][1]]; | |
1234 | 578 |
2979 | 579 if (left[1] == -1 || left[2] == -1){ |
2452 | 580 av_log(h->s.avctx, AV_LOG_ERROR, "weird prediction\n"); |
2979 | 581 return -1; |
2452 | 582 } |
1319 | 583 } |
2979 | 584 } else { /* mb_type == 33, DC_128_PRED block type */ |
1319 | 585 for (i=0; i < 4; i++) { |
2979 | 586 memset (&h->intra4x4_pred_mode_cache[scan8[0] + 8*i], DC_PRED, 4); |
1319 | 587 } |
1234 | 588 } |
589 | |
590 write_back_intra_pred_mode (h); | |
1319 | 591 |
592 if (mb_type == 8) { | |
593 check_intra4x4_pred_mode (h); | |
594 | |
595 h->top_samples_available = (s->mb_y == 0) ? 0x33FF : 0xFFFF; | |
596 h->left_samples_available = (s->mb_x == 0) ? 0x5F5F : 0xFFFF; | |
597 } else { | |
598 for (i=0; i < 4; i++) { | |
599 memset (&h->intra4x4_pred_mode_cache[scan8[0] + 8*i], DC_128_PRED, 4); | |
600 } | |
601 | |
602 h->top_samples_available = 0x33FF; | |
603 h->left_samples_available = 0x5F5F; | |
604 } | |
1234 | 605 |
606 mb_type = MB_TYPE_INTRA4x4; | |
2979 | 607 } else { /* INTRA16x16 */ |
1234 | 608 dir = i_mb_type_info[mb_type - 8].pred_mode; |
609 dir = (dir >> 1) ^ 3*(dir & 1) ^ 1; | |
610 | |
2452 | 611 if ((h->intra16x16_pred_mode = check_intra_pred_mode (h, dir)) == -1){ |
612 av_log(h->s.avctx, AV_LOG_ERROR, "check_intra_pred_mode = -1\n"); | |
1234 | 613 return -1; |
2452 | 614 } |
1234 | 615 |
616 cbp = i_mb_type_info[mb_type - 8].cbp; | |
617 mb_type = MB_TYPE_INTRA16x16; | |
618 } | |
619 | |
620 if (!IS_INTER(mb_type) && s->pict_type != I_TYPE) { | |
621 for (i=0; i < 4; i++) { | |
622 memset (s->current_picture.motion_val[0][b_xy + i*h->b_stride], 0, 4*2*sizeof(int16_t)); | |
623 } | |
1319 | 624 if (s->pict_type == B_TYPE) { |
625 for (i=0; i < 4; i++) { | |
2979 | 626 memset (s->current_picture.motion_val[1][b_xy + i*h->b_stride], 0, 4*2*sizeof(int16_t)); |
1319 | 627 } |
628 } | |
1234 | 629 } |
630 if (!IS_INTRA4x4(mb_type)) { | |
631 memset (h->intra4x4_pred_mode[mb_xy], DC_PRED, 8); | |
632 } | |
1330
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
633 if (!IS_SKIP(mb_type) || s->pict_type == B_TYPE) { |
1252 | 634 memset (h->non_zero_count_cache + 8, 0, 4*9*sizeof(uint8_t)); |
635 s->dsp.clear_blocks(h->mb); | |
1234 | 636 } |
637 | |
1330
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
638 if (!IS_INTRA16x16(mb_type) && (!IS_SKIP(mb_type) || s->pict_type == B_TYPE)) { |
2452 | 639 if ((vlc = svq3_get_ue_golomb (&s->gb)) >= 48){ |
640 av_log(h->s.avctx, AV_LOG_ERROR, "cbp_vlc=%d\n", vlc); | |
1330
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
641 return -1; |
2452 | 642 } |
1330
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
643 |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
644 cbp = IS_INTRA(mb_type) ? golomb_to_intra4x4_cbp[vlc] : golomb_to_inter_cbp[vlc]; |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
645 } |
1234 | 646 if (IS_INTRA16x16(mb_type) || (s->pict_type != I_TYPE && s->adaptive_quant && cbp)) { |
647 s->qscale += svq3_get_se_golomb (&s->gb); | |
648 | |
2452 | 649 if (s->qscale > 31){ |
650 av_log(h->s.avctx, AV_LOG_ERROR, "qscale:%d\n", s->qscale); | |
1234 | 651 return -1; |
2452 | 652 } |
1234 | 653 } |
654 if (IS_INTRA16x16(mb_type)) { | |
2452 | 655 if (svq3_decode_block (&s->gb, h->mb, 0, 0)){ |
656 av_log(h->s.avctx, AV_LOG_ERROR, "error while decoding intra luma dc\n"); | |
1234 | 657 return -1; |
2452 | 658 } |
1234 | 659 } |
660 | |
1330
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
661 if (cbp) { |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
662 const int index = IS_INTRA16x16(mb_type) ? 1 : 0; |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
663 const int type = ((s->qscale < 24 && IS_INTRA4x4(mb_type)) ? 2 : 1); |
1234 | 664 |
665 for (i=0; i < 4; i++) { | |
666 if ((cbp & (1 << i))) { | |
2979 | 667 for (j=0; j < 4; j++) { |
668 k = index ? ((j&1) + 2*(i&1) + 2*(j&2) + 4*(i&2)) : (4*i + j); | |
669 h->non_zero_count_cache[ scan8[k] ] = 1; | |
1234 | 670 |
2979 | 671 if (svq3_decode_block (&s->gb, &h->mb[16*k], index, type)){ |
2452 | 672 av_log(h->s.avctx, AV_LOG_ERROR, "error while decoding block\n"); |
2979 | 673 return -1; |
2452 | 674 } |
2979 | 675 } |
1234 | 676 } |
677 } | |
678 | |
679 if ((cbp & 0x30)) { | |
680 for (i=0; i < 2; ++i) { | |
2979 | 681 if (svq3_decode_block (&s->gb, &h->mb[16*(16 + 4*i)], 0, 3)){ |
2452 | 682 av_log(h->s.avctx, AV_LOG_ERROR, "error while decoding chroma dc block\n"); |
2979 | 683 return -1; |
2452 | 684 } |
1234 | 685 } |
686 | |
687 if ((cbp & 0x20)) { | |
2979 | 688 for (i=0; i < 8; i++) { |
689 h->non_zero_count_cache[ scan8[16+i] ] = 1; | |
1234 | 690 |
2979 | 691 if (svq3_decode_block (&s->gb, &h->mb[16*(16 + i)], 1, 1)){ |
2452 | 692 av_log(h->s.avctx, AV_LOG_ERROR, "error while decoding chroma ac block\n"); |
2979 | 693 return -1; |
2452 | 694 } |
2979 | 695 } |
1234 | 696 } |
697 } | |
698 } | |
699 | |
700 s->current_picture.mb_type[mb_xy] = mb_type; | |
701 | |
702 if (IS_INTRA(mb_type)) { | |
703 h->chroma_pred_mode = check_intra_pred_mode (h, DC_PRED8x8); | |
704 } | |
705 | |
706 return 0; | |
707 } | |
708 | |
1319 | 709 static int svq3_decode_slice_header (H264Context *h) { |
710 MpegEncContext *const s = (MpegEncContext *) h; | |
711 const int mb_xy = s->mb_x + s->mb_y*s->mb_stride; | |
712 int i, header; | |
713 | |
714 header = get_bits (&s->gb, 8); | |
715 | |
716 if (((header & 0x9F) != 1 && (header & 0x9F) != 2) || (header & 0x60) == 0) { | |
717 /* TODO: what? */ | |
1598
932d306bf1dc
av_log() patch by (Michel Bardiaux <mbardiaux at peaktime dot be>)
michael
parents:
1330
diff
changeset
|
718 av_log(h->s.avctx, AV_LOG_ERROR, "unsupported slice header (%02X)\n", header); |
1319 | 719 return -1; |
720 } else { | |
721 int length = (header >> 5) & 3; | |
722 | |
2614
5e24800ab329
various fixes related to the non alt_bitstream_reader
michael
parents:
2453
diff
changeset
|
723 h->next_slice_index = get_bits_count(&s->gb) + 8*show_bits (&s->gb, 8*length) + 8*length; |
1319 | 724 |
1754
bdf3927bf8c5
closed gop support & flags2 as all bits in flags are used
michael
parents:
1699
diff
changeset
|
725 if (h->next_slice_index > s->gb.size_in_bits){ |
bdf3927bf8c5
closed gop support & flags2 as all bits in flags are used
michael
parents:
1699
diff
changeset
|
726 av_log(h->s.avctx, AV_LOG_ERROR, "slice after bitstream end\n"); |
1319 | 727 return -1; |
1754
bdf3927bf8c5
closed gop support & flags2 as all bits in flags are used
michael
parents:
1699
diff
changeset
|
728 } |
1319 | 729 |
730 s->gb.size_in_bits = h->next_slice_index - 8*(length - 1); | |
2614
5e24800ab329
various fixes related to the non alt_bitstream_reader
michael
parents:
2453
diff
changeset
|
731 skip_bits(&s->gb, 8); |
1319 | 732 |
733 if (length > 0) { | |
2614
5e24800ab329
various fixes related to the non alt_bitstream_reader
michael
parents:
2453
diff
changeset
|
734 memcpy ((uint8_t *) &s->gb.buffer[get_bits_count(&s->gb) >> 3], |
1319 | 735 &s->gb.buffer[s->gb.size_in_bits >> 3], (length - 1)); |
736 } | |
737 } | |
738 | |
1754
bdf3927bf8c5
closed gop support & flags2 as all bits in flags are used
michael
parents:
1699
diff
changeset
|
739 if ((i = svq3_get_ue_golomb (&s->gb)) == INVALID_VLC || i >= 3){ |
bdf3927bf8c5
closed gop support & flags2 as all bits in flags are used
michael
parents:
1699
diff
changeset
|
740 av_log(h->s.avctx, AV_LOG_ERROR, "illegal slice type %d \n", i); |
1319 | 741 return -1; |
1754
bdf3927bf8c5
closed gop support & flags2 as all bits in flags are used
michael
parents:
1699
diff
changeset
|
742 } |
1319 | 743 |
744 h->slice_type = golomb_to_pict_type[i]; | |
745 | |
746 if ((header & 0x9F) == 2) { | |
747 i = (s->mb_num < 64) ? 6 : (1 + av_log2 (s->mb_num - 1)); | |
748 s->mb_skip_run = get_bits (&s->gb, i) - (s->mb_x + (s->mb_y * s->mb_width)); | |
749 } else { | |
750 get_bits1 (&s->gb); | |
751 s->mb_skip_run = 0; | |
752 } | |
753 | |
754 h->slice_num = get_bits (&s->gb, 8); | |
755 s->qscale = get_bits (&s->gb, 5); | |
756 s->adaptive_quant = get_bits1 (&s->gb); | |
757 | |
758 /* unknown fields */ | |
759 get_bits1 (&s->gb); | |
760 | |
761 if (h->unknown_svq3_flag) { | |
762 get_bits1 (&s->gb); | |
763 } | |
764 | |
765 get_bits1 (&s->gb); | |
766 get_bits (&s->gb, 2); | |
767 | |
768 while (get_bits1 (&s->gb)) { | |
769 get_bits (&s->gb, 8); | |
770 } | |
771 | |
772 /* reset intra predictors and invalidate motion vector references */ | |
773 if (s->mb_x > 0) { | |
774 memset (h->intra4x4_pred_mode[mb_xy - 1], -1, 4*sizeof(int8_t)); | |
775 memset (h->intra4x4_pred_mode[mb_xy - s->mb_x], -1, 8*sizeof(int8_t)*s->mb_x); | |
776 } | |
777 if (s->mb_y > 0) { | |
778 memset (h->intra4x4_pred_mode[mb_xy - s->mb_stride], -1, 8*sizeof(int8_t)*(s->mb_width - s->mb_x)); | |
779 | |
780 if (s->mb_x > 0) { | |
781 h->intra4x4_pred_mode[mb_xy - s->mb_stride - 1][3] = -1; | |
782 } | |
783 } | |
784 | |
785 return 0; | |
786 } | |
787 | |
1234 | 788 static int svq3_decode_frame (AVCodecContext *avctx, |
2979 | 789 void *data, int *data_size, |
790 uint8_t *buf, int buf_size) { | |
1234 | 791 MpegEncContext *const s = avctx->priv_data; |
792 H264Context *const h = avctx->priv_data; | |
1319 | 793 int m, mb_type; |
1877
869256817a91
smarten up the SVQ3 extradata decoder without changing the external API
melanson
parents:
1754
diff
changeset
|
794 unsigned char *extradata; |
1878
838c18d1e7fc
be smarter about the variable size of the extra SVQ3 data
melanson
parents:
1877
diff
changeset
|
795 unsigned int size; |
1319 | 796 |
1234 | 797 s->flags = avctx->flags; |
1754
bdf3927bf8c5
closed gop support & flags2 as all bits in flags are used
michael
parents:
1699
diff
changeset
|
798 s->flags2 = avctx->flags2; |
1699 | 799 s->unrestricted_mv = 1; |
800 | |
1234 | 801 if (!s->context_initialized) { |
1268 | 802 s->width = avctx->width; |
803 s->height = avctx->height; | |
1234 | 804 h->pred4x4[DIAG_DOWN_LEFT_PRED] = pred4x4_down_left_svq3_c; |
805 h->pred16x16[PLANE_PRED8x8] = pred16x16_plane_svq3_c; | |
806 h->halfpel_flag = 1; | |
807 h->thirdpel_flag = 1; | |
1319 | 808 h->unknown_svq3_flag = 0; |
1234 | 809 h->chroma_qp = 4; |
810 | |
811 if (MPV_common_init (s) < 0) | |
812 return -1; | |
813 | |
1268 | 814 h->b_stride = 4*s->mb_width; |
815 | |
1234 | 816 alloc_tables (h); |
1319 | 817 |
1877
869256817a91
smarten up the SVQ3 extradata decoder without changing the external API
melanson
parents:
1754
diff
changeset
|
818 /* prowl for the "SEQH" marker in the extradata */ |
869256817a91
smarten up the SVQ3 extradata decoder without changing the external API
melanson
parents:
1754
diff
changeset
|
819 extradata = (unsigned char *)avctx->extradata; |
869256817a91
smarten up the SVQ3 extradata decoder without changing the external API
melanson
parents:
1754
diff
changeset
|
820 for (m = 0; m < avctx->extradata_size; m++) { |
869256817a91
smarten up the SVQ3 extradata decoder without changing the external API
melanson
parents:
1754
diff
changeset
|
821 if (!memcmp (extradata, "SEQH", 4)) |
869256817a91
smarten up the SVQ3 extradata decoder without changing the external API
melanson
parents:
1754
diff
changeset
|
822 break; |
869256817a91
smarten up the SVQ3 extradata decoder without changing the external API
melanson
parents:
1754
diff
changeset
|
823 extradata++; |
869256817a91
smarten up the SVQ3 extradata decoder without changing the external API
melanson
parents:
1754
diff
changeset
|
824 } |
1878
838c18d1e7fc
be smarter about the variable size of the extra SVQ3 data
melanson
parents:
1877
diff
changeset
|
825 |
1877
869256817a91
smarten up the SVQ3 extradata decoder without changing the external API
melanson
parents:
1754
diff
changeset
|
826 /* if a match was found, parse the extra data */ |
869256817a91
smarten up the SVQ3 extradata decoder without changing the external API
melanson
parents:
1754
diff
changeset
|
827 if (!memcmp (extradata, "SEQH", 4)) { |
1319 | 828 |
829 GetBitContext gb; | |
1234 | 830 |
1883
12f260596631
don't fetch the size of SEQH unless SEQH was found
melanson
parents:
1881
diff
changeset
|
831 size = BE_32(&extradata[4]); |
2966
564788471dd4
bitstream related fixes from [PATCH] from DivX, Part 9: bitstream crashes by (Steve Lhomme | slhomme divxcorp com)
michael
parents:
2792
diff
changeset
|
832 init_get_bits (&gb, extradata + 8, size*8); |
1234 | 833 |
1319 | 834 /* 'frame size code' and optional 'width, height' */ |
835 if (get_bits (&gb, 3) == 7) { | |
2979 | 836 get_bits (&gb, 12); |
837 get_bits (&gb, 12); | |
1234 | 838 } |
839 | |
1319 | 840 h->halfpel_flag = get_bits1 (&gb); |
841 h->thirdpel_flag = get_bits1 (&gb); | |
842 | |
843 /* unknown fields */ | |
844 get_bits1 (&gb); | |
845 get_bits1 (&gb); | |
846 get_bits1 (&gb); | |
847 get_bits1 (&gb); | |
848 | |
849 s->low_delay = get_bits1 (&gb); | |
850 | |
851 /* unknown field */ | |
852 get_bits1 (&gb); | |
853 | |
854 while (get_bits1 (&gb)) { | |
2979 | 855 get_bits (&gb, 8); |
1319 | 856 } |
857 | |
858 h->unknown_svq3_flag = get_bits1 (&gb); | |
859 avctx->has_b_frames = !s->low_delay; | |
1234 | 860 } |
861 } | |
862 | |
1319 | 863 /* special case for last picture */ |
864 if (buf_size == 0) { | |
865 if (s->next_picture_ptr && !s->low_delay) { | |
866 *(AVFrame *) data = *(AVFrame *) &s->next_picture; | |
867 *data_size = sizeof(AVFrame); | |
1234 | 868 } |
1319 | 869 return 0; |
1234 | 870 } |
871 | |
1319 | 872 init_get_bits (&s->gb, buf, 8*buf_size); |
1234 | 873 |
1319 | 874 s->mb_x = s->mb_y = 0; |
875 | |
876 if (svq3_decode_slice_header (h)) | |
1234 | 877 return -1; |
878 | |
1319 | 879 s->pict_type = h->slice_type; |
880 s->picture_number = h->slice_num; | |
1234 | 881 |
1250 | 882 if(avctx->debug&FF_DEBUG_PICT_INFO){ |
2967 | 883 av_log(h->s.avctx, AV_LOG_DEBUG, "%c hpel:%d, tpel:%d aqp:%d qp:%d\n", |
1264 | 884 av_get_pict_type_char(s->pict_type), h->halfpel_flag, h->thirdpel_flag, |
1250 | 885 s->adaptive_quant, s->qscale |
886 ); | |
887 } | |
1234 | 888 |
1319 | 889 /* for hurry_up==5 */ |
890 s->current_picture.pict_type = s->pict_type; | |
891 s->current_picture.key_frame = (s->pict_type == I_TYPE); | |
892 | |
893 /* skip b frames if we dont have reference frames */ | |
894 if (s->last_picture_ptr == NULL && s->pict_type == B_TYPE) return 0; | |
895 /* skip b frames if we are in a hurry */ | |
896 if (avctx->hurry_up && s->pict_type == B_TYPE) return 0; | |
897 /* skip everything if we are in a hurry >= 5 */ | |
898 if (avctx->hurry_up >= 5) return 0; | |
2792 | 899 if( (avctx->skip_frame >= AVDISCARD_NONREF && s->pict_type==B_TYPE) |
900 ||(avctx->skip_frame >= AVDISCARD_NONKEY && s->pict_type!=I_TYPE) | |
901 || avctx->skip_frame >= AVDISCARD_ALL) | |
902 return 0; | |
1319 | 903 |
904 if (s->next_p_frame_damaged) { | |
905 if (s->pict_type == B_TYPE) | |
906 return 0; | |
907 else | |
908 s->next_p_frame_damaged = 0; | |
909 } | |
1234 | 910 |
911 frame_start (h); | |
912 | |
1319 | 913 if (s->pict_type == B_TYPE) { |
914 h->frame_num_offset = (h->slice_num - h->prev_frame_num); | |
915 | |
916 if (h->frame_num_offset < 0) { | |
917 h->frame_num_offset += 256; | |
918 } | |
919 if (h->frame_num_offset == 0 || h->frame_num_offset >= h->prev_frame_num_offset) { | |
1598
932d306bf1dc
av_log() patch by (Michel Bardiaux <mbardiaux at peaktime dot be>)
michael
parents:
1330
diff
changeset
|
920 av_log(h->s.avctx, AV_LOG_ERROR, "error in B-frame picture id\n"); |
1319 | 921 return -1; |
922 } | |
923 } else { | |
924 h->prev_frame_num = h->frame_num; | |
925 h->frame_num = h->slice_num; | |
926 h->prev_frame_num_offset = (h->frame_num - h->prev_frame_num); | |
927 | |
928 if (h->prev_frame_num_offset < 0) { | |
929 h->prev_frame_num_offset += 256; | |
930 } | |
931 } | |
932 | |
933 for(m=0; m<2; m++){ | |
934 int i; | |
935 for(i=0; i<4; i++){ | |
936 int j; | |
937 for(j=-1; j<4; j++) | |
2979 | 938 h->ref_cache[m][scan8[0] + 8*i + j]= 1; |
1319 | 939 h->ref_cache[m][scan8[0] + 8*i + j]= PART_NOT_AVAILABLE; |
940 } | |
1252 | 941 } |
2967 | 942 |
1234 | 943 for (s->mb_y=0; s->mb_y < s->mb_height; s->mb_y++) { |
944 for (s->mb_x=0; s->mb_x < s->mb_width; s->mb_x++) { | |
1319 | 945 |
2614
5e24800ab329
various fixes related to the non alt_bitstream_reader
michael
parents:
2453
diff
changeset
|
946 if ( (get_bits_count(&s->gb) + 7) >= s->gb.size_in_bits && |
2979 | 947 ((get_bits_count(&s->gb) & 7) == 0 || show_bits (&s->gb, (-get_bits_count(&s->gb) & 7)) == 0)) { |
1319 | 948 |
2979 | 949 skip_bits(&s->gb, h->next_slice_index - get_bits_count(&s->gb)); |
950 s->gb.size_in_bits = 8*buf_size; | |
1319 | 951 |
2979 | 952 if (svq3_decode_slice_header (h)) |
953 return -1; | |
1319 | 954 |
2979 | 955 /* TODO: support s->mb_skip_run */ |
1319 | 956 } |
957 | |
958 mb_type = svq3_get_ue_golomb (&s->gb); | |
1234 | 959 |
960 if (s->pict_type == I_TYPE) { | |
2979 | 961 mb_type += 8; |
1319 | 962 } else if (s->pict_type == B_TYPE && mb_type >= 4) { |
2979 | 963 mb_type += 4; |
1234 | 964 } |
1319 | 965 if (mb_type > 33 || svq3_decode_mb (h, mb_type)) { |
2979 | 966 av_log(h->s.avctx, AV_LOG_ERROR, "error while decoding MB %d %d\n", s->mb_x, s->mb_y); |
967 return -1; | |
1234 | 968 } |
969 | |
970 if (mb_type != 0) { | |
2979 | 971 hl_decode_mb (h); |
1234 | 972 } |
1330
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
973 |
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
974 if (s->pict_type != B_TYPE && !s->low_delay) { |
2979 | 975 s->current_picture.mb_type[s->mb_x + s->mb_y*s->mb_stride] = |
976 (s->pict_type == P_TYPE && mb_type < 8) ? (mb_type - 1) : -1; | |
1330
c05c381a9c47
- fix PLANE_PRED8x8 prediction (H/V are swapped, this is correct!)
tmmm
parents:
1321
diff
changeset
|
977 } |
1234 | 978 } |
1268 | 979 |
980 ff_draw_horiz_band(s, 16*s->mb_y, 16); | |
1234 | 981 } |
982 | |
1319 | 983 MPV_frame_end(s); |
984 | |
985 if (s->pict_type == B_TYPE || s->low_delay) { | |
986 *(AVFrame *) data = *(AVFrame *) &s->current_picture; | |
987 } else { | |
988 *(AVFrame *) data = *(AVFrame *) &s->last_picture; | |
989 } | |
1234 | 990 |
1319 | 991 avctx->frame_number = s->picture_number - 1; |
992 | |
993 /* dont output the last pic after seeking */ | |
994 if (s->last_picture_ptr || s->low_delay) { | |
995 *data_size = sizeof(AVFrame); | |
996 } | |
997 | |
1234 | 998 return buf_size; |
999 } | |
1000 | |
1001 | |
1002 AVCodec svq3_decoder = { | |
1003 "svq3", | |
1004 CODEC_TYPE_VIDEO, | |
1005 CODEC_ID_SVQ3, | |
1006 sizeof(H264Context), | |
1007 decode_init, | |
1008 NULL, | |
1009 decode_end, | |
1010 svq3_decode_frame, | |
2453 | 1011 CODEC_CAP_DRAW_HORIZ_BAND | CODEC_CAP_DR1 | CODEC_CAP_DELAY, |
1234 | 1012 }; |