include/ast_jpeg_decoder.hpp

#pragma once

#include <aspeed/JTABLES.H>
#include <ast_video_types.hpp>
#include <array>
#include <cassert>
#include <cstdint>
#include <cstring>
#include <iostream>
#include <vector>

namespace ast_video {

struct ColorCache {
  ColorCache()
      : color{0x008080, 0xFF8080, 0x808080, 0xC08080}, index{0, 1, 2, 3} {}

  unsigned long color[4];
  unsigned char index[4];
  unsigned char bitMapBits{};
};

struct RGB {
  unsigned char b;
  unsigned char g;
  unsigned char r;
  unsigned char reserved;
};

enum class JpgBlock {
  JPEG_NO_SKIP_CODE = 0x00,
  JPEG_SKIP_CODE = 0x08,

  JPEG_PASS2_CODE = 0x02,
  JPEG_SKIP_PASS2_CODE = 0x0A,

  LOW_JPEG_NO_SKIP_CODE = 0x04,
  LOW_JPEG_SKIP_CODE = 0x0C,

  VQ_NO_SKIP_1_COLOR_CODE = 0x05,
  VQ_SKIP_1_COLOR_CODE = 0x0D,

  VQ_NO_SKIP_2_COLOR_CODE = 0x06,
  VQ_SKIP_2_COLOR_CODE = 0x0E,

  VQ_NO_SKIP_4_COLOR_CODE = 0x07,
  VQ_SKIP_4_COLOR_CODE = 0x0F,

  FRAME_END_CODE = 0x09,

};

class AstJpegDecoder {
 public:
  AstJpegDecoder() {
    // TODO(ed) figure out how to init this in the constructor
    yuvBuffer.resize(1920 * 1200);
    outBuffer.resize(1920 * 1200);
    for (auto &r : outBuffer) {
      r.r = 0x00;
      r.g = 0x00;
      r.b = 0x00;
      r.reserved = 0xAA;
    }

    int qfactor = 16;

    scalefactor = qfactor;
    scalefactoruv = qfactor;
    advancescalefactor = 16;
    advancescalefactoruv = 16;
    initJpgTable();
  }

  void loadQuantTable(std::array<long, 64> &quant_table) {
    float scalefactorF[8] = {1.0f, 1.387039845f, 1.306562965f, 1.175875602f,
                             1.0f, 0.785694958f, 0.541196100f, 0.275899379f};
    uint8_t j, row, col;
    std::array<uint8_t, 64> tempQT{};

    // Load quantization coefficients from JPG file, scale them for DCT and
    // reorder
    // from zig-zag order
    switch (ySelector) {
      case 0:
        stdLuminanceQt = tbl000Y;
        break;
      case 1:
        stdLuminanceQt = tbl014Y;
        break;
      case 2:
        stdLuminanceQt = tbl029Y;
        break;
      case 3:
        stdLuminanceQt = tbl043Y;
        break;
      case 4:
        stdLuminanceQt = tbl057Y;
        break;
      case 5:
        stdLuminanceQt = tbl071Y;
        break;
      case 6:
        stdLuminanceQt = tbl086Y;
        break;
      case 7:
        stdLuminanceQt = tbl100Y;
        break;
    }
    setQuantTable(stdLuminanceQt, static_cast<uint8_t>(scalefactor), tempQT);

    for (j = 0; j <= 63; j++) {
      quant_table[j] = tempQT[zigzag[j]];
    }
    j = 0;
    for (row = 0; row <= 7; row++) {
      for (col = 0; col <= 7; col++) {
        quant_table[j] = static_cast<long>(
            (quant_table[j] * scalefactorF[row] * scalefactorF[col]) * 65536);
        j++;
      }
    }
    bytePos += 64;
  }

  void loadQuantTableCb(std::array<long, 64> &quant_table) {
    float scalefactor[8] = {1.0f, 1.387039845f, 1.306562965f, 1.175875602f,
                            1.0f, 0.785694958f, 0.541196100f, 0.275899379f};
    uint8_t j, row, col;
    std::array<uint8_t, 64> tempQT{};

    // Load quantization coefficients from JPG file, scale them for DCT and
    // reorder from zig-zag order
    if (mapping == 0) {
      switch (uvSelector) {
        case 0:
          stdChrominanceQt = tbl000Y;
          break;
        case 1:
          stdChrominanceQt = tbl014Y;
          break;
        case 2:
          stdChrominanceQt = tbl029Y;
          break;
        case 3:
          stdChrominanceQt = tbl043Y;
          break;
        case 4:
          stdChrominanceQt = tbl057Y;
          break;
        case 5:
          stdChrominanceQt = tbl071Y;
          break;
        case 6:
          stdChrominanceQt = tbl086Y;
          break;
        case 7:
          stdChrominanceQt = tbl100Y;
          break;
      }
    } else {
      switch (uvSelector) {
        case 0:
          stdChrominanceQt = tbl000Uv;
          break;
        case 1:
          stdChrominanceQt = tbl014Uv;
          break;
        case 2:
          stdChrominanceQt = tbl029Uv;
          break;
        case 3:
          stdChrominanceQt = tbl043Uv;
          break;
        case 4:
          stdChrominanceQt = tbl057Uv;
          break;
        case 5:
          stdChrominanceQt = tbl071Uv;
          break;
        case 6:
          stdChrominanceQt = tbl086Uv;
          break;
        case 7:
          stdChrominanceQt = tbl100Uv;
          break;
      }
    }
    setQuantTable(stdChrominanceQt, static_cast<uint8_t>(scalefactoruv),
                  tempQT);

    for (j = 0; j <= 63; j++) {
      quant_table[j] = tempQT[zigzag[j]];
    }
    j = 0;
    for (row = 0; row <= 7; row++) {
      for (col = 0; col <= 7; col++) {
        quant_table[j] = static_cast<long>(
            (quant_table[j] * scalefactor[row] * scalefactor[col]) * 65536);
        j++;
      }
    }
    bytePos += 64;
  }
  //  Note: Added for Dual_JPEG
  void loadAdvanceQuantTable(std::array<long, 64> &quant_table) {
    float scalefactor[8] = {1.0f, 1.387039845f, 1.306562965f, 1.175875602f,
                            1.0f, 0.785694958f, 0.541196100f, 0.275899379f};
    uint8_t j, row, col;
    std::array<uint8_t, 64> tempQT{};

    // Load quantization coefficients from JPG file, scale them for DCT and
    // reorder
    // from zig-zag order
    switch (advanceSelector) {
      case 0:
        stdLuminanceQt = tbl000Y;
        break;
      case 1:
        stdLuminanceQt = tbl014Y;
        break;
      case 2:
        stdLuminanceQt = tbl029Y;
        break;
      case 3:
        stdLuminanceQt = tbl043Y;
        break;
      case 4:
        stdLuminanceQt = tbl057Y;
        break;
      case 5:
        stdLuminanceQt = tbl071Y;
        break;
      case 6:
        stdLuminanceQt = tbl086Y;
        break;
      case 7:
        stdLuminanceQt = tbl100Y;
        break;
    }
    //  Note: pass ADVANCE SCALE FACTOR to sub-function in Dual-JPEG
    setQuantTable(stdLuminanceQt, static_cast<uint8_t>(advancescalefactor),
                  tempQT);

    for (j = 0; j <= 63; j++) {
      quant_table[j] = tempQT[zigzag[j]];
    }
    j = 0;
    for (row = 0; row <= 7; row++) {
      for (col = 0; col <= 7; col++) {
        quant_table[j] = static_cast<long>(
            (quant_table[j] * scalefactor[row] * scalefactor[col]) * 65536);
        j++;
      }
    }
    bytePos += 64;
  }

