source: SHVCSoftware/branches/SHM-4.0-dev/source/Lib/TLibEncoder/TEncGOP.cpp @ 492

Last change on this file since 492 was 489, checked in by seregin, 12 years ago

fix setting layerId and refLayerId for upsampling
  • Property svn:eol-style set to native
File size: 140.2 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     TEncGOP.cpp
35    \brief    GOP encoder class
36*/
37
38#include <list>
39#include <algorithm>
40#include <functional>
41
42#include "TEncTop.h"
43#include "TEncGOP.h"
44#include "TEncAnalyze.h"
45#include "libmd5/MD5.h"
46#include "TLibCommon/SEI.h"
47#include "TLibCommon/NAL.h"
48#include "NALwrite.h"
49#include <time.h>
50#include <math.h>
51
52using namespace std;
53//! \ingroup TLibEncoder
54//! \{
55
56// ====================================================================================================================
57// Constructor / destructor / initialization / destroy
58// ====================================================================================================================
59Int getLSB(Int poc, Int maxLSB)
60{
61  if (poc >= 0)
62  {
63    return poc % maxLSB;
64  }
65  else
66  {
67    return (maxLSB - ((-poc) % maxLSB)) % maxLSB;
68  }
69}
70
71TEncGOP::TEncGOP()
72{
73  m_iLastIDR            = 0;
74  m_iGopSize            = 0;
75  m_iNumPicCoded        = 0; //Niko
76  m_bFirst              = true;
77 
78  m_pcCfg               = NULL;
79  m_pcSliceEncoder      = NULL;
80  m_pcListPic           = NULL;
81 
82  m_pcEntropyCoder      = NULL;
83  m_pcCavlcCoder        = NULL;
84  m_pcSbacCoder         = NULL;
85  m_pcBinCABAC          = NULL;
86 
87  m_bSeqFirst           = true;
88 
89  m_bRefreshPending     = 0;
90  m_pocCRA            = 0;
91  m_numLongTermRefPicSPS = 0;
92  ::memset(m_ltRefPicPocLsbSps, 0, sizeof(m_ltRefPicPocLsbSps));
93  ::memset(m_ltRefPicUsedByCurrPicFlag, 0, sizeof(m_ltRefPicUsedByCurrPicFlag));
94  m_cpbRemovalDelay   = 0;
95  m_lastBPSEI         = 0;
96  xResetNonNestedSEIPresentFlags();
97  xResetNestedSEIPresentFlags();
98#if SVC_UPSAMPLING
99  m_pcPredSearch        = NULL;
100#endif
101  return;
102}
103
104TEncGOP::~TEncGOP()
105{
106}
107
108/** Create list to contain pointers to LCU start addresses of slice.
109 */
110#if SVC_EXTENSION
111Void  TEncGOP::create( UInt layerId )
112{
113  m_bLongtermTestPictureHasBeenCoded = 0;
114  m_bLongtermTestPictureHasBeenCoded2 = 0;
115  m_layerId = layerId;
116}
117#else
118Void  TEncGOP::create()
119{
120  m_bLongtermTestPictureHasBeenCoded = 0;
121  m_bLongtermTestPictureHasBeenCoded2 = 0;
122}
123#endif
124
125Void  TEncGOP::destroy()
126{
127}
128
129Void TEncGOP::init ( TEncTop* pcTEncTop )
130{
131  m_pcEncTop     = pcTEncTop;
132  m_pcCfg                = pcTEncTop;
133  m_pcSliceEncoder       = pcTEncTop->getSliceEncoder();
134  m_pcListPic            = pcTEncTop->getListPic(); 
135 
136  m_pcEntropyCoder       = pcTEncTop->getEntropyCoder();
137  m_pcCavlcCoder         = pcTEncTop->getCavlcCoder();
138  m_pcSbacCoder          = pcTEncTop->getSbacCoder();
139  m_pcBinCABAC           = pcTEncTop->getBinCABAC();
140  m_pcLoopFilter         = pcTEncTop->getLoopFilter();
141  m_pcBitCounter         = pcTEncTop->getBitCounter();
142 
143  //--Adaptive Loop filter
144  m_pcSAO                = pcTEncTop->getSAO();
145  m_pcRateCtrl           = pcTEncTop->getRateCtrl();
146  m_lastBPSEI          = 0;
147  m_totalCoded         = 0;
148
149#if SVC_EXTENSION
150  m_ppcTEncTop           = pcTEncTop->getLayerEnc();
151#endif
152#if SVC_UPSAMPLING
153  m_pcPredSearch         = pcTEncTop->getPredSearch();                       ///< encoder search class
154#endif
155}
156
157SEIActiveParameterSets* TEncGOP::xCreateSEIActiveParameterSets (TComSPS *sps)
158{
159  SEIActiveParameterSets *seiActiveParameterSets = new SEIActiveParameterSets(); 
160  seiActiveParameterSets->activeVPSId = m_pcCfg->getVPS()->getVPSId(); 
161  seiActiveParameterSets->m_fullRandomAccessFlag = false;
162  seiActiveParameterSets->m_noParamSetUpdateFlag = false;
163  seiActiveParameterSets->numSpsIdsMinus1 = 0;
164  seiActiveParameterSets->activeSeqParamSetId.resize(seiActiveParameterSets->numSpsIdsMinus1 + 1); 
165  seiActiveParameterSets->activeSeqParamSetId[0] = sps->getSPSId();
166  return seiActiveParameterSets;
167}
168
169#if M0043_LAYERS_PRESENT_SEI
170SEILayersPresent* TEncGOP::xCreateSEILayersPresent ()
171{
172  UInt i = 0;
173  SEILayersPresent *seiLayersPresent = new SEILayersPresent(); 
174  seiLayersPresent->m_activeVpsId = m_pcCfg->getVPS()->getVPSId(); 
175  seiLayersPresent->m_vpsMaxLayers = m_pcCfg->getVPS()->getMaxLayers();
176  for ( ; i < seiLayersPresent->m_vpsMaxLayers; i++)
177  {
178    seiLayersPresent->m_layerPresentFlag[i] = true; 
179  }
180  for ( ; i < MAX_LAYERS; i++)
181  {
182    seiLayersPresent->m_layerPresentFlag[i] = false; 
183  }
184  return seiLayersPresent;
185}
186#endif
187
188SEIFramePacking* TEncGOP::xCreateSEIFramePacking()
189{
190  SEIFramePacking *seiFramePacking = new SEIFramePacking();
191  seiFramePacking->m_arrangementId = m_pcCfg->getFramePackingArrangementSEIId();
192  seiFramePacking->m_arrangementCancelFlag = 0;
193  seiFramePacking->m_arrangementType = m_pcCfg->getFramePackingArrangementSEIType();
194  assert((seiFramePacking->m_arrangementType > 2) && (seiFramePacking->m_arrangementType < 6) );
195  seiFramePacking->m_quincunxSamplingFlag = m_pcCfg->getFramePackingArrangementSEIQuincunx();
196  seiFramePacking->m_contentInterpretationType = m_pcCfg->getFramePackingArrangementSEIInterpretation();
197  seiFramePacking->m_spatialFlippingFlag = 0;
198  seiFramePacking->m_frame0FlippedFlag = 0;
199  seiFramePacking->m_fieldViewsFlag = (seiFramePacking->m_arrangementType == 2);
200  seiFramePacking->m_currentFrameIsFrame0Flag = ((seiFramePacking->m_arrangementType == 5) && m_iNumPicCoded&1);
201  seiFramePacking->m_frame0SelfContainedFlag = 0;
202  seiFramePacking->m_frame1SelfContainedFlag = 0;
203  seiFramePacking->m_frame0GridPositionX = 0;
204  seiFramePacking->m_frame0GridPositionY = 0;
205  seiFramePacking->m_frame1GridPositionX = 0;
206  seiFramePacking->m_frame1GridPositionY = 0;
207  seiFramePacking->m_arrangementReservedByte = 0;
208  seiFramePacking->m_arrangementPersistenceFlag = true;
209  seiFramePacking->m_upsampledAspectRatio = 0;
210  return seiFramePacking;
211}
212
213SEIDisplayOrientation* TEncGOP::xCreateSEIDisplayOrientation()
214{
215  SEIDisplayOrientation *seiDisplayOrientation = new SEIDisplayOrientation();
216  seiDisplayOrientation->cancelFlag = false;
217  seiDisplayOrientation->horFlip = false;
218  seiDisplayOrientation->verFlip = false;
219  seiDisplayOrientation->anticlockwiseRotation = m_pcCfg->getDisplayOrientationSEIAngle();
220  return seiDisplayOrientation;
221}
222
223SEIToneMappingInfo*  TEncGOP::xCreateSEIToneMappingInfo()
224{
225  SEIToneMappingInfo *seiToneMappingInfo = new SEIToneMappingInfo();
226  seiToneMappingInfo->m_toneMapId = m_pcCfg->getTMISEIToneMapId();
227  seiToneMappingInfo->m_toneMapCancelFlag = m_pcCfg->getTMISEIToneMapCancelFlag();
228  seiToneMappingInfo->m_toneMapPersistenceFlag = m_pcCfg->getTMISEIToneMapPersistenceFlag();
229
230  seiToneMappingInfo->m_codedDataBitDepth = m_pcCfg->getTMISEICodedDataBitDepth();
231  assert(seiToneMappingInfo->m_codedDataBitDepth >= 8 && seiToneMappingInfo->m_codedDataBitDepth <= 14);
232  seiToneMappingInfo->m_targetBitDepth = m_pcCfg->getTMISEITargetBitDepth();
233  assert( seiToneMappingInfo->m_targetBitDepth >= 1 && seiToneMappingInfo->m_targetBitDepth <= 17 );
234  seiToneMappingInfo->m_modelId = m_pcCfg->getTMISEIModelID();
235  assert(seiToneMappingInfo->m_modelId >=0 &&seiToneMappingInfo->m_modelId<=4);
236
237  switch( seiToneMappingInfo->m_modelId)
238  {
239  case 0:
240    {
241      seiToneMappingInfo->m_minValue = m_pcCfg->getTMISEIMinValue();
242      seiToneMappingInfo->m_maxValue = m_pcCfg->getTMISEIMaxValue();
243      break;
244    }
245  case 1:
246    {
247      seiToneMappingInfo->m_sigmoidMidpoint = m_pcCfg->getTMISEISigmoidMidpoint();
248      seiToneMappingInfo->m_sigmoidWidth = m_pcCfg->getTMISEISigmoidWidth();
249      break;
250    }
251  case 2:
252    {
253      UInt num = 1u<<(seiToneMappingInfo->m_targetBitDepth);
254      seiToneMappingInfo->m_startOfCodedInterval.resize(num);
255      Int* ptmp = m_pcCfg->getTMISEIStartOfCodedInterva();
256      if(ptmp)
257      {
258        for(int i=0; i<num;i++)
259        {
260          seiToneMappingInfo->m_startOfCodedInterval[i] = ptmp[i];
261        }
262      }
263      break;
264    }
265  case 3:
266    {
267      seiToneMappingInfo->m_numPivots = m_pcCfg->getTMISEINumPivots();
268      seiToneMappingInfo->m_codedPivotValue.resize(seiToneMappingInfo->m_numPivots);
269      seiToneMappingInfo->m_targetPivotValue.resize(seiToneMappingInfo->m_numPivots);
270      Int* ptmpcoded = m_pcCfg->getTMISEICodedPivotValue();
271      Int* ptmptarget = m_pcCfg->getTMISEITargetPivotValue();
272      if(ptmpcoded&&ptmptarget)
273      {
274        for(int i=0; i<(seiToneMappingInfo->m_numPivots);i++)
275        {
276          seiToneMappingInfo->m_codedPivotValue[i]=ptmpcoded[i];
277          seiToneMappingInfo->m_targetPivotValue[i]=ptmptarget[i];
278         }
279       }
280       break;
281     }
282  case 4:
283     {
284       seiToneMappingInfo->m_cameraIsoSpeedIdc = m_pcCfg->getTMISEICameraIsoSpeedIdc();
285       seiToneMappingInfo->m_cameraIsoSpeedValue = m_pcCfg->getTMISEICameraIsoSpeedValue();
286       assert( seiToneMappingInfo->m_cameraIsoSpeedValue !=0 );
287       seiToneMappingInfo->m_exposureCompensationValueSignFlag = m_pcCfg->getTMISEIExposureCompensationValueSignFlag();
288       seiToneMappingInfo->m_exposureCompensationValueNumerator = m_pcCfg->getTMISEIExposureCompensationValueNumerator();
289       seiToneMappingInfo->m_exposureCompensationValueDenomIdc = m_pcCfg->getTMISEIExposureCompensationValueDenomIdc();
290       seiToneMappingInfo->m_refScreenLuminanceWhite = m_pcCfg->getTMISEIRefScreenLuminanceWhite();
291       seiToneMappingInfo->m_extendedRangeWhiteLevel = m_pcCfg->getTMISEIExtendedRangeWhiteLevel();
292       assert( seiToneMappingInfo->m_extendedRangeWhiteLevel >= 100 );
293       seiToneMappingInfo->m_nominalBlackLevelLumaCodeValue = m_pcCfg->getTMISEINominalBlackLevelLumaCodeValue();
294       seiToneMappingInfo->m_nominalWhiteLevelLumaCodeValue = m_pcCfg->getTMISEINominalWhiteLevelLumaCodeValue();
295       assert( seiToneMappingInfo->m_nominalWhiteLevelLumaCodeValue > seiToneMappingInfo->m_nominalBlackLevelLumaCodeValue );
296       seiToneMappingInfo->m_extendedWhiteLevelLumaCodeValue = m_pcCfg->getTMISEIExtendedWhiteLevelLumaCodeValue();
297       assert( seiToneMappingInfo->m_extendedWhiteLevelLumaCodeValue >= seiToneMappingInfo->m_nominalWhiteLevelLumaCodeValue );
298       break;
299    }
300  default:
301    {
302      assert(!"Undefined SEIToneMapModelId");
303      break;
304    }
305  }
306  return seiToneMappingInfo;
307}
308
309#if N0383_IL_CONSTRAINED_TILE_SETS_SEI
310SEIInterLayerConstrainedTileSets* TEncGOP::xCreateSEIInterLayerConstrainedTileSets()
311{
312  SEIInterLayerConstrainedTileSets *seiInterLayerConstrainedTileSets = new SEIInterLayerConstrainedTileSets();
313  seiInterLayerConstrainedTileSets->m_ilAllTilesExactSampleValueMatchFlag = false;
314  seiInterLayerConstrainedTileSets->m_ilOneTilePerTileSetFlag = false;
315  if (!seiInterLayerConstrainedTileSets->m_ilOneTilePerTileSetFlag)
316  {
317    seiInterLayerConstrainedTileSets->m_ilNumSetsInMessageMinus1 = m_pcCfg->getIlNumSetsInMessage() - 1;
318    if (seiInterLayerConstrainedTileSets->m_ilNumSetsInMessageMinus1)
319    {
320      seiInterLayerConstrainedTileSets->m_skippedTileSetPresentFlag = m_pcCfg->getSkippedTileSetPresentFlag();
321    }
322    else
323    {
324      seiInterLayerConstrainedTileSets->m_skippedTileSetPresentFlag = false;
325    }
326    seiInterLayerConstrainedTileSets->m_ilNumSetsInMessageMinus1 += seiInterLayerConstrainedTileSets->m_skippedTileSetPresentFlag ? 1 : 0;
327    for (UInt i = 0; i < m_pcCfg->getIlNumSetsInMessage(); i++)
328    {
329      seiInterLayerConstrainedTileSets->m_ilctsId[i] = i;
330      seiInterLayerConstrainedTileSets->m_ilNumTileRectsInSetMinus1[i] = 0;
331      for( UInt j = 0; j <= seiInterLayerConstrainedTileSets->m_ilNumTileRectsInSetMinus1[i]; j++)
332      {
333        seiInterLayerConstrainedTileSets->m_ilTopLeftTileIndex[i][j]     = m_pcCfg->getTopLeftTileIndex(i);
334        seiInterLayerConstrainedTileSets->m_ilBottomRightTileIndex[i][j] = m_pcCfg->getBottomRightTileIndex(i);
335      }
336      seiInterLayerConstrainedTileSets->m_ilcIdc[i] = m_pcCfg->getIlcIdc(i);
337      if (seiInterLayerConstrainedTileSets->m_ilAllTilesExactSampleValueMatchFlag)
338      {
339        seiInterLayerConstrainedTileSets->m_ilExactSampleValueMatchFlag[i] = false;
340      }
341    }
342  }
343
344  return seiInterLayerConstrainedTileSets;
345}
346#endif
347
348Void TEncGOP::xCreateLeadingSEIMessages (/*SEIMessages seiMessages,*/ AccessUnit &accessUnit, TComSPS *sps)
349{
350  OutputNALUnit nalu(NAL_UNIT_PREFIX_SEI);
351
352  if(m_pcCfg->getActiveParameterSetsSEIEnabled())
353  {
354    SEIActiveParameterSets *sei = xCreateSEIActiveParameterSets (sps);
355
356    //nalu = NALUnit(NAL_UNIT_SEI);
357    m_pcEntropyCoder->setBitstream(&nalu.m_Bitstream);
358    m_seiWriter.writeSEImessage(nalu.m_Bitstream, *sei, sps); 
359    writeRBSPTrailingBits(nalu.m_Bitstream);
360    accessUnit.push_back(new NALUnitEBSP(nalu));
361    delete sei;
362    m_activeParameterSetSEIPresentInAU = true;
363  }
364
365#if M0043_LAYERS_PRESENT_SEI
366  if(m_pcCfg->getLayersPresentSEIEnabled())
367  {
368    SEILayersPresent *sei = xCreateSEILayersPresent ();
369    m_pcEntropyCoder->setBitstream(&nalu.m_Bitstream);
370    m_seiWriter.writeSEImessage(nalu.m_Bitstream, *sei, sps); 
371    writeRBSPTrailingBits(nalu.m_Bitstream);
372    accessUnit.push_back(new NALUnitEBSP(nalu));
373    delete sei;
374  }
375#endif
376
377  if(m_pcCfg->getFramePackingArrangementSEIEnabled())
378  {
379    SEIFramePacking *sei = xCreateSEIFramePacking ();
380
381    nalu = NALUnit(NAL_UNIT_PREFIX_SEI);
382    m_pcEntropyCoder->setBitstream(&nalu.m_Bitstream);
383    m_seiWriter.writeSEImessage(nalu.m_Bitstream, *sei, sps);
384    writeRBSPTrailingBits(nalu.m_Bitstream);
385    accessUnit.push_back(new NALUnitEBSP(nalu));
386    delete sei;
387  }
388  if (m_pcCfg->getDisplayOrientationSEIAngle())
389  {
390    SEIDisplayOrientation *sei = xCreateSEIDisplayOrientation();
391
392    nalu = NALUnit(NAL_UNIT_PREFIX_SEI); 
393    m_pcEntropyCoder->setBitstream(&nalu.m_Bitstream);
394    m_seiWriter.writeSEImessage(nalu.m_Bitstream, *sei, sps); 
395    writeRBSPTrailingBits(nalu.m_Bitstream);
396    accessUnit.push_back(new NALUnitEBSP(nalu));
397    delete sei;
398  }
399  if(m_pcCfg->getToneMappingInfoSEIEnabled())
400  {
401    SEIToneMappingInfo *sei = xCreateSEIToneMappingInfo ();
402     
403    nalu = NALUnit(NAL_UNIT_PREFIX_SEI); 
404    m_pcEntropyCoder->setBitstream(&nalu.m_Bitstream);
405    m_seiWriter.writeSEImessage(nalu.m_Bitstream, *sei, sps); 
406    writeRBSPTrailingBits(nalu.m_Bitstream);
407    accessUnit.push_back(new NALUnitEBSP(nalu));
408    delete sei;
409  }
410#if N0383_IL_CONSTRAINED_TILE_SETS_SEI
411  if(m_pcCfg->getInterLayerConstrainedTileSetsSEIEnabled())
412  {
413    SEIInterLayerConstrainedTileSets *sei = xCreateSEIInterLayerConstrainedTileSets ();
414
415    nalu = NALUnit(NAL_UNIT_PREFIX_SEI, 0, m_pcCfg->getNumLayer()-1); // For highest layer
416    m_pcEntropyCoder->setBitstream(&nalu.m_Bitstream);
417    m_seiWriter.writeSEImessage(nalu.m_Bitstream, *sei, sps); 
418    writeRBSPTrailingBits(nalu.m_Bitstream);
419    accessUnit.push_back(new NALUnitEBSP(nalu));
420    delete sei;
421  }
422#endif
423}
424
425// ====================================================================================================================
426// Public member functions
427// ====================================================================================================================
428#if SVC_EXTENSION
429Void TEncGOP::compressGOP( Int iPicIdInGOP, Int iPOCLast, Int iNumPicRcvd, TComList<TComPic*>& rcListPic, TComList<TComPicYuv*>& rcListPicYuvRecOut, std::list<AccessUnit>& accessUnitsInGOP, Bool isField, Bool isTff)
430#else
431Void TEncGOP::compressGOP( Int iPOCLast, Int iNumPicRcvd, TComList<TComPic*>& rcListPic, TComList<TComPicYuv*>& rcListPicYuvRecOut, std::list<AccessUnit>& accessUnitsInGOP, Bool isField, Bool isTff)
432#endif
433{
434  TComPic*        pcPic;
435  TComPicYuv*     pcPicYuvRecOut;
436  TComSlice*      pcSlice;
437  TComOutputBitstream  *pcBitstreamRedirect;
438  pcBitstreamRedirect = new TComOutputBitstream;
439  AccessUnit::iterator  itLocationToPushSliceHeaderNALU; // used to store location where NALU containing slice header is to be inserted
440  UInt                  uiOneBitstreamPerSliceLength = 0;
441  TEncSbac* pcSbacCoders = NULL;
442  TComOutputBitstream* pcSubstreamsOut = NULL;
443
444  xInitGOP( iPOCLast, iNumPicRcvd, rcListPic, rcListPicYuvRecOut, isField );
445
446  m_iNumPicCoded = 0;
447  SEIPictureTiming pictureTimingSEI;
448  Bool writeSOP = m_pcCfg->getSOPDescriptionSEIEnabled();
449  // Initialize Scalable Nesting SEI with single layer values
450  SEIScalableNesting scalableNestingSEI;
451  scalableNestingSEI.m_bitStreamSubsetFlag           = 1;      // If the nested SEI messages are picture buffereing SEI mesages, picure timing SEI messages or sub-picture timing SEI messages, bitstream_subset_flag shall be equal to 1
452  scalableNestingSEI.m_nestingOpFlag                 = 0;
453  scalableNestingSEI.m_nestingNumOpsMinus1           = 0;      //nesting_num_ops_minus1
454  scalableNestingSEI.m_allLayersFlag                 = 0;
455  scalableNestingSEI.m_nestingNoOpMaxTemporalIdPlus1 = 6 + 1;  //nesting_no_op_max_temporal_id_plus1
456  scalableNestingSEI.m_nestingNumLayersMinus1        = 1 - 1;  //nesting_num_layers_minus1
457  scalableNestingSEI.m_nestingLayerId[0]             = 0;
458  scalableNestingSEI.m_callerOwnsSEIs                = true;
459  Int picSptDpbOutputDuDelay = 0;
460  UInt *accumBitsDU = NULL;
461  UInt *accumNalsDU = NULL;
462  SEIDecodingUnitInfo decodingUnitInfoSEI;
463#if SVC_EXTENSION
464  for ( Int iGOPid=iPicIdInGOP; iGOPid < iPicIdInGOP+1; iGOPid++ )
465#else
466  for ( Int iGOPid=0; iGOPid < m_iGopSize; iGOPid++ )
467#endif
468  {
469    UInt uiColDir = 1;
470    //-- For time output for each slice
471    long iBeforeTime = clock();
472
473    //select uiColDir
474    Int iCloseLeft=1, iCloseRight=-1;
475    for(Int i = 0; i<m_pcCfg->getGOPEntry(iGOPid).m_numRefPics; i++) 
476    {
477      Int iRef = m_pcCfg->getGOPEntry(iGOPid).m_referencePics[i];
478      if(iRef>0&&(iRef<iCloseRight||iCloseRight==-1))
479      {
480        iCloseRight=iRef;
481      }
482      else if(iRef<0&&(iRef>iCloseLeft||iCloseLeft==1))
483      {
484        iCloseLeft=iRef;
485      }
486    }
487    if(iCloseRight>-1)
488    {
489      iCloseRight=iCloseRight+m_pcCfg->getGOPEntry(iGOPid).m_POC-1;
490    }
491    if(iCloseLeft<1) 
492    {
493      iCloseLeft=iCloseLeft+m_pcCfg->getGOPEntry(iGOPid).m_POC-1;
494      while(iCloseLeft<0)
495      {
496        iCloseLeft+=m_iGopSize;
497      }
498    }
499    Int iLeftQP=0, iRightQP=0;
500    for(Int i=0; i<m_iGopSize; i++)
501    {
502      if(m_pcCfg->getGOPEntry(i).m_POC==(iCloseLeft%m_iGopSize)+1)
503      {
504        iLeftQP= m_pcCfg->getGOPEntry(i).m_QPOffset;
505      }
506      if (m_pcCfg->getGOPEntry(i).m_POC==(iCloseRight%m_iGopSize)+1)
507      {
508        iRightQP=m_pcCfg->getGOPEntry(i).m_QPOffset;
509      }
510    }
511    if(iCloseRight>-1&&iRightQP<iLeftQP)
512    {
513      uiColDir=0;
514    }
515
516    /////////////////////////////////////////////////////////////////////////////////////////////////// Initial to start encoding
517    Int iTimeOffset;
518    Int pocCurr;
519   
520    if(iPOCLast == 0) //case first frame or first top field
521    {
522      pocCurr=0;
523      iTimeOffset = 1;
524    }
525    else if(iPOCLast == 1 && isField) //case first bottom field, just like the first frame, the poc computation is not right anymore, we set the right value
526    {
527      pocCurr = 1;
528      iTimeOffset = 1;
529    }
530    else
531    {
532      pocCurr = iPOCLast - iNumPicRcvd + m_pcCfg->getGOPEntry(iGOPid).m_POC - isField;
533      iTimeOffset = m_pcCfg->getGOPEntry(iGOPid).m_POC;
534    }
535
536    if(pocCurr>=m_pcCfg->getFramesToBeEncoded())
537    {
538      continue;
539    }
540
541#if M0040_ADAPTIVE_RESOLUTION_CHANGE
542    if (m_pcEncTop->getAdaptiveResolutionChange() > 0 && ((m_layerId == 1 && pocCurr < m_pcEncTop->getAdaptiveResolutionChange()) ||
543                                                          (m_layerId == 0 && pocCurr > m_pcEncTop->getAdaptiveResolutionChange())) )
544    {
545      continue;
546    }
547#endif
548
549    if( getNalUnitType(pocCurr, m_iLastIDR) == NAL_UNIT_CODED_SLICE_IDR_W_RADL || getNalUnitType(pocCurr, m_iLastIDR) == NAL_UNIT_CODED_SLICE_IDR_N_LP )
550    {
551      m_iLastIDR = pocCurr;
552    }       
553    // start a new access unit: create an entry in the list of output access units
554    accessUnitsInGOP.push_back(AccessUnit());
555    AccessUnit& accessUnit = accessUnitsInGOP.back();
556    xGetBuffer( rcListPic, rcListPicYuvRecOut, iNumPicRcvd, iTimeOffset, pcPic, pcPicYuvRecOut, pocCurr, isField );
557
558    //  Slice data initialization
559    pcPic->clearSliceBuffer();
560    assert(pcPic->getNumAllocatedSlice() == 1);
561    m_pcSliceEncoder->setSliceIdx(0);
562    pcPic->setCurrSliceIdx(0);
563#if SVC_EXTENSION
564    pcPic->setLayerId( m_layerId );
565    m_pcSliceEncoder->initEncSlice ( pcPic, iPOCLast, pocCurr, iNumPicRcvd, iGOPid, pcSlice, m_pcEncTop->getSPS(), m_pcEncTop->getPPS(), m_pcEncTop->getVPS(), isField );
566#else
567    m_pcSliceEncoder->initEncSlice ( pcPic, iPOCLast, pocCurr, iNumPicRcvd, iGOPid, pcSlice, m_pcEncTop->getSPS(), m_pcEncTop->getPPS(), isField );
568#endif
569
570    //Set Frame/Field coding
571    pcSlice->getPic()->setField(isField);
572
573#if POC_RESET_FLAG
574    if( !pcSlice->getPocResetFlag() ) // For picture that are not reset, we should adjust the value of POC calculated from the configuration files.
575    {
576      // Subtract POC adjustment value until now.
577      pcSlice->setPOC( pcSlice->getPOC() - m_pcEncTop->getPocAdjustmentValue() );
578    }
579    else
580    {
581      // Check if this is the first slice in the picture
582      // In the encoder, the POC values are copied along with copySliceInfo, so we only need
583      // to do this for the first slice.
584      Int pocAdjustValue = pcSlice->getPOC() - m_pcEncTop->getPocAdjustmentValue();
585      if( pcSlice->getSliceIdx() == 0 )
586      {
587        TComList<TComPic*>::iterator  iterPic = rcListPic.begin(); 
588
589        // Iterate through all picture in DPB
590        while( iterPic != rcListPic.end() )
591        {             
592          TComPic *dpbPic = *iterPic;
593          if( dpbPic->getPOC() == pocCurr )
594          {
595            if( dpbPic->getReconMark() )
596            {
597              assert( !( dpbPic->getSlice(0)->isReferenced() ) && !( dpbPic->getOutputMark() ) );
598            }
599          }
600          // Check if the picture pointed to by iterPic is either used for reference or
601          // needed for output, are in the same layer, and not the current picture.
602          if( /* ( ( dpbPic->getSlice(0)->isReferenced() ) || ( dpbPic->getOutputMark() ) )
603              && */ ( dpbPic->getLayerId() == pcSlice->getLayerId() )
604              && ( dpbPic->getReconMark() ) 
605            )
606          {
607            for(Int i = dpbPic->getNumAllocatedSlice()-1; i >= 0; i--)
608            {
609              TComSlice *slice = dpbPic->getSlice(i);
610              TComReferencePictureSet *rps = slice->getRPS();
611              slice->setPOC( dpbPic->getSlice(i)->getPOC() - pocAdjustValue );
612
613              // Also adjust the POC value stored in the RPS of each such slice
614              for(Int j = rps->getNumberOfPictures(); j >= 0; j--)
615              {
616                rps->setPOC( j, rps->getPOC(j) - pocAdjustValue );
617              }
618              // Also adjust the value of refPOC
619              for(Int k = 0; k < 2; k++)  // For List 0 and List 1
620              {
621                RefPicList list = (k == 1) ? REF_PIC_LIST_1 : REF_PIC_LIST_0;
622                for(Int j = 0; j < slice->getNumRefIdx(list); j++)
623                {
624                  slice->setRefPOC( slice->getRefPOC(list, j) - pocAdjustValue, list, j);
625                }
626              }
627            }
628          }
629          iterPic++;
630        }
631        m_pcEncTop->setPocAdjustmentValue( m_pcEncTop->getPocAdjustmentValue() + pocAdjustValue );
632      }
633      pcSlice->setPocValueBeforeReset( pcSlice->getPOC() - m_pcEncTop->getPocAdjustmentValue() + pocAdjustValue );
634      pcSlice->setPOC( 0 );
635    }
636#endif
637#if M0040_ADAPTIVE_RESOLUTION_CHANGE
638    if (m_pcEncTop->getAdaptiveResolutionChange() > 0 && m_layerId == 1 && pocCurr > m_pcEncTop->getAdaptiveResolutionChange())
639    {
640      pcSlice->setActiveNumILRRefIdx(0);
641      pcSlice->setInterLayerPredEnabledFlag(false);
642#if M0457_COL_PICTURE_SIGNALING
643      pcSlice->setMFMEnabledFlag(false);
644#if !REMOVE_COL_PICTURE_SIGNALING
645      pcSlice->setAltColIndicationFlag(false);
646#endif
647#endif
648    }
649#endif
650
651#if M0457_IL_SAMPLE_PRED_ONLY_FLAG
652    pcSlice->setNumSamplePredRefLayers( m_pcEncTop->getNumSamplePredRefLayers() );
653    pcSlice->setInterLayerSamplePredOnlyFlag( 0 );
654    if( pcSlice->getNumSamplePredRefLayers() > 0 && pcSlice->getActiveNumILRRefIdx() > 0 )
655    {
656      if( m_pcEncTop->getIlSampleOnlyPred() > 0 )
657      {
658        pcSlice->setInterLayerSamplePredOnlyFlag( true );
659      }
660    }
661#endif
662
663#if IL_SL_SIGNALLING_N0371
664    m_pcEncTop->getScalingList()->setLayerId( m_layerId );
665#endif
666
667    pcSlice->setLastIDR(m_iLastIDR);
668    pcSlice->setSliceIdx(0);
669    //set default slice level flag to the same as SPS level flag
670    pcSlice->setLFCrossSliceBoundaryFlag(  pcSlice->getPPS()->getLoopFilterAcrossSlicesEnabledFlag()  );
671    pcSlice->setScalingList ( m_pcEncTop->getScalingList()  );
672    if(m_pcEncTop->getUseScalingListId() == SCALING_LIST_OFF)
673    {
674#if IL_SL_SIGNALLING_N0371
675      m_pcEncTop->getTrQuant()->setFlatScalingList( m_layerId );
676#else
677      m_pcEncTop->getTrQuant()->setFlatScalingList();
678#endif
679      m_pcEncTop->getTrQuant()->setUseScalingList(false);
680      m_pcEncTop->getSPS()->setScalingListPresentFlag(false);
681      m_pcEncTop->getPPS()->setScalingListPresentFlag(false);
682    }
683    else if(m_pcEncTop->getUseScalingListId() == SCALING_LIST_DEFAULT)
684    {
685#if IL_SL_SIGNALLING_N0371
686      pcSlice->getScalingList()->setLayerId( m_layerId );
687#endif
688
689#if IL_SL_SIGNALLING_N0371
690      pcSlice->setDefaultScalingList ( m_layerId );
691#else
692      pcSlice->setDefaultScalingList ();
693#endif
694
695      m_pcEncTop->getSPS()->setScalingListPresentFlag(false);
696      m_pcEncTop->getPPS()->setScalingListPresentFlag(false);
697      m_pcEncTop->getTrQuant()->setScalingList(pcSlice->getScalingList());
698      m_pcEncTop->getTrQuant()->setUseScalingList(true);
699    }
700    else if(m_pcEncTop->getUseScalingListId() == SCALING_LIST_FILE_READ)
701    {
702#if IL_SL_SIGNALLING_N0371
703      pcSlice->getScalingList()->setLayerId( m_layerId );
704#endif
705
706      if(pcSlice->getScalingList()->xParseScalingList(m_pcCfg->getScalingListFile()))
707      {
708#if IL_SL_SIGNALLING_N0371
709        pcSlice->setDefaultScalingList ( m_layerId );
710#else
711        pcSlice->setDefaultScalingList ();
712#endif
713      }
714#if IL_SL_SIGNALLING_N0371
715      pcSlice->getScalingList()->checkDcOfMatrix( m_layerId );
716#else
717      pcSlice->getScalingList()->checkDcOfMatrix();
718#endif
719      m_pcEncTop->getSPS()->setScalingListPresentFlag(pcSlice->checkDefaultScalingList());
720
721#if IL_SL_SIGNALLING_N0371
722      if( m_layerId > 0 )
723      {
724        m_pcEncTop->getSPS()->setPredScalingListFlag  (true); 
725        m_pcEncTop->getSPS()->setScalingListRefLayerId( 0 ); 
726      }
727#endif
728
729      m_pcEncTop->getPPS()->setScalingListPresentFlag(false);
730
731#if IL_SL_SIGNALLING_N0371
732      if( m_layerId > 0 )
733      {
734        m_pcEncTop->getPPS()->setPredScalingListFlag  (false); 
735        m_pcEncTop->getPPS()->setScalingListRefLayerId( 0   ); 
736      }
737#endif
738
739      m_pcEncTop->getTrQuant()->setScalingList(pcSlice->getScalingList());
740      m_pcEncTop->getTrQuant()->setUseScalingList(true);
741    }
742    else
743    {
744      printf("error : ScalingList == %d no support\n",m_pcEncTop->getUseScalingListId());
745      assert(0);
746    }
747
748    if(pcSlice->getSliceType()==B_SLICE&&m_pcCfg->getGOPEntry(iGOPid).m_sliceType=='P')
749    {
750      pcSlice->setSliceType(P_SLICE);
751    }
752    // Set the nal unit type
753    pcSlice->setNalUnitType(getNalUnitType(pocCurr, m_iLastIDR));
754#if SVC_EXTENSION
755    if (m_layerId > 0)
756    {
757    Int interLayerPredLayerIdcTmp[MAX_VPS_LAYER_ID_PLUS1];
758    Int activeNumILRRefIdxTmp = 0;
759
760      for( Int i = 0; i < pcSlice->getActiveNumILRRefIdx(); i++ )
761      {
762        UInt refLayerIdc = pcSlice->getInterLayerPredLayerIdc(i);
763#if VPS_EXTN_DIRECT_REF_LAYERS
764        TComList<TComPic*> *cListPic = m_ppcTEncTop[m_layerId]->getRefLayerEnc(refLayerIdc)->getListPic();
765#else
766        TComList<TComPic*> *cListPic = m_ppcTEncTop[m_layerId-1]->getListPic();
767#endif
768        pcSlice->setBaseColPic( *cListPic, refLayerIdc );
769
770        // Apply temporal layer restriction to inter-layer prediction
771#if O0225_MAX_TID_FOR_REF_LAYERS
772        Int maxTidIlRefPicsPlus1 = m_pcEncTop->getVPS()->getMaxTidIlRefPicsPlus1(pcSlice->getBaseColPic(refLayerIdc)->getSlice(0)->getLayerId(),m_layerId);
773#else
774        Int maxTidIlRefPicsPlus1 = m_pcEncTop->getVPS()->getMaxTidIlRefPicsPlus1(pcSlice->getBaseColPic(refLayerIdc)->getSlice(0)->getLayerId());
775#endif
776        if( ((Int)(pcSlice->getBaseColPic(refLayerIdc)->getSlice(0)->getTLayer())<=maxTidIlRefPicsPlus1-1) || (maxTidIlRefPicsPlus1==0 && pcSlice->getBaseColPic(refLayerIdc)->getSlice(0)->getRapPicFlag()) )
777        {
778          interLayerPredLayerIdcTmp[activeNumILRRefIdxTmp++] = refLayerIdc; // add picture to the list of valid inter-layer pictures
779        }
780        else
781        {
782          continue; // ILP is not valid due to temporal layer restriction
783        }
784
785        const Window &scalEL = m_pcEncTop->getScaledRefLayerWindow(refLayerIdc);
786
787        Int widthBL   = pcSlice->getBaseColPic(refLayerIdc)->getPicYuvRec()->getWidth();
788        Int heightBL  = pcSlice->getBaseColPic(refLayerIdc)->getPicYuvRec()->getHeight();
789
790        Int widthEL   = pcPic->getPicYuvRec()->getWidth()  - scalEL.getWindowLeftOffset() - scalEL.getWindowRightOffset();
791        Int heightEL  = pcPic->getPicYuvRec()->getHeight() - scalEL.getWindowTopOffset()  - scalEL.getWindowBottomOffset();
792
793        g_mvScalingFactor[refLayerIdc][0] = widthEL  == widthBL  ? 4096 : Clip3(-4096, 4095, ((widthEL  << 8) + (widthBL  >> 1)) / widthBL);
794        g_mvScalingFactor[refLayerIdc][1] = heightEL == heightBL ? 4096 : Clip3(-4096, 4095, ((heightEL << 8) + (heightBL >> 1)) / heightBL);
795
796        g_posScalingFactor[refLayerIdc][0] = ((widthBL  << 16) + (widthEL  >> 1)) / widthEL;
797        g_posScalingFactor[refLayerIdc][1] = ((heightBL << 16) + (heightEL >> 1)) / heightEL;
798
799#if SVC_UPSAMPLING
800        if( pcPic->isSpatialEnhLayer(refLayerIdc))
801        {
802#if O0215_PHASE_ALIGNMENT
803#if O0194_JOINT_US_BITSHIFT
804          m_pcPredSearch->upsampleBasePic( pcSlice, refLayerIdc, pcPic->getFullPelBaseRec(refLayerIdc), pcSlice->getBaseColPic(refLayerIdc)->getPicYuvRec(), pcPic->getPicYuvRec(), pcSlice->getSPS()->getScaledRefLayerWindow(refLayerIdc), pcSlice->getVPS()->getPhaseAlignFlag() );
805#else
806          m_pcPredSearch->upsampleBasePic( refLayerIdc, pcPic->getFullPelBaseRec(refLayerIdc), pcSlice->getBaseColPic(refLayerIdc)->getPicYuvRec(), pcPic->getPicYuvRec(), pcSlice->getSPS()->getScaledRefLayerWindow(refLayerIdc), pcSlice->getVPS()->getPhaseAlignFlag() );
807#endif
808#else
809#if O0194_JOINT_US_BITSHIFT
810          m_pcPredSearch->upsampleBasePic( pcSlice, refLayerIdc, pcPic->getFullPelBaseRec(refLayerIdc), pcSlice->getBaseColPic(refLayerIdc)->getPicYuvRec(), pcPic->getPicYuvRec(), pcSlice->getSPS()->getScaledRefLayerWindow(refLayerIdc) );
811#else
812          m_pcPredSearch->upsampleBasePic( refLayerIdc, pcPic->getFullPelBaseRec(refLayerIdc), pcSlice->getBaseColPic(refLayerIdc)->getPicYuvRec(), pcPic->getPicYuvRec(), pcSlice->getSPS()->getScaledRefLayerWindow(refLayerIdc) );
813#endif
814#endif
815        }
816        else
817        {
818          pcPic->setFullPelBaseRec( refLayerIdc, pcSlice->getBaseColPic(refLayerIdc)->getPicYuvRec() );
819        }
820        pcSlice->setFullPelBaseRec ( refLayerIdc, pcPic->getFullPelBaseRec(refLayerIdc) );
821#endif
822      }
823
824      // Update the list of active inter-layer pictures
825      for ( Int i = 0; i < activeNumILRRefIdxTmp; i++)
826      {
827        pcSlice->setInterLayerPredLayerIdc( interLayerPredLayerIdcTmp[i], i );
828      }
829      pcSlice->setActiveNumILRRefIdx( activeNumILRRefIdxTmp );
830      if ( pcSlice->getActiveNumILRRefIdx() == 0 )
831      {
832        // No valid inter-layer pictures -> disable inter-layer prediction
833        pcSlice->setInterLayerPredEnabledFlag(false);
834      }
835     
836      if( pocCurr % m_pcCfg->getIntraPeriod() == 0 )
837      {
838#if N0147_IRAP_ALIGN_FLAG
839        if(pcSlice->getVPS()->getCrossLayerIrapAlignFlag())
840        {
841          TComList<TComPic*> *cListPic = m_ppcTEncTop[m_layerId]->getRefLayerEnc(0)->getListPic();
842          TComPic* picLayer0 = pcSlice->getRefPic(*cListPic, pcSlice->getPOC() );
843          if(picLayer0)
844          {
845            pcSlice->setNalUnitType(picLayer0->getSlice(0)->getNalUnitType());
846          }
847          else
848          {
849            pcSlice->setNalUnitType(NAL_UNIT_CODED_SLICE_CRA);
850          }
851        }
852        else
853#endif
854        pcSlice->setNalUnitType(NAL_UNIT_CODED_SLICE_CRA);
855
856#if IDR_ALIGNMENT
857        TComList<TComPic*> *cListPic = m_ppcTEncTop[m_layerId]->getRefLayerEnc(0)->getListPic();
858        TComPic* picLayer0 = pcSlice->getRefPic(*cListPic, pcSlice->getPOC() );
859        if( picLayer0->getSlice(0)->getNalUnitType() == NAL_UNIT_CODED_SLICE_IDR_W_RADL || picLayer0->getSlice(0)->getNalUnitType() == NAL_UNIT_CODED_SLICE_IDR_N_LP )
860        {
861          pcSlice->setNalUnitType(picLayer0->getSlice(0)->getNalUnitType());
862        }
863        else
864        {
865          pcSlice->setNalUnitType(NAL_UNIT_CODED_SLICE_CRA);
866        }
867#endif
868      }
869     
870      if( pcSlice->getActiveNumILRRefIdx() == 0 && pcSlice->getNalUnitType() >= NAL_UNIT_CODED_SLICE_BLA_W_LP && pcSlice->getNalUnitType() <= NAL_UNIT_CODED_SLICE_CRA )
871      {
872        pcSlice->setSliceType(I_SLICE);
873      }
874      else if( !m_pcEncTop->getElRapSliceTypeB() )
875      {
876        if( (pcSlice->getNalUnitType() >= NAL_UNIT_CODED_SLICE_BLA_W_LP) &&
877           (pcSlice->getNalUnitType() <= NAL_UNIT_CODED_SLICE_CRA) &&
878           pcSlice->getSliceType() == B_SLICE )
879        {
880          pcSlice->setSliceType(P_SLICE);
881        }
882      }
883    }
884#endif //#if SVC_EXTENSION
885    if(pcSlice->getTemporalLayerNonReferenceFlag())
886    {
887      if (pcSlice->getNalUnitType() == NAL_UNIT_CODED_SLICE_TRAIL_R &&
888          !(m_iGopSize == 1 && pcSlice->getSliceType() == I_SLICE))
889        // Add this condition to avoid POC issues with encoder_intra_main.cfg configuration (see #1127 in bug tracker)
890      {
891        pcSlice->setNalUnitType(NAL_UNIT_CODED_SLICE_TRAIL_N);
892    }
893      if(pcSlice->getNalUnitType()==NAL_UNIT_CODED_SLICE_RADL_R)
894      {
895        pcSlice->setNalUnitType(NAL_UNIT_CODED_SLICE_RADL_N);
896      }
897      if(pcSlice->getNalUnitType()==NAL_UNIT_CODED_SLICE_RASL_R)
898      {
899        pcSlice->setNalUnitType(NAL_UNIT_CODED_SLICE_RASL_N);
900      }
901    }
902
903    // Do decoding refresh marking if any
904    pcSlice->decodingRefreshMarking(m_pocCRA, m_bRefreshPending, rcListPic);
905    m_pcEncTop->selectReferencePictureSet(pcSlice, pocCurr, iGOPid);
906    pcSlice->getRPS()->setNumberOfLongtermPictures(0);
907
908#if FIX1071
909    if ((pcSlice->checkThatAllRefPicsAreAvailable(rcListPic, pcSlice->getRPS(), false) != 0) || (pcSlice->isIRAP()))
910    {
911      pcSlice->createExplicitReferencePictureSetFromReference(rcListPic, pcSlice->getRPS(), pcSlice->isIRAP());
912    }
913#else
914    if(pcSlice->checkThatAllRefPicsAreAvailable(rcListPic, pcSlice->getRPS(), false) != 0)
915    {
916      pcSlice->createExplicitReferencePictureSetFromReference(rcListPic, pcSlice->getRPS());
917    }
918#endif
919    pcSlice->applyReferencePictureSet(rcListPic, pcSlice->getRPS());
920
921    if(pcSlice->getTLayer() > 0)
922    {
923      if(pcSlice->isTemporalLayerSwitchingPoint(rcListPic) || pcSlice->getSPS()->getTemporalIdNestingFlag())
924      {
925#if ALIGN_TSA_STSA_PICS
926        if( pcSlice->getLayerId() > 0 )
927        {
928          Bool oneRefLayerTSA = false, oneRefLayerNotTSA = false;
929          for( Int i = 0; i < pcSlice->getLayerId(); i++)
930          {
931            TComList<TComPic *> *cListPic = m_ppcTEncTop[i]->getListPic();
932            TComPic *lowerLayerPic = pcSlice->getRefPic(*cListPic, pcSlice->getPOC());
933            if( lowerLayerPic && pcSlice->getVPS()->getDirectDependencyFlag(pcSlice->getLayerId(), i) )
934            {
935              if( ( lowerLayerPic->getSlice(0)->getNalUnitType() == NAL_UNIT_CODED_SLICE_TSA_N ) ||
936                  ( lowerLayerPic->getSlice(0)->getNalUnitType() == NAL_UNIT_CODED_SLICE_TLA_R ) 
937                )
938              {
939                if(pcSlice->getTemporalLayerNonReferenceFlag() )
940                {
941                  pcSlice->setNalUnitType(NAL_UNIT_CODED_SLICE_TSA_N);
942                }
943                else
944                {
945                  pcSlice->setNalUnitType(NAL_UNIT_CODED_SLICE_TLA_R );
946                }
947                oneRefLayerTSA = true;
948              }
949              else
950              {
951                oneRefLayerNotTSA = true;
952              }
953            }
954          }
955          assert( !( oneRefLayerNotTSA && oneRefLayerTSA ) ); // Only one variable should be true - failure of this assert means
956                                                                // that two independent reference layers that are not dependent on
957                                                                // each other, but are reference for current layer have inconsistency
958          if( oneRefLayerNotTSA /*&& !oneRefLayerTSA*/ )          // No reference layer is TSA - set current as TRAIL
959          {
960            if(pcSlice->getTemporalLayerNonReferenceFlag() )
961            {
962              pcSlice->setNalUnitType( NAL_UNIT_CODED_SLICE_TRAIL_N );
963            }
964            else
965            {
966              pcSlice->setNalUnitType( NAL_UNIT_CODED_SLICE_TRAIL_R );
967            }
968          }
969          else  // This means there is no reference layer picture for current picture in this AU
970          {
971            if(pcSlice->getTemporalLayerNonReferenceFlag() )
972            {
973              pcSlice->setNalUnitType(NAL_UNIT_CODED_SLICE_TSA_N);
974            }
975            else
976            {
977              pcSlice->setNalUnitType(NAL_UNIT_CODED_SLICE_TLA_R );
978            }
979          }
980        }
981#else
982        if(pcSlice->getTemporalLayerNonReferenceFlag())
983        {
984          pcSlice->setNalUnitType(NAL_UNIT_CODED_SLICE_TSA_N);
985        }
986        else
987        {
988          pcSlice->setNalUnitType(NAL_UNIT_CODED_SLICE_TLA_R);
989        }
990#endif
991      }
992      else if(pcSlice->isStepwiseTemporalLayerSwitchingPointCandidate(rcListPic))
993      {
994        Bool isSTSA=true;
995        for(Int ii=iGOPid+1;(ii<m_pcCfg->getGOPSize() && isSTSA==true);ii++)
996        {
997          Int lTid= m_pcCfg->getGOPEntry(ii).m_temporalId;
998          if(lTid==pcSlice->getTLayer()) 
999          {
1000            TComReferencePictureSet* nRPS = pcSlice->getSPS()->getRPSList()->getReferencePictureSet(ii);
1001            for(Int jj=0;jj<nRPS->getNumberOfPictures();jj++)
1002            {
1003              if(nRPS->getUsed(jj)) 
1004              {
1005                Int tPoc=m_pcCfg->getGOPEntry(ii).m_POC+nRPS->getDeltaPOC(jj);
1006                Int kk=0;
1007                for(kk=0;kk<m_pcCfg->getGOPSize();kk++)
1008                {
1009                  if(m_pcCfg->getGOPEntry(kk).m_POC==tPoc)
1010                    break;
1011                }
1012                Int tTid=m_pcCfg->getGOPEntry(kk).m_temporalId;
1013                if(tTid >= pcSlice->getTLayer())
1014                {
1015                  isSTSA=false;
1016                  break;
1017                }
1018              }
1019            }
1020          }
1021        }
1022        if(isSTSA==true)
1023        {   
1024#if ALIGN_TSA_STSA_PICS
1025          if( pcSlice->getLayerId() > 0 )
1026          {
1027            Bool oneRefLayerSTSA = false, oneRefLayerNotSTSA = false;
1028            for( Int i = 0; i < pcSlice->getLayerId(); i++)
1029            {
1030              TComList<TComPic *> *cListPic = m_ppcTEncTop[i]->getListPic();
1031              TComPic *lowerLayerPic = pcSlice->getRefPic(*cListPic, pcSlice->getPOC());
1032              if( lowerLayerPic && pcSlice->getVPS()->getDirectDependencyFlag(pcSlice->getLayerId(), i) )
1033              {
1034                if( ( lowerLayerPic->getSlice(0)->getNalUnitType() == NAL_UNIT_CODED_SLICE_STSA_N ) ||
1035                    ( lowerLayerPic->getSlice(0)->getNalUnitType() == NAL_UNIT_CODED_SLICE_STSA_R ) 
1036                  )
1037                {
1038                  if(pcSlice->getTemporalLayerNonReferenceFlag() )
1039                  {
1040                    pcSlice->setNalUnitType(NAL_UNIT_CODED_SLICE_STSA_N);
1041                  }
1042                  else
1043                  {
1044                    pcSlice->setNalUnitType(NAL_UNIT_CODED_SLICE_STSA_R );
1045                  }
1046                  oneRefLayerSTSA = true;
1047                }
1048                else
1049                {
1050                  oneRefLayerNotSTSA = true;
1051                }
1052              }
1053            }
1054            assert( !( oneRefLayerNotSTSA && oneRefLayerSTSA ) ); // Only one variable should be true - failure of this assert means
1055                                                                  // that two independent reference layers that are not dependent on
1056                                                                  // each other, but are reference for current layer have inconsistency
1057            if( oneRefLayerNotSTSA /*&& !oneRefLayerSTSA*/ )          // No reference layer is STSA - set current as TRAIL
1058            {
1059              if(pcSlice->getTemporalLayerNonReferenceFlag() )
1060              {
1061                pcSlice->setNalUnitType( NAL_UNIT_CODED_SLICE_TRAIL_N );
1062              }
1063              else
1064              {
1065                pcSlice->setNalUnitType( NAL_UNIT_CODED_SLICE_TRAIL_R );
1066              }
1067            }
1068            else  // This means there is no reference layer picture for current picture in this AU
1069            {
1070              if(pcSlice->getTemporalLayerNonReferenceFlag() )
1071              {
1072                pcSlice->setNalUnitType(NAL_UNIT_CODED_SLICE_STSA_N);
1073              }
1074              else
1075              {
1076                pcSlice->setNalUnitType(NAL_UNIT_CODED_SLICE_STSA_R );
1077              }
1078            }
1079          }
1080#else
1081          if(pcSlice->getTemporalLayerNonReferenceFlag())
1082          {
1083            pcSlice->setNalUnitType(NAL_UNIT_CODED_SLICE_STSA_N);
1084          }
1085          else
1086          {
1087            pcSlice->setNalUnitType(NAL_UNIT_CODED_SLICE_STSA_R);
1088          }
1089#endif
1090        }
1091      }
1092    }
1093    arrangeLongtermPicturesInRPS(pcSlice, rcListPic);
1094    TComRefPicListModification* refPicListModification = pcSlice->getRefPicListModification();
1095    refPicListModification->setRefPicListModificationFlagL0(0);
1096    refPicListModification->setRefPicListModificationFlagL1(0);
1097    pcSlice->setNumRefIdx(REF_PIC_LIST_0,min(m_pcCfg->getGOPEntry(iGOPid).m_numRefPicsActive,pcSlice->getRPS()->getNumberOfPictures()));
1098    pcSlice->setNumRefIdx(REF_PIC_LIST_1,min(m_pcCfg->getGOPEntry(iGOPid).m_numRefPicsActive,pcSlice->getRPS()->getNumberOfPictures()));
1099
1100#if SVC_EXTENSION
1101    if( m_layerId > 0 && pcSlice->getActiveNumILRRefIdx() )
1102    {
1103#if RESTR_CHK
1104#if POC_RESET_FLAG
1105      if ( pocCurr > 0          && pcSlice->isRADL() && pcPic->getSlice(0)->getBaseColPic(pcPic->getSlice(0)->getInterLayerPredLayerIdc(0))->getSlice(0)->isRASL())
1106#else
1107      if (pcSlice->getPOC()>0  && pcSlice->isRADL() && pcPic->getSlice(0)->getBaseColPic(pcPic->getSlice(0)->getInterLayerPredLayerIdc(0))->getSlice(0)->isRASL())
1108#endif
1109      {
1110#if JCTVC_M0458_INTERLAYER_RPS_SIG
1111        pcSlice->setActiveNumILRRefIdx(0);
1112        pcSlice->setInterLayerPredEnabledFlag(0);
1113#else
1114        pcSlice->setNumILRRefIdx(0);
1115#endif
1116      }
1117#endif
1118#if JCTVC_M0458_INTERLAYER_RPS_SIG
1119      if( pcSlice->getNalUnitType() >= NAL_UNIT_CODED_SLICE_BLA_W_LP && pcSlice->getNalUnitType() <= NAL_UNIT_CODED_SLICE_CRA )
1120      {
1121        pcSlice->setNumRefIdx(REF_PIC_LIST_0, pcSlice->getActiveNumILRRefIdx());
1122        pcSlice->setNumRefIdx(REF_PIC_LIST_1, pcSlice->getActiveNumILRRefIdx());
1123      }
1124      else
1125      {
1126        pcSlice->setNumRefIdx(REF_PIC_LIST_0, pcSlice->getNumRefIdx(REF_PIC_LIST_0)+pcSlice->getActiveNumILRRefIdx());
1127        pcSlice->setNumRefIdx(REF_PIC_LIST_1, pcSlice->getNumRefIdx(REF_PIC_LIST_1)+pcSlice->getActiveNumILRRefIdx());
1128      }
1129#else
1130      if( pcSlice->getNalUnitType() >= NAL_UNIT_CODED_SLICE_BLA_W_LP && pcSlice->getNalUnitType() <= NAL_UNIT_CODED_SLICE_CRA )
1131      {
1132        pcSlice->setNumRefIdx(REF_PIC_LIST_0, pcSlice->getNumILRRefIdx());
1133        pcSlice->setNumRefIdx(REF_PIC_LIST_1, pcSlice->getNumILRRefIdx());
1134      }
1135      else
1136      {
1137        pcSlice->setNumRefIdx(REF_PIC_LIST_0, pcSlice->getNumRefIdx(REF_PIC_LIST_0)+pcSlice->getNumILRRefIdx());
1138        pcSlice->setNumRefIdx(REF_PIC_LIST_1, pcSlice->getNumRefIdx(REF_PIC_LIST_1)+pcSlice->getNumILRRefIdx());
1139      }
1140#endif
1141    }
1142#endif //SVC_EXTENSION
1143
1144#if ADAPTIVE_QP_SELECTION
1145    pcSlice->setTrQuant( m_pcEncTop->getTrQuant() );
1146#endif     
1147
1148#if SVC_EXTENSION
1149#if M0457_COL_PICTURE_SIGNALING && !REMOVE_COL_PICTURE_SIGNALING
1150    if ( pcSlice->getSliceType() == B_SLICE && !(pcSlice->getActiveNumILRRefIdx() > 0 && m_pcEncTop->getNumMotionPredRefLayers() > 0) )
1151#else
1152    if( pcSlice->getSliceType() == B_SLICE )
1153#endif
1154    {
1155      pcSlice->setColFromL0Flag(1-uiColDir);
1156    }
1157
1158    //  Set reference list
1159    if(m_layerId ==  0 || ( m_layerId > 0 && pcSlice->getActiveNumILRRefIdx() == 0 ) )
1160    {
1161      pcSlice->setRefPicList( rcListPic);
1162    }
1163
1164    if( m_layerId > 0 && pcSlice->getActiveNumILRRefIdx() )
1165    {
1166      pcSlice->setILRPic( m_pcEncTop->getIlpList() );
1167#if REF_IDX_MFM
1168#if M0457_COL_PICTURE_SIGNALING
1169      if( pcSlice->getMFMEnabledFlag() )
1170#else
1171      if( pcSlice->getSPS()->getMFMEnabledFlag() )
1172#endif
1173      {
1174        pcSlice->setRefPOCListILP(m_pcEncTop->getIlpList(), pcSlice->getBaseColPic());
1175#if M0457_COL_PICTURE_SIGNALING && !REMOVE_COL_PICTURE_SIGNALING
1176        pcSlice->setMotionPredIlp(getMotionPredIlp(pcSlice));
1177#endif
1178      }
1179#else
1180      //  Set reference list
1181      pcSlice->setRefPicList ( rcListPic );
1182#endif //SVC_EXTENSION
1183      pcSlice->setRefPicListModificationSvc();
1184      pcSlice->setRefPicList( rcListPic, false, m_pcEncTop->getIlpList());
1185
1186#if REF_IDX_MFM
1187#if M0457_COL_PICTURE_SIGNALING
1188#if REMOVE_COL_PICTURE_SIGNALING
1189      if( pcSlice->getMFMEnabledFlag() )
1190#else
1191      if( pcSlice->getMFMEnabledFlag() && !(pcSlice->getActiveNumILRRefIdx() > 0 && m_pcEncTop->getNumMotionPredRefLayers() > 0) )
1192#endif
1193#else
1194      if( pcSlice->getSPS()->getMFMEnabledFlag() )
1195#endif
1196      {
1197        Bool found         = false;
1198        UInt ColFromL0Flag = pcSlice->getColFromL0Flag();
1199        UInt ColRefIdx     = pcSlice->getColRefIdx();
1200
1201        for(Int colIdx = 0; colIdx < pcSlice->getNumRefIdx( RefPicList(1 - ColFromL0Flag) ); colIdx++) 
1202        { 
1203          if( pcSlice->getRefPic( RefPicList(1 - ColFromL0Flag), colIdx)->isILR(m_layerId) 
1204#if MFM_ENCCONSTRAINT
1205            && pcSlice->getBaseColPic( *m_ppcTEncTop[pcSlice->getRefPic( RefPicList(1 - ColFromL0Flag), colIdx)->getLayerId()]->getListPic() )->checkSameRefInfo() == true 
1206#endif
1207            ) 
1208          { 
1209            ColRefIdx = colIdx; 
1210            found = true;
1211            break; 
1212          }
1213        }
1214
1215        if( found == false )
1216        {
1217          ColFromL0Flag = 1 - ColFromL0Flag;
1218          for(Int colIdx = 0; colIdx < pcSlice->getNumRefIdx( RefPicList(1 - ColFromL0Flag) ); colIdx++) 
1219          { 
1220            if( pcSlice->getRefPic( RefPicList(1 - ColFromL0Flag), colIdx)->isILR(m_layerId) 
1221#if MFM_ENCCONSTRAINT
1222              && pcSlice->getBaseColPic( *m_ppcTEncTop[pcSlice->getRefPic( RefPicList(1 - ColFromL0Flag), colIdx)->getLayerId()]->getListPic() )->checkSameRefInfo() == true 
1223#endif
1224              ) 
1225            { 
1226              ColRefIdx = colIdx; 
1227              found = true; 
1228              break; 
1229            } 
1230          }
1231        }
1232
1233        if(found == true)
1234        {
1235          pcSlice->setColFromL0Flag(ColFromL0Flag);
1236          pcSlice->setColRefIdx(ColRefIdx);
1237        }
1238      }
1239#endif
1240    }
1241#else //SVC_EXTENSION
1242    //  Set reference list
1243    pcSlice->setRefPicList ( rcListPic );
1244#endif //#if SVC_EXTENSION
1245
1246    //  Slice info. refinement
1247    if ( (pcSlice->getSliceType() == B_SLICE) && (pcSlice->getNumRefIdx(REF_PIC_LIST_1) == 0) )
1248    {
1249      pcSlice->setSliceType ( P_SLICE );
1250    }
1251
1252#if M0457_IL_SAMPLE_PRED_ONLY_FLAG
1253    if (pcSlice->getSliceType() == B_SLICE && m_pcEncTop->getIlSampleOnlyPred() == 0)
1254#else
1255    if (pcSlice->getSliceType() == B_SLICE)
1256#endif
1257    {
1258#if !SVC_EXTENSION
1259      pcSlice->setColFromL0Flag(1-uiColDir);
1260#endif
1261      Bool bLowDelay = true;
1262      Int  iCurrPOC  = pcSlice->getPOC();
1263      Int iRefIdx = 0;
1264
1265      for (iRefIdx = 0; iRefIdx < pcSlice->getNumRefIdx(REF_PIC_LIST_0) && bLowDelay; iRefIdx++)
1266      {
1267        if ( pcSlice->getRefPic(REF_PIC_LIST_0, iRefIdx)->getPOC() > iCurrPOC )
1268        {
1269          bLowDelay = false;
1270        }
1271      }
1272      for (iRefIdx = 0; iRefIdx < pcSlice->getNumRefIdx(REF_PIC_LIST_1) && bLowDelay; iRefIdx++)
1273      {
1274        if ( pcSlice->getRefPic(REF_PIC_LIST_1, iRefIdx)->getPOC() > iCurrPOC )
1275        {
1276          bLowDelay = false;
1277        }
1278      }
1279
1280      pcSlice->setCheckLDC(bLowDelay); 
1281    }
1282    else
1283    {
1284      pcSlice->setCheckLDC(true); 
1285    }
1286
1287    uiColDir = 1-uiColDir;
1288
1289    //-------------------------------------------------------------
1290    pcSlice->setRefPOCList();
1291
1292    pcSlice->setList1IdxToList0Idx();
1293
1294    if (m_pcEncTop->getTMVPModeId() == 2)
1295    {
1296      if (iGOPid == 0) // first picture in SOP (i.e. forward B)
1297      {
1298        pcSlice->setEnableTMVPFlag(0);
1299      }
1300      else
1301      {
1302        // Note: pcSlice->getColFromL0Flag() is assumed to be always 0 and getcolRefIdx() is always 0.
1303        pcSlice->setEnableTMVPFlag(1);
1304      }
1305      pcSlice->getSPS()->setTMVPFlagsPresent(1);
1306    }
1307    else if (m_pcEncTop->getTMVPModeId() == 1)
1308    {
1309      pcSlice->getSPS()->setTMVPFlagsPresent(1);
1310#if SVC_EXTENSION
1311      if( pcSlice->getIdrPicFlag() )
1312      {
1313        pcSlice->setEnableTMVPFlag(0);
1314      }
1315      else
1316#endif
1317      pcSlice->setEnableTMVPFlag(1);
1318    }
1319    else
1320    {
1321      pcSlice->getSPS()->setTMVPFlagsPresent(0);
1322      pcSlice->setEnableTMVPFlag(0);
1323    }
1324    /////////////////////////////////////////////////////////////////////////////////////////////////// Compress a slice
1325    //  Slice compression
1326    if (m_pcCfg->getUseASR())
1327    {
1328      m_pcSliceEncoder->setSearchRange(pcSlice);
1329    }
1330
1331    Bool bGPBcheck=false;
1332    if ( pcSlice->getSliceType() == B_SLICE)
1333    {
1334      if ( pcSlice->getNumRefIdx(RefPicList( 0 ) ) == pcSlice->getNumRefIdx(RefPicList( 1 ) ) )
1335      {
1336        bGPBcheck=true;
1337        Int i;
1338        for ( i=0; i < pcSlice->getNumRefIdx(RefPicList( 1 ) ); i++ )
1339        {
1340          if ( pcSlice->getRefPOC(RefPicList(1), i) != pcSlice->getRefPOC(RefPicList(0), i) ) 
1341          {
1342            bGPBcheck=false;
1343            break;
1344          }
1345        }
1346      }
1347    }
1348    if(bGPBcheck)
1349    {
1350      pcSlice->setMvdL1ZeroFlag(true);
1351    }
1352    else
1353    {
1354      pcSlice->setMvdL1ZeroFlag(false);
1355    }
1356    pcPic->getSlice(pcSlice->getSliceIdx())->setMvdL1ZeroFlag(pcSlice->getMvdL1ZeroFlag());
1357
1358#if RATE_CONTROL_LAMBDA_DOMAIN
1359    Double lambda            = 0.0;
1360    Int actualHeadBits       = 0;
1361    Int actualTotalBits      = 0;
1362    Int estimatedBits        = 0;
1363    Int tmpBitsBeforeWriting = 0;
1364    if ( m_pcCfg->getUseRateCtrl() )
1365    {
1366      Int frameLevel = m_pcRateCtrl->getRCSeq()->getGOPID2Level( iGOPid );
1367      if ( pcPic->getSlice(0)->getSliceType() == I_SLICE )
1368      {
1369        frameLevel = 0;
1370      }
1371      m_pcRateCtrl->initRCPic( frameLevel );
1372      estimatedBits = m_pcRateCtrl->getRCPic()->getTargetBits();
1373
1374      Int sliceQP = m_pcCfg->getInitialQP();
1375#if POC_RESET_FLAG
1376      if ( ( pocCurr == 0 && m_pcCfg->getInitialQP() > 0 ) || ( frameLevel == 0 && m_pcCfg->getForceIntraQP() ) ) // QP is specified
1377#else
1378      if ( ( pcSlice->getPOC() == 0 && m_pcCfg->getInitialQP() > 0 ) || ( frameLevel == 0 && m_pcCfg->getForceIntraQP() ) ) // QP is specified
1379#endif
1380      {
1381        Int    NumberBFrames = ( m_pcCfg->getGOPSize() - 1 );
1382        Double dLambda_scale = 1.0 - Clip3( 0.0, 0.5, 0.05*(Double)NumberBFrames );
1383        Double dQPFactor     = 0.57*dLambda_scale;
1384        Int    SHIFT_QP      = 12;
1385        Int    bitdepth_luma_qp_scale = 0;
1386        Double qp_temp = (Double) sliceQP + bitdepth_luma_qp_scale - SHIFT_QP;
1387        lambda = dQPFactor*pow( 2.0, qp_temp/3.0 );
1388      }
1389      else if ( frameLevel == 0 )   // intra case, but use the model
1390      {
1391#if RATE_CONTROL_INTRA
1392        m_pcSliceEncoder->calCostSliceI(pcPic);
1393#endif
1394        if ( m_pcCfg->getIntraPeriod() != 1 )   // do not refine allocated bits for all intra case
1395        {
1396          Int bits = m_pcRateCtrl->getRCSeq()->getLeftAverageBits();
1397#if RATE_CONTROL_INTRA
1398          bits = m_pcRateCtrl->getRCPic()->getRefineBitsForIntra( bits );
1399#else
1400          bits = m_pcRateCtrl->getRCSeq()->getRefineBitsForIntra( bits );
1401#endif
1402          if ( bits < 200 )
1403          {
1404            bits = 200;
1405          }
1406          m_pcRateCtrl->getRCPic()->setTargetBits( bits );
1407        }
1408
1409        list<TEncRCPic*> listPreviousPicture = m_pcRateCtrl->getPicList();
1410#if RATE_CONTROL_INTRA
1411        m_pcRateCtrl->getRCPic()->getLCUInitTargetBits();
1412        lambda  = m_pcRateCtrl->getRCPic()->estimatePicLambda( listPreviousPicture, pcSlice->getSliceType());
1413#else
1414        lambda  = m_pcRateCtrl->getRCPic()->estimatePicLambda( listPreviousPicture );
1415#endif
1416        sliceQP = m_pcRateCtrl->getRCPic()->estimatePicQP( lambda, listPreviousPicture );
1417      }
1418      else    // normal case
1419      {
1420        list<TEncRCPic*> listPreviousPicture = m_pcRateCtrl->getPicList();
1421#if RATE_CONTROL_INTRA
1422        lambda  = m_pcRateCtrl->getRCPic()->estimatePicLambda( listPreviousPicture, pcSlice->getSliceType());
1423#else
1424        lambda  = m_pcRateCtrl->getRCPic()->estimatePicLambda( listPreviousPicture );
1425#endif
1426        sliceQP = m_pcRateCtrl->getRCPic()->estimatePicQP( lambda, listPreviousPicture );
1427      }
1428
1429#if REPN_FORMAT_IN_VPS
1430      sliceQP = Clip3( -pcSlice->getQpBDOffsetY(), MAX_QP, sliceQP );
1431#else
1432      sliceQP = Clip3( -pcSlice->getSPS()->getQpBDOffsetY(), MAX_QP, sliceQP );
1433#endif
1434      m_pcRateCtrl->getRCPic()->setPicEstQP( sliceQP );
1435
1436      m_pcSliceEncoder->resetQP( pcPic, sliceQP, lambda );
1437    }
1438#endif
1439
1440    UInt uiNumSlices = 1;
1441
1442    UInt uiInternalAddress = pcPic->getNumPartInCU()-4;
1443    UInt uiExternalAddress = pcPic->getPicSym()->getNumberOfCUsInFrame()-1;
1444    UInt uiPosX = ( uiExternalAddress % pcPic->getFrameWidthInCU() ) * g_uiMaxCUWidth+ g_auiRasterToPelX[ g_auiZscanToRaster[uiInternalAddress] ];
1445    UInt uiPosY = ( uiExternalAddress / pcPic->getFrameWidthInCU() ) * g_uiMaxCUHeight+ g_auiRasterToPelY[ g_auiZscanToRaster[uiInternalAddress] ];
1446#if REPN_FORMAT_IN_VPS
1447    UInt uiWidth = pcSlice->getPicWidthInLumaSamples();
1448    UInt uiHeight = pcSlice->getPicHeightInLumaSamples();
1449#else
1450    UInt uiWidth = pcSlice->getSPS()->getPicWidthInLumaSamples();
1451    UInt uiHeight = pcSlice->getSPS()->getPicHeightInLumaSamples();
1452#endif
1453    while(uiPosX>=uiWidth||uiPosY>=uiHeight) 
1454    {
1455      uiInternalAddress--;
1456      uiPosX = ( uiExternalAddress % pcPic->getFrameWidthInCU() ) * g_uiMaxCUWidth+ g_auiRasterToPelX[ g_auiZscanToRaster[uiInternalAddress] ];
1457      uiPosY = ( uiExternalAddress / pcPic->getFrameWidthInCU() ) * g_uiMaxCUHeight+ g_auiRasterToPelY[ g_auiZscanToRaster[uiInternalAddress] ];
1458    }
1459    uiInternalAddress++;
1460    if(uiInternalAddress==pcPic->getNumPartInCU()) 
1461    {
1462      uiInternalAddress = 0;
1463      uiExternalAddress++;
1464    }
1465    UInt uiRealEndAddress = uiExternalAddress*pcPic->getNumPartInCU()+uiInternalAddress;
1466
1467    UInt uiCummulativeTileWidth;
1468    UInt uiCummulativeTileHeight;
1469    Int  p, j;
1470    UInt uiEncCUAddr;
1471
1472    //set NumColumnsMinus1 and NumRowsMinus1
1473    pcPic->getPicSym()->setNumColumnsMinus1( pcSlice->getPPS()->getNumColumnsMinus1() );
1474    pcPic->getPicSym()->setNumRowsMinus1( pcSlice->getPPS()->getNumRowsMinus1() );
1475
1476    //create the TComTileArray
1477    pcPic->getPicSym()->xCreateTComTileArray();
1478
1479    if( pcSlice->getPPS()->getUniformSpacingFlag() == 1 )
1480    {
1481      //set the width for each tile
1482      for(j=0; j < pcPic->getPicSym()->getNumRowsMinus1()+1; j++)
1483      {
1484        for(p=0; p < pcPic->getPicSym()->getNumColumnsMinus1()+1; p++)
1485        {
1486          pcPic->getPicSym()->getTComTile( j * (pcPic->getPicSym()->getNumColumnsMinus1()+1) + p )->
1487            setTileWidth( (p+1)*pcPic->getPicSym()->getFrameWidthInCU()/(pcPic->getPicSym()->getNumColumnsMinus1()+1) 
1488            - (p*pcPic->getPicSym()->getFrameWidthInCU())/(pcPic->getPicSym()->getNumColumnsMinus1()+1) );
1489        }
1490      }
1491
1492      //set the height for each tile
1493      for(j=0; j < pcPic->getPicSym()->getNumColumnsMinus1()+1; j++)
1494      {
1495        for(p=0; p < pcPic->getPicSym()->getNumRowsMinus1()+1; p++)
1496        {
1497          pcPic->getPicSym()->getTComTile( p * (pcPic->getPicSym()->getNumColumnsMinus1()+1) + j )->
1498            setTileHeight( (p+1)*pcPic->getPicSym()->getFrameHeightInCU()/(pcPic->getPicSym()->getNumRowsMinus1()+1) 
1499            - (p*pcPic->getPicSym()->getFrameHeightInCU())/(pcPic->getPicSym()->getNumRowsMinus1()+1) );   
1500        }
1501      }
1502    }
1503    else
1504    {
1505      //set the width for each tile
1506      for(j=0; j < pcPic->getPicSym()->getNumRowsMinus1()+1; j++)
1507      {
1508        uiCummulativeTileWidth = 0;
1509        for(p=0; p < pcPic->getPicSym()->getNumColumnsMinus1(); p++)
1510        {
1511          pcPic->getPicSym()->getTComTile( j * (pcPic->getPicSym()->getNumColumnsMinus1()+1) + p )->setTileWidth( pcSlice->getPPS()->getColumnWidth(p) );
1512          uiCummulativeTileWidth += pcSlice->getPPS()->getColumnWidth(p);
1513        }
1514        pcPic->getPicSym()->getTComTile(j * (pcPic->getPicSym()->getNumColumnsMinus1()+1) + p)->setTileWidth( pcPic->getPicSym()->getFrameWidthInCU()-uiCummulativeTileWidth );
1515      }
1516
1517      //set the height for each tile
1518      for(j=0; j < pcPic->getPicSym()->getNumColumnsMinus1()+1; j++)
1519      {
1520        uiCummulativeTileHeight = 0;
1521        for(p=0; p < pcPic->getPicSym()->getNumRowsMinus1(); p++)
1522        {
1523          pcPic->getPicSym()->getTComTile( p * (pcPic->getPicSym()->getNumColumnsMinus1()+1) + j )->setTileHeight( pcSlice->getPPS()->getRowHeight(p) );
1524          uiCummulativeTileHeight += pcSlice->getPPS()->getRowHeight(p);
1525        }
1526        pcPic->getPicSym()->getTComTile(p * (pcPic->getPicSym()->getNumColumnsMinus1()+1) + j)->setTileHeight( pcPic->getPicSym()->getFrameHeightInCU()-uiCummulativeTileHeight );
1527      }
1528    }
1529    //intialize each tile of the current picture
1530    pcPic->getPicSym()->xInitTiles();
1531
1532#if N0383_IL_CONSTRAINED_TILE_SETS_SEI
1533    if (m_pcCfg->getInterLayerConstrainedTileSetsSEIEnabled())
1534    {
1535      xBuildTileSetsMap(pcPic->getPicSym());
1536    }
1537#endif
1538
1539    // Allocate some coders, now we know how many tiles there are.
1540    Int iNumSubstreams = pcSlice->getPPS()->getNumSubstreams();
1541
1542    //generate the Coding Order Map and Inverse Coding Order Map
1543    for(p=0, uiEncCUAddr=0; p<pcPic->getPicSym()->getNumberOfCUsInFrame(); p++, uiEncCUAddr = pcPic->getPicSym()->xCalculateNxtCUAddr(uiEncCUAddr))
1544    {
1545      pcPic->getPicSym()->setCUOrderMap(p, uiEncCUAddr);
1546      pcPic->getPicSym()->setInverseCUOrderMap(uiEncCUAddr, p);
1547    }
1548    pcPic->getPicSym()->setCUOrderMap(pcPic->getPicSym()->getNumberOfCUsInFrame(), pcPic->getPicSym()->getNumberOfCUsInFrame());   
1549    pcPic->getPicSym()->setInverseCUOrderMap(pcPic->getPicSym()->getNumberOfCUsInFrame(), pcPic->getPicSym()->getNumberOfCUsInFrame());
1550
1551    // Allocate some coders, now we know how many tiles there are.
1552    m_pcEncTop->createWPPCoders(iNumSubstreams);
1553    pcSbacCoders = m_pcEncTop->getSbacCoders();
1554    pcSubstreamsOut = new TComOutputBitstream[iNumSubstreams];
1555
1556    UInt startCUAddrSliceIdx = 0; // used to index "m_uiStoredStartCUAddrForEncodingSlice" containing locations of slice boundaries
1557    UInt startCUAddrSlice    = 0; // used to keep track of current slice's starting CU addr.
1558    pcSlice->setSliceCurStartCUAddr( startCUAddrSlice ); // Setting "start CU addr" for current slice
1559    m_storedStartCUAddrForEncodingSlice.clear();
1560
1561    UInt startCUAddrSliceSegmentIdx = 0; // used to index "m_uiStoredStartCUAddrForEntropyEncodingSlice" containing locations of slice boundaries
1562    UInt startCUAddrSliceSegment    = 0; // used to keep track of current Dependent slice's starting CU addr.
1563    pcSlice->setSliceSegmentCurStartCUAddr( startCUAddrSliceSegment ); // Setting "start CU addr" for current Dependent slice
1564
1565    m_storedStartCUAddrForEncodingSliceSegment.clear();
1566    UInt nextCUAddr = 0;
1567    m_storedStartCUAddrForEncodingSlice.push_back (nextCUAddr);
1568    startCUAddrSliceIdx++;
1569    m_storedStartCUAddrForEncodingSliceSegment.push_back(nextCUAddr);
1570    startCUAddrSliceSegmentIdx++;
1571#if AVC_BASE
1572    if( m_layerId == 0 && m_pcEncTop->getVPS()->getAvcBaseLayerFlag() )
1573    {
1574      pcPic->getPicYuvOrg()->copyToPic( pcPic->getPicYuvRec() );
1575#if AVC_SYNTAX
1576      pcPic->readBLSyntax( m_ppcTEncTop[0]->getBLSyntaxFile(), SYNTAX_BYTES );
1577#endif
1578      return;
1579    }
1580#endif
1581
1582    while(nextCUAddr<uiRealEndAddress) // determine slice boundaries
1583    {
1584      pcSlice->setNextSlice       ( false );
1585      pcSlice->setNextSliceSegment( false );
1586      assert(pcPic->getNumAllocatedSlice() == startCUAddrSliceIdx);
1587      m_pcSliceEncoder->precompressSlice( pcPic );
1588      m_pcSliceEncoder->compressSlice   ( pcPic );
1589
1590      Bool bNoBinBitConstraintViolated = (!pcSlice->isNextSlice() && !pcSlice->isNextSliceSegment());
1591      if (pcSlice->isNextSlice() || (bNoBinBitConstraintViolated && m_pcCfg->getSliceMode()==FIXED_NUMBER_OF_LCU))
1592      {
1593        startCUAddrSlice = pcSlice->getSliceCurEndCUAddr();
1594        // Reconstruction slice
1595        m_storedStartCUAddrForEncodingSlice.push_back(startCUAddrSlice);
1596        startCUAddrSliceIdx++;
1597        // Dependent slice
1598        if (startCUAddrSliceSegmentIdx>0 && m_storedStartCUAddrForEncodingSliceSegment[startCUAddrSliceSegmentIdx-1] != startCUAddrSlice)
1599        {
1600          m_storedStartCUAddrForEncodingSliceSegment.push_back(startCUAddrSlice);
1601          startCUAddrSliceSegmentIdx++;
1602        }
1603
1604        if (startCUAddrSlice < uiRealEndAddress)
1605        {
1606          pcPic->allocateNewSlice();         
1607          pcPic->setCurrSliceIdx                  ( startCUAddrSliceIdx-1 );
1608          m_pcSliceEncoder->setSliceIdx           ( startCUAddrSliceIdx-1 );
1609          pcSlice = pcPic->getSlice               ( startCUAddrSliceIdx-1 );
1610          pcSlice->copySliceInfo                  ( pcPic->getSlice(0)      );
1611          pcSlice->setSliceIdx                    ( startCUAddrSliceIdx-1 );
1612          pcSlice->setSliceCurStartCUAddr         ( startCUAddrSlice      );
1613          pcSlice->setSliceSegmentCurStartCUAddr  ( startCUAddrSlice      );
1614          pcSlice->setSliceBits(0);
1615#if SVC_EXTENSION
1616          // copy reference list modification info from the first slice, assuming that this information is the same across all slices in the picture
1617          memcpy( pcSlice->getRefPicListModification(), pcPic->getSlice(0)->getRefPicListModification(), sizeof(TComRefPicListModification) );
1618#endif
1619          uiNumSlices ++;
1620        }
1621      }
1622      else if (pcSlice->isNextSliceSegment() || (bNoBinBitConstraintViolated && m_pcCfg->getSliceSegmentMode()==FIXED_NUMBER_OF_LCU))
1623      {
1624        startCUAddrSliceSegment                                                     = pcSlice->getSliceSegmentCurEndCUAddr();
1625        m_storedStartCUAddrForEncodingSliceSegment.push_back(startCUAddrSliceSegment);
1626        startCUAddrSliceSegmentIdx++;
1627        pcSlice->setSliceSegmentCurStartCUAddr( startCUAddrSliceSegment );
1628      }
1629      else
1630      {
1631        startCUAddrSlice                                                            = pcSlice->getSliceCurEndCUAddr();
1632        startCUAddrSliceSegment                                                     = pcSlice->getSliceSegmentCurEndCUAddr();
1633      }       
1634
1635      nextCUAddr = (startCUAddrSlice > startCUAddrSliceSegment) ? startCUAddrSlice : startCUAddrSliceSegment;
1636    }
1637    m_storedStartCUAddrForEncodingSlice.push_back( pcSlice->getSliceCurEndCUAddr());
1638    startCUAddrSliceIdx++;
1639    m_storedStartCUAddrForEncodingSliceSegment.push_back(pcSlice->getSliceCurEndCUAddr());
1640    startCUAddrSliceSegmentIdx++;
1641
1642    pcSlice = pcPic->getSlice(0);
1643
1644    // SAO parameter estimation using non-deblocked pixels for LCU bottom and right boundary areas
1645    if( m_pcCfg->getSaoLcuBasedOptimization() && m_pcCfg->getSaoLcuBoundary() )
1646    {
1647      m_pcSAO->resetStats();
1648      m_pcSAO->calcSaoStatsCu_BeforeDblk( pcPic );
1649    }
1650
1651    //-- Loop filter
1652    Bool bLFCrossTileBoundary = pcSlice->getPPS()->getLoopFilterAcrossTilesEnabledFlag();
1653    m_pcLoopFilter->setCfg(bLFCrossTileBoundary);
1654    if ( m_pcCfg->getDeblockingFilterMetric() )
1655    {
1656      dblMetric(pcPic, uiNumSlices);
1657    }
1658    m_pcLoopFilter->loopFilterPic( pcPic );
1659
1660    pcSlice = pcPic->getSlice(0);
1661    if(pcSlice->getSPS()->getUseSAO())
1662    {
1663      std::vector<Bool> LFCrossSliceBoundaryFlag;
1664      for(Int s=0; s< uiNumSlices; s++)
1665      {
1666        LFCrossSliceBoundaryFlag.push_back(  ((uiNumSlices==1)?true:pcPic->getSlice(s)->getLFCrossSliceBoundaryFlag()) );
1667      }
1668      m_storedStartCUAddrForEncodingSlice.resize(uiNumSlices+1);
1669      pcPic->createNonDBFilterInfo(m_storedStartCUAddrForEncodingSlice, 0, &LFCrossSliceBoundaryFlag ,pcPic->getPicSym()->getNumTiles() ,bLFCrossTileBoundary);
1670    }
1671
1672
1673    pcSlice = pcPic->getSlice(0);
1674
1675    if(pcSlice->getSPS()->getUseSAO())
1676    {
1677      m_pcSAO->createPicSaoInfo(pcPic);
1678    }
1679
1680    /////////////////////////////////////////////////////////////////////////////////////////////////// File writing
1681    // Set entropy coder
1682    m_pcEntropyCoder->setEntropyCoder   ( m_pcCavlcCoder, pcSlice );
1683
1684    /* write various header sets. */
1685    if ( m_bSeqFirst )
1686    {
1687#if SVC_EXTENSION
1688#if VPS_NUH_LAYER_ID
1689      OutputNALUnit nalu(NAL_UNIT_VPS, 0, 0        ); // The value of nuh_layer_id of VPS NAL unit shall be equal to 0.
1690#else
1691      OutputNALUnit nalu(NAL_UNIT_VPS, 0, m_layerId);
1692#endif
1693#if AVC_BASE
1694      if( ( m_layerId == 1 && m_pcEncTop->getVPS()->getAvcBaseLayerFlag() ) || ( m_layerId == 0 && !m_pcEncTop->getVPS()->getAvcBaseLayerFlag() ) )
1695#else
1696      if( m_layerId == 0 )
1697#endif
1698      {
1699#else
1700      OutputNALUnit nalu(NAL_UNIT_VPS);
1701#endif
1702#if VPS_EXTN_OFFSET_CALC
1703      OutputNALUnit tempNalu(NAL_UNIT_VPS, 0, 0        ); // The value of nuh_layer_id of VPS NAL unit shall be equal to 0.
1704      m_pcEntropyCoder->setBitstream(&tempNalu.m_Bitstream);
1705      m_pcEntropyCoder->encodeVPS(m_pcEncTop->getVPS());  // Use to calculate the VPS extension offset
1706#endif
1707      m_pcEntropyCoder->setBitstream(&nalu.m_Bitstream);
1708      m_pcEntropyCoder->encodeVPS(m_pcEncTop->getVPS());
1709      writeRBSPTrailingBits(nalu.m_Bitstream);
1710      accessUnit.push_back(new NALUnitEBSP(nalu));
1711#if RATE_CONTROL_LAMBDA_DOMAIN
1712      actualTotalBits += UInt(accessUnit.back()->m_nalUnitData.str().size()) * 8;
1713#endif
1714#if SVC_EXTENSION
1715      }
1716#endif
1717
1718#if SVC_EXTENSION
1719      nalu = NALUnit(NAL_UNIT_SPS, 0, m_layerId);
1720#if IL_SL_SIGNALLING_N0371
1721      pcSlice->getSPS()->setVPS( pcSlice->getVPS() );
1722#endif
1723#else
1724      nalu = NALUnit(NAL_UNIT_SPS);
1725#endif
1726      m_pcEntropyCoder->setBitstream(&nalu.m_Bitstream);
1727      if (m_bSeqFirst)
1728      {
1729        pcSlice->getSPS()->setNumLongTermRefPicSPS(m_numLongTermRefPicSPS);
1730        for (Int k = 0; k < m_numLongTermRefPicSPS; k++)
1731        {
1732          pcSlice->getSPS()->setLtRefPicPocLsbSps(k, m_ltRefPicPocLsbSps[k]);
1733          pcSlice->getSPS()->setUsedByCurrPicLtSPSFlag(k, m_ltRefPicUsedByCurrPicFlag[k]);
1734        }
1735      }
1736      if( m_pcCfg->getPictureTimingSEIEnabled() || m_pcCfg->getDecodingUnitInfoSEIEnabled() )
1737      {
1738        UInt maxCU = m_pcCfg->getSliceArgument() >> ( pcSlice->getSPS()->getMaxCUDepth() << 1);
1739        UInt numDU = ( m_pcCfg->getSliceMode() == 1 ) ? ( pcPic->getNumCUsInFrame() / maxCU ) : ( 0 );
1740        if( pcPic->getNumCUsInFrame() % maxCU != 0 || numDU == 0 )
1741        {
1742          numDU ++;
1743        }
1744        pcSlice->getSPS()->getVuiParameters()->getHrdParameters()->setNumDU( numDU );
1745        pcSlice->getSPS()->setHrdParameters( m_pcCfg->getFrameRate(), numDU, m_pcCfg->getTargetBitrate(), ( m_pcCfg->getIntraPeriod() > 0 ) );
1746      }
1747      if( m_pcCfg->getBufferingPeriodSEIEnabled() || m_pcCfg->getPictureTimingSEIEnabled() || m_pcCfg->getDecodingUnitInfoSEIEnabled() )
1748      {
1749        pcSlice->getSPS()->getVuiParameters()->setHrdParametersPresentFlag( true );
1750      }
1751      m_pcEntropyCoder->encodeSPS(pcSlice->getSPS());
1752      writeRBSPTrailingBits(nalu.m_Bitstream);
1753      accessUnit.push_back(new NALUnitEBSP(nalu));
1754#if RATE_CONTROL_LAMBDA_DOMAIN
1755      actualTotalBits += UInt(accessUnit.back()->m_nalUnitData.str().size()) * 8;
1756#endif
1757
1758#if SVC_EXTENSION
1759      nalu = NALUnit(NAL_UNIT_PPS, 0, m_layerId);
1760#else
1761      nalu = NALUnit(NAL_UNIT_PPS);
1762#endif
1763      m_pcEntropyCoder->setBitstream(&nalu.m_Bitstream);
1764      m_pcEntropyCoder->encodePPS(pcSlice->getPPS());
1765      writeRBSPTrailingBits(nalu.m_Bitstream);
1766      accessUnit.push_back(new NALUnitEBSP(nalu));
1767#if RATE_CONTROL_LAMBDA_DOMAIN
1768      actualTotalBits += UInt(accessUnit.back()->m_nalUnitData.str().size()) * 8;
1769#endif
1770
1771      xCreateLeadingSEIMessages(accessUnit, pcSlice->getSPS());
1772
1773      m_bSeqFirst = false;
1774    }
1775
1776    if (writeSOP) // write SOP description SEI (if enabled) at the beginning of GOP
1777    {
1778      Int SOPcurrPOC = pocCurr;
1779
1780      OutputNALUnit nalu(NAL_UNIT_PREFIX_SEI);
1781      m_pcEntropyCoder->setEntropyCoder(m_pcCavlcCoder, pcSlice);
1782      m_pcEntropyCoder->setBitstream(&nalu.m_Bitstream);
1783
1784      SEISOPDescription SOPDescriptionSEI;
1785      SOPDescriptionSEI.m_sopSeqParameterSetId = pcSlice->getSPS()->getSPSId();
1786
1787      UInt i = 0;
1788      UInt prevEntryId = iGOPid;
1789      for (j = iGOPid; j < m_iGopSize; j++)
1790      {
1791        Int deltaPOC = m_pcCfg->getGOPEntry(j).m_POC - m_pcCfg->getGOPEntry(prevEntryId).m_POC;
1792        if ((SOPcurrPOC + deltaPOC) < m_pcCfg->getFramesToBeEncoded())
1793        {
1794          SOPcurrPOC += deltaPOC;
1795          SOPDescriptionSEI.m_sopDescVclNaluType[i] = getNalUnitType(SOPcurrPOC, m_iLastIDR);
1796          SOPDescriptionSEI.m_sopDescTemporalId[i] = m_pcCfg->getGOPEntry(j).m_temporalId;
1797          SOPDescriptionSEI.m_sopDescStRpsIdx[i] = m_pcEncTop->getReferencePictureSetIdxForSOP(pcSlice, SOPcurrPOC, j);
1798          SOPDescriptionSEI.m_sopDescPocDelta[i] = deltaPOC;
1799
1800          prevEntryId = j;
1801          i++;
1802        }
1803      }
1804
1805      SOPDescriptionSEI.m_numPicsInSopMinus1 = i - 1;
1806
1807      m_seiWriter.writeSEImessage( nalu.m_Bitstream, SOPDescriptionSEI, pcSlice->getSPS());
1808      writeRBSPTrailingBits(nalu.m_Bitstream);
1809      accessUnit.push_back(new NALUnitEBSP(nalu));
1810
1811      writeSOP = false;
1812    }
1813
1814    if( ( m_pcCfg->getPictureTimingSEIEnabled() || m_pcCfg->getDecodingUnitInfoSEIEnabled() ) &&
1815        ( pcSlice->getSPS()->getVuiParametersPresentFlag() ) &&
1816        ( ( pcSlice->getSPS()->getVuiParameters()->getHrdParameters()->getNalHrdParametersPresentFlag() ) 
1817       || ( pcSlice->getSPS()->getVuiParameters()->getHrdParameters()->getVclHrdParametersPresentFlag() ) ) )
1818    {
1819      if( pcSlice->getSPS()->getVuiParameters()->getHrdParameters()->getSubPicCpbParamsPresentFlag() )
1820      {
1821        UInt numDU = pcSlice->getSPS()->getVuiParameters()->getHrdParameters()->getNumDU();
1822        pictureTimingSEI.m_numDecodingUnitsMinus1     = ( numDU - 1 );
1823        pictureTimingSEI.m_duCommonCpbRemovalDelayFlag = false;
1824
1825        if( pictureTimingSEI.m_numNalusInDuMinus1 == NULL )
1826        {
1827          pictureTimingSEI.m_numNalusInDuMinus1       = new UInt[ numDU ];
1828        }
1829        if( pictureTimingSEI.m_duCpbRemovalDelayMinus1  == NULL )
1830        {
1831          pictureTimingSEI.m_duCpbRemovalDelayMinus1  = new UInt[ numDU ];
1832        }
1833        if( accumBitsDU == NULL )
1834        {
1835          accumBitsDU                                  = new UInt[ numDU ];
1836        }
1837        if( accumNalsDU == NULL )
1838        {
1839          accumNalsDU                                  = new UInt[ numDU ];
1840        }
1841      }
1842      pictureTimingSEI.m_auCpbRemovalDelay = std::min<Int>(std::max<Int>(1, m_totalCoded - m_lastBPSEI), static_cast<Int>(pow(2, static_cast<double>(pcSlice->getSPS()->getVuiParameters()->getHrdParameters()->getCpbRemovalDelayLengthMinus1()+1)))); // Syntax element signalled as minus, hence the .
1843#if POC_RESET_FLAG
1844      pictureTimingSEI.m_picDpbOutputDelay = pcSlice->getSPS()->getNumReorderPics(0) + pocCurr - m_totalCoded;
1845#else
1846      pictureTimingSEI.m_picDpbOutputDelay = pcSlice->getSPS()->getNumReorderPics(0) + pcSlice->getPOC() - m_totalCoded;
1847#endif
1848      Int factor = pcSlice->getSPS()->getVuiParameters()->getHrdParameters()->getTickDivisorMinus2() + 2;
1849      pictureTimingSEI.m_picDpbOutputDuDelay = factor * pictureTimingSEI.m_picDpbOutputDelay;
1850      if( m_pcCfg->getDecodingUnitInfoSEIEnabled() )
1851      {
1852        picSptDpbOutputDuDelay = factor * pictureTimingSEI.m_picDpbOutputDelay;
1853      }
1854    }
1855
1856    if( ( m_pcCfg->getBufferingPeriodSEIEnabled() ) && ( pcSlice->getSliceType() == I_SLICE ) &&
1857        ( pcSlice->getSPS()->getVuiParametersPresentFlag() ) && 
1858        ( ( pcSlice->getSPS()->getVuiParameters()->getHrdParameters()->getNalHrdParametersPresentFlag() ) 
1859       || ( pcSlice->getSPS()->getVuiParameters()->getHrdParameters()->getVclHrdParametersPresentFlag() ) ) )
1860    {
1861      OutputNALUnit nalu(NAL_UNIT_PREFIX_SEI);
1862      m_pcEntropyCoder->setEntropyCoder(m_pcCavlcCoder, pcSlice);
1863      m_pcEntropyCoder->setBitstream(&nalu.m_Bitstream);
1864
1865      SEIBufferingPeriod sei_buffering_period;
1866     
1867      UInt uiInitialCpbRemovalDelay = (90000/2);                      // 0.5 sec
1868      sei_buffering_period.m_initialCpbRemovalDelay      [0][0]     = uiInitialCpbRemovalDelay;
1869      sei_buffering_period.m_initialCpbRemovalDelayOffset[0][0]     = uiInitialCpbRemovalDelay;
1870      sei_buffering_period.m_initialCpbRemovalDelay      [0][1]     = uiInitialCpbRemovalDelay;
1871      sei_buffering_period.m_initialCpbRemovalDelayOffset[0][1]     = uiInitialCpbRemovalDelay;
1872
1873      Double dTmp = (Double)pcSlice->getSPS()->getVuiParameters()->getTimingInfo()->getNumUnitsInTick() / (Double)pcSlice->getSPS()->getVuiParameters()->getTimingInfo()->getTimeScale();
1874
1875      UInt uiTmp = (UInt)( dTmp * 90000.0 ); 
1876      uiInitialCpbRemovalDelay -= uiTmp;
1877      uiInitialCpbRemovalDelay -= uiTmp / ( pcSlice->getSPS()->getVuiParameters()->getHrdParameters()->getTickDivisorMinus2() + 2 );
1878      sei_buffering_period.m_initialAltCpbRemovalDelay      [0][0]  = uiInitialCpbRemovalDelay;
1879      sei_buffering_period.m_initialAltCpbRemovalDelayOffset[0][0]  = uiInitialCpbRemovalDelay;
1880      sei_buffering_period.m_initialAltCpbRemovalDelay      [0][1]  = uiInitialCpbRemovalDelay;
1881      sei_buffering_period.m_initialAltCpbRemovalDelayOffset[0][1]  = uiInitialCpbRemovalDelay;
1882
1883      sei_buffering_period.m_rapCpbParamsPresentFlag              = 0;
1884      //for the concatenation, it can be set to one during splicing.
1885      sei_buffering_period.m_concatenationFlag = 0;
1886      //since the temporal layer HRD is not ready, we assumed it is fixed
1887      sei_buffering_period.m_auCpbRemovalDelayDelta = 1;
1888      sei_buffering_period.m_cpbDelayOffset = 0;
1889      sei_buffering_period.m_dpbDelayOffset = 0;
1890
1891      m_seiWriter.writeSEImessage( nalu.m_Bitstream, sei_buffering_period, pcSlice->getSPS());
1892      writeRBSPTrailingBits(nalu.m_Bitstream);
1893      {
1894      UInt seiPositionInAu = xGetFirstSeiLocation(accessUnit);
1895      UInt offsetPosition = m_activeParameterSetSEIPresentInAU;   // Insert BP SEI after APS SEI
1896      AccessUnit::iterator it;
1897      for(j = 0, it = accessUnit.begin(); j < seiPositionInAu + offsetPosition; j++)
1898      {
1899        it++;
1900      }
1901      accessUnit.insert(it, new NALUnitEBSP(nalu));
1902      m_bufferingPeriodSEIPresentInAU = true;
1903      }
1904
1905      if (m_pcCfg->getScalableNestingSEIEnabled())
1906      {
1907        OutputNALUnit naluTmp(NAL_UNIT_PREFIX_SEI);
1908        m_pcEntropyCoder->setEntropyCoder(m_pcCavlcCoder, pcSlice);
1909        m_pcEntropyCoder->setBitstream(&naluTmp.m_Bitstream);
1910        scalableNestingSEI.m_nestedSEIs.clear();
1911        scalableNestingSEI.m_nestedSEIs.push_back(&sei_buffering_period);
1912        m_seiWriter.writeSEImessage( naluTmp.m_Bitstream, scalableNestingSEI, pcSlice->getSPS());
1913        writeRBSPTrailingBits(naluTmp.m_Bitstream);
1914        UInt seiPositionInAu = xGetFirstSeiLocation(accessUnit);
1915        UInt offsetPosition = m_activeParameterSetSEIPresentInAU + m_bufferingPeriodSEIPresentInAU + m_pictureTimingSEIPresentInAU;   // Insert BP SEI after non-nested APS, BP and PT SEIs
1916        AccessUnit::iterator it;
1917        for(j = 0, it = accessUnit.begin(); j < seiPositionInAu + offsetPosition; j++)
1918        {
1919          it++;
1920        }
1921        accessUnit.insert(it, new NALUnitEBSP(naluTmp));
1922        m_nestedBufferingPeriodSEIPresentInAU = true;
1923      }
1924
1925      m_lastBPSEI = m_totalCoded;
1926      m_cpbRemovalDelay = 0;
1927    }
1928    m_cpbRemovalDelay ++;
1929    if( ( m_pcEncTop->getRecoveryPointSEIEnabled() ) && ( pcSlice->getSliceType() == I_SLICE ) )
1930    {
1931      if( m_pcEncTop->getGradualDecodingRefreshInfoEnabled() && !pcSlice->getRapPicFlag() )
1932      {
1933        // Gradual decoding refresh SEI
1934        OutputNALUnit nalu(NAL_UNIT_PREFIX_SEI);
1935        m_pcEntropyCoder->setEntropyCoder(m_pcCavlcCoder, pcSlice);
1936        m_pcEntropyCoder->setBitstream(&nalu.m_Bitstream);
1937
1938        SEIGradualDecodingRefreshInfo seiGradualDecodingRefreshInfo;
1939        seiGradualDecodingRefreshInfo.m_gdrForegroundFlag = true; // Indicating all "foreground"
1940
1941        m_seiWriter.writeSEImessage( nalu.m_Bitstream, seiGradualDecodingRefreshInfo, pcSlice->getSPS() );
1942        writeRBSPTrailingBits(nalu.m_Bitstream);
1943        accessUnit.push_back(new NALUnitEBSP(nalu));
1944      }
1945    // Recovery point SEI
1946      OutputNALUnit nalu(NAL_UNIT_PREFIX_SEI);
1947      m_pcEntropyCoder->setEntropyCoder(m_pcCavlcCoder, pcSlice);
1948      m_pcEntropyCoder->setBitstream(&nalu.m_Bitstream);
1949
1950      SEIRecoveryPoint sei_recovery_point;
1951      sei_recovery_point.m_recoveryPocCnt    = 0;
1952#if POC_RESET_FLAG
1953      sei_recovery_point.m_exactMatchingFlag = ( pocCurr == 0 ) ? (true) : (false);
1954#else
1955      sei_recovery_point.m_exactMatchingFlag = ( pcSlice->getPOC() == 0 ) ? (true) : (false);
1956#endif
1957      sei_recovery_point.m_brokenLinkFlag    = false;
1958
1959      m_seiWriter.writeSEImessage( nalu.m_Bitstream, sei_recovery_point, pcSlice->getSPS() );
1960      writeRBSPTrailingBits(nalu.m_Bitstream);
1961      accessUnit.push_back(new NALUnitEBSP(nalu));
1962    }
1963
1964    /* use the main bitstream buffer for storing the marshalled picture */
1965    m_pcEntropyCoder->setBitstream(NULL);
1966
1967    startCUAddrSliceIdx = 0;
1968    startCUAddrSlice    = 0; 
1969
1970    startCUAddrSliceSegmentIdx = 0;
1971    startCUAddrSliceSegment    = 0; 
1972    nextCUAddr                 = 0;
1973    pcSlice = pcPic->getSlice(startCUAddrSliceIdx);
1974
1975    Int processingState = (pcSlice->getSPS()->getUseSAO())?(EXECUTE_INLOOPFILTER):(ENCODE_SLICE);
1976    Bool skippedSlice=false;
1977    while (nextCUAddr < uiRealEndAddress) // Iterate over all slices
1978    {
1979      switch(processingState)
1980      {
1981      case ENCODE_SLICE:
1982        {
1983          pcSlice->setNextSlice       ( false );
1984          pcSlice->setNextSliceSegment( false );
1985          if (nextCUAddr == m_storedStartCUAddrForEncodingSlice[startCUAddrSliceIdx])
1986          {
1987            pcSlice = pcPic->getSlice(startCUAddrSliceIdx);
1988            if(startCUAddrSliceIdx > 0 && pcSlice->getSliceType()!= I_SLICE)
1989            {
1990              pcSlice->checkColRefIdx(startCUAddrSliceIdx, pcPic);
1991            }
1992            pcPic->setCurrSliceIdx(startCUAddrSliceIdx);
1993            m_pcSliceEncoder->setSliceIdx(startCUAddrSliceIdx);
1994            assert(startCUAddrSliceIdx == pcSlice->getSliceIdx());
1995            // Reconstruction slice
1996            pcSlice->setSliceCurStartCUAddr( nextCUAddr );  // to be used in encodeSlice() + context restriction
1997            pcSlice->setSliceCurEndCUAddr  ( m_storedStartCUAddrForEncodingSlice[startCUAddrSliceIdx+1 ] );
1998            // Dependent slice
1999            pcSlice->setSliceSegmentCurStartCUAddr( nextCUAddr );  // to be used in encodeSlice() + context restriction
2000            pcSlice->setSliceSegmentCurEndCUAddr  ( m_storedStartCUAddrForEncodingSliceSegment[startCUAddrSliceSegmentIdx+1 ] );
2001
2002            pcSlice->setNextSlice       ( true );
2003
2004            startCUAddrSliceIdx++;
2005            startCUAddrSliceSegmentIdx++;
2006          } 
2007          else if (nextCUAddr == m_storedStartCUAddrForEncodingSliceSegment[startCUAddrSliceSegmentIdx])
2008          {
2009            // Dependent slice
2010            pcSlice->setSliceSegmentCurStartCUAddr( nextCUAddr );  // to be used in encodeSlice() + context restriction
2011            pcSlice->setSliceSegmentCurEndCUAddr  ( m_storedStartCUAddrForEncodingSliceSegment[startCUAddrSliceSegmentIdx+1 ] );
2012
2013            pcSlice->setNextSliceSegment( true );
2014
2015            startCUAddrSliceSegmentIdx++;
2016          }
2017#if SVC_EXTENSION && M0457_COL_PICTURE_SIGNALING
2018          pcSlice->setNumMotionPredRefLayers(m_pcEncTop->getNumMotionPredRefLayers());
2019#endif
2020          pcSlice->setRPS(pcPic->getSlice(0)->getRPS());
2021          pcSlice->setRPSidx(pcPic->getSlice(0)->getRPSidx());
2022          UInt uiDummyStartCUAddr;
2023          UInt uiDummyBoundingCUAddr;
2024          m_pcSliceEncoder->xDetermineStartAndBoundingCUAddr(uiDummyStartCUAddr,uiDummyBoundingCUAddr,pcPic,true);
2025
2026          uiInternalAddress = pcPic->getPicSym()->getPicSCUAddr(pcSlice->getSliceSegmentCurEndCUAddr()-1) % pcPic->getNumPartInCU();
2027          uiExternalAddress = pcPic->getPicSym()->getPicSCUAddr(pcSlice->getSliceSegmentCurEndCUAddr()-1) / pcPic->getNumPartInCU();
2028          uiPosX = ( uiExternalAddress % pcPic->getFrameWidthInCU() ) * g_uiMaxCUWidth+ g_auiRasterToPelX[ g_auiZscanToRaster[uiInternalAddress] ];
2029          uiPosY = ( uiExternalAddress / pcPic->getFrameWidthInCU() ) * g_uiMaxCUHeight+ g_auiRasterToPelY[ g_auiZscanToRaster[uiInternalAddress] ];
2030
2031#if REPN_FORMAT_IN_VPS
2032          uiWidth = pcSlice->getPicWidthInLumaSamples();
2033          uiHeight = pcSlice->getPicHeightInLumaSamples();
2034#else
2035          uiWidth = pcSlice->getSPS()->getPicWidthInLumaSamples();
2036          uiHeight = pcSlice->getSPS()->getPicHeightInLumaSamples();
2037#endif
2038          while(uiPosX>=uiWidth||uiPosY>=uiHeight)
2039          {
2040            uiInternalAddress--;
2041            uiPosX = ( uiExternalAddress % pcPic->getFrameWidthInCU() ) * g_uiMaxCUWidth+ g_auiRasterToPelX[ g_auiZscanToRaster[uiInternalAddress] ];
2042            uiPosY = ( uiExternalAddress / pcPic->getFrameWidthInCU() ) * g_uiMaxCUHeight+ g_auiRasterToPelY[ g_auiZscanToRaster[uiInternalAddress] ];
2043          }
2044          uiInternalAddress++;
2045          if(uiInternalAddress==pcPic->getNumPartInCU())
2046          {
2047            uiInternalAddress = 0;
2048            uiExternalAddress = pcPic->getPicSym()->getCUOrderMap(pcPic->getPicSym()->getInverseCUOrderMap(uiExternalAddress)+1);
2049          }
2050          UInt endAddress = pcPic->getPicSym()->getPicSCUEncOrder(uiExternalAddress*pcPic->getNumPartInCU()+uiInternalAddress);
2051          if(endAddress<=pcSlice->getSliceSegmentCurStartCUAddr()) 
2052          {
2053            UInt boundingAddrSlice, boundingAddrSliceSegment;
2054            boundingAddrSlice          = m_storedStartCUAddrForEncodingSlice[startCUAddrSliceIdx];         
2055            boundingAddrSliceSegment = m_storedStartCUAddrForEncodingSliceSegment[startCUAddrSliceSegmentIdx];         
2056            nextCUAddr               = min(boundingAddrSlice, boundingAddrSliceSegment);
2057            if(pcSlice->isNextSlice())
2058            {
2059              skippedSlice=true;
2060            }
2061            continue;
2062          }
2063          if(skippedSlice) 
2064          {
2065            pcSlice->setNextSlice       ( true );
2066            pcSlice->setNextSliceSegment( false );
2067          }
2068          skippedSlice=false;
2069          pcSlice->allocSubstreamSizes( iNumSubstreams );
2070          for ( UInt ui = 0 ; ui < iNumSubstreams; ui++ )
2071          {
2072            pcSubstreamsOut[ui].clear();
2073          }
2074
2075          m_pcEntropyCoder->setEntropyCoder   ( m_pcCavlcCoder, pcSlice );
2076          m_pcEntropyCoder->resetEntropy      ();
2077          /* start slice NALunit */
2078#if SVC_EXTENSION
2079          OutputNALUnit nalu( pcSlice->getNalUnitType(), pcSlice->getTLayer(), m_layerId );
2080#else
2081          OutputNALUnit nalu( pcSlice->getNalUnitType(), pcSlice->getTLayer() );
2082#endif
2083          Bool sliceSegment = (!pcSlice->isNextSlice());
2084          if (!sliceSegment)
2085          {
2086            uiOneBitstreamPerSliceLength = 0; // start of a new slice
2087          }
2088          m_pcEntropyCoder->setBitstream(&nalu.m_Bitstream);
2089#if RATE_CONTROL_LAMBDA_DOMAIN
2090          tmpBitsBeforeWriting = m_pcEntropyCoder->getNumberOfWrittenBits();
2091#endif
2092
2093          m_pcEntropyCoder->encodeSliceHeader(pcSlice);
2094
2095#if RATE_CONTROL_LAMBDA_DOMAIN
2096          actualHeadBits += ( m_pcEntropyCoder->getNumberOfWrittenBits() - tmpBitsBeforeWriting );
2097#endif
2098
2099          // is it needed?
2100          {
2101            if (!sliceSegment)
2102            {
2103              pcBitstreamRedirect->writeAlignOne();
2104            }
2105            else
2106            {
2107              // We've not completed our slice header info yet, do the alignment later.
2108            }
2109            m_pcSbacCoder->init( (TEncBinIf*)m_pcBinCABAC );
2110            m_pcEntropyCoder->setEntropyCoder ( m_pcSbacCoder, pcSlice );
2111            m_pcEntropyCoder->resetEntropy    ();
2112            for ( UInt ui = 0 ; ui < pcSlice->getPPS()->getNumSubstreams() ; ui++ )
2113            {
2114              m_pcEntropyCoder->setEntropyCoder ( &pcSbacCoders[ui], pcSlice );
2115              m_pcEntropyCoder->resetEntropy    ();
2116            }
2117          }
2118
2119          if(pcSlice->isNextSlice())
2120          {
2121            // set entropy coder for writing
2122            m_pcSbacCoder->init( (TEncBinIf*)m_pcBinCABAC );
2123            {
2124              for ( UInt ui = 0 ; ui < pcSlice->getPPS()->getNumSubstreams() ; ui++ )
2125              {
2126                m_pcEntropyCoder->setEntropyCoder ( &pcSbacCoders[ui], pcSlice );
2127                m_pcEntropyCoder->resetEntropy    ();
2128              }
2129              pcSbacCoders[0].load(m_pcSbacCoder);
2130              m_pcEntropyCoder->setEntropyCoder ( &pcSbacCoders[0], pcSlice );  //ALF is written in substream #0 with CABAC coder #0 (see ALF param encoding below)
2131            }
2132            m_pcEntropyCoder->resetEntropy    ();
2133            // File writing
2134            if (!sliceSegment)
2135            {
2136              m_pcEntropyCoder->setBitstream(pcBitstreamRedirect);
2137            }
2138            else
2139            {
2140              m_pcEntropyCoder->setBitstream(&nalu.m_Bitstream);
2141            }
2142            // for now, override the TILES_DECODER setting in order to write substreams.
2143            m_pcEntropyCoder->setBitstream    ( &pcSubstreamsOut[0] );
2144
2145          }
2146          pcSlice->setFinalized(true);
2147
2148          m_pcSbacCoder->load( &pcSbacCoders[0] );
2149
2150          pcSlice->setTileOffstForMultES( uiOneBitstreamPerSliceLength );
2151            pcSlice->setTileLocationCount ( 0 );
2152          m_pcSliceEncoder->encodeSlice(pcPic, pcSubstreamsOut);
2153
2154          {
2155            // Construct the final bitstream by flushing and concatenating substreams.
2156            // The final bitstream is either nalu.m_Bitstream or pcBitstreamRedirect;
2157            UInt* puiSubstreamSizes = pcSlice->getSubstreamSizes();
2158            UInt uiTotalCodedSize = 0; // for padding calcs.
2159            UInt uiNumSubstreamsPerTile = iNumSubstreams;
2160            if (iNumSubstreams > 1)
2161            {
2162              uiNumSubstreamsPerTile /= pcPic->getPicSym()->getNumTiles();
2163            }
2164            for ( UInt ui = 0 ; ui < iNumSubstreams; ui++ )
2165            {
2166              // Flush all substreams -- this includes empty ones.
2167              // Terminating bit and flush.
2168              m_pcEntropyCoder->setEntropyCoder   ( &pcSbacCoders[ui], pcSlice );
2169              m_pcEntropyCoder->setBitstream      (  &pcSubstreamsOut[ui] );
2170              m_pcEntropyCoder->encodeTerminatingBit( 1 );
2171              m_pcEntropyCoder->encodeSliceFinish();
2172
2173              pcSubstreamsOut[ui].writeByteAlignment();   // Byte-alignment in slice_data() at end of sub-stream
2174              // Byte alignment is necessary between tiles when tiles are independent.
2175              uiTotalCodedSize += pcSubstreamsOut[ui].getNumberOfWrittenBits();
2176
2177              Bool bNextSubstreamInNewTile = ((ui+1) < iNumSubstreams)&& ((ui+1)%uiNumSubstreamsPerTile == 0);
2178              if (bNextSubstreamInNewTile)
2179              {
2180                pcSlice->setTileLocation(ui/uiNumSubstreamsPerTile, pcSlice->getTileOffstForMultES()+(uiTotalCodedSize>>3));
2181              }
2182              if (ui+1 < pcSlice->getPPS()->getNumSubstreams())
2183              {
2184                puiSubstreamSizes[ui] = pcSubstreamsOut[ui].getNumberOfWrittenBits() + (pcSubstreamsOut[ui].countStartCodeEmulations()<<3);
2185              }
2186            }
2187
2188            // Complete the slice header info.
2189            m_pcEntropyCoder->setEntropyCoder   ( m_pcCavlcCoder, pcSlice );
2190            m_pcEntropyCoder->setBitstream(&nalu.m_Bitstream);
2191            m_pcEntropyCoder->encodeTilesWPPEntryPoint( pcSlice );
2192
2193            // Substreams...
2194            TComOutputBitstream *pcOut = pcBitstreamRedirect;
2195          Int offs = 0;
2196          Int nss = pcSlice->getPPS()->getNumSubstreams();
2197          if (pcSlice->getPPS()->getEntropyCodingSyncEnabledFlag())
2198          {
2199            // 1st line present for WPP.
2200            offs = pcSlice->getSliceSegmentCurStartCUAddr()/pcSlice->getPic()->getNumPartInCU()/pcSlice->getPic()->getFrameWidthInCU();
2201            nss  = pcSlice->getNumEntryPointOffsets()+1;
2202          }
2203          for ( UInt ui = 0 ; ui < nss; ui++ )
2204          {
2205            pcOut->addSubstream(&pcSubstreamsOut[ui+offs]);
2206            }
2207          }
2208
2209          UInt boundingAddrSlice, boundingAddrSliceSegment;
2210          boundingAddrSlice        = m_storedStartCUAddrForEncodingSlice[startCUAddrSliceIdx];         
2211          boundingAddrSliceSegment = m_storedStartCUAddrForEncodingSliceSegment[startCUAddrSliceSegmentIdx];         
2212          nextCUAddr               = min(boundingAddrSlice, boundingAddrSliceSegment);
2213          // If current NALU is the first NALU of slice (containing slice header) and more NALUs exist (due to multiple dependent slices) then buffer it.
2214          // If current NALU is the last NALU of slice and a NALU was buffered, then (a) Write current NALU (b) Update an write buffered NALU at approproate location in NALU list.
2215          Bool bNALUAlignedWrittenToList    = false; // used to ensure current NALU is not written more than once to the NALU list.
2216          xAttachSliceDataToNalUnit(nalu, pcBitstreamRedirect);
2217          accessUnit.push_back(new NALUnitEBSP(nalu));
2218#if RATE_CONTROL_LAMBDA_DOMAIN
2219          actualTotalBits += UInt(accessUnit.back()->m_nalUnitData.str().size()) * 8;
2220#endif
2221          bNALUAlignedWrittenToList = true; 
2222          uiOneBitstreamPerSliceLength += nalu.m_Bitstream.getNumberOfWrittenBits(); // length of bitstream after byte-alignment
2223
2224          if (!bNALUAlignedWrittenToList)
2225          {
2226            {
2227              nalu.m_Bitstream.writeAlignZero();
2228            }
2229            accessUnit.push_back(new NALUnitEBSP(nalu));
2230            uiOneBitstreamPerSliceLength += nalu.m_Bitstream.getNumberOfWrittenBits() + 24; // length of bitstream after byte-alignment + 3 byte startcode 0x000001
2231          }
2232
2233          if( ( m_pcCfg->getPictureTimingSEIEnabled() || m_pcCfg->getDecodingUnitInfoSEIEnabled() ) &&
2234              ( pcSlice->getSPS()->getVuiParametersPresentFlag() ) &&
2235              ( ( pcSlice->getSPS()->getVuiParameters()->getHrdParameters()->getNalHrdParametersPresentFlag() ) 
2236             || ( pcSlice->getSPS()->getVuiParameters()->getHrdParameters()->getVclHrdParametersPresentFlag() ) ) &&
2237              ( pcSlice->getSPS()->getVuiParameters()->getHrdParameters()->getSubPicCpbParamsPresentFlag() ) )
2238          {
2239              UInt numNalus = 0;
2240            UInt numRBSPBytes = 0;
2241            for (AccessUnit::const_iterator it = accessUnit.begin(); it != accessUnit.end(); it++)
2242            {
2243              UInt numRBSPBytes_nal = UInt((*it)->m_nalUnitData.str().size());
2244              if ((*it)->m_nalUnitType != NAL_UNIT_PREFIX_SEI && (*it)->m_nalUnitType != NAL_UNIT_SUFFIX_SEI)
2245              {
2246                numRBSPBytes += numRBSPBytes_nal;
2247                numNalus ++;
2248              }
2249            }
2250            accumBitsDU[ pcSlice->getSliceIdx() ] = ( numRBSPBytes << 3 );
2251            accumNalsDU[ pcSlice->getSliceIdx() ] = numNalus;   // SEI not counted for bit count; hence shouldn't be counted for # of NALUs - only for consistency
2252          }
2253          processingState = ENCODE_SLICE;
2254          }
2255          break;
2256        case EXECUTE_INLOOPFILTER:
2257          {
2258            // set entropy coder for RD
2259            m_pcEntropyCoder->setEntropyCoder ( m_pcSbacCoder, pcSlice );
2260            if ( pcSlice->getSPS()->getUseSAO() )
2261            {
2262              m_pcEntropyCoder->resetEntropy();
2263              m_pcEntropyCoder->setBitstream( m_pcBitCounter );
2264              m_pcSAO->startSaoEnc(pcPic, m_pcEntropyCoder, m_pcEncTop->getRDSbacCoder(), m_pcEncTop->getRDGoOnSbacCoder());
2265              SAOParam& cSaoParam = *pcSlice->getPic()->getPicSym()->getSaoParam();
2266
2267#if SAO_CHROMA_LAMBDA
2268#if SAO_ENCODING_CHOICE
2269              m_pcSAO->SAOProcess(&cSaoParam, pcPic->getSlice(0)->getLambdaLuma(), pcPic->getSlice(0)->getLambdaChroma(), pcPic->getSlice(0)->getDepth());
2270#else
2271              m_pcSAO->SAOProcess(&cSaoParam, pcPic->getSlice(0)->getLambdaLuma(), pcPic->getSlice(0)->getLambdaChroma());
2272#endif
2273#else
2274              m_pcSAO->SAOProcess(&cSaoParam, pcPic->getSlice(0)->getLambda());
2275#endif
2276              m_pcSAO->endSaoEnc();
2277              m_pcSAO->PCMLFDisableProcess(pcPic);
2278            }
2279#if SAO_RDO
2280            m_pcEntropyCoder->setEntropyCoder ( m_pcCavlcCoder, pcSlice );
2281#endif
2282            processingState = ENCODE_SLICE;
2283
2284            for(Int s=0; s< uiNumSlices; s++)
2285            {
2286              if (pcSlice->getSPS()->getUseSAO())
2287              {
2288                pcPic->getSlice(s)->setSaoEnabledFlag((pcSlice->getPic()->getPicSym()->getSaoParam()->bSaoFlag[0]==1)?true:false);
2289              }
2290            }
2291          }
2292          break;
2293        default:
2294          {
2295            printf("Not a supported encoding state\n");
2296            assert(0);
2297            exit(-1);
2298          }
2299        }
2300      } // end iteration over slices
2301
2302      if(pcSlice->getSPS()->getUseSAO())
2303      {
2304        if(pcSlice->getSPS()->getUseSAO())
2305        {
2306          m_pcSAO->destroyPicSaoInfo();
2307        }
2308        pcPic->destroyNonDBFilterInfo();
2309      }
2310
2311      pcPic->compressMotion(); 
2312     
2313      //-- For time output for each slice
2314      Double dEncTime = (Double)(clock()-iBeforeTime) / CLOCKS_PER_SEC;
2315
2316      const Char* digestStr = NULL;
2317      if (m_pcCfg->getDecodedPictureHashSEIEnabled())
2318      {
2319        /* calculate MD5sum for entire reconstructed picture */
2320        SEIDecodedPictureHash sei_recon_picture_digest;
2321        if(m_pcCfg->getDecodedPictureHashSEIEnabled() == 1)
2322        {
2323          sei_recon_picture_digest.method = SEIDecodedPictureHash::MD5;
2324          calcMD5(*pcPic->getPicYuvRec(), sei_recon_picture_digest.digest);
2325          digestStr = digestToString(sei_recon_picture_digest.digest, 16);
2326        }
2327        else if(m_pcCfg->getDecodedPictureHashSEIEnabled() == 2)
2328        {
2329          sei_recon_picture_digest.method = SEIDecodedPictureHash::CRC;
2330          calcCRC(*pcPic->getPicYuvRec(), sei_recon_picture_digest.digest);
2331          digestStr = digestToString(sei_recon_picture_digest.digest, 2);
2332        }
2333        else if(m_pcCfg->getDecodedPictureHashSEIEnabled() == 3)
2334        {
2335          sei_recon_picture_digest.method = SEIDecodedPictureHash::CHECKSUM;
2336          calcChecksum(*pcPic->getPicYuvRec(), sei_recon_picture_digest.digest);
2337          digestStr = digestToString(sei_recon_picture_digest.digest, 4);
2338        }
2339#if SVC_EXTENSION
2340        OutputNALUnit nalu(NAL_UNIT_SUFFIX_SEI, pcSlice->getTLayer(), m_layerId);
2341#else
2342        OutputNALUnit nalu(NAL_UNIT_SUFFIX_SEI, pcSlice->getTLayer());
2343#endif
2344
2345        /* write the SEI messages */
2346        m_pcEntropyCoder->setEntropyCoder(m_pcCavlcCoder, pcSlice);
2347        m_seiWriter.writeSEImessage(nalu.m_Bitstream, sei_recon_picture_digest, pcSlice->getSPS());
2348        writeRBSPTrailingBits(nalu.m_Bitstream);
2349
2350        accessUnit.insert(accessUnit.end(), new NALUnitEBSP(nalu));
2351      }
2352      if (m_pcCfg->getTemporalLevel0IndexSEIEnabled())
2353      {
2354        SEITemporalLevel0Index sei_temporal_level0_index;
2355        if (pcSlice->getRapPicFlag())
2356        {
2357          m_tl0Idx = 0;
2358          m_rapIdx = (m_rapIdx + 1) & 0xFF;
2359        }
2360        else
2361        {
2362          m_tl0Idx = (m_tl0Idx + (pcSlice->getTLayer() ? 0 : 1)) & 0xFF;
2363        }
2364        sei_temporal_level0_index.tl0Idx = m_tl0Idx;
2365        sei_temporal_level0_index.rapIdx = m_rapIdx;
2366
2367        OutputNALUnit nalu(NAL_UNIT_PREFIX_SEI); 
2368
2369        /* write the SEI messages */
2370        m_pcEntropyCoder->setEntropyCoder(m_pcCavlcCoder, pcSlice);
2371        m_seiWriter.writeSEImessage(nalu.m_Bitstream, sei_temporal_level0_index, pcSlice->getSPS());
2372        writeRBSPTrailingBits(nalu.m_Bitstream);
2373
2374        /* insert the SEI message NALUnit before any Slice NALUnits */
2375        AccessUnit::iterator it = find_if(accessUnit.begin(), accessUnit.end(), mem_fun(&NALUnit::isSlice));
2376        accessUnit.insert(it, new NALUnitEBSP(nalu));
2377      }
2378
2379      xCalculateAddPSNR( pcPic, pcPic->getPicYuvRec(), accessUnit, dEncTime );
2380
2381      //In case of field coding, compute the interlaced PSNR for both fields
2382      if (isField && ((!pcPic->isTopField() && isTff) || (pcPic->isTopField() && !isTff)))
2383      {
2384        //get complementary top field
2385        TComPic* pcPicTop;
2386        TComList<TComPic*>::iterator   iterPic = rcListPic.begin();
2387        while ((*iterPic)->getPOC() != pcPic->getPOC()-1)
2388        {
2389          iterPic ++;
2390        }
2391        pcPicTop = *(iterPic);
2392        xCalculateInterlacedAddPSNR(pcPicTop, pcPic, pcPicTop->getPicYuvRec(), pcPic->getPicYuvRec(), accessUnit, dEncTime );
2393      }
2394
2395      if (digestStr)
2396      {
2397        if(m_pcCfg->getDecodedPictureHashSEIEnabled() == 1)
2398        {
2399          printf(" [MD5:%s]", digestStr);
2400        }
2401        else if(m_pcCfg->getDecodedPictureHashSEIEnabled() == 2)
2402        {
2403          printf(" [CRC:%s]", digestStr);
2404        }
2405        else if(m_pcCfg->getDecodedPictureHashSEIEnabled() == 3)
2406        {
2407          printf(" [Checksum:%s]", digestStr);
2408        }
2409      }
2410#if RATE_CONTROL_LAMBDA_DOMAIN
2411      if ( m_pcCfg->getUseRateCtrl() )
2412      {
2413#if !M0036_RC_IMPROVEMENT
2414        Double effectivePercentage = m_pcRateCtrl->getRCPic()->getEffectivePercentage();
2415#endif
2416        Double avgQP     = m_pcRateCtrl->getRCPic()->calAverageQP();
2417        Double avgLambda = m_pcRateCtrl->getRCPic()->calAverageLambda();
2418        if ( avgLambda < 0.0 )
2419        {
2420          avgLambda = lambda;
2421        }
2422#if M0036_RC_IMPROVEMENT
2423#if RATE_CONTROL_INTRA
2424        m_pcRateCtrl->getRCPic()->updateAfterPicture( actualHeadBits, actualTotalBits, avgQP, avgLambda, pcSlice->getSliceType());
2425#else
2426        m_pcRateCtrl->getRCPic()->updateAfterPicture( actualHeadBits, actualTotalBits, avgQP, avgLambda );
2427#endif
2428#else
2429        m_pcRateCtrl->getRCPic()->updateAfterPicture( actualHeadBits, actualTotalBits, avgQP, avgLambda, effectivePercentage );
2430#endif
2431        m_pcRateCtrl->getRCPic()->addToPictureLsit( m_pcRateCtrl->getPicList() );
2432
2433        m_pcRateCtrl->getRCSeq()->updateAfterPic( actualTotalBits );
2434        if ( pcSlice->getSliceType() != I_SLICE )
2435        {
2436          m_pcRateCtrl->getRCGOP()->updateAfterPicture( actualTotalBits );
2437        }
2438        else    // for intra picture, the estimated bits are used to update the current status in the GOP
2439        {
2440          m_pcRateCtrl->getRCGOP()->updateAfterPicture( estimatedBits );
2441        }
2442      }
2443#else
2444      if(m_pcCfg->getUseRateCtrl())
2445      {
2446        UInt  frameBits = m_vRVM_RP[m_vRVM_RP.size()-1];
2447        m_pcRateCtrl->updataRCFrameStatus((Int)frameBits, pcSlice->getSliceType());
2448      }
2449#endif
2450
2451      if( ( m_pcCfg->getPictureTimingSEIEnabled() || m_pcCfg->getDecodingUnitInfoSEIEnabled() ) &&
2452          ( pcSlice->getSPS()->getVuiParametersPresentFlag() ) &&
2453          ( ( pcSlice->getSPS()->getVuiParameters()->getHrdParameters()->getNalHrdParametersPresentFlag() ) 
2454         || ( pcSlice->getSPS()->getVuiParameters()->getHrdParameters()->getVclHrdParametersPresentFlag() ) ) )
2455      {
2456        TComVUI *vui = pcSlice->getSPS()->getVuiParameters();
2457        TComHRD *hrd = vui->getHrdParameters();
2458
2459        if( hrd->getSubPicCpbParamsPresentFlag() )
2460        {
2461          Int i;
2462          UInt64 ui64Tmp;
2463          UInt uiPrev = 0;
2464          UInt numDU = ( pictureTimingSEI.m_numDecodingUnitsMinus1 + 1 );
2465          UInt *pCRD = &pictureTimingSEI.m_duCpbRemovalDelayMinus1[0];
2466          UInt maxDiff = ( hrd->getTickDivisorMinus2() + 2 ) - 1;
2467
2468          for( i = 0; i < numDU; i ++ )
2469          {
2470            pictureTimingSEI.m_numNalusInDuMinus1[ i ]       = ( i == 0 ) ? ( accumNalsDU[ i ] - 1 ) : ( accumNalsDU[ i ] - accumNalsDU[ i - 1] - 1 );
2471          }
2472
2473          if( numDU == 1 )
2474          {
2475            pCRD[ 0 ] = 0; /* don't care */
2476          }
2477          else
2478          {
2479            pCRD[ numDU - 1 ] = 0;/* by definition */
2480            UInt tmp = 0;
2481            UInt accum = 0;
2482
2483            for( i = ( numDU - 2 ); i >= 0; i -- )
2484            {
2485              ui64Tmp = ( ( ( accumBitsDU[ numDU - 1 ]  - accumBitsDU[ i ] ) * ( vui->getTimingInfo()->getTimeScale() / vui->getTimingInfo()->getNumUnitsInTick() ) * ( hrd->getTickDivisorMinus2() + 2 ) ) / ( m_pcCfg->getTargetBitrate() ) );
2486              if( (UInt)ui64Tmp > maxDiff )
2487              {
2488                tmp ++;
2489              }
2490            }
2491            uiPrev = 0;
2492
2493            UInt flag = 0;
2494            for( i = ( numDU - 2 ); i >= 0; i -- )
2495            {
2496              flag = 0;
2497              ui64Tmp = ( ( ( accumBitsDU[ numDU - 1 ]  - accumBitsDU[ i ] ) * ( vui->getTimingInfo()->getTimeScale() / vui->getTimingInfo()->getNumUnitsInTick() ) * ( hrd->getTickDivisorMinus2() + 2 ) ) / ( m_pcCfg->getTargetBitrate() ) );
2498
2499              if( (UInt)ui64Tmp > maxDiff )
2500              {
2501                if(uiPrev >= maxDiff - tmp)
2502                {
2503                  ui64Tmp = uiPrev + 1;
2504                  flag = 1;
2505                }
2506                else                            ui64Tmp = maxDiff - tmp + 1;
2507              }
2508              pCRD[ i ] = (UInt)ui64Tmp - uiPrev - 1;
2509              if( (Int)pCRD[ i ] < 0 )
2510              {
2511                pCRD[ i ] = 0;
2512              }
2513              else if (tmp > 0 && flag == 1) 
2514              {
2515                tmp --;
2516              }
2517              accum += pCRD[ i ] + 1;
2518              uiPrev = accum;
2519            }
2520          }
2521        }
2522        if( m_pcCfg->getPictureTimingSEIEnabled() )
2523        {
2524          {
2525            OutputNALUnit nalu(NAL_UNIT_PREFIX_SEI, pcSlice->getTLayer());
2526          m_pcEntropyCoder->setEntropyCoder(m_pcCavlcCoder, pcSlice);
2527          pictureTimingSEI.m_picStruct = (isField && pcSlice->getPic()->isTopField())? 1 : isField? 2 : 0;
2528          m_seiWriter.writeSEImessage(nalu.m_Bitstream, pictureTimingSEI, pcSlice->getSPS());
2529          writeRBSPTrailingBits(nalu.m_Bitstream);
2530          UInt seiPositionInAu = xGetFirstSeiLocation(accessUnit);
2531          UInt offsetPosition = m_activeParameterSetSEIPresentInAU
2532                                    + m_bufferingPeriodSEIPresentInAU;    // Insert PT SEI after APS and BP SEI
2533          AccessUnit::iterator it;
2534          for(j = 0, it = accessUnit.begin(); j < seiPositionInAu + offsetPosition; j++)
2535          {
2536            it++;
2537          }
2538          accessUnit.insert(it, new NALUnitEBSP(nalu));
2539          m_pictureTimingSEIPresentInAU = true;
2540        }
2541          if ( m_pcCfg->getScalableNestingSEIEnabled() ) // put picture timing SEI into scalable nesting SEI
2542          {
2543            OutputNALUnit nalu(NAL_UNIT_PREFIX_SEI, pcSlice->getTLayer());
2544            m_pcEntropyCoder->setEntropyCoder(m_pcCavlcCoder, pcSlice);
2545            scalableNestingSEI.m_nestedSEIs.clear();
2546            scalableNestingSEI.m_nestedSEIs.push_back(&pictureTimingSEI);
2547            m_seiWriter.writeSEImessage(nalu.m_Bitstream, scalableNestingSEI, pcSlice->getSPS());
2548            writeRBSPTrailingBits(nalu.m_Bitstream);
2549            UInt seiPositionInAu = xGetFirstSeiLocation(accessUnit);
2550            UInt offsetPosition = m_activeParameterSetSEIPresentInAU
2551              + m_bufferingPeriodSEIPresentInAU + m_pictureTimingSEIPresentInAU + m_nestedBufferingPeriodSEIPresentInAU;    // Insert PT SEI after APS and BP SEI
2552            AccessUnit::iterator it;
2553            for(j = 0, it = accessUnit.begin(); j < seiPositionInAu + offsetPosition; j++)
2554            {
2555              it++;
2556            }
2557            accessUnit.insert(it, new NALUnitEBSP(nalu));
2558            m_nestedPictureTimingSEIPresentInAU = true;
2559          }
2560        }
2561        if( m_pcCfg->getDecodingUnitInfoSEIEnabled() && hrd->getSubPicCpbParamsPresentFlag() )
2562        {             
2563          m_pcEntropyCoder->setEntropyCoder(m_pcCavlcCoder, pcSlice);
2564          for( Int i = 0; i < ( pictureTimingSEI.m_numDecodingUnitsMinus1 + 1 ); i ++ )
2565          {
2566            OutputNALUnit nalu(NAL_UNIT_PREFIX_SEI, pcSlice->getTLayer());
2567
2568            SEIDecodingUnitInfo tempSEI;
2569            tempSEI.m_decodingUnitIdx = i;
2570            tempSEI.m_duSptCpbRemovalDelay = pictureTimingSEI.m_duCpbRemovalDelayMinus1[i] + 1;
2571            tempSEI.m_dpbOutputDuDelayPresentFlag = false;
2572            tempSEI.m_picSptDpbOutputDuDelay = picSptDpbOutputDuDelay;
2573
2574            AccessUnit::iterator it;
2575            // Insert the first one in the right location, before the first slice
2576            if(i == 0)
2577            {
2578              // Insert before the first slice.
2579              m_seiWriter.writeSEImessage(nalu.m_Bitstream, tempSEI, pcSlice->getSPS());
2580              writeRBSPTrailingBits(nalu.m_Bitstream);
2581
2582              UInt seiPositionInAu = xGetFirstSeiLocation(accessUnit);
2583              UInt offsetPosition = m_activeParameterSetSEIPresentInAU
2584                                    + m_bufferingPeriodSEIPresentInAU
2585                                    + m_pictureTimingSEIPresentInAU;  // Insert DU info SEI after APS, BP and PT SEI
2586              for(j = 0, it = accessUnit.begin(); j < seiPositionInAu + offsetPosition; j++)
2587              {
2588                it++;
2589              }
2590              accessUnit.insert(it, new NALUnitEBSP(nalu));
2591            }
2592            else
2593            {
2594              Int ctr;
2595              // For the second decoding unit onwards we know how many NALUs are present
2596              for (ctr = 0, it = accessUnit.begin(); it != accessUnit.end(); it++)
2597              {           
2598                if(ctr == accumNalsDU[ i - 1 ])
2599                {
2600                  // Insert before the first slice.
2601                  m_seiWriter.writeSEImessage(nalu.m_Bitstream, tempSEI, pcSlice->getSPS());
2602                  writeRBSPTrailingBits(nalu.m_Bitstream);
2603
2604                  accessUnit.insert(it, new NALUnitEBSP(nalu));
2605                  break;
2606                }
2607                if ((*it)->m_nalUnitType != NAL_UNIT_PREFIX_SEI && (*it)->m_nalUnitType != NAL_UNIT_SUFFIX_SEI)
2608                {
2609                  ctr++;
2610                }
2611              }
2612            }           
2613          }
2614        }
2615      }
2616      xResetNonNestedSEIPresentFlags();
2617      xResetNestedSEIPresentFlags();
2618      pcPic->getPicYuvRec()->copyToPic(pcPicYuvRecOut);
2619
2620#if M0040_ADAPTIVE_RESOLUTION_CHANGE
2621      pcPicYuvRecOut->setReconstructed(true);
2622#endif
2623
2624      pcPic->setReconMark   ( true );
2625      m_bFirst = false;
2626      m_iNumPicCoded++;
2627      m_totalCoded ++;
2628      /* logging: insert a newline at end of picture period */
2629      printf("\n");
2630      fflush(stdout);
2631
2632      delete[] pcSubstreamsOut;
2633  }
2634#if !RATE_CONTROL_LAMBDA_DOMAIN
2635  if(m_pcCfg->getUseRateCtrl())
2636  {
2637    m_pcRateCtrl->updateRCGOPStatus();
2638  }
2639#endif
2640  delete pcBitstreamRedirect;
2641
2642  if( accumBitsDU != NULL) delete accumBitsDU;
2643  if( accumNalsDU != NULL) delete accumNalsDU;
2644
2645#if SVC_EXTENSION
2646  assert ( m_iNumPicCoded <= 1 || (isField && iPOCLast == 1) );
2647#else
2648  assert ( (m_iNumPicCoded == iNumPicRcvd) || (isField && iPOCLast == 1) );
2649#endif
2650}
2651
2652#if !SVC_EXTENSION
2653Void TEncGOP::printOutSummary(UInt uiNumAllPicCoded, Bool isField)
2654{
2655  assert (uiNumAllPicCoded == m_gcAnalyzeAll.getNumPic());
2656 
2657   
2658  //--CFG_KDY
2659  if(isField)
2660  {
2661    m_gcAnalyzeAll.setFrmRate( m_pcCfg->getFrameRate() * 2);
2662    m_gcAnalyzeI.setFrmRate( m_pcCfg->getFrameRate() * 2);
2663    m_gcAnalyzeP.setFrmRate( m_pcCfg->getFrameRate() * 2);
2664    m_gcAnalyzeB.setFrmRate( m_pcCfg->getFrameRate() * 2);
2665  }
2666   else
2667  {
2668    m_gcAnalyzeAll.setFrmRate( m_pcCfg->getFrameRate() );
2669    m_gcAnalyzeI.setFrmRate( m_pcCfg->getFrameRate() );
2670    m_gcAnalyzeP.setFrmRate( m_pcCfg->getFrameRate() );
2671    m_gcAnalyzeB.setFrmRate( m_pcCfg->getFrameRate() );
2672  }
2673 
2674  //-- all
2675  printf( "\n\nSUMMARY --------------------------------------------------------\n" );
2676  m_gcAnalyzeAll.printOut('a');
2677 
2678  printf( "\n\nI Slices--------------------------------------------------------\n" );
2679  m_gcAnalyzeI.printOut('i');
2680 
2681  printf( "\n\nP Slices--------------------------------------------------------\n" );
2682  m_gcAnalyzeP.printOut('p');
2683 
2684  printf( "\n\nB Slices--------------------------------------------------------\n" );
2685  m_gcAnalyzeB.printOut('b');
2686 
2687#if _SUMMARY_OUT_
2688  m_gcAnalyzeAll.printSummaryOut();
2689#endif
2690#if _SUMMARY_PIC_
2691  m_gcAnalyzeI.printSummary('I');
2692  m_gcAnalyzeP.printSummary('P');
2693  m_gcAnalyzeB.printSummary('B');
2694#endif
2695
2696  if(isField)
2697  {
2698    //-- interlaced summary
2699    m_gcAnalyzeAll_in.setFrmRate( m_pcCfg->getFrameRate());
2700    printf( "\n\nSUMMARY INTERLACED ---------------------------------------------\n" );
2701    m_gcAnalyzeAll_in.printOutInterlaced('a',  m_gcAnalyzeAll.getBits());
2702   
2703#if _SUMMARY_OUT_
2704    m_gcAnalyzeAll_in.printSummaryOutInterlaced();
2705#endif
2706  }
2707
2708  printf("\nRVM: %.3lf\n" , xCalculateRVM());
2709}
2710#endif
2711
2712Void TEncGOP::preLoopFilterPicAll( TComPic* pcPic, UInt64& ruiDist, UInt64& ruiBits )
2713{
2714  TComSlice* pcSlice = pcPic->getSlice(pcPic->getCurrSliceIdx());
2715  Bool bCalcDist = false;
2716  m_pcLoopFilter->setCfg(m_pcCfg->getLFCrossTileBoundaryFlag());
2717  m_pcLoopFilter->loopFilterPic( pcPic );
2718 
2719  m_pcEntropyCoder->setEntropyCoder ( m_pcEncTop->getRDGoOnSbacCoder(), pcSlice );
2720  m_pcEntropyCoder->resetEntropy    ();
2721  m_pcEntropyCoder->setBitstream    ( m_pcBitCounter );
2722  pcSlice = pcPic->getSlice(0);
2723  if(pcSlice->getSPS()->getUseSAO())
2724  {
2725    std::vector<Bool> LFCrossSliceBoundaryFlag(1, true);
2726    std::vector<Int>  sliceStartAddress;
2727    sliceStartAddress.push_back(0);
2728    sliceStartAddress.push_back(pcPic->getNumCUsInFrame()* pcPic->getNumPartInCU());
2729    pcPic->createNonDBFilterInfo(sliceStartAddress, 0, &LFCrossSliceBoundaryFlag);
2730  }
2731 
2732  if( pcSlice->getSPS()->getUseSAO())
2733  {
2734    pcPic->destroyNonDBFilterInfo();
2735  }
2736 
2737  m_pcEntropyCoder->resetEntropy    ();
2738  ruiBits += m_pcEntropyCoder->getNumberOfWrittenBits();
2739 
2740  if (!bCalcDist)
2741    ruiDist = xFindDistortionFrame(pcPic->getPicYuvOrg(), pcPic->getPicYuvRec());
2742}
2743
2744// ====================================================================================================================
2745// Protected member functions
2746// ====================================================================================================================
2747
2748
2749Void TEncGOP::xInitGOP( Int iPOCLast, Int iNumPicRcvd, TComList<TComPic*>& rcListPic, TComList<TComPicYuv*>& rcListPicYuvRecOut, Bool isField )
2750{
2751  assert( iNumPicRcvd > 0 );
2752  //  Exception for the first frames
2753  if ( ( isField && (iPOCLast == 0 || iPOCLast == 1) ) || (!isField  && (iPOCLast == 0))  )
2754  {
2755    m_iGopSize    = 1;
2756  }
2757  else
2758  {
2759    m_iGopSize    = m_pcCfg->getGOPSize();
2760  }
2761  assert (m_iGopSize > 0);
2762 
2763  return;
2764}
2765
2766Void TEncGOP::xInitGOP( Int iPOCLast, Int iNumPicRcvd, TComList<TComPic*>& rcListPic, TComList<TComPicYuv*>& rcListPicYuvRecOut )
2767{
2768  assert( iNumPicRcvd > 0 );
2769  //  Exception for the first frame
2770  if ( iPOCLast == 0 )
2771  {
2772    m_iGopSize    = 1;
2773  }
2774  else
2775    m_iGopSize    = m_pcCfg->getGOPSize();
2776 
2777  assert (m_iGopSize > 0); 
2778
2779  return;
2780}
2781
2782Void TEncGOP::xGetBuffer( TComList<TComPic*>&      rcListPic,
2783                         TComList<TComPicYuv*>&    rcListPicYuvRecOut,
2784                         Int                       iNumPicRcvd,
2785                         Int                       iTimeOffset,
2786                         TComPic*&                 rpcPic,
2787                         TComPicYuv*&              rpcPicYuvRecOut,
2788                         Int                       pocCurr,
2789                         Bool                      isField)
2790{
2791  Int i;
2792  //  Rec. output
2793  TComList<TComPicYuv*>::iterator     iterPicYuvRec = rcListPicYuvRecOut.end();
2794
2795  if (isField)
2796  {
2797    for ( i = 0; i < ( (pocCurr == 0 ) || (pocCurr == 1 ) ? (iNumPicRcvd - iTimeOffset + 1) : (iNumPicRcvd - iTimeOffset + 2) ); i++ )
2798    {
2799      iterPicYuvRec--;
2800    }
2801  }
2802  else
2803  {
2804    for ( i = 0; i < (iNumPicRcvd - iTimeOffset + 1); i++ )
2805    {
2806      iterPicYuvRec--;
2807    }
2808  }
2809 
2810  if (isField)
2811  {
2812    if(pocCurr == 1)
2813    {
2814      iterPicYuvRec++;
2815    }
2816  }
2817  rpcPicYuvRecOut = *(iterPicYuvRec);
2818 
2819  //  Current pic.
2820  TComList<TComPic*>::iterator        iterPic       = rcListPic.begin();
2821  while (iterPic != rcListPic.end())
2822  {
2823    rpcPic = *(iterPic);
2824    rpcPic->setCurrSliceIdx(0);
2825    if (rpcPic->getPOC() == pocCurr)
2826    {
2827      break;
2828    }
2829    iterPic++;
2830  }
2831 
2832  assert( rpcPic != NULL );
2833  assert( rpcPic->getPOC() == pocCurr );
2834 
2835  return;
2836}
2837
2838UInt64 TEncGOP::xFindDistortionFrame (TComPicYuv* pcPic0, TComPicYuv* pcPic1)
2839{
2840  Int     x, y;
2841  Pel*  pSrc0   = pcPic0 ->getLumaAddr();
2842  Pel*  pSrc1   = pcPic1 ->getLumaAddr();
2843  UInt  uiShift = 2 * DISTORTION_PRECISION_ADJUSTMENT(g_bitDepthY-8);
2844  Int   iTemp;
2845 
2846  Int   iStride = pcPic0->getStride();
2847  Int   iWidth  = pcPic0->getWidth();
2848  Int   iHeight = pcPic0->getHeight();
2849 
2850  UInt64  uiTotalDiff = 0;
2851 
2852  for( y = 0; y < iHeight; y++ )
2853  {
2854    for( x = 0; x < iWidth; x++ )
2855    {
2856      iTemp = pSrc0[x] - pSrc1[x]; uiTotalDiff += (iTemp*iTemp) >> uiShift;
2857    }
2858    pSrc0 += iStride;
2859    pSrc1 += iStride;
2860  }
2861 
2862  uiShift = 2 * DISTORTION_PRECISION_ADJUSTMENT(g_bitDepthC-8);
2863  iHeight >>= 1;
2864  iWidth  >>= 1;
2865  iStride >>= 1;
2866 
2867  pSrc0  = pcPic0->getCbAddr();
2868  pSrc1  = pcPic1->getCbAddr();
2869 
2870  for( y = 0; y < iHeight; y++ )
2871  {
2872    for( x = 0; x < iWidth; x++ )
2873    {
2874      iTemp = pSrc0[x] - pSrc1[x]; uiTotalDiff += (iTemp*iTemp) >> uiShift;
2875    }
2876    pSrc0 += iStride;
2877    pSrc1 += iStride;
2878  }
2879 
2880  pSrc0  = pcPic0->getCrAddr();
2881  pSrc1  = pcPic1->getCrAddr();
2882 
2883  for( y = 0; y < iHeight; y++ )
2884  {
2885    for( x = 0; x < iWidth; x++ )
2886    {
2887      iTemp = pSrc0[x] - pSrc1[x]; uiTotalDiff += (iTemp*iTemp) >> uiShift;
2888    }
2889    pSrc0 += iStride;
2890    pSrc1 += iStride;
2891  }
2892 
2893  return uiTotalDiff;
2894}
2895
2896#if VERBOSE_RATE
2897static const Char* nalUnitTypeToString(NalUnitType type)
2898{
2899  switch (type)
2900  {
2901    case NAL_UNIT_CODED_SLICE_TRAIL_R:    return "TRAIL_R";
2902    case NAL_UNIT_CODED_SLICE_TRAIL_N:    return "TRAIL_N";
2903    case NAL_UNIT_CODED_SLICE_TLA_R:      return "TLA_R";
2904    case NAL_UNIT_CODED_SLICE_TSA_N:      return "TSA_N";
2905    case NAL_UNIT_CODED_SLICE_STSA_R:     return "STSA_R";
2906    case NAL_UNIT_CODED_SLICE_STSA_N:     return "STSA_N";
2907    case NAL_UNIT_CODED_SLICE_BLA_W_LP:   return "BLA_W_LP";
2908    case NAL_UNIT_CODED_SLICE_BLA_W_RADL: return "BLA_W_RADL";
2909    case NAL_UNIT_CODED_SLICE_BLA_N_LP:   return "BLA_N_LP";
2910    case NAL_UNIT_CODED_SLICE_IDR_W_RADL: return "IDR_W_RADL";
2911    case NAL_UNIT_CODED_SLICE_IDR_N_LP:   return "IDR_N_LP";
2912    case NAL_UNIT_CODED_SLICE_CRA:        return "CRA";
2913    case NAL_UNIT_CODED_SLICE_RADL_R:     return "RADL_R";
2914    case NAL_UNIT_CODED_SLICE_RASL_R:     return "RASL_R";
2915    case NAL_UNIT_VPS:                    return "VPS";
2916    case NAL_UNIT_SPS:                    return "SPS";
2917    case NAL_UNIT_PPS:                    return "PPS";
2918    case NAL_UNIT_ACCESS_UNIT_DELIMITER:  return "AUD";
2919    case NAL_UNIT_EOS:                    return "EOS";
2920    case NAL_UNIT_EOB:                    return "EOB";
2921    case NAL_UNIT_FILLER_DATA:            return "FILLER";
2922    case NAL_UNIT_PREFIX_SEI:             return "SEI";
2923    case NAL_UNIT_SUFFIX_SEI:             return "SEI";
2924    default:                              return "UNK";
2925  }
2926}
2927#endif
2928
2929Void TEncGOP::xCalculateAddPSNR( TComPic* pcPic, TComPicYuv* pcPicD, const AccessUnit& accessUnit, Double dEncTime )
2930{
2931  Int     x, y;
2932  UInt64 uiSSDY  = 0;
2933  UInt64 uiSSDU  = 0;
2934  UInt64 uiSSDV  = 0;
2935 
2936  Double  dYPSNR  = 0.0;
2937  Double  dUPSNR  = 0.0;
2938  Double  dVPSNR  = 0.0;
2939 
2940  //===== calculate PSNR =====
2941  Pel*  pOrg    = pcPic ->getPicYuvOrg()->getLumaAddr();
2942  Pel*  pRec    = pcPicD->getLumaAddr();
2943  Int   iStride = pcPicD->getStride();
2944 
2945  Int   iWidth;
2946  Int   iHeight;
2947 
2948  iWidth  = pcPicD->getWidth () - m_pcEncTop->getPad(0);
2949  iHeight = pcPicD->getHeight() - m_pcEncTop->getPad(1);
2950 
2951  Int   iSize   = iWidth*iHeight;
2952 
2953  for( y = 0; y < iHeight; y++ )
2954  {
2955    for( x = 0; x < iWidth; x++ )
2956    {
2957      Int iDiff = (Int)( pOrg[x] - pRec[x] );
2958      uiSSDY   += iDiff * iDiff;
2959    }
2960    pOrg += iStride;
2961    pRec += iStride;
2962  }
2963 
2964  iHeight >>= 1;
2965  iWidth  >>= 1;
2966  iStride >>= 1;
2967  pOrg  = pcPic ->getPicYuvOrg()->getCbAddr();
2968  pRec  = pcPicD->getCbAddr();
2969 
2970  for( y = 0; y < iHeight; y++ )
2971  {
2972    for( x = 0; x < iWidth; x++ )
2973    {
2974      Int iDiff = (Int)( pOrg[x] - pRec[x] );
2975      uiSSDU   += iDiff * iDiff;
2976    }
2977    pOrg += iStride;
2978    pRec += iStride;
2979  }
2980 
2981  pOrg  = pcPic ->getPicYuvOrg()->getCrAddr();
2982  pRec  = pcPicD->getCrAddr();
2983 
2984  for( y = 0; y < iHeight; y++ )
2985  {
2986    for( x = 0; x < iWidth; x++ )
2987    {
2988      Int iDiff = (Int)( pOrg[x] - pRec[x] );
2989      uiSSDV   += iDiff * iDiff;
2990    }
2991    pOrg += iStride;
2992    pRec += iStride;
2993  }
2994 
2995  Int maxvalY = 255 << (g_bitDepthY-8);
2996  Int maxvalC = 255 << (g_bitDepthC-8);
2997  Double fRefValueY = (Double) maxvalY * maxvalY * iSize;
2998  Double fRefValueC = (Double) maxvalC * maxvalC * iSize / 4.0;
2999  dYPSNR            = ( uiSSDY ? 10.0 * log10( fRefValueY / (Double)uiSSDY ) : 99.99 );
3000  dUPSNR            = ( uiSSDU ? 10.0 * log10( fRefValueC / (Double)uiSSDU ) : 99.99 );
3001  dVPSNR            = ( uiSSDV ? 10.0 * log10( fRefValueC / (Double)uiSSDV ) : 99.99 );
3002
3003  /* calculate the size of the access unit, excluding:
3004   *  - any AnnexB contributions (start_code_prefix, zero_byte, etc.,)
3005   *  - SEI NAL units
3006   */
3007  UInt numRBSPBytes = 0;
3008  for (AccessUnit::const_iterator it = accessUnit.begin(); it != accessUnit.end(); it++)
3009  {
3010    UInt numRBSPBytes_nal = UInt((*it)->m_nalUnitData.str().size());
3011#if VERBOSE_RATE
3012    printf("*** %6s numBytesInNALunit: %u\n", nalUnitTypeToString((*it)->m_nalUnitType), numRBSPBytes_nal);
3013#endif
3014    if ((*it)->m_nalUnitType != NAL_UNIT_PREFIX_SEI && (*it)->m_nalUnitType != NAL_UNIT_SUFFIX_SEI)
3015    {
3016      numRBSPBytes += numRBSPBytes_nal;
3017    }
3018  }
3019
3020  UInt uibits = numRBSPBytes * 8;
3021  m_vRVM_RP.push_back( uibits );
3022
3023  //===== add PSNR =====
3024#if SVC_EXTENSION
3025  m_gcAnalyzeAll[m_layerId].addResult (dYPSNR, dUPSNR, dVPSNR, (Double)uibits);
3026  TComSlice*  pcSlice = pcPic->getSlice(0);
3027  if (pcSlice->isIntra())
3028  {
3029    m_gcAnalyzeI[m_layerId].addResult (dYPSNR, dUPSNR, dVPSNR, (Double)uibits);
3030  }
3031  if (pcSlice->isInterP())
3032  {
3033    m_gcAnalyzeP[m_layerId].addResult (dYPSNR, dUPSNR, dVPSNR, (Double)uibits);
3034  }
3035  if (pcSlice->isInterB())
3036  {
3037    m_gcAnalyzeB[m_layerId].addResult (dYPSNR, dUPSNR, dVPSNR, (Double)uibits);
3038  }
3039#else
3040  m_gcAnalyzeAll.addResult (dYPSNR, dUPSNR, dVPSNR, (Double)uibits);
3041  TComSlice*  pcSlice = pcPic->getSlice(0);
3042  if (pcSlice->isIntra())
3043  {
3044    m_gcAnalyzeI.addResult (dYPSNR, dUPSNR, dVPSNR, (Double)uibits);
3045  }
3046  if (pcSlice->isInterP())
3047  {
3048    m_gcAnalyzeP.addResult (dYPSNR, dUPSNR, dVPSNR, (Double)uibits);
3049  }
3050  if (pcSlice->isInterB())
3051  {
3052    m_gcAnalyzeB.addResult (dYPSNR, dUPSNR, dVPSNR, (Double)uibits);
3053  }
3054#endif
3055
3056  Char c = (pcSlice->isIntra() ? 'I' : pcSlice->isInterP() ? 'P' : 'B');
3057  if (!pcSlice->isReferenced()) c += 32;
3058
3059#if SVC_EXTENSION
3060#if ADAPTIVE_QP_SELECTION
3061  printf("POC %4d LId: %1d TId: %1d ( %c-SLICE, nQP %d QP %d ) %10d bits",
3062         pcSlice->getPOC(),
3063         pcSlice->getLayerId(),
3064         pcSlice->getTLayer(),
3065         c,
3066         pcSlice->getSliceQpBase(),
3067         pcSlice->getSliceQp(),
3068         uibits );
3069#else
3070  printf("POC %4d LId: %1d TId: %1d ( %c-SLICE, QP %d ) %10d bits",
3071         pcSlice->getPOC()-pcSlice->getLastIDR(),
3072         pcSlice->getLayerId(),
3073         pcSlice->getTLayer(),
3074         c,
3075         pcSlice->getSliceQp(),
3076         uibits );
3077#endif
3078#else
3079#if ADAPTIVE_QP_SELECTION
3080  printf("POC %4d TId: %1d ( %c-SLICE, nQP %d QP %d ) %10d bits",
3081         pcSlice->getPOC(),
3082         pcSlice->getTLayer(),
3083         c,
3084         pcSlice->getSliceQpBase(),
3085         pcSlice->getSliceQp(),
3086         uibits );
3087#else
3088  printf("POC %4d TId: %1d ( %c-SLICE, QP %d ) %10d bits",
3089         pcSlice->getPOC()-pcSlice->getLastIDR(),
3090         pcSlice->getTLayer(),
3091         c,
3092         pcSlice->getSliceQp(),
3093         uibits );
3094#endif
3095#endif
3096
3097  printf(" [Y %6.4lf dB    U %6.4lf dB    V %6.4lf dB]", dYPSNR, dUPSNR, dVPSNR );
3098  printf(" [ET %5.0f ]", dEncTime );
3099 
3100  for (Int iRefList = 0; iRefList < 2; iRefList++)
3101  {
3102    printf(" [L%d ", iRefList);
3103    for (Int iRefIndex = 0; iRefIndex < pcSlice->getNumRefIdx(RefPicList(iRefList)); iRefIndex++)
3104    {
3105#if SVC_EXTENSION
3106#if VPS_EXTN_DIRECT_REF_LAYERS
3107      if( pcSlice->getRefPic(RefPicList(iRefList), iRefIndex)->isILR(m_layerId) )
3108      {
3109        printf( "%d(%d)", pcSlice->getRefPOC(RefPicList(iRefList), iRefIndex)-pcSlice->getLastIDR(), pcSlice->getRefPic(RefPicList(iRefList), iRefIndex)->getLayerId() );
3110      }
3111      else
3112      {
3113        printf ("%d", pcSlice->getRefPOC(RefPicList(iRefList), iRefIndex)-pcSlice->getLastIDR());
3114      }
3115#endif
3116      if( pcSlice->getEnableTMVPFlag() && iRefList == 1 - pcSlice->getColFromL0Flag() && iRefIndex == pcSlice->getColRefIdx() )
3117      {
3118        printf( "c" );
3119      }
3120
3121      printf( " " );
3122#else
3123      printf ("%d ", pcSlice->getRefPOC(RefPicList(iRefList), iRefIndex)-pcSlice->getLastIDR());
3124#endif
3125    }
3126    printf("]");
3127  }
3128}
3129
3130
3131Void reinterlace(Pel* top, Pel* bottom, Pel* dst, UInt stride, UInt width, UInt height, Bool isTff)
3132{
3133 
3134  for (Int y = 0; y < height; y++)
3135  {
3136    for (Int x = 0; x < width; x++)
3137    {
3138      dst[x] = isTff ? (UChar) top[x] : (UChar) bottom[x];
3139      dst[stride+x] = isTff ? (UChar) bottom[x] : (UChar) top[x];
3140    }
3141    top += stride;
3142    bottom += stride;
3143    dst += stride*2;
3144  }
3145}
3146
3147Void TEncGOP::xCalculateInterlacedAddPSNR( TComPic* pcPicOrgTop, TComPic* pcPicOrgBottom, TComPicYuv* pcPicRecTop, TComPicYuv* pcPicRecBottom, const AccessUnit& accessUnit, Double dEncTime )
3148{
3149  Int     x, y;
3150 
3151  UInt64 uiSSDY_in  = 0;
3152  UInt64 uiSSDU_in  = 0;
3153  UInt64 uiSSDV_in  = 0;
3154 
3155  Double  dYPSNR_in  = 0.0;
3156  Double  dUPSNR_in  = 0.0;
3157  Double  dVPSNR_in  = 0.0;
3158 
3159  /*------ INTERLACED PSNR -----------*/
3160 
3161  /* Luma */
3162 
3163  Pel*  pOrgTop = pcPicOrgTop->getPicYuvOrg()->getLumaAddr();
3164  Pel*  pOrgBottom = pcPicOrgBottom->getPicYuvOrg()->getLumaAddr();
3165  Pel*  pRecTop = pcPicRecTop->getLumaAddr();
3166  Pel*  pRecBottom = pcPicRecBottom->getLumaAddr();
3167 
3168  Int   iWidth;
3169  Int   iHeight;
3170  Int iStride;
3171 
3172  iWidth  = pcPicOrgTop->getPicYuvOrg()->getWidth () - m_pcEncTop->getPad(0);
3173  iHeight = pcPicOrgTop->getPicYuvOrg()->getHeight() - m_pcEncTop->getPad(1);
3174  iStride = pcPicOrgTop->getPicYuvOrg()->getStride();
3175  Int   iSize   = iWidth*iHeight;
3176  bool isTff = pcPicOrgTop->isTopField();
3177 
3178  TComPicYuv* pcOrgInterlaced = new TComPicYuv;
3179#if AUXILIARY_PICTURES
3180  pcOrgInterlaced->create( iWidth, iHeight << 1, pcPicOrgTop->getChromaFormat(), g_uiMaxCUWidth, g_uiMaxCUHeight, g_uiMaxCUDepth );
3181#else
3182  pcOrgInterlaced->create( iWidth, iHeight << 1, g_uiMaxCUWidth, g_uiMaxCUHeight, g_uiMaxCUDepth );
3183#endif
3184 
3185  TComPicYuv* pcRecInterlaced = new TComPicYuv;
3186#if AUXILIARY_PICTURES
3187  pcRecInterlaced->create( iWidth, iHeight << 1, pcPicOrgTop->getChromaFormat(), g_uiMaxCUWidth, g_uiMaxCUHeight, g_uiMaxCUDepth );
3188#else
3189  pcRecInterlaced->create( iWidth, iHeight << 1, g_uiMaxCUWidth, g_uiMaxCUHeight, g_uiMaxCUDepth );
3190#endif
3191 
3192  Pel* pOrgInterlaced = pcOrgInterlaced->getLumaAddr();
3193  Pel* pRecInterlaced = pcRecInterlaced->getLumaAddr();
3194 
3195  //=== Interlace fields ====
3196  reinterlace(pOrgTop, pOrgBottom, pOrgInterlaced, iStride, iWidth, iHeight, isTff);
3197  reinterlace(pRecTop, pRecBottom, pRecInterlaced, iStride, iWidth, iHeight, isTff);
3198 
3199  //===== calculate PSNR =====
3200  for( y = 0; y < iHeight << 1; y++ )
3201  {
3202    for( x = 0; x < iWidth; x++ )
3203    {
3204      Int iDiff = (Int)( pOrgInterlaced[x] - pRecInterlaced[x] );
3205      uiSSDY_in   += iDiff * iDiff;
3206    }
3207    pOrgInterlaced += iStride;
3208    pRecInterlaced += iStride;
3209  }
3210 
3211  /*Chroma*/
3212 
3213  iHeight >>= 1;
3214  iWidth  >>= 1;
3215  iStride >>= 1;
3216 
3217  pOrgTop = pcPicOrgTop->getPicYuvOrg()->getCbAddr();
3218  pOrgBottom = pcPicOrgBottom->getPicYuvOrg()->getCbAddr();
3219  pRecTop = pcPicRecTop->getCbAddr();
3220  pRecBottom = pcPicRecBottom->getCbAddr();
3221  pOrgInterlaced = pcOrgInterlaced->getCbAddr();
3222  pRecInterlaced = pcRecInterlaced->getCbAddr();
3223 
3224  //=== Interlace fields ====
3225  reinterlace(pOrgTop, pOrgBottom, pOrgInterlaced, iStride, iWidth, iHeight, isTff);
3226  reinterlace(pRecTop, pRecBottom, pRecInterlaced, iStride, iWidth, iHeight, isTff);
3227 
3228  //===== calculate PSNR =====
3229  for( y = 0; y < iHeight << 1; y++ )
3230  {
3231    for( x = 0; x < iWidth; x++ )
3232    {
3233      Int iDiff = (Int)( pOrgInterlaced[x] - pRecInterlaced[x] );
3234      uiSSDU_in   += iDiff * iDiff;
3235    }
3236    pOrgInterlaced += iStride;
3237    pRecInterlaced += iStride;
3238  }
3239 
3240  pOrgTop = pcPicOrgTop->getPicYuvOrg()->getCrAddr();
3241  pOrgBottom = pcPicOrgBottom->getPicYuvOrg()->getCrAddr();
3242  pRecTop = pcPicRecTop->getCrAddr();
3243  pRecBottom = pcPicRecBottom->getCrAddr();
3244  pOrgInterlaced = pcOrgInterlaced->getCrAddr();
3245  pRecInterlaced = pcRecInterlaced->getCrAddr();
3246 
3247  //=== Interlace fields ====
3248  reinterlace(pOrgTop, pOrgBottom, pOrgInterlaced, iStride, iWidth, iHeight, isTff);
3249  reinterlace(pRecTop, pRecBottom, pRecInterlaced, iStride, iWidth, iHeight, isTff);
3250 
3251  //===== calculate PSNR =====
3252  for( y = 0; y < iHeight << 1; y++ )
3253  {
3254    for( x = 0; x < iWidth; x++ )
3255    {
3256      Int iDiff = (Int)( pOrgInterlaced[x] - pRecInterlaced[x] );
3257      uiSSDV_in   += iDiff * iDiff;
3258    }
3259    pOrgInterlaced += iStride;
3260    pRecInterlaced += iStride;
3261  }
3262 
3263  Int maxvalY = 255 << (g_bitDepthY-8);
3264  Int maxvalC = 255 << (g_bitDepthC-8);
3265  Double fRefValueY = (Double) maxvalY * maxvalY * iSize*2;
3266  Double fRefValueC = (Double) maxvalC * maxvalC * iSize*2 / 4.0;
3267  dYPSNR_in            = ( uiSSDY_in ? 10.0 * log10( fRefValueY / (Double)uiSSDY_in ) : 99.99 );
3268  dUPSNR_in            = ( uiSSDU_in ? 10.0 * log10( fRefValueC / (Double)uiSSDU_in ) : 99.99 );
3269  dVPSNR_in            = ( uiSSDV_in ? 10.0 * log10( fRefValueC / (Double)uiSSDV_in ) : 99.99 );
3270 
3271  /* calculate the size of the access unit, excluding:
3272   *  - any AnnexB contributions (start_code_prefix, zero_byte, etc.,)
3273   *  - SEI NAL units
3274   */
3275  UInt numRBSPBytes = 0;
3276  for (AccessUnit::const_iterator it = accessUnit.begin(); it != accessUnit.end(); it++)
3277  {
3278    UInt numRBSPBytes_nal = UInt((*it)->m_nalUnitData.str().size());
3279   
3280    if ((*it)->m_nalUnitType != NAL_UNIT_PREFIX_SEI && (*it)->m_nalUnitType != NAL_UNIT_SUFFIX_SEI)
3281      numRBSPBytes += numRBSPBytes_nal;
3282  }
3283 
3284  UInt uibits = numRBSPBytes * 8 ;
3285 
3286  //===== add PSNR =====
3287  m_gcAnalyzeAll_in.addResult (dYPSNR_in, dUPSNR_in, dVPSNR_in, (Double)uibits);
3288 
3289  printf("\n                                      Interlaced frame %d: [Y %6.4lf dB    U %6.4lf dB    V %6.4lf dB]", pcPicOrgBottom->getPOC()/2 , dYPSNR_in, dUPSNR_in, dVPSNR_in );
3290 
3291  pcOrgInterlaced->destroy();
3292  delete pcOrgInterlaced;
3293  pcRecInterlaced->destroy();
3294  delete pcRecInterlaced;
3295}
3296
3297/** Function for deciding the nal_unit_type.
3298 * \param pocCurr POC of the current picture
3299 * \returns the nal unit type of the picture
3300 * This function checks the configuration and returns the appropriate nal_unit_type for the picture.
3301 */
3302NalUnitType TEncGOP::getNalUnitType(Int pocCurr, Int lastIDR)
3303{
3304  if (pocCurr == 0)
3305  {
3306    return NAL_UNIT_CODED_SLICE_IDR_W_RADL;
3307  }
3308  if (pocCurr % m_pcCfg->getIntraPeriod() == 0)
3309  {
3310    if (m_pcCfg->getDecodingRefreshType() == 1)
3311    {
3312      return NAL_UNIT_CODED_SLICE_CRA;
3313    }
3314    else if (m_pcCfg->getDecodingRefreshType() == 2)
3315    {
3316      return NAL_UNIT_CODED_SLICE_IDR_W_RADL;
3317    }
3318  }
3319  if(m_pocCRA>0)
3320  {
3321    if(pocCurr<m_pocCRA)
3322    {
3323      // All leading pictures are being marked as TFD pictures here since current encoder uses all
3324      // reference pictures while encoding leading pictures. An encoder can ensure that a leading
3325      // picture can be still decodable when random accessing to a CRA/CRANT/BLA/BLANT picture by
3326      // controlling the reference pictures used for encoding that leading picture. Such a leading
3327      // picture need not be marked as a TFD picture.
3328      return NAL_UNIT_CODED_SLICE_RASL_R;
3329    }
3330  }
3331  if (lastIDR>0)
3332  {
3333    if (pocCurr < lastIDR)
3334    {
3335      return NAL_UNIT_CODED_SLICE_RADL_R;
3336    }
3337  }
3338  return NAL_UNIT_CODED_SLICE_TRAIL_R;
3339}
3340
3341Double TEncGOP::xCalculateRVM()
3342{
3343  Double dRVM = 0;
3344 
3345  if( m_pcCfg->getGOPSize() == 1 && m_pcCfg->getIntraPeriod() != 1 && m_pcCfg->getFramesToBeEncoded() > RVM_VCEGAM10_M * 2 )
3346  {
3347    // calculate RVM only for lowdelay configurations
3348    std::vector<Double> vRL , vB;
3349    size_t N = m_vRVM_RP.size();
3350    vRL.resize( N );
3351    vB.resize( N );
3352   
3353    Int i;
3354    Double dRavg = 0 , dBavg = 0;
3355    vB[RVM_VCEGAM10_M] = 0;
3356    for( i = RVM_VCEGAM10_M + 1 ; i < N - RVM_VCEGAM10_M + 1 ; i++ )
3357    {
3358      vRL[i] = 0;
3359      for( Int j = i - RVM_VCEGAM10_M ; j <= i + RVM_VCEGAM10_M - 1 ; j++ )
3360        vRL[i] += m_vRVM_RP[j];
3361      vRL[i] /= ( 2 * RVM_VCEGAM10_M );
3362      vB[i] = vB[i-1] + m_vRVM_RP[i] - vRL[i];
3363      dRavg += m_vRVM_RP[i];
3364      dBavg += vB[i];
3365    }
3366   
3367    dRavg /= ( N - 2 * RVM_VCEGAM10_M );
3368    dBavg /= ( N - 2 * RVM_VCEGAM10_M );
3369   
3370    Double dSigamB = 0;
3371    for( i = RVM_VCEGAM10_M + 1 ; i < N - RVM_VCEGAM10_M + 1 ; i++ )
3372    {
3373      Double tmp = vB[i] - dBavg;
3374      dSigamB += tmp * tmp;
3375    }
3376    dSigamB = sqrt( dSigamB / ( N - 2 * RVM_VCEGAM10_M ) );
3377   
3378    Double f = sqrt( 12.0 * ( RVM_VCEGAM10_M - 1 ) / ( RVM_VCEGAM10_M + 1 ) );
3379   
3380    dRVM = dSigamB / dRavg * f;
3381  }
3382 
3383  return( dRVM );
3384}
3385
3386/** Attaches the input bitstream to the stream in the output NAL unit
3387    Updates rNalu to contain concatenated bitstream. rpcBitstreamRedirect is cleared at the end of this function call.
3388 *  \param codedSliceData contains the coded slice data (bitstream) to be concatenated to rNalu
3389 *  \param rNalu          target NAL unit
3390 */
3391Void TEncGOP::xAttachSliceDataToNalUnit (OutputNALUnit& rNalu, TComOutputBitstream*& codedSliceData)
3392{
3393  // Byte-align
3394  rNalu.m_Bitstream.writeByteAlignment();   // Slice header byte-alignment
3395
3396  // Perform bitstream concatenation
3397  if (codedSliceData->getNumberOfWrittenBits() > 0)
3398    {
3399    rNalu.m_Bitstream.addSubstream(codedSliceData);
3400  }
3401
3402  m_pcEntropyCoder->setBitstream(&rNalu.m_Bitstream);
3403
3404  codedSliceData->clear();
3405}
3406
3407// Function will arrange the long-term pictures in the decreasing order of poc_lsb_lt,
3408// and among the pictures with the same lsb, it arranges them in increasing delta_poc_msb_cycle_lt value
3409Void TEncGOP::arrangeLongtermPicturesInRPS(TComSlice *pcSlice, TComList<TComPic*>& rcListPic)
3410{
3411  TComReferencePictureSet *rps = pcSlice->getRPS();
3412  if(!rps->getNumberOfLongtermPictures())
3413  {
3414    return;
3415  }
3416
3417  // Arrange long-term reference pictures in the correct order of LSB and MSB,
3418  // and assign values for pocLSBLT and MSB present flag
3419  Int longtermPicsPoc[MAX_NUM_REF_PICS], longtermPicsLSB[MAX_NUM_REF_PICS], indices[MAX_NUM_REF_PICS];
3420  Int longtermPicsMSB[MAX_NUM_REF_PICS];
3421  Bool mSBPresentFlag[MAX_NUM_REF_PICS];
3422  ::memset(longtermPicsPoc, 0, sizeof(longtermPicsPoc));    // Store POC values of LTRP
3423  ::memset(longtermPicsLSB, 0, sizeof(longtermPicsLSB));    // Store POC LSB values of LTRP
3424  ::memset(longtermPicsMSB, 0, sizeof(longtermPicsMSB));    // Store POC LSB values of LTRP
3425  ::memset(indices        , 0, sizeof(indices));            // Indices to aid in tracking sorted LTRPs
3426  ::memset(mSBPresentFlag , 0, sizeof(mSBPresentFlag));     // Indicate if MSB needs to be present
3427
3428  // Get the long-term reference pictures
3429  Int offset = rps->getNumberOfNegativePictures() + rps->getNumberOfPositivePictures();
3430  Int i, ctr = 0;
3431  Int maxPicOrderCntLSB = 1 << pcSlice->getSPS()->getBitsForPOC();
3432  for(i = rps->getNumberOfPictures() - 1; i >= offset; i--, ctr++)
3433  {
3434    longtermPicsPoc[ctr] = rps->getPOC(i);                                  // LTRP POC
3435    longtermPicsLSB[ctr] = getLSB(longtermPicsPoc[ctr], maxPicOrderCntLSB); // LTRP POC LSB
3436    indices[ctr]      = i; 
3437    longtermPicsMSB[ctr] = longtermPicsPoc[ctr] - longtermPicsLSB[ctr];
3438  }
3439  Int numLongPics = rps->getNumberOfLongtermPictures();
3440  assert(ctr == numLongPics);
3441
3442  // Arrange pictures in decreasing order of MSB;
3443  for(i = 0; i < numLongPics; i++)
3444  {
3445    for(Int j = 0; j < numLongPics - 1; j++)
3446    {
3447      if(longtermPicsMSB[j] < longtermPicsMSB[j+1])
3448      {
3449        std::swap(longtermPicsPoc[j], longtermPicsPoc[j+1]);
3450        std::swap(longtermPicsLSB[j], longtermPicsLSB[j+1]);
3451        std::swap(longtermPicsMSB[j], longtermPicsMSB[j+1]);
3452        std::swap(indices[j]        , indices[j+1]        );
3453      }
3454    }
3455  }
3456
3457  for(i = 0; i < numLongPics; i++)
3458  {
3459    // Check if MSB present flag should be enabled.
3460    // Check if the buffer contains any pictures that have the same LSB.
3461    TComList<TComPic*>::iterator  iterPic = rcListPic.begin(); 
3462    TComPic*                      pcPic;
3463    while ( iterPic != rcListPic.end() )
3464    {
3465      pcPic = *iterPic;
3466      if( (getLSB(pcPic->getPOC(), maxPicOrderCntLSB) == longtermPicsLSB[i])   &&     // Same LSB
3467                                      (pcPic->getSlice(0)->isReferenced())     &&    // Reference picture
3468                                        (pcPic->getPOC() != longtermPicsPoc[i])    )  // Not the LTRP itself
3469      {
3470        mSBPresentFlag[i] = true;
3471        break;
3472      }
3473      iterPic++;     
3474    }
3475  }
3476
3477  // tempArray for usedByCurr flag
3478  Bool tempArray[MAX_NUM_REF_PICS]; ::memset(tempArray, 0, sizeof(tempArray));
3479  for(i = 0; i < numLongPics; i++)
3480  {
3481    tempArray[i] = rps->getUsed(indices[i]);
3482  }
3483  // Now write the final values;
3484  ctr = 0;
3485  Int currMSB = 0, currLSB = 0;
3486  // currPicPoc = currMSB + currLSB
3487  currLSB = getLSB(pcSlice->getPOC(), maxPicOrderCntLSB); 
3488  currMSB = pcSlice->getPOC() - currLSB;
3489
3490  for(i = rps->getNumberOfPictures() - 1; i >= offset; i--, ctr++)
3491  {
3492    rps->setPOC                   (i, longtermPicsPoc[ctr]);
3493    rps->setDeltaPOC              (i, - pcSlice->getPOC() + longtermPicsPoc[ctr]);
3494    rps->setUsed                  (i, tempArray[ctr]);
3495    rps->setPocLSBLT              (i, longtermPicsLSB[ctr]);
3496    rps->setDeltaPocMSBCycleLT    (i, (currMSB - (longtermPicsPoc[ctr] - longtermPicsLSB[ctr])) / maxPicOrderCntLSB);
3497    rps->setDeltaPocMSBPresentFlag(i, mSBPresentFlag[ctr]);     
3498
3499    assert(rps->getDeltaPocMSBCycleLT(i) >= 0);   // Non-negative value
3500  }
3501  for(i = rps->getNumberOfPictures() - 1, ctr = 1; i >= offset; i--, ctr++)
3502  {
3503    for(Int j = rps->getNumberOfPictures() - 1 - ctr; j >= offset; j--)
3504    {
3505      // Here at the encoder we know that we have set the full POC value for the LTRPs, hence we
3506      // don't have to check the MSB present flag values for this constraint.
3507      assert( rps->getPOC(i) != rps->getPOC(j) ); // If assert fails, LTRP entry repeated in RPS!!!
3508    }
3509  }
3510}
3511
3512/** Function for finding the position to insert the first of APS and non-nested BP, PT, DU info SEI messages.
3513 * \param accessUnit Access Unit of the current picture
3514 * This function finds the position to insert the first of APS and non-nested BP, PT, DU info SEI messages.
3515 */
3516Int TEncGOP::xGetFirstSeiLocation(AccessUnit &accessUnit)
3517{
3518  // Find the location of the first SEI message
3519  AccessUnit::iterator it;
3520  Int seiStartPos = 0;
3521  for(it = accessUnit.begin(); it != accessUnit.end(); it++, seiStartPos++)
3522  {
3523     if ((*it)->isSei() || (*it)->isVcl())
3524     {
3525       break;
3526     }               
3527  }
3528//  assert(it != accessUnit.end());  // Triggers with some legit configurations
3529  return seiStartPos;
3530}
3531
3532#if N0383_IL_CONSTRAINED_TILE_SETS_SEI
3533Void TEncGOP::xBuildTileSetsMap(TComPicSym* picSym)
3534{
3535  Int numCUs = picSym->getFrameWidthInCU() * picSym->getFrameHeightInCU();
3536
3537  for (Int i = 0; i < numCUs; i++)
3538  {
3539    picSym->setTileSetIdxMap(i, -1, 0, false);
3540  }
3541
3542  for (Int i = 0; i < m_pcCfg->getIlNumSetsInMessage(); i++)
3543  {
3544    TComTile* topLeftTile     = picSym->getTComTile(m_pcCfg->getTopLeftTileIndex(i));
3545    TComTile* bottomRightTile = picSym->getTComTile(m_pcCfg->getBottomRightTileIndex(i));
3546    Int tileSetLeftEdgePosInCU = topLeftTile->getRightEdgePosInCU() - topLeftTile->getTileWidth() + 1;
3547    Int tileSetRightEdgePosInCU = bottomRightTile->getRightEdgePosInCU();
3548    Int tileSetTopEdgePosInCU = topLeftTile->getBottomEdgePosInCU() - topLeftTile->getTileHeight() + 1;
3549    Int tileSetBottomEdgePosInCU = bottomRightTile->getBottomEdgePosInCU();
3550    assert(tileSetLeftEdgePosInCU < tileSetRightEdgePosInCU && tileSetTopEdgePosInCU < tileSetBottomEdgePosInCU);
3551    for (Int j = tileSetTopEdgePosInCU; j <= tileSetBottomEdgePosInCU; j++)
3552    {
3553      for (Int k = tileSetLeftEdgePosInCU; k <= tileSetRightEdgePosInCU; k++)
3554      {
3555        picSym->setTileSetIdxMap(j * picSym->getFrameWidthInCU() + k, i, m_pcCfg->getIlcIdc(i), false);
3556      }
3557    }
3558  }
3559 
3560  if (m_pcCfg->getSkippedTileSetPresentFlag())
3561  {
3562    Int skippedTileSetIdx = m_pcCfg->getIlNumSetsInMessage();
3563    for (Int i = 0; i < numCUs; i++)
3564    {
3565      if (picSym->getTileSetIdxMap(i) < 0)
3566      {
3567        picSym->setTileSetIdxMap(i, skippedTileSetIdx, 0, true);
3568      }
3569    }
3570  }
3571}
3572#endif
3573
3574Void TEncGOP::dblMetric( TComPic* pcPic, UInt uiNumSlices )
3575{
3576  TComPicYuv* pcPicYuvRec = pcPic->getPicYuvRec();
3577  Pel* Rec    = pcPicYuvRec->getLumaAddr( 0 );
3578  Pel* tempRec = Rec;
3579  Int  stride = pcPicYuvRec->getStride();
3580  UInt log2maxTB = pcPic->getSlice(0)->getSPS()->getQuadtreeTULog2MaxSize();
3581  UInt maxTBsize = (1<<log2maxTB);
3582  const UInt minBlockArtSize = 8;
3583  const UInt picWidth = pcPicYuvRec->getWidth();
3584  const UInt picHeight = pcPicYuvRec->getHeight();
3585  const UInt noCol = (picWidth>>log2maxTB);
3586  const UInt noRows = (picHeight>>log2maxTB);
3587  assert(noCol > 1);
3588  assert(noRows > 1);
3589  UInt64 *colSAD = (UInt64*)malloc(noCol*sizeof(UInt64));
3590  UInt64 *rowSAD = (UInt64*)malloc(noRows*sizeof(UInt64));
3591  UInt colIdx = 0;
3592  UInt rowIdx = 0;
3593  Pel p0, p1, p2, q0, q1, q2;
3594 
3595  Int qp = pcPic->getSlice(0)->getSliceQp();
3596  Int bitdepthScale = 1 << (g_bitDepthY-8);
3597  Int beta = TComLoopFilter::getBeta( qp ) * bitdepthScale;
3598  const Int thr2 = (beta>>2);
3599  const Int thr1 = 2*bitdepthScale;
3600  UInt a = 0;
3601 
3602  memset(colSAD, 0, noCol*sizeof(UInt64));
3603  memset(rowSAD, 0, noRows*sizeof(UInt64));
3604 
3605  if (maxTBsize > minBlockArtSize)
3606  {
3607    // Analyze vertical artifact edges
3608    for(Int c = maxTBsize; c < picWidth; c += maxTBsize)
3609    {
3610      for(Int r = 0; r < picHeight; r++)
3611      {
3612        p2 = Rec[c-3];
3613        p1 = Rec[c-2];
3614        p0 = Rec[c-1];
3615        q0 = Rec[c];
3616        q1 = Rec[c+1];
3617        q2 = Rec[c+2];
3618        a = ((abs(p2-(p1<<1)+p0)+abs(q0-(q1<<1)+q2))<<1);
3619        if ( thr1 < a && a < thr2)
3620        {
3621          colSAD[colIdx] += abs(p0 - q0);
3622        }
3623        Rec += stride;
3624      }
3625      colIdx++;
3626      Rec = tempRec;
3627    }
3628   
3629    // Analyze horizontal artifact edges
3630    for(Int r = maxTBsize; r < picHeight; r += maxTBsize)
3631    {
3632      for(Int c = 0; c < picWidth; c++)
3633      {
3634        p2 = Rec[c + (r-3)*stride];
3635        p1 = Rec[c + (r-2)*stride];
3636        p0 = Rec[c + (r-1)*stride];
3637        q0 = Rec[c + r*stride];
3638        q1 = Rec[c + (r+1)*stride];
3639        q2 = Rec[c + (r+2)*stride];
3640        a = ((abs(p2-(p1<<1)+p0)+abs(q0-(q1<<1)+q2))<<1);
3641        if (thr1 < a && a < thr2)
3642        {
3643          rowSAD[rowIdx] += abs(p0 - q0);
3644        }
3645      }
3646      rowIdx++;
3647    }
3648  }
3649 
3650  UInt64 colSADsum = 0;
3651  UInt64 rowSADsum = 0;
3652  for(Int c = 0; c < noCol-1; c++)
3653  {
3654    colSADsum += colSAD[c];
3655  }
3656  for(Int r = 0; r < noRows-1; r++)
3657  {
3658    rowSADsum += rowSAD[r];
3659  }
3660 
3661  colSADsum <<= 10;
3662  rowSADsum <<= 10;
3663  colSADsum /= (noCol-1);
3664  colSADsum /= picHeight;
3665  rowSADsum /= (noRows-1);
3666  rowSADsum /= picWidth;
3667 
3668  UInt64 avgSAD = ((colSADsum + rowSADsum)>>1);
3669  avgSAD >>= (g_bitDepthY-8);
3670 
3671  if ( avgSAD > 2048 )
3672  {
3673    avgSAD >>= 9;
3674    Int offset = Clip3(2,6,(Int)avgSAD);
3675    for (Int i=0; i<uiNumSlices; i++)
3676    {
3677      pcPic->getSlice(i)->setDeblockingFilterOverrideFlag(true);
3678      pcPic->getSlice(i)->setDeblockingFilterDisable(false);
3679      pcPic->getSlice(i)->setDeblockingFilterBetaOffsetDiv2( offset );
3680      pcPic->getSlice(i)->setDeblockingFilterTcOffsetDiv2( offset );
3681    }
3682  }
3683  else
3684  {
3685    for (Int i=0; i<uiNumSlices; i++)
3686    {
3687      pcPic->getSlice(i)->setDeblockingFilterOverrideFlag(false);
3688      pcPic->getSlice(i)->setDeblockingFilterDisable(        pcPic->getSlice(i)->getPPS()->getPicDisableDeblockingFilterFlag() );
3689      pcPic->getSlice(i)->setDeblockingFilterBetaOffsetDiv2( pcPic->getSlice(i)->getPPS()->getDeblockingFilterBetaOffsetDiv2() );
3690      pcPic->getSlice(i)->setDeblockingFilterTcOffsetDiv2(   pcPic->getSlice(i)->getPPS()->getDeblockingFilterTcOffsetDiv2()   );
3691    }
3692  }
3693 
3694  free(colSAD);
3695  free(rowSAD);
3696}
3697
3698#if M0457_COL_PICTURE_SIGNALING && !REMOVE_COL_PICTURE_SIGNALING
3699TComPic* TEncGOP::getMotionPredIlp(TComSlice* pcSlice)
3700{
3701  TComPic* ilpPic = NULL;
3702  Int activeMotionPredReflayerIdx = 0;
3703
3704  for( Int i = 0; i < pcSlice->getActiveNumILRRefIdx(); i++ )
3705  {
3706    UInt refLayerIdc = pcSlice->getInterLayerPredLayerIdc(i);
3707    if( m_pcEncTop->getMotionPredEnabledFlag( pcSlice->getVPS()->getRefLayerId( m_layerId, refLayerIdc ) ) )
3708    {
3709      if (activeMotionPredReflayerIdx == pcSlice->getColRefLayerIdx())
3710      {
3711        ilpPic = m_pcEncTop->getIlpList()[refLayerIdc];
3712        break;
3713      }
3714      else
3715      {
3716        activeMotionPredReflayerIdx++;
3717      }
3718    }
3719  }
3720
3721  assert(ilpPic != NULL);
3722
3723  return ilpPic;
3724}
3725#endif
3726
3727//! \}
Note: See TracBrowser for help on using the repository browser.