source: 3DVCSoftware/branches/HTM-DEV-0.3-dev1/source/Lib/TLibCommon/TComRom.cpp @ 459

Last change on this file since 459 was 459, checked in by hhi, 11 years ago

Integation of depth intra methods in macro H_3D_DIM, including:

  • DMM coding modes in H_3D_DIM_DMM.
  • RBC coding mode in H_3D_DIM_RBC.
  • Property svn:eol-style set to native
File size: 26.8 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#endif
372// ====================================================================================================================
373// Scanning order & context model mapping
374// ====================================================================================================================
375
376// scanning order table
377UInt* g_auiSigLastScan[ 3 ][ MAX_CU_DEPTH ];
378
379const UInt g_sigLastScan8x8[ 3 ][ 4 ] =
380{
381  {0, 2, 1, 3},
382  {0, 1, 2, 3},
383  {0, 2, 1, 3}
384};
385UInt g_sigLastScanCG32x32[ 64 ];
386
387const UInt g_uiMinInGroup[ 10 ] = {0,1,2,3,4,6,8,12,16,24};
388const 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};
389
390// Rice parameters for absolute transform levels
391const UInt g_auiGoRiceRange[5] =
392{
393  7, 14, 26, 46, 78
394};
395
396const UInt g_auiGoRicePrefixLen[5] =
397{
398  8, 7, 6, 5, 4
399};
400
401Void initSigLastScan(UInt* pBuffD, UInt* pBuffH, UInt* pBuffV, Int iWidth, Int iHeight)
402{
403  const UInt  uiNumScanPos  = UInt( iWidth * iWidth );
404  UInt        uiNextScanPos = 0;
405
406  if( iWidth < 16 )
407  {
408  UInt* pBuffTemp = pBuffD;
409  if( iWidth == 8 )
410  {
411    pBuffTemp = g_sigLastScanCG32x32;
412  }
413  for( UInt uiScanLine = 0; uiNextScanPos < uiNumScanPos; uiScanLine++ )
414  {
415    Int    iPrimDim  = Int( uiScanLine );
416    Int    iScndDim  = 0;
417    while( iPrimDim >= iWidth )
418    {
419      iScndDim++;
420      iPrimDim--;
421    }
422    while( iPrimDim >= 0 && iScndDim < iWidth )
423    {
424      pBuffTemp[ uiNextScanPos ] = iPrimDim * iWidth + iScndDim ;
425      uiNextScanPos++;
426      iScndDim++;
427      iPrimDim--;
428    }
429  }
430  }
431  if( iWidth > 4 )
432  {
433    UInt uiNumBlkSide = iWidth >> 2;
434    UInt uiNumBlks    = uiNumBlkSide * uiNumBlkSide;
435    UInt log2Blk      = g_aucConvertToBit[ uiNumBlkSide ] + 1;
436
437    for( UInt uiBlk = 0; uiBlk < uiNumBlks; uiBlk++ )
438    {
439      uiNextScanPos   = 0;
440      UInt initBlkPos = g_auiSigLastScan[ SCAN_DIAG ][ log2Blk ][ uiBlk ];
441      if( iWidth == 32 )
442      {
443        initBlkPos = g_sigLastScanCG32x32[ uiBlk ];
444      }
445      UInt offsetY    = initBlkPos / uiNumBlkSide;
446      UInt offsetX    = initBlkPos - offsetY * uiNumBlkSide;
447      UInt offsetD    = 4 * ( offsetX + offsetY * iWidth );
448      UInt offsetScan = 16 * uiBlk;
449      for( UInt uiScanLine = 0; uiNextScanPos < 16; uiScanLine++ )
450      {
451        Int    iPrimDim  = Int( uiScanLine );
452        Int    iScndDim  = 0;
453        while( iPrimDim >= 4 )
454        {
455          iScndDim++;
456          iPrimDim--;
457        }
458        while( iPrimDim >= 0 && iScndDim < 4 )
459        {
460          pBuffD[ uiNextScanPos + offsetScan ] = iPrimDim * iWidth + iScndDim + offsetD;
461          uiNextScanPos++;
462          iScndDim++;
463          iPrimDim--;
464        }
465      }
466    }
467  }
468 
469  UInt uiCnt = 0;
470  if( iWidth > 2 )
471  {
472    UInt numBlkSide = iWidth >> 2;
473    for(Int blkY=0; blkY < numBlkSide; blkY++)
474    {
475      for(Int blkX=0; blkX < numBlkSide; blkX++)
476      {
477        UInt offset    = blkY * 4 * iWidth + blkX * 4;
478        for(Int y=0; y < 4; y++)
479        {
480          for(Int x=0; x < 4; x++)
481          {
482            pBuffH[uiCnt] = y*iWidth + x + offset;
483            uiCnt ++;
484          }
485        }
486      }
487    }
488
489    uiCnt = 0;
490    for(Int blkX=0; blkX < numBlkSide; blkX++)
491    {
492      for(Int blkY=0; blkY < numBlkSide; blkY++)
493      {
494        UInt offset    = blkY * 4 * iWidth + blkX * 4;
495        for(Int x=0; x < 4; x++)
496        {
497          for(Int y=0; y < 4; y++)
498          {
499            pBuffV[uiCnt] = y*iWidth + x + offset;
500            uiCnt ++;
501          }
502        }
503      }
504    }
505  }
506  else
507  {
508  for(Int iY=0; iY < iHeight; iY++)
509  {
510    for(Int iX=0; iX < iWidth; iX++)
511    {
512      pBuffH[uiCnt] = iY*iWidth + iX;
513      uiCnt ++;
514    }
515  }
516
517  uiCnt = 0;
518  for(Int iX=0; iX < iWidth; iX++)
519  {
520    for(Int iY=0; iY < iHeight; iY++)
521    {
522      pBuffV[uiCnt] = iY*iWidth + iX;
523      uiCnt ++;
524    }
525  }   
526  }
527}
528
529Int g_quantTSDefault4x4[16] =
530{
531  16,16,16,16,
532  16,16,16,16,
533  16,16,16,16,
534  16,16,16,16
535};
536
537Int g_quantIntraDefault8x8[64] =
538{
539  16,16,16,16,17,18,21,24,
540  16,16,16,16,17,19,22,25,
541  16,16,17,18,20,22,25,29,
542  16,16,18,21,24,27,31,36,
543  17,17,20,24,30,35,41,47,
544  18,19,22,27,35,44,54,65,
545  21,22,25,31,41,54,70,88,
546  24,25,29,36,47,65,88,115
547};
548
549Int g_quantInterDefault8x8[64] =
550{
551  16,16,16,16,17,18,20,24,
552  16,16,16,17,18,20,24,25,
553  16,16,17,18,20,24,25,28,
554  16,17,18,20,24,25,28,33,
555  17,18,20,24,25,28,33,41,
556  18,20,24,25,28,33,41,54,
557  20,24,25,28,33,41,54,71,
558  24,25,28,33,41,54,71,91
559};
560UInt g_scalingListSize   [4] = {16,64,256,1024}; 
561UInt g_scalingListSizeX  [4] = { 4, 8, 16,  32};
562UInt g_scalingListNum[SCALING_LIST_SIZE_NUM]={6,6,6,2};
563Int  g_eTTable[4] = {0,3,1,2};
564
565#if H_3D_DIM_DMM
566std::vector< std::vector<TComWedgelet>  > g_dmmWedgeLists;
567std::vector< std::vector<TComWedgeRef>  > g_dmmWedgeRefLists;
568std::vector< std::vector<TComWedgeNode> > g_dmmWedgeNodeLists;
569std::vector< std::vector< std::vector<UInt> > > g_aauiWdgLstM3;
570
571Void initWedgeLists( Bool initRefinements )
572{
573  if( !g_dmmWedgeLists.empty() ) return;
574
575  for( UInt ui = g_aucConvertToBit[DIM_MIN_SIZE]; ui < (g_aucConvertToBit[DIM_MAX_SIZE]+1); ui++ )
576  {
577    UInt uiWedgeBlockSize = ((UInt)DIM_MIN_SIZE)<<ui;
578    std::vector<TComWedgelet> acWedgeList;
579    std::vector<TComWedgeRef> acWedgeRefList;
580    createWedgeList( uiWedgeBlockSize, uiWedgeBlockSize, acWedgeList, acWedgeRefList, g_dmmWedgeResolution[ui] );
581    g_dmmWedgeLists.push_back( acWedgeList );
582    g_dmmWedgeRefLists.push_back( acWedgeRefList );
583
584    // create WedgeNodeList
585    std::vector<TComWedgeNode> acWedgeNodeList;
586    for( UInt uiPos = 0; uiPos < acWedgeList.size(); uiPos++ )
587    {
588      if( acWedgeList[uiPos].getIsCoarse() )
589      {
590        TComWedgeNode cWedgeNode;
591        cWedgeNode.setPatternIdx( uiPos );
592
593        if( initRefinements )
594        {
595          UInt uiRefPos = 0;
596          for( Int iOffS = -1; iOffS <= 1; iOffS++ )
597          {
598            for( Int iOffE = -1; iOffE <= 1; iOffE++ )
599            {
600              if( iOffS == 0 && iOffE == 0 ) { continue; }
601
602              Int iSx = (Int)acWedgeList[uiPos].getStartX();
603              Int iSy = (Int)acWedgeList[uiPos].getStartY();
604              Int iEx = (Int)acWedgeList[uiPos].getEndX();
605              Int iEy = (Int)acWedgeList[uiPos].getEndY();
606
607              switch( acWedgeList[uiPos].getOri() )
608              {
609              case( 0 ): { iSx += iOffS; iEy += iOffE; } break;
610              case( 1 ): { iSy += iOffS; iEx -= iOffE; } break;
611              case( 2 ): { iSx -= iOffS; iEy -= iOffE; } break;
612              case( 3 ): { iSy -= iOffS; iEx += iOffE; } break;
613              case( 4 ): { iSx += iOffS; iEx += iOffE; } break;
614              case( 5 ): { iSy += iOffS; iEy += iOffE; } break;
615              default: assert( 0 );
616              }
617
618              for( UInt k = 0; k < acWedgeRefList.size(); k++ )
619              {
620                if( iSx == (Int)acWedgeRefList[k].getStartX() && 
621                    iSy == (Int)acWedgeRefList[k].getStartY() && 
622                    iEx == (Int)acWedgeRefList[k].getEndX()   && 
623                    iEy == (Int)acWedgeRefList[k].getEndY()      )
624                {
625                  if( acWedgeRefList[k].getRefIdx() != cWedgeNode.getPatternIdx() )
626                  {
627                    Bool bNew = true;
628                    for( UInt m = 0; m < uiRefPos; m++ ) { if( acWedgeRefList[k].getRefIdx() == cWedgeNode.getRefineIdx( m ) ) { bNew = false; break; } }
629
630                    if( bNew ) 
631                    {
632                      cWedgeNode.setRefineIdx( acWedgeRefList[k].getRefIdx(), uiRefPos );
633                      uiRefPos++;
634                      break;
635                    }
636                  }
637                }
638              }
639            }
640          }
641        }
642        acWedgeNodeList.push_back( cWedgeNode );
643      }
644    }
645    g_dmmWedgeNodeLists.push_back( acWedgeNodeList );
646  }
647  return;
648}
649
650Void createWedgeList( UInt uiWidth, UInt uiHeight, std::vector<TComWedgelet> &racWedgeList, std::vector<TComWedgeRef> &racWedgeRefList, WedgeResolution eWedgeRes )
651{
652  assert( uiWidth == uiHeight );
653
654  UChar    uhStartX = 0,    uhStartY = 0,    uhEndX = 0,    uhEndY = 0;
655  Int   iStepStartX = 0, iStepStartY = 0, iStepEndX = 0, iStepEndY = 0;
656
657  UInt uiBlockSize = 0;
658  switch( eWedgeRes )
659  {
660  case( DOUBLE_PEL ): { uiBlockSize = (uiWidth>>1); break; }
661  case(   FULL_PEL ): { uiBlockSize =  uiWidth;     break; }
662  case(   HALF_PEL ): { uiBlockSize = (uiWidth<<1); break; }
663  }
664
665  TComWedgelet cTempWedgelet( uiWidth, uiHeight );
666  for( UInt uiOri = 0; uiOri < 6; uiOri++ )
667  {
668    // init the edge line parameters for each of the 6 wedgelet types
669    switch( uiOri )
670    {
671    case( 0 ): {  uhStartX = 0;               uhStartY = 0;               uhEndX = 0;               uhEndY = 0;               iStepStartX = +1; iStepStartY =  0; iStepEndX =  0; iStepEndY = +1; break; }
672    case( 1 ): {  uhStartX = (uiBlockSize-1); uhStartY = 0;               uhEndX = (uiBlockSize-1); uhEndY = 0;               iStepStartX =  0; iStepStartY = +1; iStepEndX = -1; iStepEndY =  0; break; }
673    case( 2 ): {  uhStartX = (uiBlockSize-1); uhStartY = (uiBlockSize-1); uhEndX = (uiBlockSize-1); uhEndY = (uiBlockSize-1); iStepStartX = -1; iStepStartY =  0; iStepEndX =  0; iStepEndY = -1; break; }
674    case( 3 ): {  uhStartX = 0;               uhStartY = (uiBlockSize-1); uhEndX = 0;               uhEndY = (uiBlockSize-1); iStepStartX =  0; iStepStartY = -1; iStepEndX = +1; iStepEndY =  0; break; }
675    case( 4 ): {  uhStartX = 0;               uhStartY = 0;               uhEndX = 0;               uhEndY = (uiBlockSize-1); iStepStartX = +1; iStepStartY =  0; iStepEndX = +1; iStepEndY =  0; break; }
676    case( 5 ): {  uhStartX = (uiBlockSize-1); uhStartY = 0;               uhEndX = 0;               uhEndY = 0;               iStepStartX =  0; iStepStartY = +1; iStepEndX =  0; iStepEndY = +1; break; }
677    }
678
679    for( Int iK = 0; iK < uiBlockSize; iK++ )
680    {
681      for( Int iL = 0; iL < uiBlockSize; iL++ )
682      {
683        cTempWedgelet.setWedgelet( uhStartX + (iK*iStepStartX) , uhStartY + (iK*iStepStartY), uhEndX + (iL*iStepEndX), uhEndY + (iL*iStepEndY), (UChar)uiOri, eWedgeRes, ((iL%2)==0 && (iK%2)==0) );
684        addWedgeletToList( cTempWedgelet, racWedgeList, racWedgeRefList );
685      }
686    }
687  }
688
689  UInt uiThrSz = DMM3_SIMPLIFY_TR;
690  std::vector< std::vector<UInt> > auiWdgListSz;
691  for( Int idxM=2; idxM<=34 ; idxM++)
692  {
693    std::vector<UInt> auiWdgList;
694    for( Int idxW=0; idxW<racWedgeList.size(); idxW++)
695    {
696      UInt uiAbsDiff = abs(idxM-(Int)racWedgeList[idxW].getAng());
697      if( uiAbsDiff <= uiThrSz )
698      {
699        auiWdgList.push_back(idxW);
700      }
701    }
702    auiWdgListSz.push_back(auiWdgList);
703  }
704  g_aauiWdgLstM3.push_back(auiWdgListSz);
705}
706
707Void addWedgeletToList( TComWedgelet cWedgelet, std::vector<TComWedgelet> &racWedgeList, std::vector<TComWedgeRef> &racWedgeRefList )
708{
709  Bool bValid = cWedgelet.checkNotPlain();
710  if( bValid )
711  {
712    for( UInt uiPos = 0; uiPos < racWedgeList.size(); uiPos++ )
713    {
714      if( cWedgelet.checkIdentical( racWedgeList[uiPos].getPattern() ) )
715      {
716        TComWedgeRef cWedgeRef;
717        cWedgeRef.setWedgeRef( cWedgelet.getStartX(), cWedgelet.getStartY(), cWedgelet.getEndX(), cWedgelet.getEndY(), uiPos );
718        racWedgeRefList.push_back( cWedgeRef );
719        bValid = false;
720        return;
721      }
722    }
723  }
724  if( bValid )
725  {
726    for( UInt uiPos = 0; uiPos < racWedgeList.size(); uiPos++ )
727    {
728      if( cWedgelet.checkInvIdentical( racWedgeList[uiPos].getPattern() ) )
729      {
730        TComWedgeRef cWedgeRef;
731        cWedgeRef.setWedgeRef( cWedgelet.getStartX(), cWedgelet.getStartY(), cWedgelet.getEndX(), cWedgelet.getEndY(), uiPos );
732        racWedgeRefList.push_back( cWedgeRef );
733        bValid = false;
734        return;
735      }
736    }
737  }
738  if( bValid )
739  {
740    cWedgelet.findClosestAngle();
741    racWedgeList.push_back( cWedgelet );
742    TComWedgeRef cWedgeRef;
743    cWedgeRef.setWedgeRef( cWedgelet.getStartX(), cWedgelet.getStartY(), cWedgelet.getEndX(), cWedgelet.getEndY(), (UInt)(racWedgeList.size()-1) );
744    racWedgeRefList.push_back( cWedgeRef );
745  }
746}
747#endif //H_3D_DIM_DMM
748//! \}
Note: See TracBrowser for help on using the repository browser.