  //  Note: Added for Dual-JPEG
  void loadAdvanceQuantTableCb(std::array<long, 64> &quant_table) {
    float scalefactor[8] = {1.0f, 1.387039845f, 1.306562965f, 1.175875602f,
                            1.0f, 0.785694958f, 0.541196100f, 0.275899379f};
    uint8_t j, row, col;
    std::array<uint8_t, 64> tempQT{};

    // Load quantization coefficients from JPG file, scale them for DCT and
    // reorder
    // from zig-zag order
    if (mapping == 1) {
      switch (advanceSelector) {
        case 0:
          stdChrominanceQt = tbl000Y;
          break;
        case 1:
          stdChrominanceQt = tbl014Y;
          break;
        case 2:
          stdChrominanceQt = tbl029Y;
          break;
        case 3:
          stdChrominanceQt = tbl043Y;
          break;
        case 4:
          stdChrominanceQt = tbl057Y;
          break;
        case 5:
          stdChrominanceQt = tbl071Y;
          break;
        case 6:
          stdChrominanceQt = tbl086Y;
          break;
        case 7:
          stdChrominanceQt = tbl100Y;
          break;
      }
    } else {
      switch (advanceSelector) {
        case 0:
          stdChrominanceQt = tbl000Uv;
          break;
        case 1:
          stdChrominanceQt = tbl014Uv;
          break;
        case 2:
          stdChrominanceQt = tbl029Uv;
          break;
        case 3:
          stdChrominanceQt = tbl043Uv;
          break;
        case 4:
          stdChrominanceQt = tbl057Uv;
          break;
        case 5:
          stdChrominanceQt = tbl071Uv;
          break;
        case 6:
          stdChrominanceQt = tbl086Uv;
          break;
        case 7:
          stdChrominanceQt = tbl100Uv;
          break;
      }
    }
    //  Note: pass ADVANCE SCALE FACTOR to sub-function in Dual-JPEG
    setQuantTable(stdChrominanceQt, static_cast<uint8_t>(advancescalefactoruv),
                  tempQT);

    for (j = 0; j <= 63; j++) {
      quant_table[j] = tempQT[zigzag[j]];
    }
    j = 0;
    for (row = 0; row <= 7; row++) {
      for (col = 0; col <= 7; col++) {
        quant_table[j] = static_cast<long>(
            (quant_table[j] * scalefactor[row] * scalefactor[col]) * 65536);
        j++;
      }
    }
    bytePos += 64;
  }

