Mercurial > libavcodec.hg
annotate msmpeg4.c @ 246:214ac77c0948 libavcodec
* simplified indexing
* ignoring error from msmpeg4_decode_block
author | kabi |
---|---|
date | Mon, 18 Feb 2002 09:34:54 +0000 |
parents | 3c2bad70a196 |
children | 75091bfc577b |
rev | line source |
---|---|
0 | 1 /* |
2 * MSMPEG4 backend for ffmpeg encoder and decoder | |
3 * Copyright (c) 2001 Gerard Lantau. | |
4 * | |
5 * This program is free software; you can redistribute it and/or modify | |
6 * it under the terms of the GNU General Public License as published by | |
7 * the Free Software Foundation; either version 2 of the License, or | |
8 * (at your option) any later version. | |
9 * | |
10 * This program is distributed in the hope that it will be useful, | |
11 * but WITHOUT ANY WARRANTY; without even the implied warranty of | |
12 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the | |
13 * GNU General Public License for more details. | |
14 * | |
15 * You should have received a copy of the GNU General Public License | |
16 * along with this program; if not, write to the Free Software | |
17 * Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. | |
18 */ | |
19 #include <stdlib.h> | |
20 #include <stdio.h> | |
21 #include "common.h" | |
22 #include "dsputil.h" | |
23 #include "mpegvideo.h" | |
24 | |
25 /* | |
26 * You can also call this codec : MPEG4 with a twist ! | |
27 * | |
28 * TODO: | |
29 * - (encoding) select best mv table (two choices) | |
30 * - (encoding) select best vlc/dc table | |
31 * - (decoding) handle slice indication | |
32 */ | |
33 //#define DEBUG | |
34 | |
35 /* motion vector table */ | |
36 typedef struct MVTable { | |
37 int n; | |
38 const UINT16 *table_mv_code; | |
39 const UINT8 *table_mv_bits; | |
40 const UINT8 *table_mvx; | |
41 const UINT8 *table_mvy; | |
42 UINT16 *table_mv_index; /* encoding: convert mv to index in table_mv */ | |
43 VLC vlc; /* decoding: vlc */ | |
44 } MVTable; | |
45 | |
46 static void msmpeg4_encode_block(MpegEncContext * s, DCTELEM * block, int n); | |
47 static int msmpeg4_decode_block(MpegEncContext * s, DCTELEM * block, | |
48 int n, int coded); | |
49 static int msmpeg4_decode_dc(MpegEncContext * s, int n, int *dir_ptr); | |
50 static int msmpeg4_decode_motion(MpegEncContext * s, | |
51 int *mx_ptr, int *my_ptr); | |
52 | |
225
ae145876789d
use multiply instead of divides for DC prediction on X86
michaelni
parents:
221
diff
changeset
|
53 extern UINT32 inverse[256]; |
ae145876789d
use multiply instead of divides for DC prediction on X86
michaelni
parents:
221
diff
changeset
|
54 |
0 | 55 #ifdef DEBUG |
56 int intra_count = 0; | |
57 int frame_count = 0; | |
58 #endif | |
59 /* XXX: move it to mpegvideo.h */ | |
60 | |
61 static int init_done = 0; | |
62 | |
63 #include "msmpeg4data.h" | |
64 | |
65 #ifdef STATS | |
66 | |
67 const char *st_names[ST_NB] = { | |
68 "unknown", | |
69 "dc", | |
70 "intra_ac", | |
71 "inter_ac", | |
72 "intra_mb", | |
73 "inter_mb", | |
74 "mv", | |
75 }; | |
76 | |
77 int st_current_index = 0; | |
78 unsigned int st_bit_counts[ST_NB]; | |
79 unsigned int st_out_bit_counts[ST_NB]; | |
80 | |
81 #define set_stat(var) st_current_index = var; | |
82 | |
83 void print_stats(void) | |
84 { | |
85 unsigned int total; | |
86 int i; | |
87 | |
88 printf("Input:\n"); | |
89 total = 0; | |
90 for(i=0;i<ST_NB;i++) | |
91 total += st_bit_counts[i]; | |
92 if (total == 0) | |
93 total = 1; | |
94 for(i=0;i<ST_NB;i++) { | |
95 printf("%-10s : %10.1f %5.1f%%\n", | |
96 st_names[i], | |
97 (double)st_bit_counts[i] / 8.0, | |
98 (double)st_bit_counts[i] * 100.0 / total); | |
99 } | |
100 printf("%-10s : %10.1f %5.1f%%\n", | |
101 "total", | |
102 (double)total / 8.0, | |
103 100.0); | |
104 | |
105 printf("Output:\n"); | |
106 total = 0; | |
107 for(i=0;i<ST_NB;i++) | |
108 total += st_out_bit_counts[i]; | |
109 if (total == 0) | |
110 total = 1; | |
111 for(i=0;i<ST_NB;i++) { | |
112 printf("%-10s : %10.1f %5.1f%%\n", | |
113 st_names[i], | |
114 (double)st_out_bit_counts[i] / 8.0, | |
115 (double)st_out_bit_counts[i] * 100.0 / total); | |
116 } | |
117 printf("%-10s : %10.1f %5.1f%%\n", | |
118 "total", | |
119 (double)total / 8.0, | |
120 100.0); | |
121 } | |
122 | |
123 #else | |
124 | |
125 #define set_stat(var) | |
126 | |
127 #endif | |
128 | |
129 /* build the table which associate a (x,y) motion vector to a vlc */ | |
130 static void init_mv_table(MVTable *tab) | |
131 { | |
132 int i, x, y; | |
133 | |
134 tab->table_mv_index = malloc(sizeof(UINT16) * 4096); | |
135 /* mark all entries as not used */ | |
136 for(i=0;i<4096;i++) | |
137 tab->table_mv_index[i] = tab->n; | |
138 | |
139 for(i=0;i<tab->n;i++) { | |
140 x = tab->table_mvx[i]; | |
141 y = tab->table_mvy[i]; | |
142 tab->table_mv_index[(x << 6) | y] = i; | |
143 } | |
144 } | |
145 | |
146 static void code012(PutBitContext *pb, int n) | |
147 { | |
148 if (n == 0) { | |
149 put_bits(pb, 1, 0); | |
150 } else { | |
151 put_bits(pb, 1, 1); | |
152 put_bits(pb, 1, (n >= 2)); | |
153 } | |
154 } | |
155 | |
156 /* write MSMPEG4 V3 compatible frame header */ | |
157 void msmpeg4_encode_picture_header(MpegEncContext * s, int picture_number) | |
158 { | |
159 int i; | |
160 | |
161 align_put_bits(&s->pb); | |
162 | |
163 put_bits(&s->pb, 2, s->pict_type - 1); | |
164 | |
165 put_bits(&s->pb, 5, s->qscale); | |
166 | |
167 s->rl_table_index = 2; | |
168 s->rl_chroma_table_index = 1; /* only for I frame */ | |
169 s->dc_table_index = 1; | |
170 s->mv_table_index = 1; /* only if P frame */ | |
171 s->use_skip_mb_code = 1; /* only if P frame */ | |
172 | |
173 if (s->pict_type == I_TYPE) { | |
174 put_bits(&s->pb, 5, 0x17); /* indicate only one "slice" */ | |
175 | |
176 code012(&s->pb, s->rl_chroma_table_index); | |
177 code012(&s->pb, s->rl_table_index); | |
178 | |
179 put_bits(&s->pb, 1, s->dc_table_index); | |
180 s->no_rounding = 1; | |
181 } else { | |
182 put_bits(&s->pb, 1, s->use_skip_mb_code); | |
183 | |
184 s->rl_chroma_table_index = s->rl_table_index; | |
185 code012(&s->pb, s->rl_table_index); | |
186 | |
187 put_bits(&s->pb, 1, s->dc_table_index); | |
188 | |
189 put_bits(&s->pb, 1, s->mv_table_index); | |
208 | 190 |
191 if(s->flipflop_rounding){ | |
192 s->no_rounding ^= 1; | |
193 }else{ | |
194 s->no_rounding = 0; | |
195 } | |
0 | 196 } |
197 | |
198 if (!init_done) { | |
199 /* init various encoding tables */ | |
200 init_done = 1; | |
201 init_mv_table(&mv_tables[0]); | |
202 init_mv_table(&mv_tables[1]); | |
203 for(i=0;i<NB_RL_TABLES;i++) | |
204 init_rl(&rl_table[i]); | |
205 } | |
206 | |
207 #ifdef DEBUG | |
208 intra_count = 0; | |
209 printf("*****frame %d:\n", frame_count++); | |
210 #endif | |
211 } | |
212 | |
208 | 213 void msmpeg4_encode_ext_header(MpegEncContext * s) |
214 { | |
215 if(s->pict_type == P_TYPE) | |
216 { | |
217 return; // P-Frames dont seem to have them and not even a 0 bit | |
218 } | |
219 else | |
220 { | |
221 s->flipflop_rounding=1; | |
222 s->bitrate= 910; | |
223 | |
224 put_bits(&s->pb, 1, 1); // ext header indicator | |
225 | |
226 put_bits(&s->pb, 4, 7); // ? | |
227 | |
228 put_bits(&s->pb, 11, s->bitrate); | |
229 | |
230 put_bits(&s->pb, 1, s->flipflop_rounding); | |
231 } | |
232 } | |
233 | |
0 | 234 /* predict coded block */ |
235 static inline int coded_block_pred(MpegEncContext * s, int n, UINT8 **coded_block_ptr) | |
236 { | |
237 int x, y, wrap, pred, a, b, c; | |
238 | |
239 x = 2 * s->mb_x + 1 + (n & 1); | |
240 y = 2 * s->mb_y + 1 + ((n & 2) >> 1); | |
241 wrap = s->mb_width * 2 + 2; | |
242 | |
243 /* B C | |
244 * A X | |
245 */ | |
246 a = s->coded_block[(x - 1) + (y) * wrap]; | |
247 b = s->coded_block[(x - 1) + (y - 1) * wrap]; | |
248 c = s->coded_block[(x) + (y - 1) * wrap]; | |
249 | |
250 if (b == c) { | |
251 pred = a; | |
252 } else { | |
253 pred = c; | |
254 } | |
255 | |
256 /* store value */ | |
257 *coded_block_ptr = &s->coded_block[(x) + (y) * wrap]; | |
258 | |
259 return pred; | |
260 } | |
261 | |
262 static void msmpeg4_encode_motion(MpegEncContext * s, | |
263 int mx, int my) | |
264 { | |
265 int code; | |
266 MVTable *mv; | |
267 | |
268 /* modulo encoding */ | |
269 /* WARNING : you cannot reach all the MVs even with the modulo | |
270 encoding. This is a somewhat strange compromise they took !!! */ | |
271 if (mx <= -64) | |
272 mx += 64; | |
273 else if (mx >= 64) | |
274 mx -= 64; | |
275 if (my <= -64) | |
276 my += 64; | |
277 else if (my >= 64) | |
278 my -= 64; | |
279 | |
280 mx += 32; | |
281 my += 32; | |
282 #if 0 | |
283 if ((unsigned)mx >= 64 || | |
284 (unsigned)my >= 64) | |
285 fprintf(stderr, "error mx=%d my=%d\n", mx, my); | |
286 #endif | |
287 mv = &mv_tables[s->mv_table_index]; | |
288 | |
289 code = mv->table_mv_index[(mx << 6) | my]; | |
290 set_stat(ST_MV); | |
291 put_bits(&s->pb, | |
292 mv->table_mv_bits[code], | |
293 mv->table_mv_code[code]); | |
294 if (code == mv->n) { | |
295 /* escape : code litterally */ | |
296 put_bits(&s->pb, 6, mx); | |
297 put_bits(&s->pb, 6, my); | |
298 } | |
299 } | |
300 | |
301 void msmpeg4_encode_mb(MpegEncContext * s, | |
302 DCTELEM block[6][64], | |
303 int motion_x, int motion_y) | |
304 { | |
305 int cbp, coded_cbp, i; | |
306 int pred_x, pred_y; | |
307 UINT8 *coded_block; | |
308 | |
309 if (!s->mb_intra) { | |
310 /* compute cbp */ | |
311 set_stat(ST_INTER_MB); | |
312 cbp = 0; | |
313 for (i = 0; i < 6; i++) { | |
314 if (s->block_last_index[i] >= 0) | |
315 cbp |= 1 << (5 - i); | |
316 } | |
317 if (s->use_skip_mb_code && (cbp | motion_x | motion_y) == 0) { | |
318 /* skip macroblock */ | |
319 put_bits(&s->pb, 1, 1); | |
320 return; | |
321 } | |
322 if (s->use_skip_mb_code) | |
323 put_bits(&s->pb, 1, 0); /* mb coded */ | |
324 | |
325 put_bits(&s->pb, | |
326 table_mb_non_intra[cbp + 64][1], | |
327 table_mb_non_intra[cbp + 64][0]); | |
328 | |
329 /* motion vector */ | |
330 h263_pred_motion(s, 0, &pred_x, &pred_y); | |
331 msmpeg4_encode_motion(s, motion_x - pred_x, | |
332 motion_y - pred_y); | |
333 } else { | |
334 /* compute cbp */ | |
335 cbp = 0; | |
336 coded_cbp = 0; | |
337 for (i = 0; i < 6; i++) { | |
338 int val, pred; | |
339 val = (s->block_last_index[i] >= 1); | |
340 cbp |= val << (5 - i); | |
341 if (i < 4) { | |
342 /* predict value for close blocks only for luma */ | |
343 pred = coded_block_pred(s, i, &coded_block); | |
344 *coded_block = val; | |
345 val = val ^ pred; | |
346 } | |
347 coded_cbp |= val << (5 - i); | |
348 } | |
349 #if 0 | |
350 if (coded_cbp) | |
351 printf("cbp=%x %x\n", cbp, coded_cbp); | |
352 #endif | |
353 | |
354 if (s->pict_type == I_TYPE) { | |
355 set_stat(ST_INTRA_MB); | |
356 put_bits(&s->pb, | |
357 table_mb_intra[coded_cbp][1], table_mb_intra[coded_cbp][0]); | |
358 } else { | |
359 if (s->use_skip_mb_code) | |
360 put_bits(&s->pb, 1, 0); /* mb coded */ | |
361 put_bits(&s->pb, | |
362 table_mb_non_intra[cbp][1], | |
363 table_mb_non_intra[cbp][0]); | |
364 } | |
365 set_stat(ST_INTRA_MB); | |
366 put_bits(&s->pb, 1, 0); /* no AC prediction yet */ | |
367 } | |
368 | |
369 for (i = 0; i < 6; i++) { | |
370 msmpeg4_encode_block(s, block[i], i); | |
371 } | |
372 } | |
373 | |
374 | |
375 /* strongly inspirated from MPEG4, but not exactly the same ! */ | |
376 void msmpeg4_dc_scale(MpegEncContext * s) | |
377 { | |
195
92f726205082
s->c_dc_scale was 7 if s->qscale==2 but should be 8 (the bug is visible in deep red areas in high bitrate clips) - patch by Michael Niedermayer <michaelni@gmx.at>
arpi_esp
parents:
186
diff
changeset
|
378 if (s->qscale < 5){ |
92f726205082
s->c_dc_scale was 7 if s->qscale==2 but should be 8 (the bug is visible in deep red areas in high bitrate clips) - patch by Michael Niedermayer <michaelni@gmx.at>
arpi_esp
parents:
186
diff
changeset
|
379 s->y_dc_scale = 8; |
92f726205082
s->c_dc_scale was 7 if s->qscale==2 but should be 8 (the bug is visible in deep red areas in high bitrate clips) - patch by Michael Niedermayer <michaelni@gmx.at>
arpi_esp
parents:
186
diff
changeset
|
380 s->c_dc_scale = 8; |
92f726205082
s->c_dc_scale was 7 if s->qscale==2 but should be 8 (the bug is visible in deep red areas in high bitrate clips) - patch by Michael Niedermayer <michaelni@gmx.at>
arpi_esp
parents:
186
diff
changeset
|
381 // s->c_dc_scale = (s->qscale + 13)>>1; |
92f726205082
s->c_dc_scale was 7 if s->qscale==2 but should be 8 (the bug is visible in deep red areas in high bitrate clips) - patch by Michael Niedermayer <michaelni@gmx.at>
arpi_esp
parents:
186
diff
changeset
|
382 }else if (s->qscale < 9){ |
92f726205082
s->c_dc_scale was 7 if s->qscale==2 but should be 8 (the bug is visible in deep red areas in high bitrate clips) - patch by Michael Niedermayer <michaelni@gmx.at>
arpi_esp
parents:
186
diff
changeset
|
383 s->y_dc_scale = 2 * s->qscale; |
92f726205082
s->c_dc_scale was 7 if s->qscale==2 but should be 8 (the bug is visible in deep red areas in high bitrate clips) - patch by Michael Niedermayer <michaelni@gmx.at>
arpi_esp
parents:
186
diff
changeset
|
384 s->c_dc_scale = (s->qscale + 13)>>1; |
92f726205082
s->c_dc_scale was 7 if s->qscale==2 but should be 8 (the bug is visible in deep red areas in high bitrate clips) - patch by Michael Niedermayer <michaelni@gmx.at>
arpi_esp
parents:
186
diff
changeset
|
385 }else{ |
92f726205082
s->c_dc_scale was 7 if s->qscale==2 but should be 8 (the bug is visible in deep red areas in high bitrate clips) - patch by Michael Niedermayer <michaelni@gmx.at>
arpi_esp
parents:
186
diff
changeset
|
386 s->y_dc_scale = s->qscale + 8; |
92f726205082
s->c_dc_scale was 7 if s->qscale==2 but should be 8 (the bug is visible in deep red areas in high bitrate clips) - patch by Michael Niedermayer <michaelni@gmx.at>
arpi_esp
parents:
186
diff
changeset
|
387 s->c_dc_scale = (s->qscale + 13)>>1; |
92f726205082
s->c_dc_scale was 7 if s->qscale==2 but should be 8 (the bug is visible in deep red areas in high bitrate clips) - patch by Michael Niedermayer <michaelni@gmx.at>
arpi_esp
parents:
186
diff
changeset
|
388 } |
92f726205082
s->c_dc_scale was 7 if s->qscale==2 but should be 8 (the bug is visible in deep red areas in high bitrate clips) - patch by Michael Niedermayer <michaelni@gmx.at>
arpi_esp
parents:
186
diff
changeset
|
389 // this differs for quant >24 from mpeg4 |
92f726205082
s->c_dc_scale was 7 if s->qscale==2 but should be 8 (the bug is visible in deep red areas in high bitrate clips) - patch by Michael Niedermayer <michaelni@gmx.at>
arpi_esp
parents:
186
diff
changeset
|
390 |
92f726205082
s->c_dc_scale was 7 if s->qscale==2 but should be 8 (the bug is visible in deep red areas in high bitrate clips) - patch by Michael Niedermayer <michaelni@gmx.at>
arpi_esp
parents:
186
diff
changeset
|
391 // if(s->qscale==13) s->c_dc_scale=14; |
92f726205082
s->c_dc_scale was 7 if s->qscale==2 but should be 8 (the bug is visible in deep red areas in high bitrate clips) - patch by Michael Niedermayer <michaelni@gmx.at>
arpi_esp
parents:
186
diff
changeset
|
392 |
92f726205082
s->c_dc_scale was 7 if s->qscale==2 but should be 8 (the bug is visible in deep red areas in high bitrate clips) - patch by Michael Niedermayer <michaelni@gmx.at>
arpi_esp
parents:
186
diff
changeset
|
393 // if(s->qscale>=6) |
92f726205082
s->c_dc_scale was 7 if s->qscale==2 but should be 8 (the bug is visible in deep red areas in high bitrate clips) - patch by Michael Niedermayer <michaelni@gmx.at>
arpi_esp
parents:
186
diff
changeset
|
394 // printf("%d", s->qscale); |
92f726205082
s->c_dc_scale was 7 if s->qscale==2 but should be 8 (the bug is visible in deep red areas in high bitrate clips) - patch by Michael Niedermayer <michaelni@gmx.at>
arpi_esp
parents:
186
diff
changeset
|
395 |
92f726205082
s->c_dc_scale was 7 if s->qscale==2 but should be 8 (the bug is visible in deep red areas in high bitrate clips) - patch by Michael Niedermayer <michaelni@gmx.at>
arpi_esp
parents:
186
diff
changeset
|
396 /* s->c_dc_scale values (found by Michael Nidermayer) |
92f726205082
s->c_dc_scale was 7 if s->qscale==2 but should be 8 (the bug is visible in deep red areas in high bitrate clips) - patch by Michael Niedermayer <michaelni@gmx.at>
arpi_esp
parents:
186
diff
changeset
|
397 qscale=2 -> 8 (yes iam sure about that) |
92f726205082
s->c_dc_scale was 7 if s->qscale==2 but should be 8 (the bug is visible in deep red areas in high bitrate clips) - patch by Michael Niedermayer <michaelni@gmx.at>
arpi_esp
parents:
186
diff
changeset
|
398 qscale=3 -> 8 |
92f726205082
s->c_dc_scale was 7 if s->qscale==2 but should be 8 (the bug is visible in deep red areas in high bitrate clips) - patch by Michael Niedermayer <michaelni@gmx.at>
arpi_esp
parents:
186
diff
changeset
|
399 qscale=4 -> 8 |
92f726205082
s->c_dc_scale was 7 if s->qscale==2 but should be 8 (the bug is visible in deep red areas in high bitrate clips) - patch by Michael Niedermayer <michaelni@gmx.at>
arpi_esp
parents:
186
diff
changeset
|
400 qscale=5 -> 9 |
92f726205082
s->c_dc_scale was 7 if s->qscale==2 but should be 8 (the bug is visible in deep red areas in high bitrate clips) - patch by Michael Niedermayer <michaelni@gmx.at>
arpi_esp
parents:
186
diff
changeset
|
401 qscale=6 -> 9 |
92f726205082
s->c_dc_scale was 7 if s->qscale==2 but should be 8 (the bug is visible in deep red areas in high bitrate clips) - patch by Michael Niedermayer <michaelni@gmx.at>
arpi_esp
parents:
186
diff
changeset
|
402 qscale=7 -> 10 |
92f726205082
s->c_dc_scale was 7 if s->qscale==2 but should be 8 (the bug is visible in deep red areas in high bitrate clips) - patch by Michael Niedermayer <michaelni@gmx.at>
arpi_esp
parents:
186
diff
changeset
|
403 qscale=8 -> 10 |
92f726205082
s->c_dc_scale was 7 if s->qscale==2 but should be 8 (the bug is visible in deep red areas in high bitrate clips) - patch by Michael Niedermayer <michaelni@gmx.at>
arpi_esp
parents:
186
diff
changeset
|
404 qscale=9 -> 11 |
92f726205082
s->c_dc_scale was 7 if s->qscale==2 but should be 8 (the bug is visible in deep red areas in high bitrate clips) - patch by Michael Niedermayer <michaelni@gmx.at>
arpi_esp
parents:
186
diff
changeset
|
405 qscale=10-> 11 |
92f726205082
s->c_dc_scale was 7 if s->qscale==2 but should be 8 (the bug is visible in deep red areas in high bitrate clips) - patch by Michael Niedermayer <michaelni@gmx.at>
arpi_esp
parents:
186
diff
changeset
|
406 */ |
0 | 407 } |
408 | |
409 /* dir = 0: left, dir = 1: top prediction */ | |
410 static int msmpeg4_pred_dc(MpegEncContext * s, int n, | |
25 | 411 INT16 **dc_val_ptr, int *dir_ptr) |
0 | 412 { |
246 | 413 int a, b, c, xy, wrap, pred, scale; |
25 | 414 INT16 *dc_val; |
0 | 415 |
416 /* find prediction */ | |
417 if (n < 4) { | |
418 wrap = s->mb_width * 2 + 2; | |
246 | 419 xy = 2 * s->mb_y + 1 + ((n & 2) >> 1); |
420 xy *= wrap; | |
421 xy += 2 * s->mb_x + 1 + (n & 1); | |
0 | 422 dc_val = s->dc_val[0]; |
423 scale = s->y_dc_scale; | |
424 } else { | |
425 wrap = s->mb_width + 2; | |
246 | 426 xy = s->mb_y + 1; |
427 xy *= wrap; | |
428 xy += s->mb_x + 1; | |
0 | 429 dc_val = s->dc_val[n - 4 + 1]; |
430 scale = s->c_dc_scale; | |
431 } | |
432 | |
433 /* B C | |
434 * A X | |
435 */ | |
246 | 436 a = dc_val[xy - 1]; |
437 b = dc_val[xy - 1 - wrap]; | |
438 c = dc_val[xy - wrap]; | |
0 | 439 |
440 /* XXX: the following solution consumes divisions, but it does not | |
441 necessitate to modify mpegvideo.c. The problem comes from the | |
442 fact they decided to store the quantized DC (which would lead | |
443 to problems if Q could vary !) */ | |
221 | 444 #if defined ARCH_X86 && !defined PIC |
204 | 445 asm volatile( |
446 "movl %3, %%eax \n\t" | |
447 "shrl $1, %%eax \n\t" | |
448 "addl %%eax, %2 \n\t" | |
449 "addl %%eax, %1 \n\t" | |
450 "addl %0, %%eax \n\t" | |
225
ae145876789d
use multiply instead of divides for DC prediction on X86
michaelni
parents:
221
diff
changeset
|
451 "mull %4 \n\t" |
ae145876789d
use multiply instead of divides for DC prediction on X86
michaelni
parents:
221
diff
changeset
|
452 "movl %%edx, %0 \n\t" |
204 | 453 "movl %1, %%eax \n\t" |
225
ae145876789d
use multiply instead of divides for DC prediction on X86
michaelni
parents:
221
diff
changeset
|
454 "mull %4 \n\t" |
ae145876789d
use multiply instead of divides for DC prediction on X86
michaelni
parents:
221
diff
changeset
|
455 "movl %%edx, %1 \n\t" |
204 | 456 "movl %2, %%eax \n\t" |
225
ae145876789d
use multiply instead of divides for DC prediction on X86
michaelni
parents:
221
diff
changeset
|
457 "mull %4 \n\t" |
ae145876789d
use multiply instead of divides for DC prediction on X86
michaelni
parents:
221
diff
changeset
|
458 "movl %%edx, %2 \n\t" |
228 | 459 : "+b" (a), "+c" (b), "+D" (c) |
460 : "g" (scale), "S" (inverse[scale]) | |
204 | 461 : "%eax", "%edx" |
462 ); | |
221 | 463 #else |
464 /* #elif defined (ARCH_ALPHA) */ | |
214
73df666cacc7
Alpha optimizations by Falk Hueffner <falk.hueffner@student.uni-tuebingen.de>
nickols_k
parents:
208
diff
changeset
|
465 /* Divisions are extremely costly on Alpha; optimize the most |
221 | 466 common case. But they are costly everywhere... |
467 */ | |
214
73df666cacc7
Alpha optimizations by Falk Hueffner <falk.hueffner@student.uni-tuebingen.de>
nickols_k
parents:
208
diff
changeset
|
468 if (scale == 8) { |
73df666cacc7
Alpha optimizations by Falk Hueffner <falk.hueffner@student.uni-tuebingen.de>
nickols_k
parents:
208
diff
changeset
|
469 a = (a + (8 >> 1)) / 8; |
73df666cacc7
Alpha optimizations by Falk Hueffner <falk.hueffner@student.uni-tuebingen.de>
nickols_k
parents:
208
diff
changeset
|
470 b = (b + (8 >> 1)) / 8; |
73df666cacc7
Alpha optimizations by Falk Hueffner <falk.hueffner@student.uni-tuebingen.de>
nickols_k
parents:
208
diff
changeset
|
471 c = (c + (8 >> 1)) / 8; |
73df666cacc7
Alpha optimizations by Falk Hueffner <falk.hueffner@student.uni-tuebingen.de>
nickols_k
parents:
208
diff
changeset
|
472 } else { |
73df666cacc7
Alpha optimizations by Falk Hueffner <falk.hueffner@student.uni-tuebingen.de>
nickols_k
parents:
208
diff
changeset
|
473 a = (a + (scale >> 1)) / scale; |
73df666cacc7
Alpha optimizations by Falk Hueffner <falk.hueffner@student.uni-tuebingen.de>
nickols_k
parents:
208
diff
changeset
|
474 b = (b + (scale >> 1)) / scale; |
73df666cacc7
Alpha optimizations by Falk Hueffner <falk.hueffner@student.uni-tuebingen.de>
nickols_k
parents:
208
diff
changeset
|
475 c = (c + (scale >> 1)) / scale; |
73df666cacc7
Alpha optimizations by Falk Hueffner <falk.hueffner@student.uni-tuebingen.de>
nickols_k
parents:
208
diff
changeset
|
476 } |
204 | 477 #endif |
0 | 478 /* XXX: WARNING: they did not choose the same test as MPEG4. This |
479 is very important ! */ | |
480 if (abs(a - b) <= abs(b - c)) { | |
481 pred = c; | |
482 *dir_ptr = 1; | |
483 } else { | |
484 pred = a; | |
485 *dir_ptr = 0; | |
486 } | |
487 | |
488 /* update predictor */ | |
246 | 489 *dc_val_ptr = &dc_val[xy]; |
0 | 490 return pred; |
491 } | |
492 | |
493 #define DC_MAX 119 | |
494 | |
495 static void msmpeg4_encode_dc(MpegEncContext * s, int level, int n, int *dir_ptr) | |
496 { | |
497 int sign, code; | |
498 int pred; | |
25 | 499 INT16 *dc_val; |
0 | 500 |
501 pred = msmpeg4_pred_dc(s, n, &dc_val, dir_ptr); | |
502 | |
503 /* update predictor */ | |
504 if (n < 4) { | |
505 *dc_val = level * s->y_dc_scale; | |
506 } else { | |
507 *dc_val = level * s->c_dc_scale; | |
508 } | |
509 | |
510 /* do the prediction */ | |
511 level -= pred; | |
512 | |
513 sign = 0; | |
514 if (level < 0) { | |
515 level = -level; | |
516 sign = 1; | |
517 } | |
518 | |
519 code = level; | |
520 if (code > DC_MAX) | |
521 code = DC_MAX; | |
522 | |
523 if (s->dc_table_index == 0) { | |
524 if (n < 4) { | |
525 put_bits(&s->pb, table0_dc_lum[code][1], table0_dc_lum[code][0]); | |
526 } else { | |
527 put_bits(&s->pb, table0_dc_chroma[code][1], table0_dc_chroma[code][0]); | |
528 } | |
529 } else { | |
530 if (n < 4) { | |
531 put_bits(&s->pb, table1_dc_lum[code][1], table1_dc_lum[code][0]); | |
532 } else { | |
533 put_bits(&s->pb, table1_dc_chroma[code][1], table1_dc_chroma[code][0]); | |
534 } | |
535 } | |
536 | |
537 if (code == DC_MAX) | |
538 put_bits(&s->pb, 8, level); | |
539 | |
540 if (level != 0) { | |
541 put_bits(&s->pb, 1, sign); | |
542 } | |
543 } | |
544 | |
545 /* Encoding of a block. Very similar to MPEG4 except for a different | |
546 escape coding (same as H263) and more vlc tables. | |
547 */ | |
548 static void msmpeg4_encode_block(MpegEncContext * s, DCTELEM * block, int n) | |
549 { | |
550 int level, run, last, i, j, last_index; | |
551 int last_non_zero, sign, slevel; | |
552 int code, run_diff, dc_pred_dir; | |
553 const RLTable *rl; | |
554 | |
555 if (s->mb_intra) { | |
556 set_stat(ST_DC); | |
557 msmpeg4_encode_dc(s, block[0], n, &dc_pred_dir); | |
558 i = 1; | |
559 if (n < 4) { | |
560 rl = &rl_table[s->rl_table_index]; | |
561 } else { | |
562 rl = &rl_table[3 + s->rl_chroma_table_index]; | |
563 } | |
564 run_diff = 0; | |
565 set_stat(ST_INTRA_AC); | |
566 } else { | |
567 i = 0; | |
568 rl = &rl_table[3 + s->rl_table_index]; | |
569 run_diff = 1; | |
570 set_stat(ST_INTER_AC); | |
571 } | |
572 | |
573 /* AC coefs */ | |
574 last_index = s->block_last_index[n]; | |
575 last_non_zero = i - 1; | |
576 for (; i <= last_index; i++) { | |
577 j = zigzag_direct[i]; | |
578 level = block[j]; | |
579 if (level) { | |
580 run = i - last_non_zero - 1; | |
581 last = (i == last_index); | |
582 sign = 0; | |
583 slevel = level; | |
584 if (level < 0) { | |
585 sign = 1; | |
586 level = -level; | |
587 } | |
588 code = get_rl_index(rl, last, run, level); | |
589 put_bits(&s->pb, rl->table_vlc[code][1], rl->table_vlc[code][0]); | |
590 if (code == rl->n) { | |
591 int level1, run1; | |
592 | |
593 level1 = level - rl->max_level[last][run]; | |
594 if (level1 < 1) | |
595 goto esc2; | |
596 code = get_rl_index(rl, last, run, level1); | |
597 if (code == rl->n) { | |
598 esc2: | |
599 put_bits(&s->pb, 1, 0); | |
600 if (level > MAX_LEVEL) | |
601 goto esc3; | |
602 run1 = run - rl->max_run[last][level] - run_diff; | |
603 if (run1 < 0) | |
604 goto esc3; | |
605 code = get_rl_index(rl, last, run1, level); | |
606 if (code == rl->n) { | |
607 esc3: | |
608 /* third escape */ | |
609 put_bits(&s->pb, 1, 0); | |
610 put_bits(&s->pb, 1, last); | |
611 put_bits(&s->pb, 6, run); | |
612 put_bits(&s->pb, 8, slevel & 0xff); | |
613 } else { | |
614 /* second escape */ | |
615 put_bits(&s->pb, 1, 1); | |
616 put_bits(&s->pb, rl->table_vlc[code][1], rl->table_vlc[code][0]); | |
617 put_bits(&s->pb, 1, sign); | |
618 } | |
619 } else { | |
620 /* first escape */ | |
621 put_bits(&s->pb, 1, 1); | |
622 put_bits(&s->pb, rl->table_vlc[code][1], rl->table_vlc[code][0]); | |
623 put_bits(&s->pb, 1, sign); | |
624 } | |
625 } else { | |
626 put_bits(&s->pb, 1, sign); | |
627 } | |
628 last_non_zero = i; | |
629 } | |
630 } | |
631 } | |
632 | |
633 /****************************************/ | |
634 /* decoding stuff */ | |
635 | |
636 static VLC mb_non_intra_vlc; | |
637 static VLC mb_intra_vlc; | |
638 static VLC dc_lum_vlc[2]; | |
639 static VLC dc_chroma_vlc[2]; | |
640 | |
641 /* init all vlc decoding tables */ | |
642 int msmpeg4_decode_init_vlc(MpegEncContext *s) | |
643 { | |
644 int i; | |
645 MVTable *mv; | |
646 | |
647 for(i=0;i<NB_RL_TABLES;i++) { | |
648 init_rl(&rl_table[i]); | |
649 init_vlc_rl(&rl_table[i]); | |
650 } | |
651 for(i=0;i<2;i++) { | |
652 mv = &mv_tables[i]; | |
653 init_vlc(&mv->vlc, 9, mv->n + 1, | |
654 mv->table_mv_bits, 1, 1, | |
655 mv->table_mv_code, 2, 2); | |
656 } | |
657 | |
658 init_vlc(&dc_lum_vlc[0], 9, 120, | |
659 &table0_dc_lum[0][1], 8, 4, | |
660 &table0_dc_lum[0][0], 8, 4); | |
661 init_vlc(&dc_chroma_vlc[0], 9, 120, | |
662 &table0_dc_chroma[0][1], 8, 4, | |
663 &table0_dc_chroma[0][0], 8, 4); | |
664 init_vlc(&dc_lum_vlc[1], 9, 120, | |
665 &table1_dc_lum[0][1], 8, 4, | |
666 &table1_dc_lum[0][0], 8, 4); | |
667 init_vlc(&dc_chroma_vlc[1], 9, 120, | |
668 &table1_dc_chroma[0][1], 8, 4, | |
669 &table1_dc_chroma[0][0], 8, 4); | |
670 | |
671 init_vlc(&mb_non_intra_vlc, 9, 128, | |
672 &table_mb_non_intra[0][1], 8, 4, | |
673 &table_mb_non_intra[0][0], 8, 4); | |
48 | 674 init_vlc(&mb_intra_vlc, 9, 64, |
0 | 675 &table_mb_intra[0][1], 4, 2, |
676 &table_mb_intra[0][0], 4, 2); | |
677 return 0; | |
678 } | |
679 | |
680 static int decode012(GetBitContext *gb) | |
681 { | |
682 int n; | |
21 | 683 n = get_bits1(gb); |
0 | 684 if (n == 0) |
685 return 0; | |
686 else | |
21 | 687 return get_bits1(gb) + 1; |
0 | 688 } |
689 | |
690 int msmpeg4_decode_picture_header(MpegEncContext * s) | |
691 { | |
692 int code; | |
693 | |
694 s->pict_type = get_bits(&s->gb, 2) + 1; | |
695 if (s->pict_type != I_TYPE && | |
696 s->pict_type != P_TYPE) | |
697 return -1; | |
698 | |
699 s->qscale = get_bits(&s->gb, 5); | |
700 | |
701 if (s->pict_type == I_TYPE) { | |
702 code = get_bits(&s->gb, 5); | |
703 /* 0x17: one slice, 0x18: three slices */ | |
704 /* XXX: implement it */ | |
200 | 705 //printf("%d %d %d\n", code, s->slice_height, s->first_slice_line); |
0 | 706 if (code < 0x17) |
707 return -1; | |
708 s->slice_height = s->mb_height / (code - 0x16); | |
709 s->rl_chroma_table_index = decode012(&s->gb); | |
710 s->rl_table_index = decode012(&s->gb); | |
711 | |
21 | 712 s->dc_table_index = get_bits1(&s->gb); |
0 | 713 s->no_rounding = 1; |
200 | 714 /* printf(" %d %d %d %d \n", |
715 s->qscale, | |
716 s->rl_chroma_table_index, | |
717 s->rl_table_index, | |
718 s->dc_table_index);*/ | |
0 | 719 } else { |
21 | 720 s->use_skip_mb_code = get_bits1(&s->gb); |
0 | 721 |
722 s->rl_table_index = decode012(&s->gb); | |
723 s->rl_chroma_table_index = s->rl_table_index; | |
724 | |
21 | 725 s->dc_table_index = get_bits1(&s->gb); |
0 | 726 |
21 | 727 s->mv_table_index = get_bits1(&s->gb); |
200 | 728 /* printf(" %d %d %d %d %d \n", |
729 s->use_skip_mb_code, | |
730 s->rl_table_index, | |
731 s->rl_chroma_table_index, | |
732 s->dc_table_index, | |
733 s->mv_table_index);*/ | |
208 | 734 if(s->flipflop_rounding){ |
735 s->no_rounding ^= 1; | |
736 }else{ | |
737 s->no_rounding = 0; | |
738 } | |
739 // printf("%d", s->no_rounding); | |
0 | 740 } |
208 | 741 |
742 | |
0 | 743 #ifdef DEBUG |
744 printf("*****frame %d:\n", frame_count++); | |
745 #endif | |
746 return 0; | |
747 } | |
748 | |
208 | 749 int msmpeg4_decode_ext_header(MpegEncContext * s, int buf_size) |
750 { | |
751 int firstBit=0; | |
752 | |
753 /* the alt_bitstream reader could read over the end so we need to check it */ | |
754 if(get_bits_count(&s->gb) < buf_size*8) firstBit= get_bits1(&s->gb); | |
755 | |
756 if(s->pict_type == P_TYPE) | |
757 { | |
758 if(firstBit) return -1; // havnt seen ext headers in P-Frames yet ;) | |
759 } | |
760 else | |
761 { | |
762 int unk; | |
763 if(!firstBit) // no header found | |
764 { | |
765 s->flipflop_rounding= 0; | |
766 s->bitrate= 0; | |
767 return 0; | |
768 } | |
769 | |
770 unk= get_bits(&s->gb, 4); | |
771 s->bitrate= get_bits(&s->gb, 11); | |
772 | |
773 // printf("%2d %4d ;; %1X\n", unk,s->bitrate, unk); | |
774 | |
775 s->flipflop_rounding= get_bits1(&s->gb); | |
776 } | |
777 | |
778 return 0; | |
779 } | |
780 | |
246 | 781 static inline void memsetw(short *tab, int val, int n) |
0 | 782 { |
783 int i; | |
784 for(i=0;i<n;i++) | |
785 tab[i] = val; | |
786 } | |
787 | |
788 int msmpeg4_decode_mb(MpegEncContext *s, | |
789 DCTELEM block[6][64]) | |
790 { | |
791 int cbp, code, i; | |
792 UINT8 *coded_val; | |
793 | |
794 /* special slice handling */ | |
795 if (s->mb_x == 0) { | |
122 | 796 if (s->slice_height && (s->mb_y % s->slice_height) == 0) { |
0 | 797 int wrap; |
798 /* reset DC pred (set previous line to 1024) */ | |
799 wrap = 2 * s->mb_width + 2; | |
800 memsetw(&s->dc_val[0][(1) + (2 * s->mb_y) * wrap], | |
801 1024, 2 * s->mb_width); | |
802 wrap = s->mb_width + 2; | |
803 memsetw(&s->dc_val[1][(1) + (s->mb_y) * wrap], | |
804 1024, s->mb_width); | |
805 memsetw(&s->dc_val[2][(1) + (s->mb_y) * wrap], | |
806 1024, s->mb_width); | |
186
cf37da86d990
fix slices when code=0x18, patch by Michael Niedermayer <michael@mplayer.dev.hu>
arpi_esp
parents:
122
diff
changeset
|
807 |
cf37da86d990
fix slices when code=0x18, patch by Michael Niedermayer <michael@mplayer.dev.hu>
arpi_esp
parents:
122
diff
changeset
|
808 /* reset AC pred (set previous line to 0) */ |
cf37da86d990
fix slices when code=0x18, patch by Michael Niedermayer <michael@mplayer.dev.hu>
arpi_esp
parents:
122
diff
changeset
|
809 wrap = s->mb_width * 2 + 2; |
cf37da86d990
fix slices when code=0x18, patch by Michael Niedermayer <michael@mplayer.dev.hu>
arpi_esp
parents:
122
diff
changeset
|
810 memsetw(s->ac_val[0][0] + (1 + (2 * s->mb_y) * wrap)*16, |
cf37da86d990
fix slices when code=0x18, patch by Michael Niedermayer <michael@mplayer.dev.hu>
arpi_esp
parents:
122
diff
changeset
|
811 0, 2 * s->mb_width*16); |
cf37da86d990
fix slices when code=0x18, patch by Michael Niedermayer <michael@mplayer.dev.hu>
arpi_esp
parents:
122
diff
changeset
|
812 wrap = s->mb_width + 2; |
cf37da86d990
fix slices when code=0x18, patch by Michael Niedermayer <michael@mplayer.dev.hu>
arpi_esp
parents:
122
diff
changeset
|
813 memsetw(s->ac_val[1][0] + (1 + (s->mb_y) * wrap)*16, |
cf37da86d990
fix slices when code=0x18, patch by Michael Niedermayer <michael@mplayer.dev.hu>
arpi_esp
parents:
122
diff
changeset
|
814 0, s->mb_width*16); |
cf37da86d990
fix slices when code=0x18, patch by Michael Niedermayer <michael@mplayer.dev.hu>
arpi_esp
parents:
122
diff
changeset
|
815 memsetw(s->ac_val[2][0] + (1 + (s->mb_y) * wrap)*16, |
cf37da86d990
fix slices when code=0x18, patch by Michael Niedermayer <michael@mplayer.dev.hu>
arpi_esp
parents:
122
diff
changeset
|
816 0, s->mb_width*16); |
0 | 817 |
818 s->first_slice_line = 1; | |
819 } else { | |
820 s->first_slice_line = 0; | |
821 } | |
822 } | |
823 | |
824 if (s->pict_type == P_TYPE) { | |
825 set_stat(ST_INTER_MB); | |
826 if (s->use_skip_mb_code) { | |
21 | 827 if (get_bits1(&s->gb)) { |
0 | 828 /* skip mb */ |
829 s->mb_intra = 0; | |
830 for(i=0;i<6;i++) | |
831 s->block_last_index[i] = -1; | |
832 s->mv_dir = MV_DIR_FORWARD; | |
833 s->mv_type = MV_TYPE_16X16; | |
834 s->mv[0][0][0] = 0; | |
835 s->mv[0][0][1] = 0; | |
7
1d3ac9654178
added skip macroblock optimization (big perf win on black regions for example)
glantau
parents:
0
diff
changeset
|
836 s->mb_skiped = 1; |
0 | 837 return 0; |
838 } | |
839 } | |
840 | |
841 code = get_vlc(&s->gb, &mb_non_intra_vlc); | |
842 if (code < 0) | |
843 return -1; | |
246 | 844 //s->mb_intra = (code & 0x40) ? 0 : 1; |
845 s->mb_intra = (~code & 0x40) >> 6; | |
0 | 846 |
847 cbp = code & 0x3f; | |
848 } else { | |
849 set_stat(ST_INTRA_MB); | |
850 s->mb_intra = 1; | |
851 code = get_vlc(&s->gb, &mb_intra_vlc); | |
852 if (code < 0) | |
853 return -1; | |
854 /* predict coded block pattern */ | |
855 cbp = 0; | |
856 for(i=0;i<6;i++) { | |
246 | 857 int val = ((code >> (5 - i)) & 1); |
0 | 858 if (i < 4) { |
246 | 859 int pred = coded_block_pred(s, i, &coded_val); |
0 | 860 val = val ^ pred; |
861 *coded_val = val; | |
862 } | |
863 cbp |= val << (5 - i); | |
864 } | |
865 } | |
866 | |
867 if (!s->mb_intra) { | |
868 int mx, my; | |
869 set_stat(ST_MV); | |
870 h263_pred_motion(s, 0, &mx, &my); | |
871 if (msmpeg4_decode_motion(s, &mx, &my) < 0) | |
872 return -1; | |
873 s->mv_dir = MV_DIR_FORWARD; | |
874 s->mv_type = MV_TYPE_16X16; | |
875 s->mv[0][0][0] = mx; | |
876 s->mv[0][0][1] = my; | |
877 } else { | |
878 set_stat(ST_INTRA_MB); | |
21 | 879 s->ac_pred = get_bits1(&s->gb); |
0 | 880 } |
881 | |
882 for (i = 0; i < 6; i++) { | |
883 if (msmpeg4_decode_block(s, block[i], i, (cbp >> (5 - i)) & 1) < 0) | |
246 | 884 { |
885 fprintf(stderr,"\nIgnoring error while decoding block: %d x %d (%d)\n", s->mb_x, s->mb_y, i); | |
886 // return -1; | |
887 } | |
0 | 888 } |
889 return 0; | |
890 } | |
891 | |
892 static int msmpeg4_decode_block(MpegEncContext * s, DCTELEM * block, | |
893 int n, int coded) | |
894 { | |
895 int code, level, i, j, last, run, run_diff; | |
896 int dc_pred_dir; | |
897 RLTable *rl; | |
898 const UINT8 *scan_table; | |
200 | 899 int qmul, qadd; |
0 | 900 |
901 if (s->mb_intra) { | |
200 | 902 qmul=1; |
903 qadd=0; | |
904 | |
0 | 905 /* DC coef */ |
906 set_stat(ST_DC); | |
907 level = msmpeg4_decode_dc(s, n, &dc_pred_dir); | |
908 if (level < 0) | |
909 return -1; | |
910 block[0] = level; | |
911 if (n < 4) { | |
912 rl = &rl_table[s->rl_table_index]; | |
913 } else { | |
914 rl = &rl_table[3 + s->rl_chroma_table_index]; | |
915 } | |
200 | 916 |
0 | 917 run_diff = 0; |
918 i = 1; | |
919 if (!coded) { | |
920 goto not_coded; | |
921 } | |
922 if (s->ac_pred) { | |
923 if (dc_pred_dir == 0) | |
924 scan_table = ff_alternate_vertical_scan; /* left */ | |
925 else | |
926 scan_table = ff_alternate_horizontal_scan; /* top */ | |
927 } else { | |
928 scan_table = zigzag_direct; | |
929 } | |
930 set_stat(ST_INTRA_AC); | |
931 } else { | |
200 | 932 qmul = s->qscale << 1; |
933 qadd = (s->qscale - 1) | 1; | |
0 | 934 i = 0; |
935 rl = &rl_table[3 + s->rl_table_index]; | |
936 run_diff = 1; | |
937 if (!coded) { | |
938 s->block_last_index[n] = i - 1; | |
939 return 0; | |
940 } | |
941 scan_table = zigzag_direct; | |
942 set_stat(ST_INTER_AC); | |
943 } | |
944 | |
945 for(;;) { | |
946 code = get_vlc(&s->gb, &rl->vlc); | |
947 if (code < 0) | |
948 return -1; | |
949 if (code == rl->n) { | |
950 /* escape */ | |
21 | 951 if (get_bits1(&s->gb) == 0) { |
952 if (get_bits1(&s->gb) == 0) { | |
0 | 953 /* third escape */ |
21 | 954 last = get_bits1(&s->gb); |
0 | 955 run = get_bits(&s->gb, 6); |
956 level = get_bits(&s->gb, 8); | |
957 level = (level << 24) >> 24; /* sign extend */ | |
246 | 958 //level = level * qmul + (level>0) * qadd - (level<=0) * qadd ; |
959 if (level>0) level= level * qmul + qadd; | |
200 | 960 else level= level * qmul - qadd; |
0 | 961 } else { |
962 /* second escape */ | |
963 code = get_vlc(&s->gb, &rl->vlc); | |
964 if (code < 0 || code >= rl->n) | |
965 return -1; | |
966 run = rl->table_run[code]; | |
201 | 967 level = rl->table_level[code]; |
0 | 968 last = code >= rl->last; |
969 run += rl->max_run[last][level] + run_diff; | |
201 | 970 level= level * qmul + qadd; |
21 | 971 if (get_bits1(&s->gb)) |
0 | 972 level = -level; |
973 } | |
974 } else { | |
975 /* first escape */ | |
976 code = get_vlc(&s->gb, &rl->vlc); | |
977 if (code < 0 || code >= rl->n) | |
978 return -1; | |
979 run = rl->table_run[code]; | |
980 level = rl->table_level[code]; | |
981 last = code >= rl->last; | |
982 level += rl->max_level[last][run]; | |
200 | 983 level= level * qmul + qadd; |
21 | 984 if (get_bits1(&s->gb)) |
0 | 985 level = -level; |
986 } | |
987 } else { | |
988 run = rl->table_run[code]; | |
200 | 989 level = rl->table_level[code] * qmul + qadd; |
0 | 990 last = code >= rl->last; |
21 | 991 if (get_bits1(&s->gb)) |
0 | 992 level = -level; |
993 } | |
994 i += run; | |
995 if (i >= 64) | |
996 return -1; | |
997 j = scan_table[i]; | |
998 block[j] = level; | |
999 i++; | |
1000 if (last) | |
1001 break; | |
1002 } | |
1003 not_coded: | |
1004 if (s->mb_intra) { | |
1005 mpeg4_pred_ac(s, block, n, dc_pred_dir); | |
1006 if (s->ac_pred) { | |
1007 i = 64; /* XXX: not optimal */ | |
1008 } | |
1009 } | |
1010 s->block_last_index[n] = i - 1; | |
1011 | |
1012 return 0; | |
1013 } | |
1014 | |
1015 static int msmpeg4_decode_dc(MpegEncContext * s, int n, int *dir_ptr) | |
1016 { | |
1017 int level, pred; | |
25 | 1018 INT16 *dc_val; |
0 | 1019 |
1020 if (n < 4) { | |
1021 level = get_vlc(&s->gb, &dc_lum_vlc[s->dc_table_index]); | |
1022 } else { | |
1023 level = get_vlc(&s->gb, &dc_chroma_vlc[s->dc_table_index]); | |
1024 } | |
1025 if (level < 0) | |
1026 return -1; | |
1027 | |
1028 if (level == DC_MAX) { | |
1029 level = get_bits(&s->gb, 8); | |
21 | 1030 if (get_bits1(&s->gb)) |
0 | 1031 level = -level; |
1032 } else if (level != 0) { | |
21 | 1033 if (get_bits1(&s->gb)) |
0 | 1034 level = -level; |
1035 } | |
1036 | |
1037 pred = msmpeg4_pred_dc(s, n, &dc_val, dir_ptr); | |
1038 level += pred; | |
1039 | |
1040 /* update predictor */ | |
1041 if (n < 4) { | |
1042 *dc_val = level * s->y_dc_scale; | |
1043 } else { | |
1044 *dc_val = level * s->c_dc_scale; | |
1045 } | |
1046 | |
1047 return level; | |
1048 } | |
1049 | |
1050 static int msmpeg4_decode_motion(MpegEncContext * s, | |
1051 int *mx_ptr, int *my_ptr) | |
1052 { | |
1053 MVTable *mv; | |
1054 int code, mx, my; | |
1055 | |
1056 mv = &mv_tables[s->mv_table_index]; | |
1057 | |
1058 code = get_vlc(&s->gb, &mv->vlc); | |
1059 if (code < 0) | |
1060 return -1; | |
1061 if (code == mv->n) { | |
1062 mx = get_bits(&s->gb, 6); | |
1063 my = get_bits(&s->gb, 6); | |
1064 } else { | |
1065 mx = mv->table_mvx[code]; | |
1066 my = mv->table_mvy[code]; | |
1067 } | |
1068 | |
1069 mx += *mx_ptr - 32; | |
1070 my += *my_ptr - 32; | |
1071 /* WARNING : they do not do exactly modulo encoding */ | |
1072 if (mx <= -64) | |
1073 mx += 64; | |
1074 else if (mx >= 64) | |
1075 mx -= 64; | |
1076 | |
1077 if (my <= -64) | |
1078 my += 64; | |
1079 else if (my >= 64) | |
1080 my -= 64; | |
1081 *mx_ptr = mx; | |
1082 *my_ptr = my; | |
1083 return 0; | |
1084 } |