rfreeimage 0.1.0

data/vendor/FreeImage/Source/LibJXR/image/encode/strenc.c

@@ -0,0 +1,2370 @@
+//*@@@+++@@@@******************************************************************
+//
+// Copyright � Microsoft Corp.
+// All rights reserved.
+// 
+// Redistribution and use in source and binary forms, with or without
+// modification, are permitted provided that the following conditions are met:
+// 
+// � Redistributions of source code must retain the above copyright notice,
+//   this list of conditions and the following disclaimer.
+// � Redistributions in binary form must reproduce the above copyright notice,
+//   this list of conditions and the following disclaimer in the documentation
+//   and/or other materials provided with the distribution.
+// 
+// THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
+// AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
+// IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
+// ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE
+// LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
+// CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
+// SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
+// INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
+// CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
+// ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
+// POSSIBILITY OF SUCH DAMAGE.
+//
+//*@@@---@@@@******************************************************************
+#include "strcodec.h"
+#include "encode.h"
+#include "strTransform.h"
+#include <math.h>
+#include "perfTimer.h"
+
+#ifdef MEM_TRACE
+#define TRACE_MALLOC    1
+#define TRACE_NEW       0
+#define TRACE_HEAP      0
+#include "memtrace.h"
+#endif
+
+#ifdef ADI_SYS_OPT
+extern char L1WW[];
+#endif
+
+#ifdef X86OPT_INLINE
+#define _FORCEINLINE __forceinline
+#else // X86OPT_INLINE
+#define _FORCEINLINE
+#endif // X86OPT_INLINE
+
+Int inputMBRow(CWMImageStrCodec *);
+
+#if defined(WMP_OPT_SSE2) || defined(WMP_OPT_CC_ENC) || defined(WMP_OPT_TRFM_ENC)
+void StrEncOpt(CWMImageStrCodec* pSC);
+#endif // OPT defined
+
+#define MINIMUM_PACKET_LENGTH 4  // as long as packet header - skipped if data is not accessed (happens only for flexbits)
+
+Void writeQuantizer(CWMIQuantizer * pQuantizer[MAX_CHANNELS], BitIOInfo * pIO, U8 cChMode, size_t cChannel, size_t iPos)
+{
+    if(cChMode > 2)
+        cChMode = 2;
+
+    if(cChannel > 1)
+        putBit16(pIO, cChMode, 2); // Channel mode
+    else
+        cChMode = 0;
+
+    putBit16(pIO, pQuantizer[0][iPos].iIndex, 8); // Y
+
+    if(cChMode == 1)  // MIXED
+        putBit16(pIO, pQuantizer[1][iPos].iIndex, 8); // UV
+    else if(cChMode > 0){ // INDEPENDENT
+        size_t i;
+
+        for(i = 1; i < cChannel; i ++)
+            putBit16(pIO, pQuantizer[i][iPos].iIndex, 8); // UV
+    }
+}
+
+// packet header: 00000000 00000000 00000001 ?????xxx
+// xxx:           000(spatial) 001(DC) 010(AD) 011(AC) 100(FL) 101-111(reserved)
+// ?????:         (iTileY * cNumOfSliceV + iTileX)
+Void writePacketHeader(BitIOInfo * pIO, U8 ptPacketType, U8 pID)
+{
+    putBit16(pIO, 0, 8);
+    putBit16(pIO, 0, 8);
+    putBit16(pIO, 1, 8);
+    putBit16(pIO, (pID << 3) + (ptPacketType & 7), 8);
+}
+
+Int writeTileHeaderDC(CWMImageStrCodec * pSC, BitIOInfo * pIO)
+{
+    size_t iTile, j = (pSC->m_pNextSC == NULL ? 1U : 2U);
+
+    for(; j > 0; j --){
+        if((pSC->m_param.uQPMode & 1) != 0){ // not DC uniform
+            CWMITile * pTile = pSC->pTile + pSC->cTileColumn;
+            size_t i;
+            
+            pTile->cChModeDC = (U8)(rand() & 3); // channel mode, just for concept proofing!
+            
+            if(pSC->cTileRow + pSC->cTileColumn == 0) // allocate DC QP info
+                for(iTile = 0; iTile <= pSC->WMISCP.cNumOfSliceMinus1V; iTile ++)
+                    if(allocateQuantizer(pSC->pTile[iTile].pQuantizerDC, pSC->m_param.cNumChannels, 1) != ICERR_OK)
+                        return ICERR_ERROR;
+            
+            for(i = 0; i < pSC->m_param.cNumChannels; i ++)
+                pTile->pQuantizerDC[i]->iIndex = (U8)((rand() & 0x2f) + 1); // QP indexes, just for concept proofing!
+            
+            formatQuantizer(pTile->pQuantizerDC, pTile->cChModeDC, pSC->m_param.cNumChannels, 0, TRUE, pSC->m_param.bScaledArith);
+
+            for(i = 0; i < pSC->m_param.cNumChannels; i ++)
+                pTile->pQuantizerDC[i]->iOffset = (pTile->pQuantizerDC[i]->iQP >> 1);
+            
+            writeQuantizer(pTile->pQuantizerDC, pIO, pTile->cChModeDC, pSC->m_param.cNumChannels, 0);
+        }
+
+        pSC = pSC->m_pNextSC;
+    }
+
+    return ICERR_OK;
+}
+
+Int writeTileHeaderLP(CWMImageStrCodec * pSC, BitIOInfo * pIO)
+{
+    size_t k = (pSC->m_pNextSC == NULL ? 1U : 2U);
+    
+    for(; k > 0; k --){
+        if(pSC->WMISCP.sbSubband != SB_DC_ONLY && (pSC->m_param.uQPMode & 2) != 0){ // not LP uniform
+            CWMITile * pTile = pSC->pTile + pSC->cTileColumn;
+            U8 i, j;
+
+            pTile->bUseDC = ((rand() & 1) == 0 ? TRUE : FALSE); // use DC quantizer?
+            putBit16(pIO, pTile->bUseDC == TRUE ? 1 : 0, 1);
+            pTile->cBitsLP = 0;
+            
+            pTile->cNumQPLP = (pTile->bUseDC == TRUE ? 1 : (U8)((rand() & 0xf) + 1)); // # of LP QPs
+            
+            if(pSC->cTileRow > 0)
+                freeQuantizer(pTile->pQuantizerLP);
+            
+            if(allocateQuantizer(pTile->pQuantizerLP, pSC->m_param.cNumChannels, pTile->cNumQPLP) != ICERR_OK)
+                return ICERR_ERROR;
+
+            if(pTile->bUseDC == TRUE)
+                useDCQuantizer(pSC, pSC->cTileColumn);
+            else{
+                putBit16(pIO, pTile->cNumQPLP - 1, 4);
+                
+                pTile->cBitsLP = dquantBits(pTile->cNumQPLP);
+                
+                for(i = 0; i < pTile->cNumQPLP; i ++){
+                    pTile->cChModeLP[i] = (U8)(rand() & 3); // channel mode, just for concept proofing!
+                    
+                    for(j = 0; j < pSC->m_param.cNumChannels; j ++)
+                        pTile->pQuantizerLP[j][i].iIndex = (U8)((rand() & 0xfe) + 1); // QP indexes, just for concept proofing!
+                    formatQuantizer(pTile->pQuantizerLP, pTile->cChModeLP[i], pSC->m_param.cNumChannels, i, TRUE, pSC->m_param.bScaledArith);
+                    writeQuantizer(pTile->pQuantizerLP, pIO, pTile->cChModeLP[i], pSC->m_param.cNumChannels, i);
+                }
+            }
+        }
+        pSC = pSC->m_pNextSC;
+    }
+
+    return ICERR_OK;
+}
+
+Int writeTileHeaderHP(CWMImageStrCodec * pSC, BitIOInfo * pIO)
+{
+    size_t k = (pSC->m_pNextSC == NULL ? 1U : 2U);
+    
+    for(; k > 0; k --){
+        if(pSC->WMISCP.sbSubband != SB_DC_ONLY && pSC->WMISCP.sbSubband != SB_NO_HIGHPASS && (pSC->m_param.uQPMode & 4) != 0){ // not HP uniform
+            CWMITile * pTile = pSC->pTile + pSC->cTileColumn;
+            U8 i, j;
+
+            pTile->bUseLP = ((rand() & 1) == 0 ? TRUE : FALSE); // use LP quantizer?
+            putBit16(pIO, pTile->bUseLP == TRUE ? 1 : 0, 1);
+            pTile->cBitsHP = 0;
+            
+            pTile->cNumQPHP = (pTile->bUseLP == TRUE ? pTile->cNumQPLP : (U8)((rand() & 0xf) + 1)); // # of LP QPs
+            
+            if(pSC->cTileRow > 0)
+                freeQuantizer(pTile->pQuantizerHP);
+            
+            if(allocateQuantizer(pTile->pQuantizerHP, pSC->m_param.cNumChannels, pTile->cNumQPHP) != ICERR_OK)
+                return ICERR_ERROR;
+            
+            if(pTile->bUseLP == TRUE)
+                useLPQuantizer(pSC, pTile->cNumQPHP, pSC->cTileColumn);
+            else{
+                putBit16(pIO, pTile->cNumQPHP - 1, 4);
+                pTile->cBitsHP = dquantBits(pTile->cNumQPHP);
+                
+                for(i = 0; i < pTile->cNumQPHP; i ++){
+                    pTile->cChModeHP[i] = (U8)(rand() & 3); // channel mode, just for concept proofing!
+                    
+                    for(j = 0; j < pSC->m_param.cNumChannels; j ++)
+                        pTile->pQuantizerHP[j][i].iIndex = (U8)((rand() & 0xfe) + 1); // QP indexes, just for concept proofing!
+                    formatQuantizer(pTile->pQuantizerHP, pTile->cChModeHP[i], pSC->m_param.cNumChannels, i, FALSE, pSC->m_param.bScaledArith);
+                    writeQuantizer(pTile->pQuantizerHP, pIO, pTile->cChModeHP[i], pSC->m_param.cNumChannels, i);
+                }
+            }
+        }
+        pSC = pSC->m_pNextSC;
+    }
+
+    return ICERR_OK;
+}
+
+Int encodeMB(CWMImageStrCodec * pSC, Int iMBX, Int iMBY)
+{
+    CCodingContext * pContext = &pSC->m_pCodingContext[pSC->cTileColumn];
+    
+    if(pSC->m_bCtxLeft && pSC->m_bCtxTop && pSC->m_bSecondary == FALSE && pSC->m_param.bTranscode == FALSE){ // write packet headers
+        U8 pID = (U8)((pSC->cTileRow * (pSC->WMISCP.cNumOfSliceMinus1V + 1) + pSC->cTileColumn) & 0x1F);
+        
+        if(pSC->WMISCP.bfBitstreamFormat == SPATIAL) {
+            writePacketHeader(pContext->m_pIODC, 0, pID);
+            if (pSC->m_param.bTrimFlexbitsFlag)
+                putBit16(pContext->m_pIODC, pContext->m_iTrimFlexBits, 4);
+            writeTileHeaderDC(pSC, pContext->m_pIODC);
+            writeTileHeaderLP(pSC, pContext->m_pIODC);
+            writeTileHeaderHP(pSC, pContext->m_pIODC);
+        }
+        else{
+            writePacketHeader(pContext->m_pIODC, 1, pID);
+            writeTileHeaderDC(pSC, pContext->m_pIODC);
+            if(pSC->cSB > 1){
+                writePacketHeader(pContext->m_pIOLP, 2, pID);
+                writeTileHeaderLP(pSC, pContext->m_pIOLP);
+            }
+            if(pSC->cSB > 2){
+                writePacketHeader(pContext->m_pIOAC, 3, pID);
+                writeTileHeaderHP(pSC, pContext->m_pIOAC);
+            }
+            if(pSC->cSB > 3) {
+                writePacketHeader(pContext->m_pIOFL, 4, pID);
+                if (pSC->m_param.bTrimFlexbitsFlag)
+                    putBit16(pContext->m_pIOFL, pContext->m_iTrimFlexBits, 4);
+            }
+        }
+    }
+    
+    if(EncodeMacroblockDC(pSC, pContext, iMBX, iMBY) != ICERR_OK)
+        return ICERR_ERROR;
+    
+    if(pSC->WMISCP.sbSubband != SB_DC_ONLY)
+        if(EncodeMacroblockLowpass(pSC, pContext, iMBX, iMBY) != ICERR_OK)
+            return ICERR_ERROR;
+
+    if(pSC->WMISCP.sbSubband != SB_DC_ONLY && pSC->WMISCP.sbSubband != SB_NO_HIGHPASS)
+        if(EncodeMacroblockHighpass(pSC, pContext, iMBX, iMBY) != ICERR_OK)
+            return ICERR_ERROR;
+    
+    if(iMBX + 1 == (int) pSC->cmbWidth && (iMBY + 1 == (int) pSC->cmbHeight || 
+        (pSC->cTileRow < pSC->WMISCP.cNumOfSliceMinus1H && iMBY == (int) pSC->WMISCP.uiTileY[pSC->cTileRow + 1] - 1)))
+    { // end of a horizontal slice
+        size_t k, l;
+
+        // get sizes of each packet and update index table
+        if (pSC->m_pNextSC == NULL || pSC->m_bSecondary) {
+            for(k = 0; k < pSC->cNumBitIO; k ++){
+                fillToByte(pSC->m_ppBitIO[k]);
+                pSC->ppWStream[k]->GetPos(pSC->ppWStream[k], &l);
+                pSC->pIndexTable[pSC->cNumBitIO * pSC->cTileRow + k] = l + getSizeWrite(pSC->m_ppBitIO[k]); // offset
+            }
+        }
+        
+        // reset coding contexts
+        if(iMBY + 1 != (int) pSC->cmbHeight){
+            for(k = 0; k <= pSC->WMISCP.cNumOfSliceMinus1V; k ++)
+                ResetCodingContextEnc(&pSC->m_pCodingContext[k]);
+        }
+    }
+
+    return ICERR_OK;
+}
+
+/*************************************************************************
+    Top level function for processing a macroblock worth of input
+*************************************************************************/
+Int processMacroblock(CWMImageStrCodec *pSC)
+{
+    Bool topORleft = (pSC->cColumn == 0 || pSC->cRow == 0);
+    ERR_CODE result = ICERR_OK;
+    size_t j, jend = (pSC->m_pNextSC != NULL);
+
+    for (j = 0; j <= jend; j++) {
+        transformMacroblock(pSC);
+        if(!topORleft){
+            getTilePos(pSC, (Int)pSC->cColumn - 1, (Int)pSC->cRow - 1);
+            if(jend){
+                pSC->m_pNextSC->cTileRow = pSC->cTileRow;
+                pSC->m_pNextSC->cTileColumn = pSC->cTileColumn;
+            }
+            if ((result = encodeMB(pSC, (Int)pSC->cColumn - 1, (Int)pSC->cRow - 1)) != ICERR_OK)
+                return result;
+        }
+
+        if (jend) {
+            pSC->m_pNextSC->cRow = pSC->cRow;
+            pSC->m_pNextSC->cColumn = pSC->cColumn;
+            pSC = pSC->m_pNextSC;
+        }
+    }
+
+    return ICERR_OK;
+}
+
+/*************************************************************************
+  forwardRGBE: forward conversion from RGBE to RGB
+*************************************************************************/
+static _FORCEINLINE PixelI forwardRGBE (PixelI RGB, PixelI E)
+{
+    PixelI iResult = 0, iAppend = 1;
+
+    if (E == 0)
+        return 0;
+
+    assert (E!=0);
+
+    E--;
+    while (((RGB & 0x80) == 0) && (E > 0)) {
+        RGB = (RGB << 1) + iAppend;
+        iAppend = 0;
+        E--;    
+    }
+
+    // result will always be one of 3 cases
+    // E  RGB       convert to
+    // 0  [0.x]      [0   x]
+    // 0  [1.x]      [1   x]
+    // e  [1.x]      [e+1 x]
+    if (E == 0) {
+        iResult = RGB;
+    }
+    else {
+        E++;
+        iResult = (RGB & 0x7f) + (E << 7);
+    }
+
+    return iResult;
+}
+
+/*************************************************************************
+  convert float-32 into float with (c, lm)!!
+*************************************************************************/
+static _FORCEINLINE PixelI float2pixel (float f, const char _c, const unsigned char _lm)
+{
+    union uif
+    {
+        I32   i;
+        float f;
+    } x;
+
+    PixelI _h, e, e1, m, s;
+
+    if (f == 0)
+    {
+        _h = 0;
+    }
+    else
+    {
+        x.f = f;
+
+        e = (x.i >> 23) & 0x000000ff;//here set e as e, not s! e includes s: [s e] 9 bits [31..23]
+        m = (x.i & 0x007fffff) | 0x800000; // actual mantissa, with normalizer
+        if (e == 0) { // denormal-land
+            m ^= 0x800000;  // actual mantissa, removing normalizer
+            e++; // actual exponent -126
+        }
+
+        e1 = e - 127 + _c;  // this is basically a division or quantization to a different exponent
+                            // note: _c cannot be greater than 127, so e1 cannot be greater than e
+        //assert (_c <= 127);
+        if (e1 <= 1) {  // denormal-land
+            if (e1 < 1)
+                m >>= (1 - e1);  // shift mantissa right to make exponent 1
+            e1 = 1;
+            if ((m & 0x800000) == 0) // if denormal, set e1 to zero else to 1
+                e1 = 0;
+        }
+        m &= 0x007fffff;
+
+        //for float-22:	    
+        _h = (e1 << _lm) + ((m + (1 << (23 - _lm - 1))) >> (23 - _lm));//take 23-bit m, shift (23-lm), get lm-bit m for float22
+        s = ((PixelI) x.i) >> 31;
+        //padding to int-32: 
+        _h = (_h ^ s) - s;	
+    }
+
+    return _h;
+}
+
+/*************************************************************************
+  convert Half-16 to internal format, only need to handle sign bit
+*************************************************************************/
+static _FORCEINLINE PixelI forwardHalf (PixelI hHalf)
+{
+    PixelI s;
+    s = hHalf >> 31;
+    hHalf = ((hHalf & 0x7fff) ^ s) - s;
+    return hHalf;
+}
+
+
+//================================================================
+// Color Conversion 
+// functions to get image data from input buffer
+// this inlcudes necessary color conversion and boundary padding
+//================================================================
+#define _CC(r, g, b) (b -= r, r += ((b + 1) >> 1) - g, g += ((r + 0) >> 1))
+#define _CC_CMYK(c, m, y, k) (y -= c, c += ((y + 1) >> 1) - m, m += (c >> 1) - k, k += ((m + 1) >> 1))
+
+//================================================================
+// BitIOInfo init/term for encoding
+const size_t MAX_MEMORY_SIZE_IN_WORDS = 64 << 20; // 1 << 20 \approx 1 million
+
+Int StrIOEncInit(CWMImageStrCodec* pSC)
+{
+    pSC->m_param.bIndexTable = !(pSC->WMISCP.bfBitstreamFormat == SPATIAL && pSC->WMISCP.cNumOfSliceMinus1H + pSC->WMISCP.cNumOfSliceMinus1V == 0);
+    if(allocateBitIOInfo(pSC) != ICERR_OK){
+        return ICERR_ERROR;
+    }
+
+    attachISWrite(pSC->pIOHeader, pSC->WMISCP.pWStream);
+
+    if(pSC->cNumBitIO > 0){
+        size_t i;
+#if defined(_WINDOWS_) || defined(UNDER_CE)  // tmpnam does not exist in VS2005 WinCE CRT
+        TCHAR szPath[MAX_PATH];
+        DWORD cSize, j, k;
+#endif
+        char * pFilename;
+
+        pSC->ppWStream = (struct WMPStream **)malloc(pSC->cNumBitIO * sizeof(struct WMPStream *));
+        if(pSC->ppWStream == NULL) return ICERR_ERROR;
+        memset(pSC->ppWStream, 0, pSC->cNumBitIO * sizeof(struct WMPStream *));
+
+        if (pSC->cmbHeight * pSC->cmbWidth * pSC->WMISCP.cChannel >= MAX_MEMORY_SIZE_IN_WORDS) {
+#ifdef _WINDOWS_
+            pSC->ppTempFile = (TCHAR **)malloc(pSC->cNumBitIO * sizeof(TCHAR *));
+            if(pSC->ppTempFile == NULL) return ICERR_ERROR;
+            memset(pSC->ppTempFile, 0, pSC->cNumBitIO * sizeof(TCHAR *)); 
+#else
+            pSC->ppTempFile = (char **)malloc(pSC->cNumBitIO * sizeof(char *));
+            if(pSC->ppTempFile == NULL) return ICERR_ERROR;
+            memset(pSC->ppTempFile, 0, pSC->cNumBitIO * sizeof(char *));
+#endif
+        }
+
+        for(i = 0; i < pSC->cNumBitIO; i ++){
+            if (pSC->cmbHeight * pSC->cmbWidth * pSC->WMISCP.cChannel >= MAX_MEMORY_SIZE_IN_WORDS) {
+#if defined(_WINDOWS_) || defined(UNDER_CE)  // tmpnam does not exist in VS2005 WinCE CRT              
+                Bool bUnicode = sizeof(TCHAR) == 2;
+                pSC->ppTempFile[i] = (TCHAR *)malloc(MAX_PATH * sizeof(TCHAR));
+                if(pSC->ppTempFile[i] == NULL) return ICERR_ERROR;
+
+                pFilename = (char *)pSC->ppTempFile[i];
+
+                cSize = GetTempPath(MAX_PATH, szPath);
+                if(cSize == 0 || cSize >= MAX_PATH)
+                    return ICERR_ERROR;
+                if(!GetTempFileName(szPath, TEXT("wdp"), 0, pSC->ppTempFile[i]))
+                    return ICERR_ERROR;
+
+                if(bUnicode){ // unicode file name
+                    for(k = j = cSize = 0; cSize < MAX_PATH; cSize ++, j += 2){
+                        if(pSC->ppTempFile[i][cSize] == '\0')
+                            break;
+                        if(pFilename[j] != '\0')
+                            pFilename[k ++] = pFilename[j];
+                        if(pFilename[j + 1] != '\0')
+                            pFilename[k ++] = pFilename[j + 1];
+                    }
+                    pFilename[cSize] = '\0';
+                }
+
+#else //DPK needs to support ANSI 
+                pSC->ppTempFile[i] = (char *)malloc(FILENAME_MAX * sizeof(char));
+                if(pSC->ppTempFile[i] == NULL) return ICERR_ERROR;
+
+                if ((pFilename = tmpnam(NULL)) == NULL)
+                    return ICERR_ERROR;                
+                strcpy(pSC->ppTempFile[i], pFilename);
+#endif
+                if(CreateWS_File(pSC->ppWStream + i, pFilename, "w+b") != ICERR_OK) return ICERR_ERROR;                
+
+            }
+            else {
+                if(CreateWS_List(pSC->ppWStream + i) != ICERR_OK) return ICERR_ERROR;
+            }
+            attachISWrite(pSC->m_ppBitIO[i], pSC->ppWStream[i]);
+        }
+    }
+
+    return ICERR_OK;
+}
+
+#define PUTBITS putBit16
+/*************************************************************************
+    Write variable length byte aligned integer
+*************************************************************************/
+static Void PutVLWordEsc(BitIOInfo* pIO, Int iEscape, size_t s)
+{
+    if (iEscape) {
+        assert(iEscape <= 0xff && iEscape > 0xfc); // fd,fe,ff are the only valid escapes
+        PUTBITS(pIO, iEscape, 8);
+    }
+    else if (s < 0xfb00) {
+        PUTBITS(pIO, (U32) s, 16);
+    }
+    else {
+        size_t t = s >> 16;
+        if ((t >> 16) == 0) {
+            PUTBITS(pIO, 0xfb, 8);
+        }
+        else {
+            t >>= 16;
+            PUTBITS(pIO, 0xfc, 8);
+            PUTBITS(pIO, (U32)(t >> 16) & 0xffff, 16);
+            PUTBITS(pIO, (U32) t & 0xffff, 16);
+        }
+        PUTBITS(pIO, (U32) t & 0xffff, 16);
+        PUTBITS(pIO, (U32) s & 0xffff, 16);
+    }
+}
+
+/*************************************************************************
+    Write index table at start (null index table)
+*************************************************************************/
+Int writeIndexTableNull(CWMImageStrCodec * pSC)
+{
+    if(pSC->cNumBitIO == 0){
+        BitIOInfo* pIO = pSC->pIOHeader;
+        fillToByte(pIO);
+
+        /* Profile / Level info */
+        PutVLWordEsc(pIO, 0, 4);    // 4 bytes
+        PUTBITS(pIO, 111, 8);       // default profile idc
+        PUTBITS(pIO, 255, 8);       // default level idc
+        PUTBITS(pIO, 1, 16);        // LAST_FLAG
+    }
+
+    return ICERR_OK;
+}
+
+/*************************************************************************
+    Write index table
+*************************************************************************/
+Int writeIndexTable(CWMImageStrCodec * pSC)
+{
+    if(pSC->cNumBitIO > 0){
+        BitIOInfo* pIO = pSC->pIOHeader;
+        size_t *pTable = pSC->pIndexTable, iSize[4] = { 0 };
+        I32 iEntry = (I32)pSC->cNumBitIO * (pSC->WMISCP.cNumOfSliceMinus1H + 1), i, k, l;
+        
+        // write index table header [0x0001] - 2 bytes
+        PUTBITS(pIO, 1, 16);
+
+        for(i = pSC->WMISCP.cNumOfSliceMinus1H; i>= 0 && pSC->bTileExtraction == FALSE; i --){
+            for(k = 0; k < (int)pSC->cNumBitIO; ){
+                for(l = 0; l < (pSC->WMISCP.bfBitstreamFormat == FREQUENCY && pSC->WMISCP.bProgressiveMode ? pSC->cSB : 1); l ++, k ++)
+                {
+                if (i > 0)
+                pTable[pSC->cNumBitIO * i + k] -= pSC->pIndexTable[pSC->cNumBitIO * (i - 1) + k]; // packet length
+                iSize[l] += pTable[pSC->cNumBitIO * i + k];
+                }
+            }
+        }
+
+        iSize[3] = iSize[2] + iSize[1] + iSize[0];
+        iSize[2] = iSize[1] + iSize[0];
+        iSize[1] = iSize[0];
+        iSize[0] = 0;
+
+        for(i = 0; i < iEntry; ){
+        for(l = 0; l < (pSC->WMISCP.bfBitstreamFormat == FREQUENCY && pSC->WMISCP.bProgressiveMode ? pSC->cSB : 1); l ++, i ++)
+        {
+            writeIS_L1(pSC, pIO);
+            PutVLWordEsc(pIO, (pTable[i] <= MINIMUM_PACKET_LENGTH) ? 0xff : 0, iSize[l]);
+            iSize[l] += (pTable[i] <= MINIMUM_PACKET_LENGTH) ? 0 : pTable[i];
+        }
+        }
+
+        writeIS_L1(pSC, pIO);
+        PutVLWordEsc(pIO, 0xff, 0); // escape to end
+        fillToByte(pIO);
+    }
+
+    return ICERR_OK;
+}
+
+Int copyTo(struct WMPStream * pSrc, struct WMPStream * pDst, size_t iBytes)
+{
+    char pData[PACKETLENGTH];
+
+    if (iBytes <= MINIMUM_PACKET_LENGTH){
+        pSrc->Read(pSrc, pData, iBytes);
+        return ICERR_OK;
+    }
+
+    while(iBytes > PACKETLENGTH){
+        pSrc->Read(pSrc, pData, PACKETLENGTH);
+        pDst->Write(pDst, pData, PACKETLENGTH);
+        iBytes -= PACKETLENGTH;
+    }
+    pSrc->Read(pSrc, pData, iBytes);
+    pDst->Write(pDst, pData, iBytes);
+
+    return ICERR_OK;
+}
+
+Int StrIOEncTerm(CWMImageStrCodec* pSC)
+{
+    BitIOInfo * pIO = pSC->pIOHeader;
+
+    fillToByte(pIO);
+
+    if(pSC->WMISCP.bVerbose){
+        U32 i, j;
+
+        printf("\n%d horizontal tiles:\n", pSC->WMISCP.cNumOfSliceMinus1H + 1);
+        for(i = 0; i <= pSC->WMISCP.cNumOfSliceMinus1H; i ++){
+            printf("    offset of tile %d in MBs: %d\n", i, pSC->WMISCP.uiTileY[i]);
+        }
+
+        printf("\n%d vertical tiles:\n", pSC->WMISCP.cNumOfSliceMinus1V + 1);
+        for(i = 0; i <= pSC->WMISCP.cNumOfSliceMinus1V; i ++){
+            printf("    offset of tile %d in MBs: %d\n", i, pSC->WMISCP.uiTileX[i]);
+        }
+
+        if(pSC->WMISCP.bfBitstreamFormat == SPATIAL){
+            printf("\nSpatial order bitstream\n");
+        }
+        else{
+            printf("\nFrequency order bitstream\n");
+        }
+
+        if(!pSC->m_param.bIndexTable){
+            printf("\nstreaming mode, no index table.\n");
+        }
+        else if(pSC->WMISCP.bfBitstreamFormat == SPATIAL){
+            for(j = 0; j <= pSC->WMISCP.cNumOfSliceMinus1H; j ++){
+                for(i = 0; i <= pSC->WMISCP.cNumOfSliceMinus1V; i ++){
+                    printf("bitstream size for tile (%d, %d): %d.\n", j, i, (int) pSC->pIndexTable[j * (pSC->WMISCP.cNumOfSliceMinus1V + 1) + i]);
+                }
+            }
+        }
+        else{
+            for(j = 0; j <= pSC->WMISCP.cNumOfSliceMinus1H; j ++){
+                for(i = 0; i <= pSC->WMISCP.cNumOfSliceMinus1V; i ++){
+                    size_t * p = &pSC->pIndexTable[(j * (pSC->WMISCP.cNumOfSliceMinus1V + 1) + i) * 4];
+                    printf("bitstream size of (DC, LP, AC, FL) for tile (%d, %d): %d %d %d %d.\n", j, i,
+                        (int) p[0], (int) p[1], (int) p[2], (int) p[3]);
+                }
+            }
+        }
+    }
+    
+    writeIndexTable(pSC); // write index table to the header
+
+    detachISWrite(pSC, pIO);
+
+    if(pSC->cNumBitIO > 0){
+        size_t i, j, k, l;
+        struct WMPStream * pDst = pSC->WMISCP.pWStream;
+        size_t * pTable = pSC->pIndexTable;
+
+        for(i = 0; i < pSC->cNumBitIO; i ++){
+            detachISWrite(pSC, pSC->m_ppBitIO[i]);
+        }
+
+        for(i = 0; i < pSC->cNumBitIO; i ++){
+            pSC->ppWStream[i]->SetPos(pSC->ppWStream[i], 0); // seek back for read
+        }
+
+        for(l = 0; l < (size_t)(pSC->WMISCP.bfBitstreamFormat == FREQUENCY && pSC->WMISCP.bProgressiveMode ? pSC->cSB : 1); l ++){
+			for(i = 0, k = l; i <= pSC->WMISCP.cNumOfSliceMinus1H; i ++){ // loop through tiles
+				for(j = 0; j <= pSC->WMISCP.cNumOfSliceMinus1V; j ++){
+
+					if(pSC->WMISCP.bfBitstreamFormat == SPATIAL)
+						copyTo(pSC->ppWStream[j], pDst, pTable[k ++]);
+					else if (!pSC->WMISCP.bProgressiveMode){
+						copyTo(pSC->ppWStream[j * pSC->cSB + 0], pDst, pTable[k ++]);
+						if(pSC->cSB > 1)
+							copyTo(pSC->ppWStream[j * pSC->cSB + 1], pDst, pTable[k ++]);
+						if(pSC->cSB > 2)
+							copyTo(pSC->ppWStream[j * pSC->cSB + 2], pDst, pTable[k ++]);
+						if(pSC->cSB > 3)
+							copyTo(pSC->ppWStream[j * pSC->cSB + 3], pDst, pTable[k ++]);
+					}
+					else{
+						copyTo(pSC->ppWStream[j * pSC->cSB + l], pDst, pTable[k]);
+						k += pSC->cSB;
+					}
+				}
+			}
+        }
+
+        if (pSC->cmbHeight * pSC->cmbWidth * pSC->WMISCP.cChannel >= MAX_MEMORY_SIZE_IN_WORDS){           
+            for(i = 0; i < pSC->cNumBitIO; i ++){
+                if(pSC->ppWStream && pSC->ppWStream[i]){
+                    if((*(pSC->ppWStream + i))->state.file.pFile){
+                        fclose((*(pSC->ppWStream + i))->state.file.pFile);
+#ifdef _WINDOWS_
+                        if(DeleteFileA((LPCSTR)pSC->ppTempFile[i]) == 0)
+                            return ICERR_ERROR;
+#else
+                        if (remove(pSC->ppTempFile[i]) == -1)
+                            return ICERR_ERROR;
+#endif
+                    }
+
+                    if (*(pSC->ppWStream + i))
+                        free(*(pSC->ppWStream + i));
+                }
+                if(pSC->ppTempFile){
+                    if(pSC->ppTempFile[i])
+                        free(pSC->ppTempFile[i]);
+                }
+            }
+
+            if(pSC->ppTempFile)
+                free(pSC->ppTempFile);
+        }
+        else{
+            for(i = 0; i < pSC->cNumBitIO; i ++){
+                if(pSC->ppWStream && pSC->ppWStream[i])
+                    pSC->ppWStream[i]->Close(pSC->ppWStream + i);
+            }
+        }
+
+        free(pSC->ppWStream);
+
+        free(pSC->m_ppBitIO);
+        free(pSC->pIndexTable);
+    }
+
+    return 0;
+}
+
+/*************************************************************************
+    Write header of image plane
+*************************************************************************/
+Int WriteImagePlaneHeader(CWMImageStrCodec * pSC)
+{
+    CWMImageInfo * pII = &pSC->WMII;
+    CWMIStrCodecParam * pSCP = &pSC->WMISCP;
+    BitIOInfo* pIO = pSC->pIOHeader;
+
+    PUTBITS(pIO, (Int) pSC->m_param.cfColorFormat, 3); // internal color format
+    PUTBITS(pIO, (Int) pSC->m_param.bScaledArith, 1); // lossless mode
+
+// subbands
+    PUTBITS(pIO, (U32)pSCP->sbSubband, 4);
+
+// color parameters
+    switch (pSC->m_param.cfColorFormat) {
+        case YUV_420:
+        case YUV_422:
+        case YUV_444:
+            PUTBITS(pIO, 0, 4);
+            PUTBITS(pIO, 0, 4);
+            break;
+        case NCOMPONENT:
+            PUTBITS(pIO, (Int) pSC->m_param.cNumChannels - 1, 4);
+            PUTBITS(pIO, 0, 4);
+            break;
+        default:
+            break;
+    }
+
+// float and 32s additional parameters
+    switch (pII->bdBitDepth) {
+        case BD_16:
+        case BD_16S:
+            PUTBITS(pIO, pSCP->nLenMantissaOrShift, 8);
+            break;
+        case BD_32:
+        case BD_32S:
+            if(pSCP->nLenMantissaOrShift == 0)
+                pSCP->nLenMantissaOrShift = 10;//default
+            PUTBITS(pIO, pSCP->nLenMantissaOrShift, 8);
+            break;
+        case BD_32F:
+            if(pSCP->nLenMantissaOrShift == 0)
+                pSCP->nLenMantissaOrShift = 13;//default
+            PUTBITS(pIO, pSCP->nLenMantissaOrShift, 8);//float conversion parameters
+            PUTBITS(pIO, pSCP->nExpBias, 8);
+            break;
+        default:
+            break;
+    }
+
+        // quantization
+    PUTBITS(pIO, (pSC->m_param.uQPMode & 1) == 1 ? 0 : 1, 1); // DC frame uniform quantization?
+    if((pSC->m_param.uQPMode & 1) == 0)
+        writeQuantizer(pSC->pTile[0].pQuantizerDC, pIO, (pSC->m_param.uQPMode >> 3) & 3, pSC->m_param.cNumChannels, 0);
+    if(pSC->WMISCP.sbSubband != SB_DC_ONLY){
+        PUTBITS(pIO, (pSC->m_param.uQPMode & 0x200) == 0 ? 1 : 0, 1); // use DC quantization?
+        if((pSC->m_param.uQPMode & 0x200) != 0){
+            PUTBITS(pIO, (pSC->m_param.uQPMode & 2) == 2 ? 0 : 1, 1); // LP frame uniform quantization?
+            if((pSC->m_param.uQPMode & 2) == 0)
+                writeQuantizer(pSC->pTile[0].pQuantizerLP, pIO, (pSC->m_param.uQPMode >> 5) & 3,  pSC->m_param.cNumChannels, 0);
+        }
+
+        if(pSC->WMISCP.sbSubband != SB_NO_HIGHPASS){
+            PUTBITS(pIO, (pSC->m_param.uQPMode & 0x400) == 0 ? 1 : 0, 1); // use LP quantization?
+            if((pSC->m_param.uQPMode & 0x400) != 0){
+                PUTBITS(pIO, (pSC->m_param.uQPMode & 4) == 4 ? 0 : 1, 1); // HP frame uniform quantization?
+                if((pSC->m_param.uQPMode & 4) == 0)
+                    writeQuantizer(pSC->pTile[0].pQuantizerHP, pIO, (pSC->m_param.uQPMode >> 7) & 3,  pSC->m_param.cNumChannels, 0);
+            }
+        }
+    }
+
+    fillToByte(pIO);  // remove this later
+    return ICERR_OK;
+}
+
+/*************************************************************************
+    Write header to buffer
+*************************************************************************/
+Int WriteWMIHeader(CWMImageStrCodec * pSC)
+{
+    CWMImageInfo * pII = &pSC->WMII;
+    CWMIStrCodecParam * pSCP = &pSC->WMISCP;
+    CCoreParameters * pCoreParam = &pSC->m_param;
+    BitIOInfo* pIO = pSC->pIOHeader;
+    U32 /*iSizeOfSize = 2,*/ i;
+    // temporary assignments / reserved words
+    // const Int HEADERSIZE = 0;
+    Bool bInscribed = FALSE;
+    Bool bAbbreviatedHeader = (((pII->cWidth + 15) / 16 > 255 || (pII->cHeight + 15) / 16 > 255) ? FALSE : TRUE);
+
+    if(pCoreParam->bTranscode == FALSE)
+        pCoreParam->cExtraPixelsTop = pCoreParam->cExtraPixelsLeft = pCoreParam->cExtraPixelsRight = pCoreParam->cExtraPixelsBottom = 0;
+
+    // num of extra boundary pixels due to compressed domain processing
+    bInscribed = (pCoreParam->cExtraPixelsTop || pCoreParam->cExtraPixelsLeft || pCoreParam->cExtraPixelsBottom || pCoreParam->cExtraPixelsRight);
+
+// 0
+    /** signature **/
+    for (i = 0; i < 8; PUTBITS(pSC->pIOHeader, gGDISignature[i++], 8));
+
+// 8
+    /** codec version and subversion **/
+    PUTBITS(pIO, CODEC_VERSION, 4);  // this should be changed to "profile" in RTM
+    if (pSC->WMISCP.bUseHardTileBoundaries)
+        PUTBITS(pIO, CODEC_SUBVERSION_NEWSCALING_HARD_TILES, 4);
+    else
+        PUTBITS(pIO, CODEC_SUBVERSION_NEWSCALING_SOFT_TILES, 4);
+
+// 9 primary parameters
+    PUTBITS(pIO, (pSCP->cNumOfSliceMinus1V || pSCP->cNumOfSliceMinus1H) ? 1 : 0, 1); // tiling present
+    PUTBITS(pIO, (Int) pSCP->bfBitstreamFormat, 1); // bitstream layout
+    PUTBITS(pIO, pII->oOrientation, 3);        // m_iRotateFlip
+    PUTBITS(pIO, pSC->m_param.bIndexTable, 1); // index table present
+    PUTBITS(pIO, pSCP->olOverlap, 2); // overlap
+
+// 10
+    PUTBITS(pIO, bAbbreviatedHeader, 1); // short words for size and tiles
+    PUTBITS(pIO, 1, 1); // long word length (use intelligence later)
+    PUTBITS(pIO, bInscribed, 1); // windowing
+    PUTBITS(pIO, pSC->m_param.bTrimFlexbitsFlag, 1); // trim flexbits flag sent
+    PUTBITS(pIO, 0, 1); // tile stretching parameters (not enabled)
+    PUTBITS(pIO, 0, 2); // reserved bits
+    PUTBITS(pIO, (Int) pSC->m_param.bAlphaChannel, 1); // alpha channel present
+
+// 11 - informational
+    PUTBITS(pIO, (Int) pII->cfColorFormat, 4); // source color format
+    if(BD_1 == pII->bdBitDepth && pSCP->bBlackWhite)
+        PUTBITS(pIO, (Int) BD_1alt, 4); // source bit depth
+    else 
+        PUTBITS(pIO, (Int) pII->bdBitDepth, 4); // source bit depth
+
+// 12 - Variable length fields
+// size
+    putBit32(pIO, (U32)(pII->cWidth - 1), bAbbreviatedHeader ? 16 : 32);
+    putBit32(pIO, (U32)(pII->cHeight - 1), bAbbreviatedHeader ? 16 : 32);
+
+// tiling
+    if (pSCP->cNumOfSliceMinus1V || pSCP->cNumOfSliceMinus1H) {
+        PUTBITS(pIO, pSCP->cNumOfSliceMinus1V, LOG_MAX_TILES); // # of vertical slices
+        PUTBITS(pIO, pSCP->cNumOfSliceMinus1H, LOG_MAX_TILES); // # of horizontal slices
+    }
+
+// tile sizes
+    for(i = 0; i < pSCP->cNumOfSliceMinus1V; i ++){ // width in MB of vertical slices, not needed for last slice!
+        PUTBITS(pIO, pSCP->uiTileX[i + 1] - pSCP->uiTileX[i], bAbbreviatedHeader ? 8 : 16);
+    }
+    for(i = 0; i < pSCP->cNumOfSliceMinus1H; i ++){ // width in MB of horizontal slices, not needed for last slice!
+        PUTBITS(pIO, pSCP->uiTileY[i + 1] - pSCP->uiTileY[i], bAbbreviatedHeader ? 8 : 16);
+    }
+
+// window due to compressed domain processing
+    if (bInscribed) {
+        PUTBITS(pIO, (U32)pCoreParam->cExtraPixelsTop, 6);
+        PUTBITS(pIO, (U32)pCoreParam->cExtraPixelsLeft, 6);
+        PUTBITS(pIO, (U32)pCoreParam->cExtraPixelsBottom, 6);
+        PUTBITS(pIO, (U32)pCoreParam->cExtraPixelsRight, 6);
+    }    
+    fillToByte(pIO);  // redundant
+
+    // write image plane headers
+    WriteImagePlaneHeader(pSC);
+
+    return ICERR_OK;
+}
+
+// streaming codec init/term
+Int StrEncInit(CWMImageStrCodec* pSC)
+{
+    COLORFORMAT cf = pSC->m_param.cfColorFormat;
+    COLORFORMAT cfE = pSC->WMII.cfColorFormat;
+    U16 iQPIndexY = 0, iQPIndexYLP = 0, iQPIndexYHP = 0;
+	U16 iQPIndexU = 0, iQPIndexULP = 0, iQPIndexUHP = 0;
+    U16 iQPIndexV = 0, iQPIndexVLP = 0, iQPIndexVHP = 0; 
+    size_t i;
+    Bool b32bit = sizeof(size_t) == 4;
+
+    /** color transcoding with resolution change **/
+    pSC->m_bUVResolutionChange = (((cfE == CF_RGB || cfE == YUV_444 || cfE == CMYK || cfE == CF_RGBE) && 
+								   (cf == YUV_422 || cf == YUV_420))
+								  || (cfE == YUV_422 && cf == YUV_420)) && !pSC->WMISCP.bYUVData;
+
+    if(pSC->m_bUVResolutionChange){
+        size_t cSize = ((cfE == YUV_422 ? 128 : 256) + (cf == YUV_420 ? 32 : 0)) * pSC->cmbWidth + 256;
+
+        if(b32bit){ // integer overlow/underflow check for 32-bit system
+            if(((pSC->cmbWidth >> 16) * ((cfE == YUV_422 ? 128 : 256) + (cf == YUV_420 ? 32 : 0))) & 0xffff0000)
+                return ICERR_ERROR;
+            if(cSize >= 0x3fffffff)
+                return ICERR_ERROR;
+        }
+        pSC->pResU = (PixelI *)malloc(cSize * sizeof(PixelI));
+        pSC->pResV = (PixelI *)malloc(cSize * sizeof(PixelI));
+        if(pSC->pResU == NULL || pSC->pResV == NULL){
+            return ICERR_ERROR;
+        }
+    }
+
+    pSC->cTileColumn = pSC->cTileRow = 0;
+
+    if(allocateTileInfo(pSC) != ICERR_OK)
+        return ICERR_ERROR;
+
+    if(pSC->m_param.bTranscode == FALSE){
+        pSC->m_param.uQPMode = 0x150;   // 101010 000
+                                        // 000    == uniform (not per tile) DC, LP, HP
+                                        // 101010 == cChMode == 2 == independent (not same) DC, LP, HP
+
+        /** lossless or Y component lossless condition: all subbands present, uniform quantization with QPIndex 1 **/
+        pSC->m_param.bScaledArith = !((pSC->m_param.uQPMode & 7) == 0 && 
+									  1 == pSC->WMISCP.uiDefaultQPIndex <= 1 && 
+									  pSC->WMISCP.sbSubband == SB_ALL && 
+									  pSC->m_bUVResolutionChange == FALSE) &&
+                                     !pSC->WMISCP.bUnscaledArith;
+        if (BD_32 == pSC->WMII.bdBitDepth || BD_32S == pSC->WMII.bdBitDepth || BD_32F == pSC->WMII.bdBitDepth) {
+            pSC->m_param.bScaledArith = FALSE;
+        }
+        pSC->m_param.uQPMode |= 0x600;  // don't use DC QP for LP, LP QP for HP
+
+        // default QPs
+        iQPIndexY = pSC->m_param.bAlphaChannel && pSC->m_param.cNumChannels == 1?
+            pSC->WMISCP.uiDefaultQPIndexAlpha : pSC->WMISCP.uiDefaultQPIndex;
+
+		// determine the U,V index
+        iQPIndexU = pSC->WMISCP.uiDefaultQPIndexU!=0? 
+			pSC->WMISCP.uiDefaultQPIndexU: iQPIndexY; 
+        iQPIndexV = pSC->WMISCP.uiDefaultQPIndexV!=0? 
+			pSC->WMISCP.uiDefaultQPIndexV: iQPIndexY; 
+
+		// determine the QPIndexYLP
+        iQPIndexYLP = pSC->m_param.bAlphaChannel && pSC->m_param.cNumChannels == 1 ?
+            pSC->WMISCP.uiDefaultQPIndexAlpha :
+            (pSC->WMISCP.uiDefaultQPIndexYLP == 0 ? 
+			 pSC->WMISCP.uiDefaultQPIndex : pSC->WMISCP.uiDefaultQPIndexYLP); // default to QPIndex if not set
+
+		// determine the QPIndexYHP
+        iQPIndexYHP = pSC->m_param.bAlphaChannel && pSC->m_param.cNumChannels == 1 ?
+            pSC->WMISCP.uiDefaultQPIndexAlpha :
+            (pSC->WMISCP.uiDefaultQPIndexYHP == 0 ? 
+			 pSC->WMISCP.uiDefaultQPIndex : pSC->WMISCP.uiDefaultQPIndexYHP); // default to QPIndex if not set
+
+		// determine the U,V LP index
+        iQPIndexULP = pSC->WMISCP.uiDefaultQPIndexULP!=0? 
+			pSC->WMISCP.uiDefaultQPIndexULP: iQPIndexU; 
+        iQPIndexVLP = pSC->WMISCP.uiDefaultQPIndexVLP!=0? 
+			pSC->WMISCP.uiDefaultQPIndexVLP: iQPIndexV; 
+
+		// determine the U,V HP index
+        iQPIndexUHP = pSC->WMISCP.uiDefaultQPIndexUHP!=0? 
+			pSC->WMISCP.uiDefaultQPIndexUHP: iQPIndexU; 
+        iQPIndexVHP = pSC->WMISCP.uiDefaultQPIndexVHP!=0? 
+			pSC->WMISCP.uiDefaultQPIndexVHP: iQPIndexV; 
+
+		// clamp the QPIndex - 0 is lossless mode
+        if(iQPIndexY < 2)
+            iQPIndexY = 0;
+        if (iQPIndexYLP < 2)
+            iQPIndexYLP = 0;
+        if (iQPIndexYHP < 2)
+            iQPIndexYHP = 0;
+		if(iQPIndexU < 2)
+            iQPIndexU = 0;
+        if (iQPIndexULP < 2)
+            iQPIndexULP = 0;
+        if (iQPIndexUHP < 2)
+            iQPIndexUHP = 0;
+		if(iQPIndexV < 2)
+            iQPIndexV = 0;
+		if (iQPIndexVLP < 2)
+            iQPIndexVLP = 0;
+		if (iQPIndexVHP < 2)
+            iQPIndexVHP = 0;
+    }
+
+    if((pSC->m_param.uQPMode & 1) == 0){ // DC frame uniform quantization
+        if(allocateQuantizer(pSC->pTile[0].pQuantizerDC, pSC->m_param.cNumChannels, 1) != ICERR_OK)
+            return ICERR_ERROR;
+        setUniformQuantizer(pSC, 0);
+        for(i = 0; i < pSC->m_param.cNumChannels; i ++)
+            if(pSC->m_param.bTranscode)
+                pSC->pTile[0].pQuantizerDC[i]->iIndex = pSC->m_param.uiQPIndexDC[i];
+            else
+                pSC->pTile[0].pQuantizerDC[i]->iIndex = pSC->m_param.uiQPIndexDC[i] = (U8)(((i == 0 ? iQPIndexY : (i == 1) ? iQPIndexU: iQPIndexV)) & 0xff);
+        formatQuantizer(pSC->pTile[0].pQuantizerDC, (pSC->m_param.uQPMode >> 3) & 3, pSC->m_param.cNumChannels, 0, TRUE, pSC->m_param.bScaledArith);
+
+        for(i = 0; i < pSC->m_param.cNumChannels; i ++)
+            pSC->pTile[0].pQuantizerDC[i]->iOffset = (pSC->pTile[0].pQuantizerDC[i]->iQP >> 1);
+    }
+
+    if(pSC->WMISCP.sbSubband != SB_DC_ONLY){
+        if((pSC->m_param.uQPMode & 2) == 0){ // LP frame uniform quantization
+            if(allocateQuantizer(pSC->pTile[0].pQuantizerLP, pSC->m_param.cNumChannels, 1) != ICERR_OK)
+                return ICERR_ERROR;
+            setUniformQuantizer(pSC, 1);
+            for(i = 0; i < pSC->m_param.cNumChannels; i ++)
+                if(pSC->m_param.bTranscode)
+                    pSC->pTile[0].pQuantizerLP[i]->iIndex = pSC->m_param.uiQPIndexLP[i];
+                else
+                    pSC->pTile[0].pQuantizerLP[i]->iIndex = pSC->m_param.uiQPIndexLP[i] = (U8)(((i == 0 ? iQPIndexYLP : (i == 1) ? iQPIndexULP: iQPIndexVLP)) & 0xff);
+            formatQuantizer(pSC->pTile[0].pQuantizerLP, (pSC->m_param.uQPMode >> 5) & 3, pSC->m_param.cNumChannels, 0, TRUE, pSC->m_param.bScaledArith);
+        }
+
+        if(pSC->WMISCP.sbSubband != SB_NO_HIGHPASS){
+            if((pSC->m_param.uQPMode & 4) == 0){ // HP frame uniform quantization
+                if(allocateQuantizer(pSC->pTile[0].pQuantizerHP, pSC->m_param.cNumChannels, 1) != ICERR_OK)
+                    return ICERR_ERROR;
+                setUniformQuantizer(pSC, 2);
+                for(i = 0; i < pSC->m_param.cNumChannels; i ++)
+                    if(pSC->m_param.bTranscode)
+                        pSC->pTile[0].pQuantizerHP[i]->iIndex = pSC->m_param.uiQPIndexHP[i];
+                    else
+                        pSC->pTile[0].pQuantizerHP[i]->iIndex = pSC->m_param.uiQPIndexHP[i] = (U8)(((i == 0 ? iQPIndexYHP : (i == 1) ? iQPIndexUHP: iQPIndexVHP)) & 0xff);
+                formatQuantizer(pSC->pTile[0].pQuantizerHP, (pSC->m_param.uQPMode >> 7) & 3, pSC->m_param.cNumChannels, 0, FALSE, pSC->m_param.bScaledArith);
+            }
+        }
+    }
+
+    if(allocatePredInfo(pSC) != ICERR_OK){
+        return ICERR_ERROR;
+    }
+
+    if(pSC->WMISCP.cNumOfSliceMinus1V >= MAX_TILES || AllocateCodingContextEnc (pSC, pSC->WMISCP.cNumOfSliceMinus1V + 1, pSC->WMISCP.uiTrimFlexBits) != ICERR_OK){
+        return ICERR_ERROR;
+    }
+    
+    if (pSC->m_bSecondary) {
+        pSC->pIOHeader = pSC->m_pNextSC->pIOHeader;
+        pSC->m_ppBitIO = pSC->m_pNextSC->m_ppBitIO;
+        pSC->cNumBitIO = pSC->m_pNextSC->cNumBitIO;
+        pSC->cSB = pSC->m_pNextSC->cSB;
+        pSC->ppWStream = pSC->m_pNextSC->ppWStream;
+        pSC->pIndexTable = pSC->m_pNextSC->pIndexTable;
+        setBitIOPointers(pSC);
+    }
+    else {
+        StrIOEncInit(pSC);
+        setBitIOPointers(pSC);
+        WriteWMIHeader(pSC);
+    }
+
+    return ICERR_OK;
+}
+
+static Int StrEncTerm(CTXSTRCODEC ctxSC)
+{
+    CWMImageStrCodec* pSC = (CWMImageStrCodec*)ctxSC;
+    size_t j, jend = (pSC->m_pNextSC != NULL);
+
+    for (j = 0; j <= jend; j++) {
+        if (sizeof(*pSC) != pSC->cbStruct) {
+            return ICERR_ERROR;
+        }
+
+        if(pSC->m_bUVResolutionChange){
+            if(pSC->pResU != NULL)
+                free(pSC->pResU);
+            if(pSC->pResV != NULL)
+                free(pSC->pResV);
+        }
+
+        freePredInfo(pSC);
+
+        if (j == 0)
+            StrIOEncTerm(pSC);
+
+        FreeCodingContextEnc(pSC);
+        
+        freeTileInfo(pSC);
+
+        pSC->WMISCP.nExpBias -= 128; // reset
+
+        pSC = pSC->m_pNextSC;
+    }
+
+    return 0;
+}
+
+U32 setUniformTiling(U32 * pTile, U32 cNumTile, U32 cNumMB)
+{
+    U32 i, j;
+
+    while((cNumMB + cNumTile - 1) / cNumTile > 65535) // too few tiles
+        cNumTile ++;
+
+    for(i = cNumTile, j = cNumMB; i > 1; i --){
+        pTile[cNumTile - i] = (j + i - 1) / i;
+        j -= pTile[cNumTile - i];
+    }
+
+    return cNumTile;
+}
+
+U32 validateTiling(U32 * pTile, U32 cNumTile, U32 cNumMB)
+{
+    U32 i, cMBs;
+
+    if(cNumTile == 0)
+        cNumTile = 1;
+    if(cNumTile > cNumMB) // too many tiles
+        cNumTile = 1;
+    if(cNumTile > MAX_TILES)
+        cNumTile = MAX_TILES;
+
+    for(i = cMBs = 0; i + 1 < cNumTile; i ++){
+        if(pTile[i] == 0 || pTile[i] > 65535){ // invalid tile setting, resetting to uniform tiling
+            cNumTile = setUniformTiling(pTile, cNumTile, cNumMB);
+            break;
+        }
+        
+        cMBs += pTile[i];
+
+        if(cMBs >= cNumMB){
+            cNumTile = i + 1;
+            break;
+        }
+    }
+
+    // last tile
+    if(cNumMB - cMBs > 65536)
+        cNumTile = setUniformTiling(pTile, cNumTile, cNumMB);
+
+    for(i = 1; i < cNumTile; i ++)
+        pTile[i] += pTile[i - 1];
+    for(i = cNumTile - 1; i > 0; i --)
+        pTile[i] = pTile[i - 1];
+    pTile[0] = 0;
+
+    return cNumTile;
+}
+
+/*************************************************************************
+  Validate and adjust input params here
+*************************************************************************/
+Int ValidateArgs(CWMImageInfo* pII, CWMIStrCodecParam *pSCP)
+{
+    int i;
+    Bool bTooNarrowTile = FALSE;
+
+    if(pII->cWidth > (1 << 28) || pII->cHeight > (1 << 28) || pII->cWidth == 0 || pII->cHeight == 0){
+        printf("Unsurpported image size!\n");
+        return ICERR_ERROR; // unsurpported image size
+    }
+
+    if (((pSCP->cfColorFormat == YUV_420) || (pSCP->cfColorFormat == YUV_422)) && (pSCP->olOverlap == OL_TWO) && ((Int)(((U32)pII->cWidth + 15) >> 4) < 2)) {
+        printf("Image width must be at least 2 MB wide for subsampled chroma and two levels of overlap!\n");
+        return ICERR_ERROR;
+    }
+
+    if(pSCP->sbSubband == SB_ISOLATED || pSCP->sbSubband >= SB_MAX) // not allowed
+        pSCP->sbSubband = SB_ALL;
+
+    if(pII->bdBitDepth == BD_5 && (pII->cfColorFormat != CF_RGB || pII->cBitsPerUnit != 16 || pII->cLeadingPadding != 0)){
+        printf("Unsupported BD_5 image format!\n");
+        return ICERR_ERROR; // BD_5 must be compact RGB!
+    }   
+    if(pII->bdBitDepth == BD_565 && (pII->cfColorFormat != CF_RGB || pII->cBitsPerUnit != 16 || pII->cLeadingPadding != 0)){
+        printf("Unsupported BD_565 image format!\n");
+        return ICERR_ERROR; // BD_5 must be compact RGB!
+    }   
+    if(pII->bdBitDepth == BD_10 && (pII->cfColorFormat != CF_RGB || pII->cBitsPerUnit != 32 || pII->cLeadingPadding != 0)){
+        printf("Unsupported BD_10 image format!\n");
+        return ICERR_ERROR; // BD_10 must be compact RGB!
+    }
+
+    if((pII->bdBitDepth == BD_5 || pII->bdBitDepth == BD_565 || pII->bdBitDepth == BD_10) && 
+        (pSCP->cfColorFormat != YUV_420 && pSCP->cfColorFormat != YUV_422 && pSCP->cfColorFormat != Y_ONLY))
+            pSCP->cfColorFormat = YUV_444;
+
+    if(BD_1 == pII->bdBitDepth){ // binary image
+        if(pII->cfColorFormat != Y_ONLY){
+            printf("BD_1 image must be black-and white!\n");
+            return ICERR_ERROR;
+        }
+        pSCP->cfColorFormat = Y_ONLY; // can only be black white
+    }
+
+    if(pSCP->bdBitDepth != BD_LONG)
+        pSCP->bdBitDepth = BD_LONG; // currently only support 32 bit internally
+
+    if(pSCP->uAlphaMode > 1 && (pII->cfColorFormat == YUV_420 || pII->cfColorFormat == YUV_422 
+								|| pII->bdBitDepth == BD_5 || pII->bdBitDepth == BD_10 
+								|| pII->bdBitDepth == BD_1))
+    {
+        printf("Alpha is not supported for this pixel format!\n");
+        return ICERR_ERROR;
+    }
+
+    if((pSCP->cfColorFormat == YUV_420 || pSCP->cfColorFormat == YUV_422) && (pII->bdBitDepth == BD_16F || pII->bdBitDepth == BD_32F || pII->cfColorFormat == CF_RGBE))
+    {
+        printf("Float or RGBE images must be encoded with YUV 444!\n");
+        return ICERR_ERROR;
+    }
+
+    // adjust tiling
+    pSCP->cNumOfSliceMinus1V = validateTiling(pSCP->uiTileX, pSCP->cNumOfSliceMinus1V + 1, (((U32)pII->cWidth + 15) >> 4)) - 1;
+    pSCP->cNumOfSliceMinus1H = validateTiling(pSCP->uiTileY, pSCP->cNumOfSliceMinus1H + 1, (((U32)pII->cHeight + 15) >> 4)) - 1;
+
+    if (pSCP->bUseHardTileBoundaries && ((pSCP->cfColorFormat == YUV_420) || (pSCP->cfColorFormat == YUV_422)) && (pSCP->olOverlap == OL_TWO)) {
+        for (i = 1; i < (int) (pSCP->cNumOfSliceMinus1H + 1); i++) {
+            if ((Int)(pSCP->uiTileY[i] - pSCP->uiTileY[i - 1]) < 2) {
+                bTooNarrowTile = TRUE;
+                break;
+            }
+        }
+        if ((Int)((((U32)pII->cWidth + 15) >> 4) - pSCP->uiTileY[pSCP->cNumOfSliceMinus1H]) < 2) 
+            bTooNarrowTile = TRUE;
+    }
+    if (bTooNarrowTile) {
+        printf("Tile width must be at least 2 MB wide for hard tiles, subsampled chroma, and two levels of overlap!\n");
+        return ICERR_ERROR;
+    }
+
+    if(pSCP->cChannel > MAX_CHANNELS)
+        return ICERR_ERROR;
+
+    /** supported color transcoding **/
+    /** ARGB, RGB => YUV_444, YUV_422, YUV_420, Y_ONLY **/
+    /** YUV_444   =>          YUV_422, YUV_420, Y_ONLY **/
+    /** YUV_422   =>                   YUV_420, Y_ONLY **/
+    /** YUV_420   =>                            Y_ONLY **/
+
+    /** unsupported color transcoding       **/
+    /** Y_ONLY, YUV_420, YUV_422 => YUV_444 **/
+    /** Y_ONLY, YUV_420          => YUV_422 **/
+    /** Y_ONLY                   => YUV_420 **/
+    if((pII->cfColorFormat == Y_ONLY &&  pSCP->cfColorFormat != Y_ONLY) || 
+        (pSCP->cfColorFormat == YUV_422 && (pII->cfColorFormat == YUV_420 || pII->cfColorFormat == Y_ONLY)) || 
+        (pSCP->cfColorFormat == YUV_444 && (pII->cfColorFormat == YUV_422 || pII->cfColorFormat == YUV_420 || pII->cfColorFormat == Y_ONLY))){
+		pSCP->cfColorFormat = pII->cfColorFormat; // force not to do color transcoding!
+    }
+    else if (pII->cfColorFormat == NCOMPONENT) {
+		pSCP->cfColorFormat = NCOMPONENT; // force not to do color transcoding!
+    }
+    if (CMYK == pII->cfColorFormat && pSCP->cfColorFormat == NCOMPONENT) 
+    {
+        pSCP->cfColorFormat = CMYK;
+    }
+
+    if(pSCP->cfColorFormat != NCOMPONENT){
+        if(pSCP->cfColorFormat == Y_ONLY)
+            pSCP->cChannel = 1;
+        else if(pSCP->cfColorFormat == CMYK)
+            pSCP->cChannel = 4;
+        else
+            pSCP->cChannel = 3;
+    }
+
+    if(pSCP->sbSubband >= SB_MAX)
+        pSCP->sbSubband = SB_ALL;
+
+
+    pII->cChromaCenteringX = 0;
+    pII->cChromaCenteringY = 0;
+
+    return ICERR_OK;
+}
+
+/*************************************************************************
+  Initialization of CWMImageStrCodec struct
+*************************************************************************/
+static Void InitializeStrEnc(CWMImageStrCodec *pSC,
+    const CWMImageInfo* pII, const CWMIStrCodecParam *pSCP)
+{
+    pSC->cbStruct = sizeof(*pSC);
+    pSC->WMII = *pII;
+    pSC->WMISCP = *pSCP;
+
+    // set nExpBias
+    if (pSC->WMISCP.nExpBias == 0)
+        pSC->WMISCP.nExpBias = 4 + 128;//default
+    pSC->WMISCP.nExpBias += 128; // rollover arithmetic
+
+    pSC->cRow = 0;
+    pSC->cColumn = 0;
+    
+    pSC->cmbWidth = (pSC->WMII.cWidth + 15) / 16;
+    pSC->cmbHeight = (pSC->WMII.cHeight + 15) / 16;
+
+    pSC->Load = inputMBRow;
+    pSC->Quantize = quantizeMacroblock;
+    pSC->ProcessTopLeft = processMacroblock;
+    pSC->ProcessTop = processMacroblock;
+    pSC->ProcessTopRight = processMacroblock;
+    pSC->ProcessLeft = processMacroblock;
+    pSC->ProcessCenter = processMacroblock;
+    pSC->ProcessRight = processMacroblock;
+    pSC->ProcessBottomLeft = processMacroblock;
+    pSC->ProcessBottom = processMacroblock;
+    pSC->ProcessBottomRight = processMacroblock;
+
+    pSC->m_pNextSC = NULL;
+    pSC->m_bSecondary = FALSE;
+}
+
+/*************************************************************************
+   Streaming API init
+*************************************************************************/
+Int ImageStrEncInit(
+    CWMImageInfo* pII,
+    CWMIStrCodecParam *pSCP,
+    CTXSTRCODEC* pctxSC)
+{
+    static size_t cbChannels[BD_MAX] = {2, 4};
+
+    size_t cbChannel = 0, cblkChroma = 0, i;
+    size_t cbMacBlockStride = 0, cbMacBlockChroma = 0, cMacBlock = 0;
+
+    CWMImageStrCodec* pSC = NULL, *pNextSC = NULL;
+    char* pb = NULL;
+    size_t cb = 0;
+    Bool b32bit = sizeof(size_t) == 4;
+
+    Int err;
+
+    if(ValidateArgs(pII, pSCP) != ICERR_OK){
+        goto ErrorExit;
+    }
+
+    //================================================
+    *pctxSC = NULL;
+
+    //================================================
+    cbChannel = cbChannels[pSCP->bdBitDepth];
+    cblkChroma = cblkChromas[pSCP->cfColorFormat];
+    cbMacBlockStride = cbChannel * 16 * 16;
+    cbMacBlockChroma = cbChannel * 16 * cblkChroma;
+    cMacBlock = (pII->cWidth + 15) / 16;
+
+    //================================================
+    cb = sizeof(*pSC) + (128 - 1) + (PACKETLENGTH * 4 - 1) + (PACKETLENGTH * 2 ) + sizeof(*pSC->pIOHeader);
+    i = cbMacBlockStride + cbMacBlockChroma * (pSCP->cChannel - 1);
+    if(b32bit) // integer overlow/underflow check for 32-bit system
+        if(((cMacBlock >> 15) * i) & 0xffff0000)
+            return ICERR_ERROR;
+    i *= cMacBlock * 2;
+    cb += i;
+
+    pb = malloc(cb);
+    if (NULL == pb)
+    {
+        goto ErrorExit;
+    }
+    memset(pb, 0, cb);
+
+    //================================================
+    pSC = (CWMImageStrCodec*)pb; pb += sizeof(*pSC);
+
+    // Set up perf timers
+    PERFTIMER_ONLY(pSC->m_fMeasurePerf = pSCP->fMeasurePerf);
+    PERFTIMER_NEW(pSC->m_fMeasurePerf, &pSC->m_ptEndToEndPerf);
+    PERFTIMER_NEW(pSC->m_fMeasurePerf, &pSC->m_ptEncDecPerf);
+    PERFTIMER_START(pSC->m_fMeasurePerf, pSC->m_ptEndToEndPerf);
+    PERFTIMER_START(pSC->m_fMeasurePerf, pSC->m_ptEncDecPerf);
+    PERFTIMER_COPYSTARTTIME(pSC->m_fMeasurePerf, pSC->m_ptEncDecPerf, pSC->m_ptEndToEndPerf);
+
+    pSC->m_param.cfColorFormat = pSCP->cfColorFormat;
+    pSC->m_param.bAlphaChannel = (pSCP->uAlphaMode == 3);
+    pSC->m_param.cNumChannels = pSCP->cChannel;
+    pSC->m_param.cExtraPixelsTop = pSC->m_param.cExtraPixelsBottom
+        = pSC->m_param.cExtraPixelsLeft = pSC->m_param.cExtraPixelsRight = 0;
+
+    pSC->cbChannel = cbChannel;
+
+    pSC->m_param.bTranscode = pSC->bTileExtraction = FALSE;
+
+    //================================================
+    InitializeStrEnc(pSC, pII, pSCP);
+
+    //================================================
+    // 2 Macro Row buffers for each channel
+    pb = ALIGNUP(pb, 128);
+    for (i = 0; i < pSC->m_param.cNumChannels; i++) {
+        pSC->a0MBbuffer[i] = (PixelI*)pb; pb += cbMacBlockStride * pSC->cmbWidth;
+        pSC->a1MBbuffer[i] = (PixelI*)pb; pb += cbMacBlockStride * pSC->cmbWidth;
+        cbMacBlockStride = cbMacBlockChroma;
+    }
+
+    //================================================
+    // lay 2 aligned IO buffers just below pIO struct
+    pb = (char*)ALIGNUP(pb, PACKETLENGTH * 4) + PACKETLENGTH * 2;
+    pSC->pIOHeader = (BitIOInfo*)pb;
+
+    //================================================
+    err = StrEncInit(pSC);
+    if (ICERR_OK != err)
+        goto ErrorExit;
+
+    // if interleaved alpha is needed
+    if (pSC->m_param.bAlphaChannel) {
+        cbMacBlockStride = cbChannel * 16 * 16;
+        // 1. allocate new pNextSC info
+        //================================================
+        cb = sizeof(*pNextSC) + (128 - 1) + cbMacBlockStride * cMacBlock * 2;
+        pb = malloc(cb);
+        if (NULL == pb)
+        {
+            goto ErrorExit;
+        }
+        memset(pb, 0, cb);
+        //================================================
+        pNextSC = (CWMImageStrCodec*)pb; pb += sizeof(*pNextSC);
+
+        // 2. initialize pNextSC
+        pNextSC->m_param.cfColorFormat = Y_ONLY;
+        pNextSC->m_param.cNumChannels = 1;
+        pNextSC->m_param.bAlphaChannel = TRUE;
+        pNextSC->cbChannel = cbChannel;
+        //================================================
+
+        // 3. initialize arrays
+        InitializeStrEnc(pNextSC, pII, pSCP);
+        //================================================
+
+        // 2 Macro Row buffers for each channel
+        pb = ALIGNUP(pb, 128);
+        pNextSC->a0MBbuffer[0] = (PixelI*)pb; pb += cbMacBlockStride * pNextSC->cmbWidth;
+        pNextSC->a1MBbuffer[0] = (PixelI*)pb; pb += cbMacBlockStride * pNextSC->cmbWidth;
+        //================================================
+        pNextSC->pIOHeader = pSC->pIOHeader;
+        //================================================
+
+        // 4. link pSC->pNextSC = pNextSC
+        pNextSC->m_pNextSC = pSC;
+        pNextSC->m_bSecondary = TRUE;
+
+        // 5. StrEncInit
+        StrEncInit(pNextSC);
+
+        // 6. Write header of image plane
+        WriteImagePlaneHeader(pNextSC);
+    }
+
+    pSC->m_pNextSC = pNextSC;
+    //================================================
+    *pctxSC = (CTXSTRCODEC)pSC;
+
+    writeIndexTableNull(pSC);
+#if defined(WMP_OPT_SSE2) || defined(WMP_OPT_CC_ENC) || defined(WMP_OPT_TRFM_ENC)
+    StrEncOpt(pSC);
+#endif // OPT defined
+
+    PERFTIMER_STOP(pSC->m_fMeasurePerf, pSC->m_ptEncDecPerf);
+    return ICERR_OK;
+
+ErrorExit:
+    return ICERR_ERROR;
+}
+
+/*************************************************************************
+   Streaming API encode
+*************************************************************************/
+Int ImageStrEncEncode(
+    CTXSTRCODEC ctxSC,
+    const CWMImageBufferInfo* pBI)
+{
+    CWMImageStrCodec* pSC = (CWMImageStrCodec*)ctxSC;
+    CWMImageStrCodec* pNextSC = pSC->m_pNextSC;
+    ImageDataProc ProcessLeft, ProcessCenter, ProcessRight;
+
+    if (sizeof(*pSC) != pSC->cbStruct)
+    {
+        return ICERR_ERROR;
+    }
+
+    //================================
+    PERFTIMER_START(pSC->m_fMeasurePerf, pSC->m_ptEncDecPerf);
+
+    pSC->WMIBI = *pBI;
+    pSC->cColumn = 0;
+    initMRPtr(pSC);
+    if (pNextSC)
+        pNextSC->WMIBI = *pBI;
+
+    if (0 == pSC->cRow) {
+        ProcessLeft = pSC->ProcessTopLeft;
+        ProcessCenter = pSC->ProcessTop;
+        ProcessRight = pSC->ProcessTopRight;
+    }
+    else {
+        ProcessLeft = pSC->ProcessLeft;
+        ProcessCenter = pSC->ProcessCenter;
+        ProcessRight = pSC->ProcessRight;
+    }
+
+    if( pSC->Load(pSC) != ICERR_OK )
+		return ICERR_ERROR;
+    if(ProcessLeft(pSC) != ICERR_OK)
+        return ICERR_ERROR;
+    advanceMRPtr(pSC);
+
+    //================================
+    for (pSC->cColumn = 1; pSC->cColumn < pSC->cmbWidth; ++pSC->cColumn) {
+        if(ProcessCenter(pSC) != ICERR_OK)
+            return ICERR_ERROR;
+        advanceMRPtr(pSC);
+    }
+
+    //================================
+    if(ProcessRight(pSC) != ICERR_OK)
+        return ICERR_ERROR;
+    if (pSC->cRow)
+        advanceOneMBRow(pSC);
+
+    ++pSC->cRow;
+    swapMRPtr(pSC);
+
+    PERFTIMER_STOP(pSC->m_fMeasurePerf, pSC->m_ptEncDecPerf);
+    return ICERR_OK;
+}
+
+/*************************************************************************
+   Streaming API term
+*************************************************************************/
+Int ImageStrEncTerm(
+    CTXSTRCODEC ctxSC)
+{
+    CWMImageStrCodec* pSC = (CWMImageStrCodec*)ctxSC;
+    // CWMImageStrCodec *pNextSC = pSC->m_pNextSC;
+
+    if (sizeof(*pSC) != pSC->cbStruct)
+    {
+        return ICERR_ERROR;
+    }
+
+    //================================
+    PERFTIMER_START(pSC->m_fMeasurePerf, pSC->m_ptEncDecPerf);
+    pSC->cColumn = 0;
+    initMRPtr(pSC);
+
+    pSC->ProcessBottomLeft(pSC);
+    advanceMRPtr(pSC);
+
+    //================================
+    for (pSC->cColumn = 1; pSC->cColumn < pSC->cmbWidth; ++pSC->cColumn) {
+        pSC->ProcessBottom(pSC);
+        advanceMRPtr(pSC);
+    }
+
+    //================================
+    pSC->ProcessBottomRight(pSC);
+
+    //================================
+    StrEncTerm(pSC);
+
+    PERFTIMER_STOP(pSC->m_fMeasurePerf, pSC->m_ptEncDecPerf);
+    PERFTIMER_STOP(pSC->m_fMeasurePerf, pSC->m_ptEndToEndPerf);
+    PERFTIMER_REPORT(pSC->m_fMeasurePerf, pSC);
+    PERFTIMER_DELETE(pSC->m_fMeasurePerf, pSC->m_ptEncDecPerf);
+    PERFTIMER_DELETE(pSC->m_fMeasurePerf, pSC->m_ptEndToEndPerf);
+
+    free(pSC);
+    return ICERR_OK;
+}
+
+// centralized UV downsampling
+#define DF_ODD ((((d1 + d2 + d3) << 2) + (d2 << 1) + d0 + d4 + 8) >> 4)
+Void downsampleUV(CWMImageStrCodec * pSC)
+{
+    const COLORFORMAT cfInt = pSC->m_param.cfColorFormat;
+    const COLORFORMAT cfExt = pSC->WMII.cfColorFormat;
+    PixelI * pSrc, * pDst;
+    PixelI d0, d1, d2, d3, d4;
+    size_t iChannel, iRow, iColumn;
+
+    for(iChannel = 1; iChannel < 3; iChannel ++){
+        if(cfExt != YUV_422){ // need to do horizontal downsampling, 444 => 422
+            const size_t cShift = (cfInt == YUV_422 ? 1 : 0);
+
+            pSrc = (iChannel == 1 ? pSC->pResU : pSC->pResV);
+            pDst = (cfInt == YUV_422 ? pSC->p1MBbuffer[iChannel] : pSrc);
+            
+            for(iRow = 0; iRow < 16; iRow ++){
+                d0 = d4 = pSrc[idxCC[iRow][2]], d1 = d3 = pSrc[idxCC[iRow][1]], d2 = pSrc[idxCC[iRow][0]]; // left boundary
+                
+                for(iColumn = 0; iColumn + 2 < pSC->cmbWidth * 16; iColumn += 2){
+                    pDst[((iColumn >> 4) << (8 - cShift)) + idxCC[iRow][(iColumn & 15) >> cShift]] = DF_ODD;
+                    d0 = d2, d1 = d3, d2 = d4;
+                    d3 = pSrc[(((iColumn + 3) >> 4) << 8) + idxCC[iRow][(iColumn + 3) & 0xf]];
+                    d4 = pSrc[(((iColumn + 4) >> 4) << 8) + idxCC[iRow][(iColumn + 4) & 0xf]];
+                }
+
+                d4 = d2; // right boundary
+                pDst[((iColumn >> 4) << (8 - cShift)) + idxCC[iRow][(iColumn & 15) >> cShift]] = DF_ODD;
+            }
+        }
+
+        if(cfInt == YUV_420){ // need to do vertical downsampling
+            const size_t cShift = (cfExt == YUV_422 ? 0 : 1);
+            PixelI * pBuf[4];
+            size_t mbOff, pxOff;
+            
+            pDst = pSC->p1MBbuffer[iChannel];
+            pSrc = (iChannel == 1 ? pSC->pResU : pSC->pResV);
+            pBuf[0] = pSrc + (pSC->cmbWidth << (cfExt == YUV_422 ? 7 : 8));
+            pBuf[1] = pBuf[0] + pSC->cmbWidth * 8, pBuf[2] = pBuf[1] + pSC->cmbWidth * 8, pBuf[3] = pBuf[2] + pSC->cmbWidth * 8;
+
+            for(iColumn = 0; iColumn < pSC->cmbWidth * 8; iColumn ++){
+                mbOff = (iColumn >> 3) << (7 + cShift);
+                pxOff = (iColumn & 7) << cShift;
+
+                if(pSC->cRow == 0) // top image boundary
+                    d0 = d4 = pSrc[mbOff + idxCC[2][pxOff]], d1 = d3 = pSrc[mbOff + idxCC[1][pxOff]], d2 = pSrc[mbOff + idxCC[0][pxOff]]; // top MB boundary
+                else{
+                    // last row of previous MB row
+                    d0 = pBuf[0][iColumn], d1 = pBuf[1][iColumn], d2 = pBuf[2][iColumn], d3 = pBuf[3][iColumn], d4 = pSrc[mbOff + idxCC[0][pxOff]];
+                    pSC->p0MBbuffer[iChannel][((iColumn >> 3) << 6) + idxCC_420[7][iColumn & 7]] = DF_ODD;
+
+                    // for first row of current MB
+                    d0 = pBuf[2][iColumn], d1 = pBuf[3][iColumn];
+                    d2 = pSrc[mbOff + idxCC[0][pxOff]], d3 = pSrc[mbOff + idxCC[1][pxOff]], d4 = pSrc[mbOff + idxCC[2][pxOff]];
+                }
+
+                for(iRow = 0; iRow < 12; iRow += 2){
+                    pDst[((iColumn >> 3) << 6) + idxCC_420[iRow >> 1][iColumn & 7]] = DF_ODD;
+                    d0 = d2, d1 = d3, d2 = d4;
+                    d3 = pSrc[mbOff + idxCC[iRow + 3][pxOff]];
+                    d4 = pSrc[mbOff + idxCC[iRow + 4][pxOff]];
+                }
+                
+                //last row of current MB
+                pDst[((iColumn >> 3) << 6) + idxCC_420[6][iColumn & 7]] = DF_ODD;
+                d0 = d2, d1 = d3, d2 = d4;
+                d3 = pSrc[mbOff + idxCC[iRow + 3][pxOff]];
+
+                if(pSC->cRow + 1 == pSC->cmbHeight){ // bottom image boundary
+                    d4 = d2;
+                    pDst[((iColumn >> 3) << 6) + idxCC_420[7][iColumn & 7]] = DF_ODD;
+                }
+                else{
+                    for(iRow = 0; iRow < 4; iRow ++)
+                        pBuf[iRow][iColumn] = pSrc[mbOff + idxCC[iRow + 12][pxOff]];
+                }
+            }
+        }
+    }
+}
+
+// centralized horizontal padding
+Void padHorizontally(CWMImageStrCodec * pSC)
+{
+    if(pSC->WMII.cWidth != pSC->cmbWidth * 16){ // horizontal padding is necessary!
+        const COLORFORMAT cfExt = pSC->WMISCP.bYUVData ?
+            pSC->m_param.cfColorFormat : pSC->WMII.cfColorFormat;
+        size_t cFullChannel = pSC->WMISCP.cChannel;
+        size_t iLast = pSC->WMII.cWidth - 1;
+        PixelI * pCh[16];
+        size_t iChannel, iColumn, iRow;
+
+        if(cfExt == YUV_420 || cfExt == YUV_422 || cfExt == Y_ONLY)
+            cFullChannel = 1;
+
+        assert(cFullChannel <= 16);
+
+        assert(pSC->WMISCP.cChannel <= 16);
+        for(iChannel = 0; iChannel < pSC->WMISCP.cChannel; iChannel ++)
+            pCh[iChannel & 15] = pSC->p1MBbuffer[iChannel & 15];
+
+        if(pSC->m_bUVResolutionChange)
+            pCh[1] = pSC->pResU, pCh[2] = pSC->pResV;
+
+        // pad full resoluton channels
+        for(iRow = 0; iRow < 16; iRow ++){
+            const size_t iPosLast = ((iLast >> 4) << 8) + idxCC[iRow][iLast & 0xf];
+            for(iColumn = iLast + 1; iColumn < pSC->cmbWidth * 16; iColumn ++){
+                const size_t iPos = ((iColumn >> 4) << 8) + idxCC[iRow][iColumn & 0xf];
+                for(iChannel = 0; iChannel < cFullChannel; iChannel ++)
+                    pCh[iChannel & 15][iPos] = pCh[iChannel & 15][iPosLast];
+            }
+        }
+
+        if(cfExt == YUV_422) // pad YUV_422 UV
+            for(iLast >>= 1, iRow = 0; iRow < 16; iRow ++){
+                const size_t iPosLast = ((iLast >> 3) << 7) + idxCC[iRow][iLast & 7];
+                for(iColumn = iLast + 1; iColumn < pSC->cmbWidth * 8; iColumn ++){
+                    const size_t iPos = ((iColumn >> 3) << 7) + idxCC[iRow][iColumn & 7];
+                    for(iChannel = 1; iChannel < 3; iChannel ++)
+                        pCh[iChannel][iPos] = pCh[iChannel][iPosLast];
+                }
+            }
+        else if(cfExt == YUV_420) // pad YUV_420 UV
+            for(iLast >>= 1, iRow = 0; iRow < 8; iRow ++){
+                const size_t iPosLast = ((iLast >> 3) << 6) + idxCC_420[iRow][iLast & 7];
+                for(iColumn = iLast + 1; iColumn < pSC->cmbWidth * 8; iColumn ++){
+                    const size_t iPos = ((iColumn >> 3) << 6) + idxCC_420[iRow][iColumn & 7];
+                    for(iChannel = 1; iChannel < 3; iChannel ++)
+                        pCh[iChannel][iPos] = pCh[iChannel][iPosLast];
+                }
+            }
+    }
+}
+
+// centralized alpha channel color conversion, small perf penalty
+Int inputMBRowAlpha(CWMImageStrCodec* pSC)
+{
+    if(pSC->m_bSecondary == FALSE && pSC->m_pNextSC != NULL){ // alpha channel is present
+        const size_t cShift = (pSC->m_pNextSC->m_param.bScaledArith ? (SHIFTZERO + QPFRACBITS) : 0);
+        const BITDEPTH_BITS bdExt = pSC->WMII.bdBitDepth;
+        const size_t iAlphaPos = pSC->WMII.cLeadingPadding + (pSC->WMII.cfColorFormat == CMYK ? 4 : 3);//only RGB and CMYK may have interleaved alpha
+        const size_t cRow = pSC->WMIBI.cLine;
+        const size_t cColumn = pSC->WMII.cWidth;
+        const U8 * pSrc0 = (U8 *)pSC->WMIBI.pv;
+        PixelI * pA = pSC->m_pNextSC->p1MBbuffer[0];
+        size_t iRow, iColumn;
+
+        for(iRow = 0; iRow < 16; iRow ++){
+            if(bdExt == BD_8){
+                const size_t cStride = (pSC->WMII.cBitsPerUnit >> 3);
+                const U8 * pSrc = pSrc0;
+                
+                for(iColumn = 0; iColumn < cColumn; iColumn ++, pSrc += cStride)
+                    pA[((iColumn >> 4) << 8) + idxCC[iRow][iColumn & 0xf]] = ((PixelI)pSrc[iAlphaPos] - (1 << 7)) << cShift;
+            }
+            else if(bdExt == BD_16){
+                const size_t cStride = (pSC->WMII.cBitsPerUnit >> 3) / sizeof(U16);
+                const U8 nLenMantissaOrShift = pSC->m_pNextSC->WMISCP.nLenMantissaOrShift;
+                const U16 * pSrc = (U16 *)pSrc0;
+
+                for(iColumn = 0; iColumn < cColumn; iColumn ++, pSrc += cStride)
+                    pA[((iColumn >> 4) << 8) + idxCC[iRow][iColumn & 0xf]] = ((((PixelI)pSrc[iAlphaPos] - (1 << 15)) >> nLenMantissaOrShift) << cShift);
+            }
+            else if(bdExt == BD_16S){
+                const size_t cStride = (pSC->WMII.cBitsPerUnit >> 3) / sizeof(I16);
+                const U8 nLenMantissaOrShift = pSC->m_pNextSC->WMISCP.nLenMantissaOrShift;
+                const I16 * pSrc = (I16 *)pSrc0;
+            
+                for(iColumn = 0; iColumn < cColumn; iColumn ++, pSrc += cStride)
+                    pA[((iColumn >> 4) << 8) + idxCC[iRow][iColumn & 0xf]] = (((PixelI)pSrc[iAlphaPos] >> nLenMantissaOrShift) << cShift);
+            }
+            else if(bdExt == BD_16F){
+                const size_t cStride = (pSC->WMII.cBitsPerUnit >> 3) / sizeof(U16);
+                const I16 * pSrc = (I16 *)pSrc0;
+
+                for(iColumn = 0; iColumn < cColumn; iColumn ++, pSrc += cStride)
+                    pA[((iColumn >> 4) << 8) + idxCC[iRow][iColumn & 0xf]] = forwardHalf (pSrc[iAlphaPos]) << cShift;
+            }
+            else if(bdExt == BD_32S){
+                const size_t cStride = (pSC->WMII.cBitsPerUnit >> 3) / sizeof(I32);
+                const U8 nLenMantissaOrShift = pSC->m_pNextSC->WMISCP.nLenMantissaOrShift;
+                const I32 * pSrc = (I32 *)pSrc0;
+            
+                for(iColumn = 0; iColumn < cColumn; iColumn ++, pSrc += cStride)
+                    pA[((iColumn >> 4) << 8) + idxCC[iRow][iColumn & 0xf]] = (((PixelI)pSrc[iAlphaPos] >> nLenMantissaOrShift) << cShift);
+            }
+            else if(bdExt == BD_32F){
+                const size_t cStride = (pSC->WMII.cBitsPerUnit >> 3) / sizeof(float);
+                const U8 nLen = pSC->m_pNextSC->WMISCP.nLenMantissaOrShift;
+                const I8 nExpBias = pSC->m_pNextSC->WMISCP.nExpBias;
+                const float * pSrc = (float *)pSrc0;
+            
+                for(iColumn = 0; iColumn < cColumn; iColumn ++, pSrc += cStride)
+                    pA[((iColumn >> 4) << 8) + idxCC[iRow][iColumn & 0xf]] = float2pixel (pSrc[iAlphaPos], nExpBias, nLen) << cShift;
+            }
+            else // not supported
+                return ICERR_ERROR;
+
+            if(iRow + 1 < cRow) // vertical padding!
+                pSrc0 += pSC->WMIBI.cbStride;
+
+            for(iColumn = cColumn; iColumn < pSC->cmbWidth * 16; iColumn ++) // horizontal padding
+                pA[((iColumn >> 4) << 8) + idxCC[iRow][iColumn & 0xf]] =  pA[(((cColumn - 1) >> 4) << 8) + idxCC[iRow][(cColumn - 1) & 0xf]];
+        }
+    }
+    
+    return ICERR_OK;
+}
+
+// input one MB row of image data from input buffer
+Int inputMBRow(CWMImageStrCodec* pSC)
+{
+    const size_t cShift = (pSC->m_param.bScaledArith ? (SHIFTZERO + QPFRACBITS) : 0);
+    const BITDEPTH_BITS bdExt = pSC->WMII.bdBitDepth;
+    COLORFORMAT cfExt = pSC->WMII.cfColorFormat;
+    const COLORFORMAT cfInt = pSC->m_param.cfColorFormat;
+    const size_t cPixelStride = (pSC->WMII.cBitsPerUnit >> 3);
+    const size_t iRowStride = 
+		(cfExt == YUV_420 || (pSC->WMISCP.bYUVData && pSC->m_param.cfColorFormat==YUV_420)) ? 2 : 1;
+    const size_t cRow = pSC->WMIBI.cLine;
+    const size_t cColumn = pSC->WMII.cWidth;
+	const size_t iB = (pSC->WMII.bRGB ? 2 : 0);
+    const size_t iR = 2 - iB;
+    const U8 * pSrc0 = (U8 *)pSC->WMIBI.pv;
+    const U8 nLen = pSC->WMISCP.nLenMantissaOrShift;
+    const I8 nExpBias = pSC->WMISCP.nExpBias;
+
+    PixelI *pY = pSC->p1MBbuffer[0], *pU = pSC->p1MBbuffer[1], *pV = pSC->p1MBbuffer[2];
+    size_t iRow, iColumn, iPos;
+
+    // guard input buffer
+    if(checkImageBuffer(pSC, cColumn, cRow) != ICERR_OK)
+        return ICERR_ERROR;
+
+    if(pSC->m_bUVResolutionChange)  // will do downsampling somewhere else!
+        pU = pSC->pResU, pV = pSC->pResV;
+    else if(cfInt == Y_ONLY) // xxx to Y_ONLY transcoding!
+        pU = pV = pY; // write pY AFTER pU and pV so Y will overwrite U&V
+
+    for(iRow = 0; iRow < 16; iRow += iRowStride){
+        if (pSC->WMISCP.bYUVData){
+            I32 * pSrc = (I32 *)pSrc0 + pSC->WMII.cLeadingPadding;
+
+            switch(pSC->m_param.cfColorFormat){
+            case Y_ONLY:
+            case YUV_444:
+            case NCOMPONENT:
+                {
+                    const size_t cChannel = pSC->m_param.cNumChannels;
+                    PixelI * pChannel[16];
+                    size_t iChannel;
+
+                    assert(cChannel <= 16);
+                    for(iChannel = 0; iChannel < cChannel; iChannel ++)
+                        pChannel[iChannel & 15] = pSC->p1MBbuffer[iChannel & 15];
+                    if(pSC->m_bUVResolutionChange)
+                        pChannel[1] = pSC->pResU, pChannel[2] = pSC->pResV;
+
+                    for(iColumn = 0; iColumn < cColumn; iColumn ++, pSrc += cChannel){
+                        iPos = ((iColumn >> 4) << 8) + idxCC[iRow][iColumn & 0xf];
+                        for(iChannel = 0; iChannel < cChannel; iChannel ++)
+                            pChannel[iChannel & 15][iPos] = (PixelI)pSrc[iChannel & 15];
+                    }
+                }
+                break;
+
+            case YUV_422:
+                for(iColumn = 0; iColumn < cColumn; iColumn += 2, pSrc += 4){
+                    if(cfInt != Y_ONLY){
+                        iPos = ((iColumn >> 4) << 7) + idxCC[iRow][(iColumn >> 1) & 7];
+                        pU[iPos] = (PixelI)pSrc[0];
+                        pV[iPos] = (PixelI)pSrc[2];
+                    }
+
+                    pY[((iColumn >> 4) << 8) + idxCC[iRow][iColumn & 15]] = (PixelI)pSrc[1];
+                    pY[(((iColumn + 1) >> 4) << 8) + idxCC[iRow][(iColumn + 1) & 15]] = (PixelI)pSrc[3];
+                }
+                break;
+
+            case YUV_420:
+                for(iColumn = 0; iColumn < cColumn; iColumn += 2, pSrc += 6){
+                    if(cfInt != Y_ONLY){
+                        iPos = ((iColumn >> 4) << 6) + idxCC_420[iRow >> 1][(iColumn >> 1) & 7];
+                        pU[iPos] = (PixelI)pSrc[4];
+                        pV[iPos] = (PixelI)pSrc[5];
+                    }
+
+                    pY[((iColumn >> 4) << 8) + idxCC[iRow][iColumn & 15]] = (PixelI)pSrc[0];
+                    pY[(((iColumn + 1) >> 4) << 8) + idxCC[iRow][(iColumn + 1) & 15]] = (PixelI)pSrc[1];
+                    pY[((iColumn >> 4) << 8) + idxCC[iRow + 1][iColumn & 15]] = (PixelI)pSrc[2];
+                    pY[(((iColumn + 1) >> 4) << 8) + idxCC[iRow + 1][(iColumn + 1) & 15]] = (PixelI)pSrc[3];
+                }
+                break;
+
+            default:
+                assert(0);
+                break;
+            }
+        }
+        else if(bdExt == BD_8){
+            const U8 * pSrc = pSrc0 + pSC->WMII.cLeadingPadding;
+            const PixelI iOffset = (128 << cShift);
+
+            switch(cfExt){
+                case CF_RGB:
+                    assert (pSC->m_bSecondary == FALSE);
+					for(iColumn = 0; iColumn < cColumn; iColumn ++, pSrc += cPixelStride){
+						PixelI r = ((PixelI)pSrc[iR]) << cShift, g = ((PixelI)pSrc[1]) << cShift, b = ((PixelI)pSrc[iB]) << cShift;   
+
+						_CC(r, g, b); // color conversion
+				   
+						iPos = ((iColumn >> 4) << 8) + idxCC[iRow][iColumn & 0xf];
+						pU[iPos] = -r, pV[iPos] = b, pY[iPos] = g - iOffset;
+					}
+                    break;
+
+                case Y_ONLY:
+                case YUV_444:
+                case NCOMPONENT:
+                {
+                    const size_t cChannel = pSC->m_param.cNumChannels;
+                    PixelI * pChannel[16];
+                    size_t iChannel;
+
+                    assert(cChannel <= 16);
+                    for(iChannel = 0; iChannel < cChannel; iChannel ++)
+                        pChannel[iChannel & 15] = pSC->p1MBbuffer[iChannel & 15];
+                    if(pSC->m_bUVResolutionChange)
+                        pChannel[1] = pSC->pResU, pChannel[2] = pSC->pResV;
+
+                    for(iColumn = 0; iColumn < cColumn; iColumn ++, pSrc += cPixelStride){
+                        iPos = ((iColumn >> 4) << 8) + idxCC[iRow][iColumn & 0xf];
+                        for(iChannel = 0; iChannel < cChannel; iChannel ++)
+                            pChannel[iChannel & 15][iPos] = (((PixelI)pSrc[iChannel & 15]) << cShift) - iOffset;
+                    }
+                    break;
+                }
+
+                case CF_RGBE:
+                    for(iColumn = 0; iColumn < cColumn; iColumn ++, pSrc += cPixelStride){
+                        PixelI iExp = (PixelI)pSrc[3];
+                        PixelI r = forwardRGBE (pSrc[0], iExp) << cShift;
+                        PixelI g = forwardRGBE (pSrc[1], iExp) << cShift;
+                        PixelI b = forwardRGBE (pSrc[2], iExp) << cShift;
+
+                        _CC(r, g, b);
+
+                        iPos = ((iColumn >> 4) << 8) + idxCC[iRow][iColumn & 0xf];
+                        pU[iPos] = -r, pV[iPos] = b, pY[iPos] = g;
+                    }
+                    break;
+
+                case CMYK:
+                {
+                    PixelI * pK = (cfInt == CMYK ? pSC->p1MBbuffer[3] : pY); // CMYK -> YUV_xxx transcoding!
+
+                    for(iColumn = 0; iColumn < cColumn; iColumn ++, pSrc += cPixelStride){
+                        PixelI c = ((PixelI)pSrc[0]) << cShift;
+                        PixelI m = ((PixelI)pSrc[1]) << cShift;
+                        PixelI y = ((PixelI)pSrc[2]) << cShift;
+                        PixelI k = ((PixelI)pSrc[3]) << cShift;
+
+                        _CC_CMYK(c, m, y, k);
+                        
+                        iPos = ((iColumn >> 4) << 8) + idxCC[iRow][iColumn & 0xf];
+                        pU[iPos] = c, pV[iPos] = -y, pK[iPos] = k, pY[iPos] = iOffset - m;
+                    }
+                    break;
+                }
+
+                case YUV_422:
+                    for(iColumn = 0; iColumn < cColumn; iColumn += 2, pSrc += cPixelStride){
+                        if(cfInt != Y_ONLY){
+                            iPos = ((iColumn >> 4) << 7) + idxCC[iRow][(iColumn >> 1) & 7];
+                            pU[iPos] = (((PixelI)pSrc[0]) << cShift) - iOffset;
+                            pV[iPos] = (((PixelI)pSrc[2]) << cShift) - iOffset;
+                        }
+
+                        pY[((iColumn >> 4) << 8) + idxCC[iRow][iColumn & 15]] = (((PixelI)pSrc[1]) << cShift) - iOffset;
+                        pY[(((iColumn + 1) >> 4) << 8) + idxCC[iRow][(iColumn + 1) & 15]] = (((PixelI)pSrc[3]) << cShift) - iOffset;
+                    }
+                    break;
+
+                case YUV_420:
+                    for(iColumn = 0; iColumn < cColumn; iColumn += 2, pSrc += cPixelStride){
+                        if(cfInt != Y_ONLY){
+                            iPos = ((iColumn >> 4) << 6) + idxCC_420[iRow >> 1][(iColumn >> 1) & 7];
+                            pU[iPos] = (((PixelI)pSrc[4]) << cShift) - iOffset;
+                            pV[iPos] = (((PixelI)pSrc[5]) << cShift) - iOffset;
+                        }
+
+                        pY[((iColumn >> 4) << 8) + idxCC[iRow][iColumn & 15]] = (((PixelI)pSrc[0]) << cShift) - iOffset;
+                        pY[(((iColumn + 1) >> 4) << 8) + idxCC[iRow][(iColumn + 1) & 15]] = (((PixelI)pSrc[1]) << cShift) - iOffset;
+                        pY[((iColumn >> 4) << 8) + idxCC[iRow + 1][iColumn & 15]] = (((PixelI)pSrc[2]) << cShift) - iOffset;
+                        pY[(((iColumn + 1) >> 4) << 8) + idxCC[iRow + 1][(iColumn + 1) & 15]] = (((PixelI)pSrc[3]) << cShift) - iOffset;
+                    }
+                    break;
+
+                default:
+                    assert(0);
+                    break;
+            }
+        }
+        else if(bdExt == BD_16){
+            const U16 * pSrc = (U16 *)pSrc0 + pSC->WMII.cLeadingPadding;
+            const size_t cStride = cPixelStride / sizeof(U16);
+            const PixelI iOffset = ((1 << 15) >> nLen) << cShift;
+
+            switch(cfExt){
+                case CF_RGB:
+                    for(iColumn = 0; iColumn < cColumn; iColumn ++, pSrc += cStride){
+                        PixelI r = ((PixelI)pSrc[0] >> nLen) << cShift, g = ((PixelI)pSrc[1] >> nLen) << cShift, b = ((PixelI)pSrc[2] >> nLen) << cShift;
+                        
+                        _CC(r, g, b); // color conversion
+
+                        iPos = ((iColumn >> 4) << 8) + idxCC[iRow][iColumn & 0xf];
+                        pU[iPos] = -r, pV[iPos] = b, pY[iPos] = g - iOffset;
+                    }
+                    break;
+
+                case Y_ONLY:
+                case YUV_444:
+                case NCOMPONENT:
+                {
+                    const size_t cChannel = pSC->WMISCP.cChannel;
+                    size_t iChannel;
+
+                    for(iColumn = 0; iColumn < cColumn; iColumn ++, pSrc += cStride){
+                        iPos = ((iColumn >> 4) << 8) + idxCC[iRow][iColumn & 0xf];
+                        for(iChannel = 0; iChannel < cChannel; iChannel ++)
+                            pSC->p1MBbuffer[iChannel][iPos] = (((PixelI)pSrc[iChannel] >> nLen) << cShift) - iOffset;
+                    }
+                    break;
+                }
+
+                case CMYK:
+                {
+                    PixelI * pK = (cfInt == CMYK ? pSC->p1MBbuffer[3] : pY); // CMYK -> YUV_xxx transcoding!
+
+                    for(iColumn = 0; iColumn < cColumn; iColumn ++, pSrc += cStride){
+                        PixelI c = ((PixelI)pSrc[0] >> nLen) << cShift;
+                        PixelI m = ((PixelI)pSrc[1] >> nLen) << cShift;
+                        PixelI y = ((PixelI)pSrc[2] >> nLen) << cShift;
+                        PixelI k = ((PixelI)pSrc[3] >> nLen) << cShift;
+
+                        _CC_CMYK(c, m, y, k);
+                        
+                        iPos = ((iColumn >> 4) << 8) + idxCC[iRow][iColumn & 0xf];
+                        pU[iPos] = c, pV[iPos] = -y, pK[iPos] = k, pY[iPos] = iOffset - m;
+                    }
+                    break;
+                }
+
+                case YUV_422:
+                    for(iColumn = 0; iColumn < cColumn; iColumn += 2, pSrc += cStride){
+                        if(cfInt != Y_ONLY){
+                            iPos = ((iColumn >> 4) << 7) + idxCC[iRow][(iColumn >> 1) & 7];
+                            pU[iPos] = (((PixelI)pSrc[0]) << cShift) - iOffset;
+                            pV[iPos] = (((PixelI)pSrc[2]) << cShift) - iOffset;
+                        }
+
+                        pY[((iColumn >> 4) << 8) + idxCC[iRow][iColumn & 15]] = (((PixelI)pSrc[1]) << cShift) - iOffset;
+                        pY[(((iColumn + 1) >> 4) << 8) + idxCC[iRow][(iColumn + 1) & 15]] = (((PixelI)pSrc[3]) << cShift) - iOffset;
+                    }
+                    break;
+
+                case YUV_420:
+                    for(iColumn = 0; iColumn < cColumn; iColumn += 2, pSrc += cStride){
+                        if(cfInt != Y_ONLY){
+                            iPos = ((iColumn >> 4) << 6) + idxCC_420[iRow >> 1][(iColumn >> 1) & 7];
+                            pU[iPos] = (((PixelI)pSrc[4]) << cShift) - iOffset;
+                            pV[iPos] = (((PixelI)pSrc[5]) << cShift) - iOffset;
+                        }
+
+                        pY[((iColumn >> 4) << 8) + idxCC[iRow][iColumn & 15]] = (((PixelI)pSrc[0]) << cShift) - iOffset;
+                        pY[(((iColumn + 1) >> 4) << 8) + idxCC[iRow][(iColumn + 1) & 15]] = (((PixelI)pSrc[1]) << cShift) - iOffset;
+                        pY[((iColumn >> 4) << 8) + idxCC[iRow + 1][iColumn & 15]] = (((PixelI)pSrc[2]) << cShift) - iOffset;
+                        pY[(((iColumn + 1) >> 4) << 8) + idxCC[iRow + 1][(iColumn + 1) & 15]] = (((PixelI)pSrc[3]) << cShift) - iOffset;
+                    }
+                    break;
+
+                default:
+                    assert(0);
+                    break;
+            }
+        }
+        else if(bdExt == BD_16S){
+            const I16 * pSrc = (I16 *)pSrc0 + pSC->WMII.cLeadingPadding;
+            const size_t cStride = cPixelStride / sizeof(I16);
+
+            switch(cfExt){
+                case CF_RGB:
+                    for(iColumn = 0; iColumn < cColumn; iColumn ++, pSrc += cStride){
+                        PixelI r = ((PixelI)pSrc[0] >> nLen) << cShift, g = ((PixelI)pSrc[1] >> nLen) << cShift, b = ((PixelI)pSrc[2] >> nLen) << cShift;
+                        
+                        _CC(r, g, b); // color conversion
+
+                        iPos = ((iColumn >> 4) << 8) + idxCC[iRow][iColumn & 0xf];
+                        pU[iPos] = -r, pV[iPos] = b, pY[iPos] = g;
+                    }
+                    break;
+
+                case Y_ONLY:
+                case YUV_444:
+                case NCOMPONENT:
+					{
+						const size_t cChannel = pSC->WMISCP.cChannel;
+						size_t iChannel;
+	
+						for(iColumn = 0; iColumn < cColumn; iColumn ++, pSrc += cStride){
+							iPos = ((iColumn >> 4) << 8) + idxCC[iRow][iColumn & 0xf];
+							for(iChannel = 0; iChannel < cChannel; iChannel ++)
+								pSC->p1MBbuffer[iChannel][iPos] = (((PixelI)pSrc[iChannel] >> nLen) << cShift);
+						}
+					}
+				    break;
+
+                case CMYK:
+					{
+						PixelI * pK = (cfInt == CMYK ? pSC->p1MBbuffer[3] : pY); // CMYK -> YUV_xxx transcoding!
+	
+						for(iColumn = 0; iColumn < cColumn; iColumn ++, pSrc += cStride){
+							PixelI c = ((PixelI)pSrc[0] >> nLen) << cShift;
+							PixelI m = ((PixelI)pSrc[1] >> nLen) << cShift;
+							PixelI y = ((PixelI)pSrc[2] >> nLen) << cShift;
+							PixelI k = ((PixelI)pSrc[3] >> nLen) << cShift;
+	
+							_CC_CMYK(c, m, y, k);
+							
+							iPos = ((iColumn >> 4) << 8) + idxCC[iRow][iColumn & 0xf];
+							pU[iPos] = c, pV[iPos] = -y, pK[iPos] = k, pY[iPos] = -m;
+						}
+					}
+					break;
+
+                default:
+                    assert(0);
+                    break;
+            }
+        }
+        else if(bdExt == BD_16F){
+            const I16 * pSrc = (I16 *)pSrc0 + pSC->WMII.cLeadingPadding;
+            const size_t cStride = cPixelStride / sizeof(U16);
+
+            switch(cfExt){
+                case CF_RGB:
+                    for(iColumn = 0; iColumn < cColumn; iColumn ++, pSrc += cStride){
+                        PixelI r = forwardHalf (pSrc[0]) << cShift;
+                        PixelI g = forwardHalf (pSrc[1]) << cShift;
+                        PixelI b = forwardHalf (pSrc[2]) << cShift;
+                        
+                        _CC(r, g, b); // color conversion
+
+                        iPos = ((iColumn >> 4) << 8) + idxCC[iRow][iColumn & 0xf];
+                        pU[iPos] = -r, pV[iPos] = b, pY[iPos] = g;
+                    }
+                    break;
+
+                case Y_ONLY:
+                case YUV_444:
+                case NCOMPONENT:
+					{
+						const size_t cChannel = pSC->WMISCP.cChannel; // check xxx => Y_ONLY transcoding!
+						size_t iChannel;
+	
+						for(iColumn = 0; iColumn < cColumn; iColumn ++, pSrc += cStride){
+							iPos = ((iColumn >> 4) << 8) + idxCC[iRow][iColumn & 0xf];
+							for(iChannel = 0; iChannel < cChannel; iChannel ++)
+								pSC->p1MBbuffer[iChannel][iPos] = forwardHalf (pSrc[iChannel]) << cShift;
+						}
+					}
+					break;
+
+                default:
+                    assert(0);
+                    break;
+            }
+        }
+        else if(bdExt == BD_32){
+            const U32 * pSrc = (U32 *)pSrc0 + pSC->WMII.cLeadingPadding;
+            const size_t cStride = cPixelStride / sizeof(U32);
+            const PixelI iOffset = ((1 << 31) >> nLen) << cShift;
+
+            switch(cfExt){
+                case CF_RGB:
+                    for(iColumn = 0; iColumn < cColumn; iColumn ++, pSrc += cStride){
+                        PixelI r = (pSrc[0] >> nLen) << cShift, g = (pSrc[1] >> nLen) << cShift, b = (pSrc[2] >> nLen) << cShift;
+                        
+                        _CC(r, g, b); // color conversion
+
+                        iPos = ((iColumn >> 4) << 8) + idxCC[iRow][iColumn & 0xf];
+                        pU[iPos] = -r, pV[iPos] = b, pY[iPos] = g - iOffset;
+                    }
+                    break;
+
+                case Y_ONLY:
+                case YUV_444:
+                case NCOMPONENT:
+                {
+                    const size_t cChannel = pSC->WMISCP.cChannel;
+                    size_t iChannel;
+
+                    for(iColumn = 0; iColumn < cColumn; iColumn ++, pSrc += cStride){
+                        iPos = ((iColumn >> 4) << 8) + idxCC[iRow][iColumn & 0xf];
+                        for(iChannel = 0; iChannel < cChannel; iChannel ++)
+                            pSC->p1MBbuffer[iChannel][iPos] = (pSrc[iChannel] >> nLen) << cShift;
+                    }
+                    break;
+                }
+
+                default:
+                    assert(0);
+                    break;
+            }
+        }
+        else if(bdExt == BD_32S){
+            const I32 * pSrc = (I32 *)pSrc0 + pSC->WMII.cLeadingPadding;
+            const size_t cStride = cPixelStride / sizeof(I32);
+
+            switch(cfExt){
+                case CF_RGB:
+                    for(iColumn = 0; iColumn < cColumn; iColumn ++, pSrc += cStride){
+                        PixelI r = (pSrc[0] >> nLen)<< cShift, g = (pSrc[1] >> nLen)<< cShift, b = (pSrc[2] >> nLen)<< cShift;
+                        
+                        _CC(r, g, b); // color conversion
+
+                        iPos = ((iColumn >> 4) << 8) + idxCC[iRow][iColumn & 0xf];
+                        pU[iPos] = -r, pV[iPos] = b, pY[iPos] = g;
+                    }
+                    break;
+
+                case Y_ONLY:
+                case YUV_444:
+                case NCOMPONENT:
+					{
+						const size_t cChannel = pSC->WMISCP.cChannel; // check xxx => Y_ONLY transcoding!
+						size_t iChannel;
+	
+						for(iColumn = 0; iColumn < cColumn; iColumn ++, pSrc += cStride){
+							iPos = ((iColumn >> 4) << 8) + idxCC[iRow][iColumn & 0xf];
+							for(iChannel = 0; iChannel < cChannel; iChannel ++)
+								pSC->p1MBbuffer[iChannel][iPos] = (pSrc[iChannel] >> nLen) << cShift;
+						}
+					}
+					break;
+
+                default:
+                    assert(0);
+                    break;
+            }
+        }
+        else if(bdExt == BD_32F){
+            const float * pSrc = (float *)pSrc0 + pSC->WMII.cLeadingPadding;
+            const size_t cStride = cPixelStride / sizeof(float);
+
+            switch(cfExt){
+                case CF_RGB:
+                    for(iColumn = 0; iColumn < cColumn; iColumn ++, pSrc += cStride){
+                        PixelI r = float2pixel (pSrc[0], nExpBias, nLen) << cShift;
+                        PixelI g = float2pixel (pSrc[1], nExpBias, nLen) << cShift;
+                        PixelI b = float2pixel (pSrc[2], nExpBias, nLen) << cShift;
+
+                        _CC(r, g, b); // color conversion
+
+                        iPos = ((iColumn >> 4) << 8) + idxCC[iRow][iColumn & 0xf];
+                        pU[iPos] = -r, pV[iPos] = b, pY[iPos] = g;
+                    }
+                    break;
+
+                case Y_ONLY:
+                case YUV_444:
+                case NCOMPONENT:
+					{
+						const size_t cChannel = pSC->WMISCP.cChannel;
+						size_t iChannel;
+	
+						for(iColumn = 0; iColumn < cColumn; iColumn ++, pSrc += cStride){
+							iPos = ((iColumn >> 4) << 8) + idxCC[iRow][iColumn & 0xf];
+							for(iChannel = 0; iChannel < cChannel; iChannel ++)
+								pSC->p1MBbuffer[iChannel][iPos] = float2pixel (pSrc[iChannel], nExpBias, nLen) << cShift;
+						}
+					}
+					break;
+                default:
+                    assert(0);
+                    break;
+            }
+        }
+        else if(bdExt == BD_5){ // RGB 555, work for both big endian and small endian!
+            const U8 * pSrc = pSrc0;
+            const PixelI iOffset = (16 << cShift);
+
+            assert(cfExt == CF_RGB);
+            
+            for(iColumn = 0; iColumn < cColumn; iColumn ++, pSrc += cPixelStride){
+                PixelI r = (PixelI)pSrc[0], g = (PixelI)pSrc[1], b = ((g >> 2) & 0x1F) << cShift;
+
+                g = ((r >> 5) + ((g & 3) << 3)) << cShift, r = (r & 0x1F) << cShift;
+
+                _CC(r, g, b); // color conversion
+                
+                iPos = ((iColumn >> 4) << 8) + idxCC[iRow][iColumn & 0xf];
+                pU[iPos] = -r, pV[iPos] = b, pY[iPos] = g - iOffset;
+            }
+        }
+        else if(bdExt == BD_565){ // RGB 555, work for both big endian and small endian!
+            const U8 * pSrc = pSrc0;
+            const PixelI iOffset = (32 << cShift);
+
+            assert(cfExt == CF_RGB);
+            
+            for(iColumn = 0; iColumn < cColumn; iColumn ++, pSrc += cPixelStride){
+                PixelI r = (PixelI)pSrc[0], g = (PixelI)pSrc[1], b = (g >> 3) << (cShift + 1);
+
+                g = ((r >> 5) + ((g & 7) << 3)) << cShift, r = (r & 0x1F) << (cShift + 1);
+
+                _CC(r, g, b); // color conversion
+                
+                iPos = ((iColumn >> 4) << 8) + idxCC[iRow][iColumn & 0xf];
+                pU[iPos] = -r, pV[iPos] = b, pY[iPos] = g - iOffset;
+            }
+        }
+        else if(bdExt == BD_10){ //RGB 101010, work for both big endian and small endian!
+            const U8 * pSrc = pSrc0;
+            const PixelI iOffset = (512 << cShift);
+
+            assert(cfExt == CF_RGB);
+            
+            for(iColumn = 0; iColumn < cColumn; iColumn ++, pSrc += cPixelStride){
+                PixelI r = (PixelI)pSrc[0], g = (PixelI)pSrc[1], b = (PixelI)pSrc[2];
+
+                r = (r + ((g & 3) << 8)) << cShift, g = ((g >> 2) + ((b & 0xF) << 6)) << cShift;
+                b = ((b >> 4) + (((PixelI)pSrc[3] & 0x3F) << 4)) << cShift;
+
+                _CC(r, g, b); // color conversion
+                
+                iPos = ((iColumn >> 4) << 8) + idxCC[iRow][iColumn & 0xf];
+                pU[iPos] = -r, pV[iPos] = b, pY[iPos] = g - iOffset;
+            }
+        }
+        else if(bdExt == BD_1){
+            assert(cfExt == Y_ONLY);
+            for(iColumn = 0; iColumn < cColumn; iColumn ++) {
+                pY[((iColumn >> 4) << 8) + idxCC[iRow][iColumn & 0xf]] = ((pSC->WMISCP.bBlackWhite + (pSrc0[iColumn >> 3] >> (7 - (iColumn & 7)))) & 1) << cShift;
+            }
+        }
+
+        if(iRow + iRowStride < cRow) // centralized vertical padding!
+            pSrc0 += pSC->WMIBI.cbStride;
+    }
+
+    padHorizontally(pSC); // centralized horizontal padding
+
+    // centralized down-sampling
+    if(pSC->m_bUVResolutionChange)
+        downsampleUV(pSC);
+
+    // centralized alpha channel handdling
+    if (pSC->WMISCP.uAlphaMode == 3)
+        if(inputMBRowAlpha(pSC) != ICERR_OK)
+            return ICERR_ERROR;
+
+    return ICERR_OK;
+}
+
+