  void idctTransform(short *coef, uint8_t *data, uint8_t nBlock) {
#define FIX_1_082392200 ((int)277) /* FIX(1.082392200) */
#define FIX_1_414213562 ((int)362) /* FIX(1.414213562) */
#define FIX_1_847759065 ((int)473) /* FIX(1.847759065) */
#define FIX_2_613125930 ((int)669) /* FIX(2.613125930) */

#define MULTIPLY(var, cons) ((int)((var) * (cons)) >> 8)

    int tmp0, tmp1, tmp2, tmp3, tmp4, tmp5, tmp6, tmp7;
    int tmp10, tmp11, tmp12, tmp13;
    int z5, z10, z11, z12, z13;
    int workspace[64]; /* buffers data between passes */

    short *inptr = coef;
    long *quantptr;
    int *wsptr = workspace;
    unsigned char *outptr;
    unsigned char *rLimit = rlimitTable + 128;
    int ctr, dcval, dctsize = 8;

    quantptr = &qt[nBlock][0];

    // Pass 1: process columns from input (inptr), store into work array(wsptr)

    for (ctr = 8; ctr > 0; ctr--) {
      /* Due to quantization, we will usually find that many of the input
       * coefficients are zero, especially the AC terms.  We can exploit this
       * by short-circuiting the IDCT calculation for any column in which all
       * the AC terms are zero.  In that case each output is equal to the
       * DC coefficient (with scale factor as needed).
       * With typical images and quantization tables, half or more of the
       * column DCT calculations can be simplified this way.
       */

      if ((inptr[dctsize * 1] | inptr[dctsize * 2] | inptr[dctsize * 3] |
           inptr[dctsize * 4] | inptr[dctsize * 5] | inptr[dctsize * 6] |
           inptr[dctsize * 7]) == 0) {
        /* AC terms all zero */
        dcval = static_cast<int>((inptr[dctsize * 0] * quantptr[dctsize * 0]) >>
                                 16);

        wsptr[dctsize * 0] = dcval;
        wsptr[dctsize * 1] = dcval;
        wsptr[dctsize * 2] = dcval;
        wsptr[dctsize * 3] = dcval;
        wsptr[dctsize * 4] = dcval;
        wsptr[dctsize * 5] = dcval;
        wsptr[dctsize * 6] = dcval;
        wsptr[dctsize * 7] = dcval;

        inptr++; /* advance pointers to next column */
        quantptr++;
        wsptr++;
        continue;
      }

      /* Even part */

      tmp0 = (inptr[dctsize * 0] * quantptr[dctsize * 0]) >> 16;
      tmp1 = (inptr[dctsize * 2] * quantptr[dctsize * 2]) >> 16;
      tmp2 = (inptr[dctsize * 4] * quantptr[dctsize * 4]) >> 16;
      tmp3 = (inptr[dctsize * 6] * quantptr[dctsize * 6]) >> 16;

      tmp10 = tmp0 + tmp2; /* phase 3 */
      tmp11 = tmp0 - tmp2;

      tmp13 = tmp1 + tmp3;                                    /* phases 5-3 */
      tmp12 = MULTIPLY(tmp1 - tmp3, FIX_1_414213562) - tmp13; /* 2*c4 */

      tmp0 = tmp10 + tmp13; /* phase 2 */
      tmp3 = tmp10 - tmp13;
      tmp1 = tmp11 + tmp12;
      tmp2 = tmp11 - tmp12;

      /* Odd part */

      tmp4 = (inptr[dctsize * 1] * quantptr[dctsize * 1]) >> 16;
      tmp5 = (inptr[dctsize * 3] * quantptr[dctsize * 3]) >> 16;
      tmp6 = (inptr[dctsize * 5] * quantptr[dctsize * 5]) >> 16;
      tmp7 = (inptr[dctsize * 7] * quantptr[dctsize * 7]) >> 16;

      z13 = tmp6 + tmp5; /* phase 6 */
      z10 = tmp6 - tmp5;
      z11 = tmp4 + tmp7;
      z12 = tmp4 - tmp7;

      tmp7 = z11 + z13;                             /* phase 5 */
      tmp11 = MULTIPLY(z11 - z13, FIX_1_414213562); /* 2*c4 */

      z5 = MULTIPLY(z10 + z12, FIX_1_847759065);    /* 2*c2 */
      tmp10 = MULTIPLY(z12, FIX_1_082392200) - z5;  /* 2*(c2-c6) */
      tmp12 = MULTIPLY(z10, -FIX_2_613125930) + z5; /* -2*(c2+c6) */

      tmp6 = tmp12 - tmp7; /* phase 2 */
      tmp5 = tmp11 - tmp6;
      tmp4 = tmp10 + tmp5;

      wsptr[dctsize * 0] = (tmp0 + tmp7);
      wsptr[dctsize * 7] = (tmp0 - tmp7);
      wsptr[dctsize * 1] = (tmp1 + tmp6);
      wsptr[dctsize * 6] = (tmp1 - tmp6);
      wsptr[dctsize * 2] = (tmp2 + tmp5);
      wsptr[dctsize * 5] = (tmp2 - tmp5);
      wsptr[dctsize * 4] = (tmp3 + tmp4);
      wsptr[dctsize * 3] = (tmp3 - tmp4);

      inptr++; /* advance pointers to next column */
      quantptr++;
      wsptr++;
    }

/* Pass 2: process rows from work array, store into output array. */
/* Note that we must descale the results by a factor of 8 == 2**3, */
/* and also undo the PASS1_BITS scaling. */

//#define RANGE_MASK 1023; //2 bits wider than legal samples
#define PASS1_BITS 0
#define IDESCALE(x, n) ((int)((x) >> (n)))

    wsptr = workspace;
    for (ctr = 0; ctr < dctsize; ctr++) {
      outptr = data + ctr * 8;

      /* Rows of zeroes can be exploited in the same way as we did with columns.
       * However, the column calculation has created many nonzero AC terms, so
       * the simplification applies less often (typically 5% to 10% of the
       * time). On machines with very fast multiplication, it's possible that
       * the test takes more time than it's worth.  In that case this section
       * may be commented out.
       */
      /* Even part */

      tmp10 = (wsptr[0] + wsptr[4]);
      tmp11 = (wsptr[0] - wsptr[4]);

      tmp13 = (wsptr[2] + wsptr[6]);
      tmp12 = MULTIPLY((int)wsptr[2] - (int)wsptr[6], FIX_1_414213562) - tmp13;

      tmp0 = tmp10 + tmp13;
      tmp3 = tmp10 - tmp13;
      tmp1 = tmp11 + tmp12;
      tmp2 = tmp11 - tmp12;

      /* Odd part */

      z13 = wsptr[5] + wsptr[3];
      z10 = wsptr[5] - wsptr[3];
      z11 = wsptr[1] + wsptr[7];
      z12 = wsptr[1] - wsptr[7];

      tmp7 = z11 + z13;                             /* phase 5 */
      tmp11 = MULTIPLY(z11 - z13, FIX_1_414213562); /* 2*c4 */

      z5 = MULTIPLY(z10 + z12, FIX_1_847759065);    /* 2*c2 */
      tmp10 = MULTIPLY(z12, FIX_1_082392200) - z5;  /* 2*(c2-c6) */
      tmp12 = MULTIPLY(z10, -FIX_2_613125930) + z5; /* -2*(c2+c6) */

      tmp6 = tmp12 - tmp7; /* phase 2 */
      tmp5 = tmp11 - tmp6;
      tmp4 = tmp10 + tmp5;

      /* Final output stage: scale down by a factor of 8 and range-limit */

      outptr[0] = rLimit[IDESCALE((tmp0 + tmp7), (PASS1_BITS + 3)) & 1023L];
      outptr[7] = rLimit[IDESCALE((tmp0 - tmp7), (PASS1_BITS + 3)) & 1023L];
      outptr[1] = rLimit[IDESCALE((tmp1 + tmp6), (PASS1_BITS + 3)) & 1023L];
      outptr[6] = rLimit[IDESCALE((tmp1 - tmp6), (PASS1_BITS + 3)) & 1023L];
      outptr[2] = rLimit[IDESCALE((tmp2 + tmp5), (PASS1_BITS + 3)) & 1023L];
      outptr[5] = rLimit[IDESCALE((tmp2 - tmp5), (PASS1_BITS + 3)) & 1023L];
      outptr[4] = rLimit[IDESCALE((tmp3 + tmp4), (PASS1_BITS + 3)) & 1023L];
      outptr[3] = rLimit[IDESCALE((tmp3 - tmp4), (PASS1_BITS + 3)) & 1023L];

      wsptr += dctsize; /* advance pointer to next row */
    }
  }
  void yuvToRgb(
      int txb, int tyb,
      unsigned char
          *pYCbCr,       // in, Y: 256 or 64 bytes; Cb: 64 bytes; Cr: 64 bytes
      struct RGB *pYUV,  // in, Y: 256 or 64 bytes; Cb: 64 bytes; Cr: 64 bytes
      unsigned char
          *pBgr  // out, BGR format, 16*16*3 = 768 bytes; or 8*8*3=192 bytes
  ) {
    int i, j, pos, m, n;
    unsigned char cb, cr, *py, *pcb, *pcr, *py420[4];
    int y;
    struct RGB *pByte;
    int nBlocksInMcu = 6;
    unsigned int pixelX, pixelY;

    pByte = reinterpret_cast<struct RGB *>(pBgr);
    if (yuvmode == YuvMode::YUV444) {
      py = pYCbCr;
      pcb = pYCbCr + 64;
      pcr = pcb + 64;

      pixelX = txb * 8;
      pixelY = tyb * 8;
      pos = (pixelY * width) + pixelX;

      for (j = 0; j < 8; j++) {
        for (i = 0; i < 8; i++) {
          m = ((j << 3) + i);
          y = py[m];
          cb = pcb[m];
          cr = pcr[m];
          n = pos + i;
          // For 2Pass. Save the YUV value
          pYUV[n].b = cb;
          pYUV[n].g = y;
          pYUV[n].r = cr;
          pByte[n].b = rlimitTable[mY[y] + mCbToB[cb]];
          pByte[n].g = rlimitTable[mY[y] + mCbToG[cb] + mCrToG[cr]];
          pByte[n].r = rlimitTable[mY[y] + mCrToR[cr]];
        }
        pos += width;
      }
    } else {
      for (i = 0; i < nBlocksInMcu - 2; i++) {
        py420[i] = pYCbCr + i * 64;
      }
      pcb = pYCbCr + (nBlocksInMcu - 2) * 64;
      pcr = pcb + 64;

      pixelX = txb * 16;
      pixelY = tyb * 16;
      pos = (pixelY * width) + pixelX;

      for (j = 0; j < 16; j++) {
        for (i = 0; i < 16; i++) {
          //	block number is ((j/8) * 2 + i/8)={0, 1, 2, 3}
          y = *(py420[(j >> 3) * 2 + (i >> 3)]++);
          m = ((j >> 1) << 3) + (i >> 1);
          cb = pcb[m];
          cr = pcr[m];
          n = pos + i;
          pByte[n].b = rlimitTable[mY[y] + mCbToB[cb]];
          pByte[n].g = rlimitTable[mY[y] + mCbToG[cb] + mCrToG[cr]];
          pByte[n].r = rlimitTable[mY[y] + mCrToR[cr]];
        }
        pos += width;
      }
    }
  }
  void yuvToBuffer(
      int txb, int tyb,
      unsigned char
          *pYCbCr,  // in, Y: 256 or 64 bytes; Cb: 64 bytes; Cr: 64 bytes
      struct RGB
          *pYUV,  // out, BGR format, 16*16*3 = 768 bytes; or 8*8*3=192 bytes
      unsigned char
          *pBgr  // out, BGR format, 16*16*3 = 768 bytes; or 8*8*3=192 bytes
  ) {
    int i, j, pos, m, n;
    unsigned char cb, cr, *py, *pcb, *pcr, *py420[4];
    int y;
    struct RGB *pByte;
    int nBlocksInMcu = 6;
    unsigned int pixelX, pixelY;

    pByte = reinterpret_cast<struct RGB *>(pBgr);
    if (yuvmode == YuvMode::YUV444) {
      py = pYCbCr;
      pcb = pYCbCr + 64;
      pcr = pcb + 64;

      pixelX = txb * 8;
      pixelY = tyb * 8;
      pos = (pixelY * width) + pixelX;

      for (j = 0; j < 8; j++) {
        for (i = 0; i < 8; i++) {
          m = ((j << 3) + i);
          n = pos + i;
          y = pYUV[n].g + (py[m] - 128);
          cb = pYUV[n].b + (pcb[m] - 128);
          cr = pYUV[n].r + (pcr[m] - 128);
          pYUV[n].b = cb;
          pYUV[n].g = y;
          pYUV[n].r = cr;
          pByte[n].b = rlimitTable[mY[y] + mCbToB[cb]];
          pByte[n].g = rlimitTable[mY[y] + mCbToG[cb] + mCrToG[cr]];
          pByte[n].r = rlimitTable[mY[y] + mCrToR[cr]];
        }
        pos += width;
      }
    } else {
      for (i = 0; i < nBlocksInMcu - 2; i++) {
        py420[i] = pYCbCr + i * 64;
      }
      pcb = pYCbCr + (nBlocksInMcu - 2) * 64;
      pcr = pcb + 64;

      pixelX = txb * 16;
      pixelY = tyb * 16;
      pos = (pixelY * width) + pixelX;

      for (j = 0; j < 16; j++) {
        for (i = 0; i < 16; i++) {
          //	block number is ((j/8) * 2 + i/8)={0, 1, 2, 3}
          y = *(py420[(j >> 3) * 2 + (i >> 3)]++);
          m = ((j >> 1) << 3) + (i >> 1);
          cb = pcb[m];
          cr = pcr[m];
          n = pos + i;
          pByte[n].b = rlimitTable[mY[y] + mCbToB[cb]];
          pByte[n].g = rlimitTable[mY[y] + mCbToG[cb] + mCrToG[cr]];
          pByte[n].r = rlimitTable[mY[y] + mCrToR[cr]];
        }
        pos += width;
      }
    }
  }
  void decompress(int txb, int tyb, char *outBuf, uint8_t QT_TableSelection) {
    unsigned char *ptr;
    unsigned char byTileYuv[768] = {};

    memset(dctCoeff, 0, 384 * 2);
    ptr = byTileYuv;
    processHuffmanDataUnit(ydcNr, yacNr, &dcy, 0);
    idctTransform(dctCoeff, ptr, QT_TableSelection);
    ptr += 64;

    if (yuvmode == YuvMode::YUV420) {
      processHuffmanDataUnit(ydcNr, yacNr, &dcy, 64);
      idctTransform(dctCoeff + 64, ptr, QT_TableSelection);
      ptr += 64;

      processHuffmanDataUnit(ydcNr, yacNr, &dcy, 128);
      idctTransform(dctCoeff + 128, ptr, QT_TableSelection);
      ptr += 64;

      processHuffmanDataUnit(ydcNr, yacNr, &dcy, 192);
      idctTransform(dctCoeff + 192, ptr, QT_TableSelection);
      ptr += 64;

      processHuffmanDataUnit(cbDcNr, cbAcNr, &dcCb, 256);
      idctTransform(dctCoeff + 256, ptr, QT_TableSelection + 1);
      ptr += 64;

      processHuffmanDataUnit(crDcNr, crAcNr, &dcCr, 320);
      idctTransform(dctCoeff + 320, ptr, QT_TableSelection + 1);
    } else {
      processHuffmanDataUnit(cbDcNr, cbAcNr, &dcCb, 64);
      idctTransform(dctCoeff + 64, ptr, QT_TableSelection + 1);
      ptr += 64;

      processHuffmanDataUnit(crDcNr, crAcNr, &dcCr, 128);
      idctTransform(dctCoeff + 128, ptr, QT_TableSelection + 1);
    }

    //    yuvToRgb (txb, tyb, byTileYuv, (unsigned char *)outBuf);
    //  yuvBuffer for YUV record
    yuvToRgb(txb, tyb, byTileYuv, yuvBuffer.data(),
             reinterpret_cast<unsigned char *>(outBuf));
  }

