source: 3DVCSoftware/branches/HTM-DEV-2.0-dev3-LG/source/Lib/TLibCommon/TComRom.cpp @ 556

Last change on this file since 556 was 531, checked in by tech, 11 years ago

Merged HTM-DEV-0.3-dev1 Rev. 520.
  • Property svn:eol-style set to native
File size: 27.7 KB
Line 
1/* The copyright in this software is being made available under the BSD
2 * License, included below. This software may be subject to other third party
3 * and contributor rights, including patent rights, and no such rights are
4 * granted under this license. 
5 *
6 * Copyright (c) 2010-2013, ITU/ISO/IEC
7 * All rights reserved.
8 *
9 * Redistribution and use in source and binary forms, with or without
10 * modification, are permitted provided that the following conditions are met:
11 *
12 *  * Redistributions of source code must retain the above copyright notice,
13 *    this list of conditions and the following disclaimer.
14 *  * Redistributions in binary form must reproduce the above copyright notice,
15 *    this list of conditions and the following disclaimer in the documentation
16 *    and/or other materials provided with the distribution.
17 *  * Neither the name of the ITU/ISO/IEC nor the names of its contributors may
18 *    be used to endorse or promote products derived from this software without
19 *    specific prior written permission.
20 *
21 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
22 * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
23 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
24 * ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS
25 * BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
26 * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
27 * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
28 * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
29 * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
30 * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF
31 * THE POSSIBILITY OF SUCH DAMAGE.
32 */
33
34/** \file     TComRom.cpp
35    \brief    global variables & functions
36*/
37
38#include "TComRom.h"
39#include <memory.h>
40#include <stdlib.h>
41#include <stdio.h>
42// ====================================================================================================================
43// Initialize / destroy functions
44// ====================================================================================================================
45
46//! \ingroup TLibCommon
47//! \{
48
49// initialize ROM variables
50Void initROM()
51{
52  Int i, c;
53 
54  // g_aucConvertToBit[ x ]: log2(x/4), if x=4 -> 0, x=8 -> 1, x=16 -> 2, ...
55  ::memset( g_aucConvertToBit,   -1, sizeof( g_aucConvertToBit ) );
56  c=0;
57  for ( i=4; i<MAX_CU_SIZE; i*=2 )
58  {
59    g_aucConvertToBit[ i ] = c;
60    c++;
61  }
62  g_aucConvertToBit[ i ] = c;
63 
64  c=2;
65  for ( i=0; i<MAX_CU_DEPTH; i++ )
66  {
67    g_auiSigLastScan[0][i] = new UInt[ c*c ];
68    g_auiSigLastScan[1][i] = new UInt[ c*c ];
69    g_auiSigLastScan[2][i] = new UInt[ c*c ];
70    initSigLastScan( g_auiSigLastScan[0][i], g_auiSigLastScan[1][i], g_auiSigLastScan[2][i], c, c);
71
72    c <<= 1;
73  } 
74}
75
76Void destroyROM()
77{
78  for (Int i=0; i<MAX_CU_DEPTH; i++ )
79  {
80    delete[] g_auiSigLastScan[0][i];
81    delete[] g_auiSigLastScan[1][i];
82    delete[] g_auiSigLastScan[2][i];
83  }
84
85#if H_3D_DIM_DMM
86  if( !g_dmmWedgeLists.empty() ) 
87  {
88    for( UInt ui = 0; ui < g_dmmWedgeLists.size(); ui++ ) { g_dmmWedgeLists[ui].clear(); }
89    g_dmmWedgeLists.clear();
90  }
91  if( !g_dmmWedgeRefLists.empty() )
92  {
93    for( UInt ui = 0; ui < g_dmmWedgeRefLists.size(); ui++ ) { g_dmmWedgeRefLists[ui].clear(); }
94    g_dmmWedgeRefLists.clear();
95  }
96
97  if( !g_dmmWedgeNodeLists.empty() )
98  {
99    for( UInt ui = 0; ui < g_dmmWedgeNodeLists.size(); ui++ ) { g_dmmWedgeNodeLists[ui].clear(); }
100    g_dmmWedgeNodeLists.clear();
101  }
102#endif
103}
104
105// ====================================================================================================================
106// Data structure related table & variable
107// ====================================================================================================================
108
109UInt g_uiMaxCUWidth  = MAX_CU_SIZE;
110UInt g_uiMaxCUHeight = MAX_CU_SIZE;
111UInt g_uiMaxCUDepth  = MAX_CU_DEPTH;
112UInt g_uiAddCUDepth  = 0;
113UInt g_auiZscanToRaster [ MAX_NUM_SPU_W*MAX_NUM_SPU_W ] = { 0, };
114UInt g_auiRasterToZscan [ MAX_NUM_SPU_W*MAX_NUM_SPU_W ] = { 0, };
115UInt g_auiRasterToPelX  [ MAX_NUM_SPU_W*MAX_NUM_SPU_W ] = { 0, };
116UInt g_auiRasterToPelY  [ MAX_NUM_SPU_W*MAX_NUM_SPU_W ] = { 0, };
117
118UInt g_auiPUOffset[8] = { 0, 8, 4, 4, 2, 10, 1, 5};
119
120Void initZscanToRaster ( Int iMaxDepth, Int iDepth, UInt uiStartVal, UInt*& rpuiCurrIdx )
121{
122  Int iStride = 1 << ( iMaxDepth - 1 );
123 
124  if ( iDepth == iMaxDepth )
125  {
126    rpuiCurrIdx[0] = uiStartVal;
127    rpuiCurrIdx++;
128  }
129  else
130  {
131    Int iStep = iStride >> iDepth;
132    initZscanToRaster( iMaxDepth, iDepth+1, uiStartVal,                     rpuiCurrIdx );
133    initZscanToRaster( iMaxDepth, iDepth+1, uiStartVal+iStep,               rpuiCurrIdx );
134    initZscanToRaster( iMaxDepth, iDepth+1, uiStartVal+iStep*iStride,       rpuiCurrIdx );
135    initZscanToRaster( iMaxDepth, iDepth+1, uiStartVal+iStep*iStride+iStep, rpuiCurrIdx );
136  }
137}
138
139Void initRasterToZscan ( UInt uiMaxCUWidth, UInt uiMaxCUHeight, UInt uiMaxDepth )
140{
141  UInt  uiMinCUWidth  = uiMaxCUWidth  >> ( uiMaxDepth - 1 );
142  UInt  uiMinCUHeight = uiMaxCUHeight >> ( uiMaxDepth - 1 );
143 
144  UInt  uiNumPartInWidth  = (UInt)uiMaxCUWidth  / uiMinCUWidth;
145  UInt  uiNumPartInHeight = (UInt)uiMaxCUHeight / uiMinCUHeight;
146 
147  for ( UInt i = 0; i < uiNumPartInWidth*uiNumPartInHeight; i++ )
148  {
149    g_auiRasterToZscan[ g_auiZscanToRaster[i] ] = i;
150  }
151}
152
153Void initRasterToPelXY ( UInt uiMaxCUWidth, UInt uiMaxCUHeight, UInt uiMaxDepth )
154{
155  UInt    i;
156 
157  UInt* uiTempX = &g_auiRasterToPelX[0];
158  UInt* uiTempY = &g_auiRasterToPelY[0];
159 
160  UInt  uiMinCUWidth  = uiMaxCUWidth  >> ( uiMaxDepth - 1 );
161  UInt  uiMinCUHeight = uiMaxCUHeight >> ( uiMaxDepth - 1 );
162 
163  UInt  uiNumPartInWidth  = uiMaxCUWidth  / uiMinCUWidth;
164  UInt  uiNumPartInHeight = uiMaxCUHeight / uiMinCUHeight;
165 
166  uiTempX[0] = 0; uiTempX++;
167  for ( i = 1; i < uiNumPartInWidth; i++ )
168  {
169    uiTempX[0] = uiTempX[-1] + uiMinCUWidth; uiTempX++;
170  }
171  for ( i = 1; i < uiNumPartInHeight; i++ )
172  {
173    memcpy(uiTempX, uiTempX-uiNumPartInWidth, sizeof(UInt)*uiNumPartInWidth);
174    uiTempX += uiNumPartInWidth;
175  }
176 
177  for ( i = 1; i < uiNumPartInWidth*uiNumPartInHeight; i++ )
178  {
179    uiTempY[i] = ( i / uiNumPartInWidth ) * uiMinCUWidth;
180  }
181};
182
183
184Int g_quantScales[6] =
185{
186  26214,23302,20560,18396,16384,14564
187};   
188
189Int g_invQuantScales[6] =
190{
191  40,45,51,57,64,72
192};
193
194const Short g_aiT4[4][4] =
195{
196  { 64, 64, 64, 64},
197  { 83, 36,-36,-83},
198  { 64,-64,-64, 64},
199  { 36,-83, 83,-36}
200};
201
202const Short g_aiT8[8][8] =
203{
204  { 64, 64, 64, 64, 64, 64, 64, 64},
205  { 89, 75, 50, 18,-18,-50,-75,-89},
206  { 83, 36,-36,-83,-83,-36, 36, 83},
207  { 75,-18,-89,-50, 50, 89, 18,-75},
208  { 64,-64,-64, 64, 64,-64,-64, 64},
209  { 50,-89, 18, 75,-75,-18, 89,-50},
210  { 36,-83, 83,-36,-36, 83,-83, 36},
211  { 18,-50, 75,-89, 89,-75, 50,-18}
212};
213
214const Short g_aiT16[16][16] =
215{
216  { 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64},
217  { 90, 87, 80, 70, 57, 43, 25,  9, -9,-25,-43,-57,-70,-80,-87,-90},
218  { 89, 75, 50, 18,-18,-50,-75,-89,-89,-75,-50,-18, 18, 50, 75, 89},
219  { 87, 57,  9,-43,-80,-90,-70,-25, 25, 70, 90, 80, 43, -9,-57,-87},
220  { 83, 36,-36,-83,-83,-36, 36, 83, 83, 36,-36,-83,-83,-36, 36, 83},
221  { 80,  9,-70,-87,-25, 57, 90, 43,-43,-90,-57, 25, 87, 70, -9,-80},
222  { 75,-18,-89,-50, 50, 89, 18,-75,-75, 18, 89, 50,-50,-89,-18, 75},
223  { 70,-43,-87,  9, 90, 25,-80,-57, 57, 80,-25,-90, -9, 87, 43,-70},
224  { 64,-64,-64, 64, 64,-64,-64, 64, 64,-64,-64, 64, 64,-64,-64, 64},
225  { 57,-80,-25, 90, -9,-87, 43, 70,-70,-43, 87,  9,-90, 25, 80,-57},
226  { 50,-89, 18, 75,-75,-18, 89,-50,-50, 89,-18,-75, 75, 18,-89, 50},
227  { 43,-90, 57, 25,-87, 70,  9,-80, 80, -9,-70, 87,-25,-57, 90,-43},
228  { 36,-83, 83,-36,-36, 83,-83, 36, 36,-83, 83,-36,-36, 83,-83, 36},
229  { 25,-70, 90,-80, 43,  9,-57, 87,-87, 57, -9,-43, 80,-90, 70,-25},
230  { 18,-50, 75,-89, 89,-75, 50,-18,-18, 50,-75, 89,-89, 75,-50, 18},
231  {  9,-25, 43,-57, 70,-80, 87,-90, 90,-87, 80,-70, 57,-43, 25, -9}
232};
233
234const Short g_aiT32[32][32] =
235{
236  { 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64},
237  { 90, 90, 88, 85, 82, 78, 73, 67, 61, 54, 46, 38, 31, 22, 13,  4, -4,-13,-22,-31,-38,-46,-54,-61,-67,-73,-78,-82,-85,-88,-90,-90},
238  { 90, 87, 80, 70, 57, 43, 25,  9, -9,-25,-43,-57,-70,-80,-87,-90,-90,-87,-80,-70,-57,-43,-25, -9,  9, 25, 43, 57, 70, 80, 87, 90},
239  { 90, 82, 67, 46, 22, -4,-31,-54,-73,-85,-90,-88,-78,-61,-38,-13, 13, 38, 61, 78, 88, 90, 85, 73, 54, 31,  4,-22,-46,-67,-82,-90},
240  { 89, 75, 50, 18,-18,-50,-75,-89,-89,-75,-50,-18, 18, 50, 75, 89, 89, 75, 50, 18,-18,-50,-75,-89,-89,-75,-50,-18, 18, 50, 75, 89},
241  { 88, 67, 31,-13,-54,-82,-90,-78,-46, -4, 38, 73, 90, 85, 61, 22,-22,-61,-85,-90,-73,-38,  4, 46, 78, 90, 82, 54, 13,-31,-67,-88},
242  { 87, 57,  9,-43,-80,-90,-70,-25, 25, 70, 90, 80, 43, -9,-57,-87,-87,-57, -9, 43, 80, 90, 70, 25,-25,-70,-90,-80,-43,  9, 57, 87},
243  { 85, 46,-13,-67,-90,-73,-22, 38, 82, 88, 54, -4,-61,-90,-78,-31, 31, 78, 90, 61,  4,-54,-88,-82,-38, 22, 73, 90, 67, 13,-46,-85},
244  { 83, 36,-36,-83,-83,-36, 36, 83, 83, 36,-36,-83,-83,-36, 36, 83, 83, 36,-36,-83,-83,-36, 36, 83, 83, 36,-36,-83,-83,-36, 36, 83},
245  { 82, 22,-54,-90,-61, 13, 78, 85, 31,-46,-90,-67,  4, 73, 88, 38,-38,-88,-73, -4, 67, 90, 46,-31,-85,-78,-13, 61, 90, 54,-22,-82},
246  { 80,  9,-70,-87,-25, 57, 90, 43,-43,-90,-57, 25, 87, 70, -9,-80,-80, -9, 70, 87, 25,-57,-90,-43, 43, 90, 57,-25,-87,-70,  9, 80},
247  { 78, -4,-82,-73, 13, 85, 67,-22,-88,-61, 31, 90, 54,-38,-90,-46, 46, 90, 38,-54,-90,-31, 61, 88, 22,-67,-85,-13, 73, 82,  4,-78},
248  { 75,-18,-89,-50, 50, 89, 18,-75,-75, 18, 89, 50,-50,-89,-18, 75, 75,-18,-89,-50, 50, 89, 18,-75,-75, 18, 89, 50,-50,-89,-18, 75},
249  { 73,-31,-90,-22, 78, 67,-38,-90,-13, 82, 61,-46,-88, -4, 85, 54,-54,-85,  4, 88, 46,-61,-82, 13, 90, 38,-67,-78, 22, 90, 31,-73},
250  { 70,-43,-87,  9, 90, 25,-80,-57, 57, 80,-25,-90, -9, 87, 43,-70,-70, 43, 87, -9,-90,-25, 80, 57,-57,-80, 25, 90,  9,-87,-43, 70},
251  { 67,-54,-78, 38, 85,-22,-90,  4, 90, 13,-88,-31, 82, 46,-73,-61, 61, 73,-46,-82, 31, 88,-13,-90, -4, 90, 22,-85,-38, 78, 54,-67},
252  { 64,-64,-64, 64, 64,-64,-64, 64, 64,-64,-64, 64, 64,-64,-64, 64, 64,-64,-64, 64, 64,-64,-64, 64, 64,-64,-64, 64, 64,-64,-64, 64},
253  { 61,-73,-46, 82, 31,-88,-13, 90, -4,-90, 22, 85,-38,-78, 54, 67,-67,-54, 78, 38,-85,-22, 90,  4,-90, 13, 88,-31,-82, 46, 73,-61},
254  { 57,-80,-25, 90, -9,-87, 43, 70,-70,-43, 87,  9,-90, 25, 80,-57,-57, 80, 25,-90,  9, 87,-43,-70, 70, 43,-87, -9, 90,-25,-80, 57},
255  { 54,-85, -4, 88,-46,-61, 82, 13,-90, 38, 67,-78,-22, 90,-31,-73, 73, 31,-90, 22, 78,-67,-38, 90,-13,-82, 61, 46,-88,  4, 85,-54},
256  { 50,-89, 18, 75,-75,-18, 89,-50,-50, 89,-18,-75, 75, 18,-89, 50, 50,-89, 18, 75,-75,-18, 89,-50,-50, 89,-18,-75, 75, 18,-89, 50},
257  { 46,-90, 38, 54,-90, 31, 61,-88, 22, 67,-85, 13, 73,-82,  4, 78,-78, -4, 82,-73,-13, 85,-67,-22, 88,-61,-31, 90,-54,-38, 90,-46},
258  { 43,-90, 57, 25,-87, 70,  9,-80, 80, -9,-70, 87,-25,-57, 90,-43,-43, 90,-57,-25, 87,-70, -9, 80,-80,  9, 70,-87, 25, 57,-90, 43},
259  { 38,-88, 73, -4,-67, 90,-46,-31, 85,-78, 13, 61,-90, 54, 22,-82, 82,-22,-54, 90,-61,-13, 78,-85, 31, 46,-90, 67,  4,-73, 88,-38},
260  { 36,-83, 83,-36,-36, 83,-83, 36, 36,-83, 83,-36,-36, 83,-83, 36, 36,-83, 83,-36,-36, 83,-83, 36, 36,-83, 83,-36,-36, 83,-83, 36},
261  { 31,-78, 90,-61,  4, 54,-88, 82,-38,-22, 73,-90, 67,-13,-46, 85,-85, 46, 13,-67, 90,-73, 22, 38,-82, 88,-54, -4, 61,-90, 78,-31},
262  { 25,-70, 90,-80, 43,  9,-57, 87,-87, 57, -9,-43, 80,-90, 70,-25,-25, 70,-90, 80,-43, -9, 57,-87, 87,-57,  9, 43,-80, 90,-70, 25},
263  { 22,-61, 85,-90, 73,-38, -4, 46,-78, 90,-82, 54,-13,-31, 67,-88, 88,-67, 31, 13,-54, 82,-90, 78,-46,  4, 38,-73, 90,-85, 61,-22},
264  { 18,-50, 75,-89, 89,-75, 50,-18,-18, 50,-75, 89,-89, 75,-50, 18, 18,-50, 75,-89, 89,-75, 50,-18,-18, 50,-75, 89,-89, 75,-50, 18},
265  { 13,-38, 61,-78, 88,-90, 85,-73, 54,-31,  4, 22,-46, 67,-82, 90,-90, 82,-67, 46,-22, -4, 31,-54, 73,-85, 90,-88, 78,-61, 38,-13},
266  {  9,-25, 43,-57, 70,-80, 87,-90, 90,-87, 80,-70, 57,-43, 25, -9, -9, 25,-43, 57,-70, 80,-87, 90,-90, 87,-80, 70,-57, 43,-25,  9},
267  {  4,-13, 22,-31, 38,-46, 54,-61, 67,-73, 78,-82, 85,-88, 90,-90, 90,-90, 88,-85, 82,-78, 73,-67, 61,-54, 46,-38, 31,-22, 13, -4}
268};
269
270const UChar g_aucChromaScale[58]=
271{
272   0, 1, 2, 3, 4, 5, 6, 7, 8, 9,10,11,12,13,14,15,16,
273  17,18,19,20,21,22,23,24,25,26,27,28,29,29,30,31,32,
274  33,33,34,34,35,35,36,36,37,37,38,39,40,41,42,43,44,
275  45,46,47,48,49,50,51
276};
277
278
279// Mode-Dependent DCT/DST
280const Short g_as_DST_MAT_4 [4][4]=
281{
282  {29,   55,    74,   84},
283  {74,   74,    0 ,  -74},
284  {84,  -29,   -74,   55},
285  {55,  -84,    74,  -29},
286};
287
288
289// ====================================================================================================================
290// ADI
291// ====================================================================================================================
292
293#if FAST_UDI_USE_MPM
294const UChar g_aucIntraModeNumFast[7] =
295{
296  3,  //   2x2
297  8,  //   4x4
298  8,  //   8x8
299  3,  //  16x16   
300  3,  //  32x32   
301  3,  //  64x64   
302  3   // 128x128 
303};
304#else // FAST_UDI_USE_MPM
305const UChar g_aucIntraModeNumFast[7] =
306{
307  3,  //   2x2
308  9,  //   4x4
309  9,  //   8x8
310  4,  //  16x16   33
311  4,  //  32x32   33
312  5,  //  64x64   33
313  4   // 128x128  33
314};
315#endif // FAST_UDI_USE_MPM
316
317// chroma
318
319const UChar g_aucConvertTxtTypeToIdx[4] = { 0, 1, 1, 2 };
320
321
322// ====================================================================================================================
323// Bit-depth
324// ====================================================================================================================
325
326Int  g_bitDepthY = 8;
327Int  g_bitDepthC = 8;
328
329UInt g_uiPCMBitDepthLuma     = 8;    // PCM bit-depth
330UInt g_uiPCMBitDepthChroma   = 8;    // PCM bit-depth
331
332// ====================================================================================================================
333// Depth coding modes
334// ====================================================================================================================
335#if H_3D_DIM_DMM
336const WedgeResolution g_dmmWedgeResolution[6] = 
337{
338  HALF_PEL,    //   4x4
339  HALF_PEL,    //   8x8
340  FULL_PEL,    //  16x16
341  DOUBLE_PEL,  //  32x32
342  DOUBLE_PEL,  //  64x64
343  DOUBLE_PEL   // 128x128
344};
345
346const UChar g_dmm1TabIdxBits[6] =
347{ //2x2   4x4   8x8 16x16 32x32 64x64
348     0,    7,   10,   11,   11,   13 };
349const UChar g_dmm3IntraTabIdxBits[6] =
350{ //2x2   4x4   8x8 16x16 32x32 64x64
351     0,    6,    9,    9,    9,    0 };
352
353extern std::vector< std::vector<TComWedgelet> >   g_dmmWedgeLists;
354extern std::vector< std::vector<TComWedgeRef> >   g_dmmWedgeRefLists;
355extern std::vector< std::vector<TComWedgeNode> >  g_dmmWedgeNodeLists;
356#endif
357
358// ====================================================================================================================
359// Misc.
360// ====================================================================================================================
361
362Char  g_aucConvertToBit  [ MAX_CU_SIZE+1 ];
363
364#if ENC_DEC_TRACE
365FILE*  g_hTrace = NULL;
366const Bool g_bEncDecTraceEnable  = true;
367const Bool g_bEncDecTraceDisable = false;
368Bool   g_HLSTraceEnable = true;
369Bool   g_bJustDoIt = false;
370UInt64 g_nSymbolCounter = 0;
371#if H_MV_ENC_DEC_TRAC
372Bool g_traceCU = true; 
373Bool g_tracePU = true; 
374Bool g_traceTU = true; 
375Bool g_disableHLSTrace = false; 
376UInt64 g_stopAtCounter   = 10803; 
377#endif
378#endif
379// ====================================================================================================================
380// Scanning order & context model mapping
381// ====================================================================================================================
382
383// scanning order table
384UInt* g_auiSigLastScan[ 3 ][ MAX_CU_DEPTH ];
385
386const UInt g_sigLastScan8x8[ 3 ][ 4 ] =
387{
388  {0, 2, 1, 3},
389  {0, 1, 2, 3},
390  {0, 2, 1, 3}
391};
392UInt g_sigLastScanCG32x32[ 64 ];
393
394const UInt g_uiMinInGroup[ 10 ] = {0,1,2,3,4,6,8,12,16,24};
395const UInt g_uiGroupIdx[ 32 ]   = {0,1,2,3,4,4,5,5,6,6,6,6,7,7,7,7,8,8,8,8,8,8,8,8,9,9,9,9,9,9,9,9};
396
397// Rice parameters for absolute transform levels
398const UInt g_auiGoRiceRange[5] =
399{
400  7, 14, 26, 46, 78
401};
402
403const UInt g_auiGoRicePrefixLen[5] =
404{
405  8, 7, 6, 5, 4
406};
407
408Void initSigLastScan(UInt* pBuffD, UInt* pBuffH, UInt* pBuffV, Int iWidth, Int iHeight)
409{
410  const UInt  uiNumScanPos  = UInt( iWidth * iWidth );
411  UInt        uiNextScanPos = 0;
412
413  if( iWidth < 16 )
414  {
415  UInt* pBuffTemp = pBuffD;
416  if( iWidth == 8 )
417  {
418    pBuffTemp = g_sigLastScanCG32x32;
419  }
420  for( UInt uiScanLine = 0; uiNextScanPos < uiNumScanPos; uiScanLine++ )
421  {
422    Int    iPrimDim  = Int( uiScanLine );
423    Int    iScndDim  = 0;
424    while( iPrimDim >= iWidth )
425    {
426      iScndDim++;
427      iPrimDim--;
428    }
429    while( iPrimDim >= 0 && iScndDim < iWidth )
430    {
431      pBuffTemp[ uiNextScanPos ] = iPrimDim * iWidth + iScndDim ;
432      uiNextScanPos++;
433      iScndDim++;
434      iPrimDim--;
435    }
436  }
437  }
438  if( iWidth > 4 )
439  {
440    UInt uiNumBlkSide = iWidth >> 2;
441    UInt uiNumBlks    = uiNumBlkSide * uiNumBlkSide;
442    UInt log2Blk      = g_aucConvertToBit[ uiNumBlkSide ] + 1;
443
444    for( UInt uiBlk = 0; uiBlk < uiNumBlks; uiBlk++ )
445    {
446      uiNextScanPos   = 0;
447      UInt initBlkPos = g_auiSigLastScan[ SCAN_DIAG ][ log2Blk ][ uiBlk ];
448      if( iWidth == 32 )
449      {
450        initBlkPos = g_sigLastScanCG32x32[ uiBlk ];
451      }
452      UInt offsetY    = initBlkPos / uiNumBlkSide;
453      UInt offsetX    = initBlkPos - offsetY * uiNumBlkSide;
454      UInt offsetD    = 4 * ( offsetX + offsetY * iWidth );
455      UInt offsetScan = 16 * uiBlk;
456      for( UInt uiScanLine = 0; uiNextScanPos < 16; uiScanLine++ )
457      {
458        Int    iPrimDim  = Int( uiScanLine );
459        Int    iScndDim  = 0;
460        while( iPrimDim >= 4 )
461        {
462          iScndDim++;
463          iPrimDim--;
464        }
465        while( iPrimDim >= 0 && iScndDim < 4 )
466        {
467          pBuffD[ uiNextScanPos + offsetScan ] = iPrimDim * iWidth + iScndDim + offsetD;
468          uiNextScanPos++;
469          iScndDim++;
470          iPrimDim--;
471        }
472      }
473    }
474  }
475 
476  UInt uiCnt = 0;
477  if( iWidth > 2 )
478  {
479    UInt numBlkSide = iWidth >> 2;
480    for(Int blkY=0; blkY < numBlkSide; blkY++)
481    {
482      for(Int blkX=0; blkX < numBlkSide; blkX++)
483      {
484        UInt offset    = blkY * 4 * iWidth + blkX * 4;
485        for(Int y=0; y < 4; y++)
486        {
487          for(Int x=0; x < 4; x++)
488          {
489            pBuffH[uiCnt] = y*iWidth + x + offset;
490            uiCnt ++;
491          }
492        }
493      }
494    }
495
496    uiCnt = 0;
497    for(Int blkX=0; blkX < numBlkSide; blkX++)
498    {
499      for(Int blkY=0; blkY < numBlkSide; blkY++)
500      {
501        UInt offset    = blkY * 4 * iWidth + blkX * 4;
502        for(Int x=0; x < 4; x++)
503        {
504          for(Int y=0; y < 4; y++)
505          {
506            pBuffV[uiCnt] = y*iWidth + x + offset;
507            uiCnt ++;
508          }
509        }
510      }
511    }
512  }
513  else
514  {
515  for(Int iY=0; iY < iHeight; iY++)
516  {
517    for(Int iX=0; iX < iWidth; iX++)
518    {
519      pBuffH[uiCnt] = iY*iWidth + iX;
520      uiCnt ++;
521    }
522  }
523
524  uiCnt = 0;
525  for(Int iX=0; iX < iWidth; iX++)
526  {
527    for(Int iY=0; iY < iHeight; iY++)
528    {
529      pBuffV[uiCnt] = iY*iWidth + iX;
530      uiCnt ++;
531    }
532  }   
533  }
534}
535
536Int g_quantTSDefault4x4[16] =
537{
538  16,16,16,16,
539  16,16,16,16,
540  16,16,16,16,
541  16,16,16,16
542};
543
544Int g_quantIntraDefault8x8[64] =
545{
546  16,16,16,16,17,18,21,24,
547  16,16,16,16,17,19,22,25,
548  16,16,17,18,20,22,25,29,
549  16,16,18,21,24,27,31,36,
550  17,17,20,24,30,35,41,47,
551  18,19,22,27,35,44,54,65,
552  21,22,25,31,41,54,70,88,
553  24,25,29,36,47,65,88,115
554};
555
556Int g_quantInterDefault8x8[64] =
557{
558  16,16,16,16,17,18,20,24,
559  16,16,16,17,18,20,24,25,
560  16,16,17,18,20,24,25,28,
561  16,17,18,20,24,25,28,33,
562  17,18,20,24,25,28,33,41,
563  18,20,24,25,28,33,41,54,
564  20,24,25,28,33,41,54,71,
565  24,25,28,33,41,54,71,91
566};
567UInt g_scalingListSize   [4] = {16,64,256,1024}; 
568UInt g_scalingListSizeX  [4] = { 4, 8, 16,  32};
569UInt g_scalingListNum[SCALING_LIST_SIZE_NUM]={6,6,6,2};
570Int  g_eTTable[4] = {0,3,1,2};
571
572#if H_MV_ENC_DEC_TRAC
573#if ENC_DEC_TRACE
574Void writeToTraceFile( Char* symbolName, Int val, Bool doIt )
575{
576  if ( ( ( g_nSymbolCounter >= COUNTER_START && g_nSymbolCounter <= COUNTER_END )|| g_bJustDoIt ) && doIt  ) 
577  {
578    if ( g_stopAtCounter == g_nSymbolCounter )
579    {
580      std::cout << "Break point here." << std::endl; 
581    }
582    fprintf( g_hTrace, "%8lld  ", g_nSymbolCounter++ );
583    fprintf( g_hTrace, "%-50s       : %d\n", symbolName, val );     
584    fflush ( g_hTrace );
585    g_nSymbolCounter++; 
586  }
587}
588
589Void writeToTraceFile( Char* symbolName, Bool doIt )
590{
591  if ( ( ( g_nSymbolCounter >= COUNTER_START && g_nSymbolCounter <= COUNTER_END )|| g_bJustDoIt ) && doIt  ) 
592  {
593    fprintf( g_hTrace, "%s", symbolName );   
594    fflush ( g_hTrace );
595    g_nSymbolCounter++; 
596  }
597}
598
599#endif
600#endif
601#if H_3D_DIM_DMM
602std::vector< std::vector<TComWedgelet>  > g_dmmWedgeLists;
603std::vector< std::vector<TComWedgeRef>  > g_dmmWedgeRefLists;
604std::vector< std::vector<TComWedgeNode> > g_dmmWedgeNodeLists;
605std::vector< std::vector< std::vector<UInt> > > g_aauiWdgLstM3;
606
607Void initWedgeLists( Bool initRefinements )
608{
609  if( !g_dmmWedgeLists.empty() ) return;
610
611  for( UInt ui = g_aucConvertToBit[DIM_MIN_SIZE]; ui < (g_aucConvertToBit[DIM_MAX_SIZE]+1); ui++ )
612  {
613    UInt uiWedgeBlockSize = ((UInt)DIM_MIN_SIZE)<<ui;
614    std::vector<TComWedgelet> acWedgeList;
615    std::vector<TComWedgeRef> acWedgeRefList;
616    createWedgeList( uiWedgeBlockSize, uiWedgeBlockSize, acWedgeList, acWedgeRefList, g_dmmWedgeResolution[ui] );
617    g_dmmWedgeLists.push_back( acWedgeList );
618    g_dmmWedgeRefLists.push_back( acWedgeRefList );
619
620    // create WedgeNodeList
621    std::vector<TComWedgeNode> acWedgeNodeList;
622    for( UInt uiPos = 0; uiPos < acWedgeList.size(); uiPos++ )
623    {
624      if( acWedgeList[uiPos].getIsCoarse() )
625      {
626        TComWedgeNode cWedgeNode;
627        cWedgeNode.setPatternIdx( uiPos );
628
629        if( initRefinements )
630        {
631          UInt uiRefPos = 0;
632          for( Int iOffS = -1; iOffS <= 1; iOffS++ )
633          {
634            for( Int iOffE = -1; iOffE <= 1; iOffE++ )
635            {
636              if( iOffS == 0 && iOffE == 0 ) { continue; }
637
638              Int iSx = (Int)acWedgeList[uiPos].getStartX();
639              Int iSy = (Int)acWedgeList[uiPos].getStartY();
640              Int iEx = (Int)acWedgeList[uiPos].getEndX();
641              Int iEy = (Int)acWedgeList[uiPos].getEndY();
642
643              switch( acWedgeList[uiPos].getOri() )
644              {
645              case( 0 ): { iSx += iOffS; iEy += iOffE; } break;
646              case( 1 ): { iSy += iOffS; iEx -= iOffE; } break;
647              case( 2 ): { iSx -= iOffS; iEy -= iOffE; } break;
648              case( 3 ): { iSy -= iOffS; iEx += iOffE; } break;
649              case( 4 ): { iSx += iOffS; iEx += iOffE; } break;
650              case( 5 ): { iSy += iOffS; iEy += iOffE; } break;
651              default: assert( 0 );
652              }
653
654              for( UInt k = 0; k < acWedgeRefList.size(); k++ )
655              {
656                if( iSx == (Int)acWedgeRefList[k].getStartX() && 
657                    iSy == (Int)acWedgeRefList[k].getStartY() && 
658                    iEx == (Int)acWedgeRefList[k].getEndX()   && 
659                    iEy == (Int)acWedgeRefList[k].getEndY()      )
660                {
661                  if( acWedgeRefList[k].getRefIdx() != cWedgeNode.getPatternIdx() )
662                  {
663                    Bool bNew = true;
664                    for( UInt m = 0; m < uiRefPos; m++ ) { if( acWedgeRefList[k].getRefIdx() == cWedgeNode.getRefineIdx( m ) ) { bNew = false; break; } }
665
666                    if( bNew ) 
667                    {
668                      cWedgeNode.setRefineIdx( acWedgeRefList[k].getRefIdx(), uiRefPos );
669                      uiRefPos++;
670                      break;
671                    }
672                  }
673                }
674              }
675            }
676          }
677        }
678        acWedgeNodeList.push_back( cWedgeNode );
679      }
680    }
681    g_dmmWedgeNodeLists.push_back( acWedgeNodeList );
682  }
683  return;
684}
685
686Void createWedgeList( UInt uiWidth, UInt uiHeight, std::vector<TComWedgelet> &racWedgeList, std::vector<TComWedgeRef> &racWedgeRefList, WedgeResolution eWedgeRes )
687{
688  assert( uiWidth == uiHeight );
689
690  UChar    uhStartX = 0,    uhStartY = 0,    uhEndX = 0,    uhEndY = 0;
691  Int   iStepStartX = 0, iStepStartY = 0, iStepEndX = 0, iStepEndY = 0;
692
693  UInt uiBlockSize = 0;
694  switch( eWedgeRes )
695  {
696  case( DOUBLE_PEL ): { uiBlockSize = (uiWidth>>1); break; }
697  case(   FULL_PEL ): { uiBlockSize =  uiWidth;     break; }
698  case(   HALF_PEL ): { uiBlockSize = (uiWidth<<1); break; }
699  }
700
701  TComWedgelet cTempWedgelet( uiWidth, uiHeight );
702  for( UInt uiOri = 0; uiOri < 6; uiOri++ )
703  {
704    // init the edge line parameters for each of the 6 wedgelet types
705    switch( uiOri )
706    {
707    case( 0 ): {  uhStartX = 0;               uhStartY = 0;               uhEndX = 0;               uhEndY = 0;               iStepStartX = +1; iStepStartY =  0; iStepEndX =  0; iStepEndY = +1; break; }
708    case( 1 ): {  uhStartX = (uiBlockSize-1); uhStartY = 0;               uhEndX = (uiBlockSize-1); uhEndY = 0;               iStepStartX =  0; iStepStartY = +1; iStepEndX = -1; iStepEndY =  0; break; }
709    case( 2 ): {  uhStartX = (uiBlockSize-1); uhStartY = (uiBlockSize-1); uhEndX = (uiBlockSize-1); uhEndY = (uiBlockSize-1); iStepStartX = -1; iStepStartY =  0; iStepEndX =  0; iStepEndY = -1; break; }
710    case( 3 ): {  uhStartX = 0;               uhStartY = (uiBlockSize-1); uhEndX = 0;               uhEndY = (uiBlockSize-1); iStepStartX =  0; iStepStartY = -1; iStepEndX = +1; iStepEndY =  0; break; }
711    case( 4 ): {  uhStartX = 0;               uhStartY = 0;               uhEndX = 0;               uhEndY = (uiBlockSize-1); iStepStartX = +1; iStepStartY =  0; iStepEndX = +1; iStepEndY =  0; break; }
712    case( 5 ): {  uhStartX = (uiBlockSize-1); uhStartY = 0;               uhEndX = 0;               uhEndY = 0;               iStepStartX =  0; iStepStartY = +1; iStepEndX =  0; iStepEndY = +1; break; }
713    }
714
715    for( Int iK = 0; iK < uiBlockSize; iK++ )
716    {
717      for( Int iL = 0; iL < uiBlockSize; iL++ )
718      {
719        cTempWedgelet.setWedgelet( uhStartX + (iK*iStepStartX) , uhStartY + (iK*iStepStartY), uhEndX + (iL*iStepEndX), uhEndY + (iL*iStepEndY), (UChar)uiOri, eWedgeRes, ((iL%2)==0 && (iK%2)==0) );
720        addWedgeletToList( cTempWedgelet, racWedgeList, racWedgeRefList );
721      }
722    }
723  }
724
725  UInt uiThrSz = DMM3_SIMPLIFY_TR;
726  std::vector< std::vector<UInt> > auiWdgListSz;
727  for( Int idxM=2; idxM<=34 ; idxM++)
728  {
729    std::vector<UInt> auiWdgList;
730    for( Int idxW=0; idxW<racWedgeList.size(); idxW++)
731    {
732      UInt uiAbsDiff = abs(idxM-(Int)racWedgeList[idxW].getAng());
733      if( uiAbsDiff <= uiThrSz )
734      {
735        auiWdgList.push_back(idxW);
736      }
737    }
738    auiWdgListSz.push_back(auiWdgList);
739  }
740  g_aauiWdgLstM3.push_back(auiWdgListSz);
741}
742
743Void addWedgeletToList( TComWedgelet cWedgelet, std::vector<TComWedgelet> &racWedgeList, std::vector<TComWedgeRef> &racWedgeRefList )
744{
745  Bool bValid = cWedgelet.checkNotPlain();
746  if( bValid )
747  {
748    for( UInt uiPos = 0; uiPos < racWedgeList.size(); uiPos++ )
749    {
750      if( cWedgelet.checkIdentical( racWedgeList[uiPos].getPattern() ) )
751      {
752        TComWedgeRef cWedgeRef;
753        cWedgeRef.setWedgeRef( cWedgelet.getStartX(), cWedgelet.getStartY(), cWedgelet.getEndX(), cWedgelet.getEndY(), uiPos );
754        racWedgeRefList.push_back( cWedgeRef );
755        bValid = false;
756        return;
757      }
758    }
759  }
760  if( bValid )
761  {
762    for( UInt uiPos = 0; uiPos < racWedgeList.size(); uiPos++ )
763    {
764      if( cWedgelet.checkInvIdentical( racWedgeList[uiPos].getPattern() ) )
765      {
766        TComWedgeRef cWedgeRef;
767        cWedgeRef.setWedgeRef( cWedgelet.getStartX(), cWedgelet.getStartY(), cWedgelet.getEndX(), cWedgelet.getEndY(), uiPos );
768        racWedgeRefList.push_back( cWedgeRef );
769        bValid = false;
770        return;
771      }
772    }
773  }
774  if( bValid )
775  {
776    cWedgelet.findClosestAngle();
777    racWedgeList.push_back( cWedgelet );
778    TComWedgeRef cWedgeRef;
779    cWedgeRef.setWedgeRef( cWedgelet.getStartX(), cWedgelet.getStartY(), cWedgelet.getEndX(), cWedgelet.getEndY(), (UInt)(racWedgeList.size()-1) );
780    racWedgeRefList.push_back( cWedgeRef );
781  }
782}
783#endif //H_3D_DIM_DMM
784
785//! \}
Note: See TracBrowser for help on using the repository browser.