Parent Directory | Revision Log
Revision 1132 - (view) (download)
1 : | edgomez | 1054 | /***************************************************************************** |
2 : | edgomez | 851 | * |
3 : | * XVID MPEG-4 VIDEO CODEC | ||
4 : | edgomez | 1054 | * - Motion Estimation related header - |
5 : | edgomez | 851 | * |
6 : | edgomez | 1054 | * Copyright(C) 2002 Christoph Lampert <gruel@web.de> |
7 : | * 2002 Michael Militzer <michael@xvid.org> | ||
8 : | * 2002-2003 Radoslaw Czyz <xvid@syskin.cjb.net> | ||
9 : | edgomez | 851 | * |
10 : | edgomez | 1054 | * This program is free software ; you can redistribute it and/or modify |
11 : | edgomez | 851 | * it under the terms of the GNU General Public License as published by |
12 : | edgomez | 1054 | * the Free Software Foundation ; either version 2 of the License, or |
13 : | edgomez | 851 | * (at your option) any later version. |
14 : | * | ||
15 : | * This program is distributed in the hope that it will be useful, | ||
16 : | edgomez | 1054 | * but WITHOUT ANY WARRANTY ; without even the implied warranty of |
17 : | edgomez | 851 | * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the |
18 : | * GNU General Public License for more details. | ||
19 : | * | ||
20 : | * You should have received a copy of the GNU General Public License | ||
21 : | edgomez | 1054 | * along with this program ; if not, write to the Free Software |
22 : | edgomez | 851 | * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA |
23 : | * | ||
24 : | syskin | 1132 | * $Id: motion_est.h,v 1.3.2.18 2003-08-28 11:06:16 syskin Exp $ |
25 : | edgomez | 851 | * |
26 : | edgomez | 1054 | ****************************************************************************/ |
27 : | edgomez | 851 | |
28 : | #ifndef _MOTION_EST_H_ | ||
29 : | #define _MOTION_EST_H_ | ||
30 : | |||
31 : | #include "../portab.h" | ||
32 : | #include "../global.h" | ||
33 : | #include "../image/reduced.h" | ||
34 : | |||
35 : | /* hard coded motion search parameters for motion_est and smp_motion_est */ | ||
36 : | |||
37 : | edgomez | 1053 | /* very large value */ |
38 : | edgomez | 851 | #define MV_MAX_ERROR (4096 * 256) |
39 : | |||
40 : | /* INTER bias for INTER/INTRA decision; mpeg4 spec suggests 2*nb */ | ||
41 : | #define MV16_INTER_BIAS 512 | ||
42 : | |||
43 : | /* vector map (vlc delta size) smoother parameters ! float !*/ | ||
44 : | #define NEIGH_TEND_16X16 10.5 | ||
45 : | #define NEIGH_TEND_8X8 40.0 | ||
46 : | #define NEIGH_8X8_BIAS 30 | ||
47 : | |||
48 : | edgomez | 1023 | #define BITS_MULT 16 |
49 : | |||
50 : | edgomez | 851 | /* Parameters which control inter/inter4v decision */ |
51 : | #define IMV16X16 2 | ||
52 : | |||
53 : | static const int lambda_vec16[32] = | ||
54 : | { 0 ,(int)(1.00235 * NEIGH_TEND_16X16 + 0.5), | ||
55 : | (int)(1.15582*NEIGH_TEND_16X16 + 0.5), (int)(1.31976*NEIGH_TEND_16X16 + 0.5), | ||
56 : | (int)(1.49591*NEIGH_TEND_16X16 + 0.5), (int)(1.68601*NEIGH_TEND_16X16 + 0.5), | ||
57 : | (int)(1.89187*NEIGH_TEND_16X16 + 0.5), (int)(2.11542*NEIGH_TEND_16X16 + 0.5), | ||
58 : | (int)(2.35878*NEIGH_TEND_16X16 + 0.5), (int)(2.62429*NEIGH_TEND_16X16 + 0.5), | ||
59 : | (int)(2.91455*NEIGH_TEND_16X16 + 0.5), (int)(3.23253*NEIGH_TEND_16X16 + 0.5), | ||
60 : | (int)(3.58158*NEIGH_TEND_16X16 + 0.5), (int)(3.96555*NEIGH_TEND_16X16 + 0.5), | ||
61 : | (int)(4.38887*NEIGH_TEND_16X16 + 0.5), (int)(4.85673*NEIGH_TEND_16X16 + 0.5), | ||
62 : | (int)(5.37519*NEIGH_TEND_16X16 + 0.5), (int)(5.95144*NEIGH_TEND_16X16 + 0.5), | ||
63 : | (int)(6.59408*NEIGH_TEND_16X16 + 0.5), (int)(7.31349*NEIGH_TEND_16X16 + 0.5), | ||
64 : | (int)(8.12242*NEIGH_TEND_16X16 + 0.5), (int)(9.03669*NEIGH_TEND_16X16 + 0.5), | ||
65 : | (int)(10.0763*NEIGH_TEND_16X16 + 0.5), (int)(11.2669*NEIGH_TEND_16X16 + 0.5), | ||
66 : | (int)(12.6426*NEIGH_TEND_16X16 + 0.5), (int)(14.2493*NEIGH_TEND_16X16 + 0.5), | ||
67 : | (int)(16.1512*NEIGH_TEND_16X16 + 0.5), (int)(18.442*NEIGH_TEND_16X16 + 0.5), | ||
68 : | (int)(21.2656*NEIGH_TEND_16X16 + 0.5), (int)(24.8580*NEIGH_TEND_16X16 + 0.5), | ||
69 : | (int)(29.6436*NEIGH_TEND_16X16 + 0.5), (int)(36.4949*NEIGH_TEND_16X16 + 0.5) }; | ||
70 : | |||
71 : | static const int lambda_vec8[32] = | ||
72 : | { 0 ,(int)(1.00235 * NEIGH_TEND_8X8 + 0.5), | ||
73 : | (int)(1.15582 + NEIGH_TEND_8X8 + 0.5), (int)(1.31976*NEIGH_TEND_8X8 + 0.5), | ||
74 : | (int)(1.49591*NEIGH_TEND_8X8 + 0.5), (int)(1.68601*NEIGH_TEND_8X8 + 0.5), | ||
75 : | (int)(1.89187*NEIGH_TEND_8X8 + 0.5), (int)(2.11542*NEIGH_TEND_8X8 + 0.5), | ||
76 : | (int)(2.35878*NEIGH_TEND_8X8 + 0.5), (int)(2.62429*NEIGH_TEND_8X8 + 0.5), | ||
77 : | (int)(2.91455*NEIGH_TEND_8X8 + 0.5), (int)(3.23253*NEIGH_TEND_8X8 + 0.5), | ||
78 : | (int)(3.58158*NEIGH_TEND_8X8 + 0.5), (int)(3.96555*NEIGH_TEND_8X8 + 0.5), | ||
79 : | (int)(4.38887*NEIGH_TEND_8X8 + 0.5), (int)(4.85673*NEIGH_TEND_8X8 + 0.5), | ||
80 : | (int)(5.37519*NEIGH_TEND_8X8 + 0.5), (int)(5.95144*NEIGH_TEND_8X8 + 0.5), | ||
81 : | (int)(6.59408*NEIGH_TEND_8X8 + 0.5), (int)(7.31349*NEIGH_TEND_8X8 + 0.5), | ||
82 : | (int)(8.12242*NEIGH_TEND_8X8 + 0.5), (int)(9.03669*NEIGH_TEND_8X8 + 0.5), | ||
83 : | (int)(10.0763*NEIGH_TEND_8X8 + 0.5), (int)(11.2669*NEIGH_TEND_8X8 + 0.5), | ||
84 : | (int)(12.6426*NEIGH_TEND_8X8 + 0.5), (int)(14.2493*NEIGH_TEND_8X8 + 0.5), | ||
85 : | (int)(16.1512*NEIGH_TEND_8X8 + 0.5), (int)(18.442*NEIGH_TEND_8X8 + 0.5), | ||
86 : | (int)(21.2656*NEIGH_TEND_8X8 + 0.5), (int)(24.8580*NEIGH_TEND_8X8 + 0.5), | ||
87 : | (int)(29.6436*NEIGH_TEND_8X8 + 0.5), (int)(36.4949*NEIGH_TEND_8X8 + 0.5) }; | ||
88 : | |||
89 : | edgomez | 1053 | /* mv.length table */ |
90 : | edgomez | 974 | static const int mvtab[64] = { |
91 : | 1, 2, 3, 4, 6, 7, 7, 7, | ||
92 : | 9, 9, 9, 10, 10, 10, 10, 10, | ||
93 : | 10, 10, 10, 10, 10, 10, 10, 10, | ||
94 : | 10, 11, 11, 11, 11, 11, 11, 12, | ||
95 : | 12, 12, 12, 12, 12, 12, 12, 12, | ||
96 : | 12, 12, 12, 12, 12, 12, 12, 12, | ||
97 : | 12, 12, 12, 12, 12, 12, 12, 12, 12 }; | ||
98 : | edgomez | 851 | |
99 : | static const int DQtab[4] = { | ||
100 : | -1, -2, 1, 2 | ||
101 : | }; | ||
102 : | |||
103 : | #define RRV_MV_SCALEDOWN(a) ( (a)>=0 ? (a+1)/2 : (a-1)/2 ) | ||
104 : | |||
105 : | typedef struct | ||
106 : | { | ||
107 : | edgomez | 1053 | /* general fields */ |
108 : | edgomez | 851 | int max_dx, min_dx, max_dy, min_dy; |
109 : | uint32_t rounding; | ||
110 : | VECTOR predMV; | ||
111 : | VECTOR * currentMV; | ||
112 : | VECTOR * currentQMV; | ||
113 : | int32_t * iMinSAD; | ||
114 : | edgomez | 1053 | const uint8_t * RefP[6]; /* N, V, H, HV, cU, cV */ |
115 : | edgomez | 851 | const uint8_t * CurU; |
116 : | const uint8_t * CurV; | ||
117 : | uint8_t * RefQ; | ||
118 : | const uint8_t * Cur; | ||
119 : | uint32_t lambda16; | ||
120 : | uint32_t lambda8; | ||
121 : | uint32_t iEdgedWidth; | ||
122 : | uint32_t iFcode; | ||
123 : | int * temp; | ||
124 : | syskin | 1132 | int * dir; |
125 : | edgomez | 851 | int qpel, qpel_precision; |
126 : | int chroma; | ||
127 : | int rrv; | ||
128 : | edgomez | 1053 | |
129 : | /* fields for interpolate and direct modes */ | ||
130 : | const uint8_t * b_RefP[6]; /* N, V, H, HV, cU, cV */ | ||
131 : | edgomez | 851 | VECTOR bpredMV; |
132 : | uint32_t bFcode; | ||
133 : | edgomez | 1053 | |
134 : | /* fields for direct mode */ | ||
135 : | edgomez | 851 | VECTOR directmvF[4]; |
136 : | VECTOR directmvB[4]; | ||
137 : | const VECTOR * referencemv; | ||
138 : | edgomez | 1053 | |
139 : | /* BITS/R-D stuff */ | ||
140 : | edgomez | 959 | int16_t * dctSpace; |
141 : | edgomez | 1022 | uint32_t iQuant; |
142 : | uint32_t quant_type; | ||
143 : | syskin | 1084 | int * cbp; |
144 : | edgomez | 851 | |
145 : | } SearchData; | ||
146 : | |||
147 : | |||
148 : | typedef void(CheckFunc)(const int x, const int y, | ||
149 : | syskin | 1132 | const SearchData * const Data, |
150 : | const int Direction); | ||
151 : | edgomez | 851 | |
152 : | /* | ||
153 : | * Calculate the min/max range | ||
154 : | * relative to the _MACROBLOCK_ position | ||
155 : | */ | ||
156 : | static void __inline | ||
157 : | get_range(int32_t * const min_dx, | ||
158 : | int32_t * const max_dx, | ||
159 : | int32_t * const min_dy, | ||
160 : | int32_t * const max_dy, | ||
161 : | const uint32_t x, | ||
162 : | const uint32_t y, | ||
163 : | syskin | 1129 | uint32_t block_sz, /* block dimension, 3(8) or 4(16) */ |
164 : | edgomez | 851 | const uint32_t width, |
165 : | const uint32_t height, | ||
166 : | const uint32_t fcode, | ||
167 : | syskin | 1129 | const int precision, /* 2 for qpel, 1 for halfpel */ |
168 : | edgomez | 851 | const int rrv) |
169 : | { | ||
170 : | syskin | 1129 | int k; |
171 : | const int search_range = 16 << fcode; | ||
172 : | edgomez | 851 | int high = search_range - 1; |
173 : | int low = -search_range; | ||
174 : | |||
175 : | if (rrv) { | ||
176 : | high = RRV_MV_SCALEUP(high); | ||
177 : | low = RRV_MV_SCALEUP(low); | ||
178 : | syskin | 1129 | block_sz++; |
179 : | edgomez | 851 | } |
180 : | |||
181 : | syskin | 1129 | k = (int)(width - (x<<block_sz))<<precision; |
182 : | edgomez | 1130 | *max_dx = MIN(high, k); |
183 : | syskin | 1129 | k = (int)(height - (y<<block_sz))<<precision; |
184 : | edgomez | 1130 | *max_dy = MIN(high, k); |
185 : | edgomez | 851 | |
186 : | syskin | 1129 | k = (-(int)((x+1)<<block_sz))<<precision; |
187 : | edgomez | 1130 | *min_dx = MAX(low, k); |
188 : | syskin | 1129 | k = (-(int)((y+1)<<block_sz))<<precision; |
189 : | edgomez | 1130 | *min_dy = MAX(low, k); |
190 : | edgomez | 851 | } |
191 : | |||
192 : | syskin | 1132 | typedef void |
193 : | MainSearchFunc(int x, int y, const SearchData * const Data, | ||
194 : | int bDirection, CheckFunc * const CheckCandidate); | ||
195 : | edgomez | 851 | |
196 : | static MainSearchFunc DiamondSearch, AdvDiamondSearch, SquareSearch; | ||
197 : | |||
198 : | static void Search8(const SearchData * const OldData, | ||
199 : | const int x, const int y, | ||
200 : | const uint32_t MotionFlags, | ||
201 : | const MBParam * const pParam, | ||
202 : | MACROBLOCK * const pMB, | ||
203 : | const MACROBLOCK * const pMBs, | ||
204 : | const int block, | ||
205 : | SearchData * const Data); | ||
206 : | |||
207 : | bool | ||
208 : | MotionEstimation(MBParam * const pParam, | ||
209 : | syskin | 884 | FRAMEINFO * const current, |
210 : | FRAMEINFO * const reference, | ||
211 : | const IMAGE * const pRefH, | ||
212 : | const IMAGE * const pRefV, | ||
213 : | const IMAGE * const pRefHV, | ||
214 : | chl | 1077 | const IMAGE * const pGMC, |
215 : | syskin | 884 | const uint32_t iLimit); |
216 : | edgomez | 851 | |
217 : | static void | ||
218 : | SearchP(const IMAGE * const pRef, | ||
219 : | const uint8_t * const pRefH, | ||
220 : | const uint8_t * const pRefV, | ||
221 : | const uint8_t * const pRefHV, | ||
222 : | const IMAGE * const pCur, | ||
223 : | const int x, | ||
224 : | const int y, | ||
225 : | const uint32_t MotionFlags, | ||
226 : | edgomez | 1022 | const uint32_t VopFlags, |
227 : | const uint32_t VolFlags, | ||
228 : | edgomez | 851 | SearchData * const Data, |
229 : | const MBParam * const pParam, | ||
230 : | const MACROBLOCK * const pMBs, | ||
231 : | const MACROBLOCK * const prevMBs, | ||
232 : | MACROBLOCK * const pMB); | ||
233 : | |||
234 : | chl | 1077 | static __inline void |
235 : | GMEanalyzeMB (const uint8_t * const pCur, | ||
236 : | const uint8_t * const pRef, | ||
237 : | const uint8_t * const pRefH, | ||
238 : | const uint8_t * const pRefV, | ||
239 : | const uint8_t * const pRefHV, | ||
240 : | const int x, | ||
241 : | const int y, | ||
242 : | const MBParam * const pParam, | ||
243 : | MACROBLOCK * const pMBs, | ||
244 : | SearchData * const Data); | ||
245 : | |||
246 : | void | ||
247 : | GMEanalysis(const MBParam * const pParam, | ||
248 : | const FRAMEINFO * const current, | ||
249 : | const FRAMEINFO * const reference, | ||
250 : | const IMAGE * const pRefH, | ||
251 : | const IMAGE * const pRefV, | ||
252 : | const IMAGE * const pRefHV); | ||
253 : | |||
254 : | |||
255 : | |||
256 : | WARPPOINTS | ||
257 : | GlobalMotionEst(MACROBLOCK * const pMBs, | ||
258 : | edgomez | 851 | const MBParam * const pParam, |
259 : | const FRAMEINFO * const current, | ||
260 : | const FRAMEINFO * const reference, | ||
261 : | const IMAGE * const pRefH, | ||
262 : | const IMAGE * const pRefV, | ||
263 : | const IMAGE * const pRefHV ); | ||
264 : | |||
265 : | chl | 1077 | |
266 : | int | ||
267 : | GlobalMotionEstRefine(WARPPOINTS *const startwp, | ||
268 : | MACROBLOCK * const pMBs, | ||
269 : | const MBParam * const pParam, | ||
270 : | const FRAMEINFO * const current, | ||
271 : | const FRAMEINFO * const reference, | ||
272 : | const IMAGE * const pCurr, | ||
273 : | const IMAGE * const pRef, | ||
274 : | const IMAGE * const pRefH, | ||
275 : | const IMAGE * const pRefV, | ||
276 : | const IMAGE * const pRefHV); | ||
277 : | |||
278 : | |||
279 : | |||
280 : | |||
281 : | int | ||
282 : | globalSAD(const WARPPOINTS *const wp, | ||
283 : | const MBParam * const pParam, | ||
284 : | const MACROBLOCK * const pMBs, | ||
285 : | const FRAMEINFO * const current, | ||
286 : | const IMAGE * const pRef, | ||
287 : | const IMAGE * const pCurr, | ||
288 : | uint8_t *const GMCblock); | ||
289 : | |||
290 : | edgomez | 851 | #define iDiamondSize 2 |
291 : | |||
292 : | static __inline uint32_t | ||
293 : | suxen_drol | 890 | MakeGoodMotionFlags(const uint32_t MotionFlags, const uint32_t VopFlags, const uint32_t VolFlags) |
294 : | edgomez | 851 | { |
295 : | uint32_t Flags = MotionFlags; | ||
296 : | |||
297 : | edgomez | 1107 | if (!(VopFlags & XVID_VOP_MODEDECISION_RD)) |
298 : | Flags &= ~(XVID_ME_QUARTERPELREFINE16_RD+XVID_ME_QUARTERPELREFINE8_RD+XVID_ME_HALFPELREFINE16_RD+XVID_ME_HALFPELREFINE8_RD+XVID_ME_EXTSEARCH_RD); | ||
299 : | edgomez | 851 | |
300 : | edgomez | 1107 | if (Flags & XVID_ME_EXTSEARCH_RD) |
301 : | Flags |= XVID_ME_HALFPELREFINE16_RD; | ||
302 : | edgomez | 851 | |
303 : | edgomez | 1107 | if (Flags & XVID_ME_EXTSEARCH_RD && MotionFlags & XVID_ME_EXTSEARCH8) |
304 : | Flags |= XVID_ME_HALFPELREFINE8_RD; | ||
305 : | edgomez | 851 | |
306 : | edgomez | 1107 | if (Flags & XVID_ME_HALFPELREFINE16_RD) |
307 : | Flags |= XVID_ME_QUARTERPELREFINE16_RD; | ||
308 : | edgomez | 851 | |
309 : | edgomez | 1107 | if (Flags & XVID_ME_HALFPELREFINE8_RD) { |
310 : | Flags |= XVID_ME_QUARTERPELREFINE8_RD; | ||
311 : | edgomez | 949 | Flags &= ~XVID_ME_HALFPELREFINE8; |
312 : | edgomez | 851 | } |
313 : | |||
314 : | edgomez | 1107 | if (Flags & XVID_ME_QUARTERPELREFINE8_RD) |
315 : | edgomez | 949 | Flags &= ~XVID_ME_QUARTERPELREFINE8; |
316 : | edgomez | 851 | |
317 : | edgomez | 949 | if (!(VolFlags & XVID_VOL_QUARTERPEL)) |
318 : | edgomez | 1107 | Flags &= ~(XVID_ME_QUARTERPELREFINE16+XVID_ME_QUARTERPELREFINE8+XVID_ME_QUARTERPELREFINE16_RD+XVID_ME_QUARTERPELREFINE8_RD); |
319 : | edgomez | 851 | |
320 : | edgomez | 949 | if (!(VopFlags & XVID_VOP_HALFPEL)) |
321 : | edgomez | 1107 | Flags &= ~(XVID_ME_EXTSEARCH16+XVID_ME_HALFPELREFINE16+XVID_ME_HALFPELREFINE8+XVID_ME_HALFPELREFINE16_RD+XVID_ME_HALFPELREFINE8_RD); |
322 : | edgomez | 851 | |
323 : | edgomez | 949 | if ((VopFlags & XVID_VOP_GREYSCALE) || (VopFlags & XVID_VOP_REDUCED)) |
324 : | edgomez | 1107 | Flags &= ~(XVID_ME_CHROMA_PVOP + XVID_ME_CHROMA_BVOP); |
325 : | edgomez | 851 | |
326 : | return Flags; | ||
327 : | } | ||
328 : | |||
329 : | edgomez | 1107 | /* RD mode decision and search */ |
330 : | edgomez | 851 | |
331 : | #include "../bitstream/zigzag.h" | ||
332 : | #include "../quant/quant_mpeg4.h" | ||
333 : | #include "../quant/quant_h263.h" | ||
334 : | #include "../bitstream/vlc_codes.h" | ||
335 : | edgomez | 1022 | #include "../dct/fdct.h" |
336 : | edgomez | 851 | |
337 : | static int | ||
338 : | syskin | 1108 | findRDinter(SearchData * const Data, |
339 : | const MACROBLOCK * const pMBs, const int x, const int y, | ||
340 : | const MBParam * const pParam, | ||
341 : | const uint32_t MotionFlags); | ||
342 : | edgomez | 851 | |
343 : | static int | ||
344 : | syskin | 1108 | findRDinter4v(const SearchData * const Data, |
345 : | MACROBLOCK * const pMB, const MACROBLOCK * const pMBs, | ||
346 : | const int x, const int y, | ||
347 : | const MBParam * const pParam, const uint32_t MotionFlags, | ||
348 : | const VECTOR * const backup); | ||
349 : | edgomez | 851 | |
350 : | static int | ||
351 : | syskin | 1108 | findRDintra(const SearchData * const Data); |
352 : | edgomez | 851 | |
353 : | syskin | 1084 | static int |
354 : | syskin | 1108 | findRDgmc(const SearchData * const Data, const IMAGE * const vGMC, const int x, const int y); |
355 : | syskin | 1084 | |
356 : | edgomez | 851 | int CodeCoeffIntra_CalcBits(const int16_t qcoeff[64], const uint16_t * zigzag); |
357 : | int CodeCoeffInter_CalcBits(const int16_t qcoeff[64], const uint16_t * zigzag); | ||
358 : | edgomez | 1023 | |
359 : | syskin | 1075 | /* one over lambda for R-D mode decision and motion search */ |
360 : | #define LAMBDA ( (int)(BITS_MULT/1.0) ) | ||
361 : | edgomez | 1023 | |
362 : | edgomez | 1022 | static __inline unsigned int |
363 : | edgomez | 1023 | Block_CalcBits( int16_t * const coeff, |
364 : | int16_t * const data, | ||
365 : | syskin | 1075 | int16_t * const dqcoeff, |
366 : | edgomez | 1023 | const uint32_t quant, const int quant_type, |
367 : | uint32_t * cbp, | ||
368 : | const int block) | ||
369 : | edgomez | 1022 | { |
370 : | int sum; | ||
371 : | edgomez | 1023 | int bits; |
372 : | int distortion = 0; | ||
373 : | int i; | ||
374 : | edgomez | 851 | |
375 : | edgomez | 1022 | fdct(data); |
376 : | |||
377 : | syskin | 1084 | if (quant_type) sum = quant_inter(coeff, data, quant); |
378 : | edgomez | 1022 | else sum = quant4_inter(coeff, data, quant); |
379 : | |||
380 : | if (sum > 0) { | ||
381 : | *cbp |= 1 << (5 - block); | ||
382 : | edgomez | 1023 | bits = BITS_MULT * CodeCoeffInter_CalcBits(coeff, scan_tables[0]); |
383 : | |||
384 : | syskin | 1084 | if (quant_type) dequant_inter(dqcoeff, coeff, quant); |
385 : | else dequant4_inter(dqcoeff, coeff, quant); | ||
386 : | edgomez | 1023 | |
387 : | syskin | 1084 | for (i = 0; i < 64; i++) |
388 : | distortion += (data[i] - dqcoeff[i])*(data[i] - dqcoeff[i]); | ||
389 : | |||
390 : | } else { | ||
391 : | bits = 0; | ||
392 : | for (i = 0; i < 64; i++) | ||
393 : | distortion += data[i]*data[i]; | ||
394 : | edgomez | 1023 | } |
395 : | |||
396 : | syskin | 1108 | return bits + (LAMBDA*distortion)/(quant*quant); |
397 : | edgomez | 1022 | } |
398 : | |||
399 : | edgomez | 1023 | static __inline unsigned int |
400 : | Block_CalcBitsIntra(int16_t * const coeff, | ||
401 : | int16_t * const data, | ||
402 : | syskin | 1075 | int16_t * const dqcoeff, |
403 : | edgomez | 1023 | const uint32_t quant, const int quant_type, |
404 : | uint32_t * cbp, | ||
405 : | const int block, | ||
406 : | int * dcpred) | ||
407 : | { | ||
408 : | int bits, i; | ||
409 : | int distortion = 0; | ||
410 : | syskin | 1071 | uint32_t iDcScaler = get_dc_scaler(quant, block < 4); |
411 : | edgomez | 1023 | int b_dc; |
412 : | |||
413 : | fdct(data); | ||
414 : | data[0] -= 1024; | ||
415 : | |||
416 : | syskin | 1084 | if (quant_type) quant_intra(coeff, data, quant, iDcScaler); |
417 : | syskin | 1075 | else quant4_intra(coeff, data, quant, iDcScaler); |
418 : | edgomez | 1023 | |
419 : | b_dc = coeff[0]; | ||
420 : | if (block < 4) { | ||
421 : | coeff[0] -= *dcpred; | ||
422 : | *dcpred = b_dc; | ||
423 : | } | ||
424 : | |||
425 : | bits = BITS_MULT*CodeCoeffIntra_CalcBits(coeff, scan_tables[0]); | ||
426 : | if (bits != 0) *cbp |= 1 << (5 - block); | ||
427 : | |||
428 : | syskin | 1071 | if (block < 4) bits += BITS_MULT*dcy_tab[coeff[0] + 255].len; |
429 : | else bits += BITS_MULT*dcc_tab[coeff[0] + 255].len; | ||
430 : | |||
431 : | edgomez | 1023 | coeff[0] = b_dc; |
432 : | syskin | 1084 | if (quant_type) dequant_intra(dqcoeff, coeff, quant, iDcScaler); |
433 : | syskin | 1075 | else dequant4_intra(dqcoeff, coeff, quant, iDcScaler); |
434 : | edgomez | 1023 | |
435 : | syskin | 1108 | for (i = 0; i < 64; i++) |
436 : | syskin | 1075 | distortion += (data[i] - dqcoeff[i])*(data[i] - dqcoeff[i]); |
437 : | edgomez | 1023 | |
438 : | syskin | 1108 | return bits + (LAMBDA*distortion)/(quant*quant); |
439 : | edgomez | 1023 | } |
440 : | |||
441 : | edgomez | 851 | #endif /* _MOTION_EST_H_ */ |
No admin address has been configured | ViewVC Help |
Powered by ViewVC 1.0.4 |