  void decompress2Pass(int txb, int tyb, char *outBuf,
                       uint8_t QT_TableSelection) {
    unsigned char *ptr;
    unsigned char byTileYuv[768];
    memset(dctCoeff, 0, 384 * 2);

    ptr = byTileYuv;
    processHuffmanDataUnit(ydcNr, yacNr, &dcy, 0);
    idctTransform(dctCoeff, ptr, QT_TableSelection);
    ptr += 64;

    processHuffmanDataUnit(cbDcNr, cbAcNr, &dcCb, 64);
    idctTransform(dctCoeff + 64, ptr, QT_TableSelection + 1);
    ptr += 64;

    processHuffmanDataUnit(crDcNr, crAcNr, &dcCr, 128);
    idctTransform(dctCoeff + 128, ptr, QT_TableSelection + 1);

    yuvToBuffer(txb, tyb, byTileYuv, yuvBuffer.data(),
                reinterpret_cast<unsigned char *>(outBuf));
    //    yuvToRgb (txb, tyb, byTileYuv, (unsigned char *)outBuf);
  }

  void vqDecompress(int txb, int tyb, char *outBuf, uint8_t QT_TableSelection,
                    struct ColorCache *VQ) {
    unsigned char *ptr, i;
    unsigned char byTileYuv[192];
    int data;

    ptr = byTileYuv;
    if (VQ->bitMapBits == 0) {
      for (i = 0; i < 64; i++) {
        ptr[0] = (VQ->color[VQ->index[0]] & 0xFF0000) >> 16;
        ptr[64] = (VQ->color[VQ->index[0]] & 0x00FF00) >> 8;
        ptr[128] = VQ->color[VQ->index[0]] & 0x0000FF;
        ptr += 1;
      }
    } else {
      for (i = 0; i < 64; i++) {
        data = static_cast<int>(lookKbits(VQ->bitMapBits));
        ptr[0] = (VQ->color[VQ->index[data]] & 0xFF0000) >> 16;
        ptr[64] = (VQ->color[VQ->index[data]] & 0x00FF00) >> 8;
        ptr[128] = VQ->color[VQ->index[data]] & 0x0000FF;
        ptr += 1;
        skipKbits(VQ->bitMapBits);
      }
    }
    //    yuvToRgb (txb, tyb, byTileYuv, (unsigned char *)outBuf);
    yuvToRgb(txb, tyb, byTileYuv, yuvBuffer.data(),
             reinterpret_cast<unsigned char *>(outBuf));
  }

  void moveBlockIndex() {
    if (yuvmode == YuvMode::YUV444) {
      txb++;
      if (txb >= static_cast<int>(width / 8)) {
        tyb++;
        if (tyb >= static_cast<int>(height / 8)) {
          tyb = 0;
        }
        txb = 0;
      }
    } else {
      txb++;
      if (txb >= static_cast<int>(width / 16)) {
        tyb++;
        if (tyb >= static_cast<int>(height / 16)) {
          tyb = 0;
        }
        txb = 0;
      }
    }
  }

  void initColorTable() {
    int i, x;
    int nScale = 1L << 16;  // equal to power(2,16)
    int nHalf = nScale >> 1;

#define FIX(x) ((int)((x)*nScale + 0.5))

    /* i is the actual input pixel value, in the range 0..MAXJSAMPLE */
    /* The Cb or Cr value we are thinking of is x = i - CENTERJSAMPLE */
    /* Cr=>r value is nearest int to 1.597656 * x */
    /* Cb=>b value is nearest int to 2.015625 * x */
    /* Cr=>g value is scaled-up -0.8125 * x */
    /* Cb=>g value is scaled-up -0.390625 * x */
    for (i = 0, x = -128; i < 256; i++, x++) {
      mCrToR[i] = (FIX(1.597656) * x + nHalf) >> 16;
      mCbToB[i] = (FIX(2.015625) * x + nHalf) >> 16;
      mCrToG[i] = (-FIX(0.8125) * x + nHalf) >> 16;
      mCbToG[i] = (-FIX(0.390625) * x + nHalf) >> 16;
    }
    for (i = 0, x = -16; i < 256; i++, x++) {
      mY[i] = (FIX(1.164) * x + nHalf) >> 16;
    }
    // For color Text Enchance Y Re-map. Recommend to disable in default
    /*
            for (i = 0; i < (VideoEngineInfo->INFData.Gamma1_Gamma2_Seperate);
       i++) {
                    temp = (double)i /
       VideoEngineInfo->INFData.Gamma1_Gamma2_Seperate;
                    temp1 = 1.0 / VideoEngineInfo->INFData.Gamma1Parameter;
                    mY[i] =
       (BYTE)(VideoEngineInfo->INFData.Gamma1_Gamma2_Seperate * pow (temp,
       temp1));
                    if (mY[i] > 255) mY[i] = 255;
            }
            for (i = (VideoEngineInfo->INFData.Gamma1_Gamma2_Seperate); i < 256;
       i++) {
                    mY[i] =
       (BYTE)((VideoEngineInfo->INFData.Gamma1_Gamma2_Seperate) + (256 -
       VideoEngineInfo->INFData.Gamma1_Gamma2_Seperate) * ( pow((double)((i -
       VideoEngineInfo->INFData.Gamma1_Gamma2_Seperate) / (256 -
       (VideoEngineInfo->INFData.Gamma1_Gamma2_Seperate))), (1.0 /
       VideoEngineInfo->INFData.Gamma2Parameter)) ));
                    if (mY[i] > 255) mY[i] = 255;
            }
    */
  }
  void loadHuffmanTable(HuffmanTable *HT, const unsigned char *nrcode,
                        const unsigned char *value,
                        const unsigned short int *Huff_code) {
    unsigned char k, j, i;
    unsigned int code, codeIndex;

    for (j = 1; j <= 16; j++) {
      HT->length[j] = nrcode[j];
    }
    for (i = 0, k = 1; k <= 16; k++) {
      for (j = 0; j < HT->length[k]; j++) {
        HT->v[wordHiLo(k, j)] = value[i];
        i++;
      }
    }

    code = 0;
    for (k = 1; k <= 16; k++) {
      HT->minorCode[k] = static_cast<unsigned short int>(code);
      for (j = 1; j <= HT->length[k]; j++) {
        code++;
      }
      HT->majorCode[k] = static_cast<unsigned short int>(code - 1);
      code *= 2;
      if (HT->length[k] == 0) {
        HT->minorCode[k] = 0xFFFF;
        HT->majorCode[k] = 0;
      }
    }

    HT->len[0] = 2;
    i = 2;

    for (codeIndex = 1; codeIndex < 65535; codeIndex++) {
      if (codeIndex < Huff_code[i]) {
        HT->len[codeIndex] = static_cast<unsigned char>(Huff_code[i + 1]);
      } else {
        i = i + 2;
        HT->len[codeIndex] = static_cast<unsigned char>(Huff_code[i + 1]);
      }
    }
  }
  void initJpgTable() {
    initColorTable();
    prepareRangeLimitTable();
    loadHuffmanTable(&htdc[0], stdDcLuminanceNrcodes, stdDcLuminanceValues,
                     dcLuminanceHuffmancode);
    loadHuffmanTable(&htac[0], stdAcLuminanceNrcodes, stdAcLuminanceValues,
                     acLuminanceHuffmancode);
    loadHuffmanTable(&htdc[1], stdDcChrominanceNrcodes, stdDcChrominanceValues,
                     dcChrominanceHuffmancode);
    loadHuffmanTable(&htac[1], stdAcChrominanceNrcodes, stdAcChrominanceValues,
                     acChrominanceHuffmancode);
  }

  void prepareRangeLimitTable()
  /* Allocate and fill in the sample_range_limit table */
  {
    int j;
    rlimitTable = reinterpret_cast<unsigned char *>(malloc(5 * 256L + 128));
    /* First segment of "simple" table: limit[x] = 0 for x < 0 */
    memset((void *)rlimitTable, 0, 256);
    rlimitTable += 256; /* allow negative subscripts of simple table */
    /* Main part of "simple" table: limit[x] = x */
    for (j = 0; j < 256; j++) {
      rlimitTable[j] = j;
    }
    /* End of simple table, rest of first half of post-IDCT table */
    for (j = 256; j < 640; j++) {
      rlimitTable[j] = 255;
    }

    /* Second half of post-IDCT table */
    memset((void *)(rlimitTable + 640), 0, 384);
    for (j = 0; j < 128; j++) {
      rlimitTable[j + 1024] = j;
    }
  }

  inline unsigned short int wordHiLo(uint8_t byte_high, uint8_t byte_low) {
    return (byte_high << 8) + byte_low;
  }

  // river
  void processHuffmanDataUnit(uint8_t DC_nr, uint8_t AC_nr,
                              signed short int *previous_DC,
                              unsigned short int position) {
    uint8_t nr = 0;
    uint8_t k;
    unsigned short int tmpHcode;
    uint8_t sizeVal, count0;
    unsigned short int *minCode;
    uint8_t *huffValues;
    uint8_t byteTemp;

    minCode = htdc[DC_nr].minorCode;
    //   maj_code=htdc[DC_nr].majorCode;
    huffValues = htdc[DC_nr].v;

    // DC
    k = htdc[DC_nr].len[static_cast<unsigned short int>(codebuf >> 16)];
    // river
    //	 tmp_Hcode=lookKbits(k);
    tmpHcode = static_cast<unsigned short int>(codebuf >> (32 - k));
    skipKbits(k);
    sizeVal =
        huffValues[wordHiLo(k, static_cast<uint8_t>(tmpHcode - minCode[k]))];
    if (sizeVal == 0) {
      dctCoeff[position + 0] = *previous_DC;
    } else {
      dctCoeff[position + 0] = *previous_DC + getKbits(sizeVal);
      *previous_DC = dctCoeff[position + 0];
    }

    // Second, AC coefficient decoding
    minCode = htac[AC_nr].minorCode;
    //   maj_code=htac[AC_nr].majorCode;
    huffValues = htac[AC_nr].v;

    nr = 1;  // AC coefficient
    do {
      k = htac[AC_nr].len[static_cast<unsigned short int>(codebuf >> 16)];
      tmpHcode = static_cast<unsigned short int>(codebuf >> (32 - k));
      skipKbits(k);

      byteTemp =
          huffValues[wordHiLo(k, static_cast<uint8_t>(tmpHcode - minCode[k]))];
      sizeVal = byteTemp & 0xF;
      count0 = byteTemp >> 4;
      if (sizeVal == 0) {
        if (count0 != 0xF) {
          break;
        }
        nr += 16;
      } else {
        nr += count0;  // skip count_0 zeroes
        dctCoeff[position + dezigzag[nr++]] = getKbits(sizeVal);
      }
    } while (nr < 64);
  }

  unsigned short int lookKbits(uint8_t k) {
    unsigned short int revcode;

    revcode = static_cast<unsigned short int>(codebuf >> (32 - k));

    return (revcode);
  }

  void skipKbits(uint8_t k) {
    unsigned long readbuf;

    if ((newbits - k) <= 0) {
      readbuf = buffer[bufferIndex];
      bufferIndex++;
      codebuf =
          (codebuf << k) | ((newbuf | (readbuf >> (newbits))) >> (32 - k));
      newbuf = readbuf << (k - newbits);
      newbits = 32 + newbits - k;
    } else {
      codebuf = (codebuf << k) | (newbuf >> (32 - k));
      newbuf = newbuf << k;
      newbits -= k;
    }
  }

  signed short int getKbits(uint8_t k) {
    signed short int signedWordvalue;

    // river
    // signed_wordvalue=lookKbits(k);
    signedWordvalue = static_cast<unsigned short int>(codebuf >> (32 - k));
    if (((1L << (k - 1)) & signedWordvalue) == 0) {
      // neg_pow2 was previously defined as the below.  It seemed silly to keep
      // a table of values around for something
      // THat's relatively easy to compute, so it was replaced with the
      // appropriate math
      // signed_wordvalue = signed_wordvalue - (0xFFFF >> (16 - k));
      std::array<signed short int, 17> negPow2 = {
          0,    -1,   -3,    -7,    -15,   -31,   -63,    -127,
          -255, -511, -1023, -2047, -4095, -8191, -16383, -32767};

      signedWordvalue = signedWordvalue + negPow2[k];
    }
    skipKbits(k);
    return signedWordvalue;
  }
  int initJpgDecoding() {
    bytePos = 0;
    loadQuantTable(qt[0]);
    loadQuantTableCb(qt[1]);
    //  Note: Added for Dual-JPEG
    loadAdvanceQuantTable(qt[2]);
    loadAdvanceQuantTableCb(qt[3]);
    return 1;
  }

  void setQuantTable(const uint8_t *basic_table, uint8_t scale_factor,
                     std::array<uint8_t, 64> &newtable)
  // Set quantization table and zigzag reorder it
  {
    uint8_t i;
    long temp;
    for (i = 0; i < 64; i++) {
      temp = (static_cast<long>(basic_table[i] * 16) / scale_factor);
      /* limit the values to the valid range */
      if (temp <= 0L) {
        temp = 1L;
      }
      if (temp > 255L) {
        temp = 255L; /* limit to baseline range if requested */
      }
      newtable[zigzag[i]] = static_cast<uint8_t>(temp);
    }
  }

  void updatereadbuf(uint32_t *codebuf, uint32_t *newbuf, int walks,
                     int *newbits, std::vector<uint32_t> &buffer) {
    unsigned long readbuf;

    if ((*newbits - walks) <= 0) {
      readbuf = buffer[bufferIndex];
      bufferIndex++;
      *codebuf = (*codebuf << walks) |
                 ((*newbuf | (readbuf >> (*newbits))) >> (32 - walks));
      *newbuf = readbuf << (walks - *newbits);
      *newbits = 32 + *newbits - walks;
    } else {
      *codebuf = (*codebuf << walks) | (*newbuf >> (32 - walks));
      *newbuf = *newbuf << walks;
      *newbits -= walks;
    }
  }

  uint32_t decode(std::vector<uint32_t> &bufferVector, unsigned long width,
                  unsigned long height, YuvMode yuvmode_in, int ySelector,
                  int uvSelector) {
    ColorCache decodeColor;
    if (width != userWidth || height != userHeight || yuvmode_in != yuvmode ||
        ySelector != ySelector || uvSelector != uvSelector) {
      yuvmode = yuvmode_in;
      ySelector = ySelector;    // 0-7
      uvSelector = uvSelector;  // 0-7
      userHeight = height;
      userWidth = width;
      width = width;
      height = height;

      // TODO(ed) Magic number section.  Document appropriately
      advanceSelector = 0;  // 0-7
      mapping = 0;          // 0 or 1

      if (yuvmode == YuvMode::YUV420) {
        if ((width % 16) != 0u) {
          width = width + 16 - (width % 16);
        }
        if ((height % 16) != 0u) {
          height = height + 16 - (height % 16);
        }
      } else {
        if ((width % 8) != 0u) {
          width = width + 8 - (width % 8);
        }
        if ((height % 8) != 0u) {
          height = height + 8 - (height % 8);
        }
      }

      initJpgDecoding();
    }
    // TODO(ed) cleanup cruft
    buffer = bufferVector.data();

    codebuf = bufferVector[0];
    newbuf = bufferVector[1];
    bufferIndex = 2;

    txb = tyb = 0;
    newbits = 32;
    dcy = dcCb = dcCr = 0;

    static const uint32_t vqHeaderMask = 0x01;
    static const uint32_t vqNoUpdateHeader = 0x00;
    static const uint32_t vqUpdateHeader = 0x01;
    static const int vqNoUpdateLength = 0x03;
    static const int vqUpdateLength = 0x1B;
    static const uint32_t vqIndexMask = 0x03;
    static const uint32_t vqColorMask = 0xFFFFFF;

    static const int blockAsT2100StartLength = 0x04;
    static const int blockAsT2100SkipLength = 20;  // S:1 H:3 X:8 Y:8

    do {
      auto blockHeader = static_cast<JpgBlock>((codebuf >> 28) & 0xFF);
      switch (blockHeader) {
        case JpgBlock::JPEG_NO_SKIP_CODE:
          updatereadbuf(&codebuf, &newbuf, blockAsT2100StartLength, &newbits,
                        bufferVector);
          decompress(txb, tyb, reinterpret_cast<char *>(outBuffer.data()), 0);
          break;
        case JpgBlock::FRAME_END_CODE:
          return 0;
          break;
        case JpgBlock::JPEG_SKIP_CODE:

          txb = (codebuf & 0x0FF00000) >> 20;
          tyb = (codebuf & 0x0FF000) >> 12;

          updatereadbuf(&codebuf, &newbuf, blockAsT2100SkipLength, &newbits,
                        bufferVector);
          decompress(txb, tyb, reinterpret_cast<char *>(outBuffer.data()), 0);
          break;
        case JpgBlock::VQ_NO_SKIP_1_COLOR_CODE:
          updatereadbuf(&codebuf, &newbuf, blockAsT2100StartLength, &newbits,
                        bufferVector);
          decodeColor.bitMapBits = 0;

          for (int i = 0; i < 1; i++) {
            decodeColor.index[i] = ((codebuf >> 29) & vqIndexMask);
            if (((codebuf >> 31) & vqHeaderMask) == vqNoUpdateHeader) {
              updatereadbuf(&codebuf, &newbuf, vqNoUpdateLength, &newbits,
                            bufferVector);
            } else {
              decodeColor.color[decodeColor.index[i]] =
                  ((codebuf >> 5) & vqColorMask);
              updatereadbuf(&codebuf, &newbuf, vqUpdateLength, &newbits,
                            bufferVector);
            }
          }
          vqDecompress(txb, tyb, reinterpret_cast<char *>(outBuffer.data()), 0,
                       &decodeColor);
          break;
        case JpgBlock::VQ_SKIP_1_COLOR_CODE:
          txb = (codebuf & 0x0FF00000) >> 20;
          tyb = (codebuf & 0x0FF000) >> 12;

          updatereadbuf(&codebuf, &newbuf, blockAsT2100SkipLength, &newbits,
                        bufferVector);
          decodeColor.bitMapBits = 0;

          for (int i = 0; i < 1; i++) {
            decodeColor.index[i] = ((codebuf >> 29) & vqIndexMask);
            if (((codebuf >> 31) & vqHeaderMask) == vqNoUpdateHeader) {
              updatereadbuf(&codebuf, &newbuf, vqNoUpdateLength, &newbits,
                            bufferVector);
            } else {
              decodeColor.color[decodeColor.index[i]] =
                  ((codebuf >> 5) & vqColorMask);
              updatereadbuf(&codebuf, &newbuf, vqUpdateLength, &newbits,
                            bufferVector);
            }
          }
          vqDecompress(txb, tyb, reinterpret_cast<char *>(outBuffer.data()), 0,
                       &decodeColor);
          break;

        case JpgBlock::VQ_NO_SKIP_2_COLOR_CODE:
          updatereadbuf(&codebuf, &newbuf, blockAsT2100StartLength, &newbits,
                        bufferVector);
          decodeColor.bitMapBits = 1;

          for (int i = 0; i < 2; i++) {
            decodeColor.index[i] = ((codebuf >> 29) & vqIndexMask);
            if (((codebuf >> 31) & vqHeaderMask) == vqNoUpdateHeader) {
              updatereadbuf(&codebuf, &newbuf, vqNoUpdateLength, &newbits,
                            bufferVector);
            } else {
              decodeColor.color[decodeColor.index[i]] =
                  ((codebuf >> 5) & vqColorMask);
              updatereadbuf(&codebuf, &newbuf, vqUpdateLength, &newbits,
                            bufferVector);
            }
          }
          vqDecompress(txb, tyb, reinterpret_cast<char *>(outBuffer.data()), 0,
                       &decodeColor);
          break;
        case JpgBlock::VQ_SKIP_2_COLOR_CODE:
          txb = (codebuf & 0x0FF00000) >> 20;
          tyb = (codebuf & 0x0FF000) >> 12;

          updatereadbuf(&codebuf, &newbuf, blockAsT2100SkipLength, &newbits,
                        bufferVector);
          decodeColor.bitMapBits = 1;

          for (int i = 0; i < 2; i++) {
            decodeColor.index[i] = ((codebuf >> 29) & vqIndexMask);
            if (((codebuf >> 31) & vqHeaderMask) == vqNoUpdateHeader) {
              updatereadbuf(&codebuf, &newbuf, vqNoUpdateLength, &newbits,
                            bufferVector);
            } else {
              decodeColor.color[decodeColor.index[i]] =
                  ((codebuf >> 5) & vqColorMask);
              updatereadbuf(&codebuf, &newbuf, vqUpdateLength, &newbits,
                            bufferVector);
            }
          }
          vqDecompress(txb, tyb, reinterpret_cast<char *>(outBuffer.data()), 0,
                       &decodeColor);

          break;
        case JpgBlock::VQ_NO_SKIP_4_COLOR_CODE:
          updatereadbuf(&codebuf, &newbuf, blockAsT2100StartLength, &newbits,
                        bufferVector);
          decodeColor.bitMapBits = 2;

          for (unsigned char &i : decodeColor.index) {
            i = ((codebuf >> 29) & vqIndexMask);
            if (((codebuf >> 31) & vqHeaderMask) == vqNoUpdateHeader) {
              updatereadbuf(&codebuf, &newbuf, vqNoUpdateLength, &newbits,
                            bufferVector);
            } else {
              decodeColor.color[i] = ((codebuf >> 5) & vqColorMask);
              updatereadbuf(&codebuf, &newbuf, vqUpdateLength, &newbits,
                            bufferVector);
            }
          }
          vqDecompress(txb, tyb, reinterpret_cast<char *>(outBuffer.data()), 0,
                       &decodeColor);

          break;

        case JpgBlock::VQ_SKIP_4_COLOR_CODE:
          txb = (codebuf & 0x0FF00000) >> 20;
          tyb = (codebuf & 0x0FF000) >> 12;

          updatereadbuf(&codebuf, &newbuf, blockAsT2100SkipLength, &newbits,
                        bufferVector);
          decodeColor.bitMapBits = 2;

          for (unsigned char &i : decodeColor.index) {
            i = ((codebuf >> 29) & vqIndexMask);
            if (((codebuf >> 31) & vqHeaderMask) == vqNoUpdateHeader) {
              updatereadbuf(&codebuf, &newbuf, vqNoUpdateLength, &newbits,
                            bufferVector);
            } else {
              decodeColor.color[i] = ((codebuf >> 5) & vqColorMask);
              updatereadbuf(&codebuf, &newbuf, vqUpdateLength, &newbits,
                            bufferVector);
            }
          }
          vqDecompress(txb, tyb, reinterpret_cast<char *>(outBuffer.data()), 0,
                       &decodeColor);

          break;
        case JpgBlock::JPEG_SKIP_PASS2_CODE:
          txb = (codebuf & 0x0FF00000) >> 20;
          tyb = (codebuf & 0x0FF000) >> 12;

          updatereadbuf(&codebuf, &newbuf, blockAsT2100SkipLength, &newbits,
                        bufferVector);
          decompress2Pass(txb, tyb, reinterpret_cast<char *>(outBuffer.data()),
                          2);

          break;
        default:
          // TODO(ed) propogate errors upstream
          return -1;
          break;
      }
      moveBlockIndex();

    } while (bufferIndex <= bufferVector.size());

    return -1;
  }

#ifdef cimg_version
  void dump_to_bitmap_file() {
    cimg_library::CImg<unsigned char> image(width, height, 1, 3);
    for (int y = 0; y < width; y++) {
      for (int x = 0; x < height; x++) {
        auto pixel = outBuffer[x + (y * width)];
        image(x, y, 0) = pixel.r;
        image(x, y, 1) = pixel.g;
        image(x, y, 2) = pixel.b;
      }
    }
    image.save("/tmp/file2.bmp");
  }
#endif

 private:
  YuvMode yuvmode{};
  // width and height are the modes your display used
  unsigned long width{};
  unsigned long height{};
  unsigned long userWidth{};
  unsigned long userHeight{};
  unsigned char ySelector{};
  int scalefactor;
  int scalefactoruv;
  int advancescalefactor;
  int advancescalefactoruv;
  int mapping{};
  unsigned char uvSelector{};
  unsigned char advanceSelector{};
  int bytePos{};  // current byte position

  // quantization tables, no more than 4 quantization tables
  std::array<std::array<long, 64>, 4> qt{};

  // DC huffman tables , no more than 4 (0..3)
  std::array<HuffmanTable, 4> htdc{};
  // AC huffman tables (0..3)
  std::array<HuffmanTable, 4> htac{};
  std::array<int, 256> mCrToR{};
  std::array<int, 256> mCbToB{};
  std::array<int, 256> mCrToG{};
  std::array<int, 256> mCbToG{};
  std::array<int, 256> mY{};
  unsigned long bufferIndex{};
  uint32_t codebuf{}, newbuf{}, readbuf{};
  const unsigned char *stdLuminanceQt{};
  const uint8_t *stdChrominanceQt{};

  signed short int dcy{}, dcCb{}, dcCr{};  // Coeficientii DC pentru Y,Cb,Cr
  signed short int dctCoeff[384]{};
  // std::vector<signed short int> dctCoeff;  // Current DCT_coefficients
  // quantization table number for Y, Cb, Cr
  uint8_t yqNr = 0, cbQNr = 1, crQNr = 1;
  // DC Huffman table number for Y,Cb, Cr
  uint8_t ydcNr = 0, cbDcNr = 1, crDcNr = 1;
  // AC Huffman table number for Y,Cb, Cr
  uint8_t yacNr = 0, cbAcNr = 1, crAcNr = 1;
  int txb = 0;
  int tyb = 0;
  int newbits{};
  uint8_t *rlimitTable{};
  std::vector<RGB> yuvBuffer;
  // TODO(ed) this shouldn't exist.  It is cruft that needs cleaning up
  uint32_t *buffer{};

 public:
  std::vector<RGB> outBuffer;
};
}  // namespace ast_video