// Copyright (c) the JPEG XL Project Authors. All rights reserved.

//

// Use of this source code is governed by a BSD-style

// license that can be found in the LICENSE file.

#include "lib/jxl/dec_frame.h"

#include <jxl/decode.h>

#include <jxl/memory_manager.h>

#include <jxl/types.h>

#include <algorithm>

#include <cstddef>

#include <cstdint>

#include <cstdlib>

#include <memory>

#include <utility>

#include <vector>

#include "lib/jxl/ac_context.h"

#include "lib/jxl/ac_strategy.h"

#include "lib/jxl/base/bits.h"

#include "lib/jxl/base/common.h"

#include "lib/jxl/base/compiler_specific.h"

#include "lib/jxl/base/data_parallel.h"

#include "lib/jxl/base/printf_macros.h"

#include "lib/jxl/base/rect.h"

#include "lib/jxl/base/span.h"

#include "lib/jxl/base/status.h"

#include "lib/jxl/chroma_from_luma.h"

#include "lib/jxl/coeff_order.h"

#include "lib/jxl/coeff_order_fwd.h"

#include "lib/jxl/common.h"  // kMaxNumPasses

#include "lib/jxl/compressed_dc.h"

#include "lib/jxl/dct_util.h"

#include "lib/jxl/dec_ans.h"

#include "lib/jxl/dec_bit_reader.h"

#include "lib/jxl/dec_cache.h"

#include "lib/jxl/dec_group.h"

#include "lib/jxl/dec_modular.h"

#include "lib/jxl/dec_noise.h"

#include "lib/jxl/dec_patch_dictionary.h"

#include "lib/jxl/entropy_coder.h"

#include "lib/jxl/epf.h"

#include "lib/jxl/fields.h"

#include "lib/jxl/frame_dimensions.h"

#include "lib/jxl/frame_header.h"

#include "lib/jxl/image_bundle.h"

#include "lib/jxl/image_metadata.h"

#include "lib/jxl/image_ops.h"

#include "lib/jxl/jpeg/jpeg_data.h"

#include "lib/jxl/loop_filter.h"

#include "lib/jxl/passes_state.h"

#include "lib/jxl/quant_weights.h"

#include "lib/jxl/quantizer.h"

#include "lib/jxl/render_pipeline/render_pipeline.h"

#include "lib/jxl/splines.h"

#include "lib/jxl/toc.h"

namespace jxl {

namespace {

Status DecodeGlobalDCInfo(BitReader* reader, bool is_jpeg,

                          PassesDecoderState* state, ThreadPool* pool) {

  JXL_RETURN_IF_ERROR(state->shared_storage.quantizer.Decode(reader));

  JXL_RETURN_IF_ERROR(DecodeBlockCtxMap(state->memory_manager(), reader,

                                        &state->shared_storage.block_ctx_map));

  JXL_RETURN_IF_ERROR(state->shared_storage.cmap.DecodeDC(reader));

  // Pre-compute info for decoding a group.

  if (is_jpeg) {

    state->shared_storage.quantizer.ClearDCMul();  // Don't dequant DC

  }

  state->shared_storage.ac_strategy.FillInvalid();

  return true;

}

}  // namespace

Status DecodeFrame(PassesDecoderState* dec_state, ThreadPool* JXL_RESTRICT pool,

                   const uint8_t* next_in, size_t avail_in,

                   FrameHeader* frame_header, ImageBundle* decoded,

                   const CodecMetadata& metadata,

                   bool use_slow_rendering_pipeline) {

  FrameDecoder frame_decoder(dec_state, metadata, pool,

                             use_slow_rendering_pipeline);

  BitReader reader(Bytes(next_in, avail_in));

  JXL_RETURN_IF_ERROR(frame_decoder.InitFrame(&reader, decoded,

                                              /*is_preview=*/false));

  JXL_RETURN_IF_ERROR(frame_decoder.InitFrameOutput());

  if (frame_header) {

    *frame_header = frame_decoder.GetFrameHeader();

  }

  JXL_RETURN_IF_ERROR(reader.AllReadsWithinBounds());

  size_t header_bytes = reader.TotalBitsConsumed() / kBitsPerByte;

  JXL_RETURN_IF_ERROR(reader.Close());

  size_t processed_bytes = header_bytes;

  Status close_ok = true;

  std::vector<std::unique_ptr<BitReader>> section_readers;

  {

    std::vector<std::unique_ptr<BitReaderScopedCloser>> section_closers;

    std::vector<FrameDecoder::SectionInfo> section_info;

    std::vector<FrameDecoder::SectionStatus> section_status;

    size_t pos = header_bytes;

    size_t index = 0;

    for (auto toc_entry : frame_decoder.Toc()) {

      JXL_RETURN_IF_ERROR(pos + toc_entry.size <= avail_in);

      auto br = make_unique<BitReader>(Bytes(next_in + pos, toc_entry.size));

      section_info.emplace_back(

          FrameDecoder::SectionInfo{br.get(), toc_entry.id, index++});

      section_closers.emplace_back(

          make_unique<BitReaderScopedCloser>(*br, close_ok));

      section_readers.emplace_back(std::move(br));

      pos += toc_entry.size;

    }

    section_status.resize(section_info.size());

    JXL_RETURN_IF_ERROR(frame_decoder.ProcessSections(

        section_info.data(), section_info.size(), section_status.data()));

    for (size_t i = 0; i < section_status.size(); i++) {

      JXL_RETURN_IF_ERROR(section_status[i] == FrameDecoder::kDone);

      processed_bytes += frame_decoder.Toc()[i].size;

    }

  }

  JXL_RETURN_IF_ERROR(close_ok);

  JXL_RETURN_IF_ERROR(frame_decoder.FinalizeFrame());

  decoded->SetDecodedBytes(processed_bytes);

  return true;

}

Status FrameDecoder::InitFrame(BitReader* JXL_RESTRICT br, ImageBundle* decoded,

                               bool is_preview) {

  decoded_ = decoded;

  JXL_ENSURE(is_finalized_);

  JxlMemoryManager* memory_manager = decoded_->memory_manager();

  // Reset the dequantization matrices to their default values.

  dec_state_->shared_storage.matrices = DequantMatrices();

  frame_header_.nonserialized_is_preview = is_preview;

  JXL_ENSURE(frame_header_.nonserialized_metadata != nullptr);

  JXL_RETURN_IF_ERROR(ReadFrameHeader(br, &frame_header_));

  frame_dim_ = frame_header_.ToFrameDimensions();

  JXL_DEBUG_V(2, "FrameHeader: %s", frame_header_.DebugString().c_str());

  const size_t num_passes = frame_header_.passes.num_passes;

  const size_t num_groups = frame_dim_.num_groups;

  // If the previous frame was not a kRegularFrame, `decoded` may have different

  // dimensions; must reset to avoid errors.

  decoded->RemoveColor();

  decoded->ClearExtraChannels();

  decoded->duration = frame_header_.animation_frame.duration;

  if (!frame_header_.nonserialized_is_preview &&

      (frame_header_.is_last || frame_header_.animation_frame.duration > 0) &&

      (frame_header_.frame_type == kRegularFrame ||

       frame_header_.frame_type == kSkipProgressive)) {

    ++dec_state_->visible_frame_index;

    dec_state_->nonvisible_frame_index = 0;

  } else {

    ++dec_state_->nonvisible_frame_index;

  }

  // Read TOC.

  const size_t toc_entries =

      NumTocEntries(num_groups, frame_dim_.num_dc_groups, num_passes);

  std::vector<uint32_t> sizes;

  std::vector<coeff_order_t> permutation;

  JXL_RETURN_IF_ERROR(

      ReadToc(memory_manager, toc_entries, br, &sizes, &permutation));

  bool have_permutation = !permutation.empty();

  toc_.resize(toc_entries);

  section_sizes_sum_ = 0;

  for (size_t i = 0; i < toc_entries; ++i) {

    toc_[i].size = sizes[i];

    size_t index = have_permutation ? permutation[i] : i;

    toc_[index].id = i;

    if (section_sizes_sum_ + toc_[i].size < section_sizes_sum_) {

      return JXL_FAILURE("group offset overflow");

    }

    section_sizes_sum_ += toc_[i].size;

  }

  if (JXL_DEBUG_V_LEVEL >= 3) {

    for (size_t i = 0; i < toc_entries; ++i) {

      JXL_DEBUG_V(3, "TOC entry %" PRIuS " size %" PRIuS " id %" PRIuS "", i,

                  toc_[i].size, toc_[i].id);

    }

  }

  JXL_ENSURE((br->TotalBitsConsumed() % kBitsPerByte) == 0);

  const size_t group_codes_begin = br->TotalBitsConsumed() / kBitsPerByte;

  JXL_ENSURE(!toc_.empty());

  // Overflow check.

  if (group_codes_begin + section_sizes_sum_ < group_codes_begin) {

    return JXL_FAILURE("Invalid group codes");

  }

  return true;

}

Status FrameDecoder::InitFrameOutput() {

  JXL_RETURN_IF_ERROR(

      InitializePassesSharedState(frame_header_, &dec_state_->shared_storage));

  JXL_RETURN_IF_ERROR(dec_state_->Init(frame_header_));

  modular_frame_decoder_.Init(frame_dim_);

  if (decoded_->IsJPEG()) {

    if (frame_header_.encoding == FrameEncoding::kModular) {

      return JXL_FAILURE("Cannot output JPEG from Modular");

    }

    jpeg::JPEGData* jpeg_data = decoded_->jpeg_data.get();

    size_t num_components = jpeg_data->components.size();

    if (num_components != 1 && num_components != 3) {

      return JXL_FAILURE("Invalid number of components");

    }

    if (frame_header_.nonserialized_metadata->m.xyb_encoded) {

      return JXL_FAILURE("Cannot decode to JPEG an XYB image");

    }

    auto jpeg_c_map = JpegOrder(ColorTransform::kYCbCr, num_components == 1);

    decoded_->jpeg_data->width = frame_dim_.xsize;

    decoded_->jpeg_data->height = frame_dim_.ysize;

    for (size_t c = 0; c < num_components; c++) {

      auto& component = jpeg_data->components[jpeg_c_map[c]];

      component.width_in_blocks =

          frame_dim_.xsize_blocks >> frame_header_.chroma_subsampling.HShift(c);

      component.height_in_blocks =

          frame_dim_.ysize_blocks >> frame_header_.chroma_subsampling.VShift(c);

      component.h_samp_factor =

          1 << frame_header_.chroma_subsampling.RawHShift(c);

      component.v_samp_factor =

          1 << frame_header_.chroma_subsampling.RawVShift(c);

      size_t num_blocks = static_cast<size_t>(component.width_in_blocks) *

                          component.height_in_blocks;

      component.coeffs.resize(num_blocks * jxl::kDCTBlockSize);

    }

  }

  // Clear the state.

  decoded_dc_global_ = false;

  decoded_ac_global_ = false;

  is_finalized_ = false;

  finalized_dc_ = false;

  num_sections_done_ = 0;

  decoded_dc_groups_.clear();

  decoded_dc_groups_.resize(frame_dim_.num_dc_groups);

  decoded_passes_per_ac_group_.clear();

  decoded_passes_per_ac_group_.resize(frame_dim_.num_groups, 0);

  processed_section_.clear();

  processed_section_.resize(toc_.size());

  allocated_ = false;

  return true;

}

Status FrameDecoder::ProcessDCGlobal(BitReader* br) {

  PassesSharedState& shared = dec_state_->shared_storage;

  JxlMemoryManager* memory_manager = shared.memory_manager;

  if (frame_header_.flags & FrameHeader::kPatches) {

    bool uses_extra_channels;

    JXL_RETURN_IF_ERROR(shared.image_features.patches.Decode(

        memory_manager, br, frame_dim_.xsize_padded, frame_dim_.ysize_padded,

        shared.metadata->m.num_extra_channels, &uses_extra_channels));

    if (frame_header_.upsampling != 1) {

      for (size_t ecups : frame_header_.extra_channel_upsampling) {

        if (ecups != frame_header_.upsampling) {

          return JXL_FAILURE(

              "Cannot use extra channels in patches if color channels are "

              "subsampled differently from extra channels");

        }

      }

    }

  } else {

    shared.image_features.patches.Clear();

  }

  shared.image_features.splines.Clear();

  if (frame_header_.flags & FrameHeader::kSplines) {

    JXL_RETURN_IF_ERROR(shared.image_features.splines.Decode(

        br, frame_dim_.xsize * frame_dim_.ysize));

  }

  if (frame_header_.flags & FrameHeader::kNoise) {

    JXL_RETURN_IF_ERROR(DecodeNoise(br, &shared.image_features.noise_params));

  }

  JXL_RETURN_IF_ERROR(dec_state_->shared_storage.matrices.DecodeDC(br));

  if (frame_header_.encoding == FrameEncoding::kVarDCT) {

    JXL_RETURN_IF_ERROR(

        jxl::DecodeGlobalDCInfo(br, decoded_->IsJPEG(), dec_state_, pool_));

  }

  // Splines' draw cache uses the color correlation map.

  if (frame_header_.flags & FrameHeader::kSplines) {

    JXL_RETURN_IF_ERROR(shared.image_features.splines.InitializeDrawCache(

        frame_dim_.xsize_upsampled, frame_dim_.ysize_upsampled,

        dec_state_->shared->cmap.base()));

  }

  Status dec_status = modular_frame_decoder_.DecodeGlobalInfo(

      br, frame_header_, /*allow_truncated_group=*/false);

  if (dec_status.IsFatalError()) return dec_status;

  if (dec_status) {

    decoded_dc_global_ = true;

  }

  return dec_status;

}

Status FrameDecoder::ProcessDCGroup(size_t dc_group_id, BitReader* br) {

  const size_t gx = dc_group_id % frame_dim_.xsize_dc_groups;

  const size_t gy = dc_group_id / frame_dim_.xsize_dc_groups;

  const LoopFilter& lf = frame_header_.loop_filter;

  if (frame_header_.encoding == FrameEncoding::kVarDCT &&

      !(frame_header_.flags & FrameHeader::kUseDcFrame)) {

    JXL_RETURN_IF_ERROR(modular_frame_decoder_.DecodeVarDCTDC(

        frame_header_, dc_group_id, br, dec_state_));

  }

  const Rect mrect(gx * frame_dim_.dc_group_dim, gy * frame_dim_.dc_group_dim,

                   frame_dim_.dc_group_dim, frame_dim_.dc_group_dim);

  JXL_RETURN_IF_ERROR(modular_frame_decoder_.DecodeGroup(

      frame_header_, mrect, br, 3, 1000,

      ModularStreamId::ModularDC(dc_group_id),

      /*zerofill=*/false, nullptr, nullptr,

      /*allow_truncated=*/false));

  if (frame_header_.encoding == FrameEncoding::kVarDCT) {

    JXL_RETURN_IF_ERROR(modular_frame_decoder_.DecodeAcMetadata(

        frame_header_, dc_group_id, br, dec_state_));

  } else if (lf.epf_iters > 0) {

    FillImage(kInvSigmaNum / lf.epf_sigma_for_modular, &dec_state_->sigma);

  }

  decoded_dc_groups_[dc_group_id] = JXL_TRUE;

  return true;

}

Status FrameDecoder::FinalizeDC() {

  // Do Adaptive DC smoothing if enabled. This *must* happen between all the

  // ProcessDCGroup and ProcessACGroup.

  JxlMemoryManager* memory_manager = dec_state_->memory_manager();

  if (frame_header_.encoding == FrameEncoding::kVarDCT &&

      !(frame_header_.flags & FrameHeader::kSkipAdaptiveDCSmoothing) &&

      !(frame_header_.flags & FrameHeader::kUseDcFrame) &&

      // Skip smoothing when reconstructing to JPEG.

      !decoded_->IsJPEG()) {

    JXL_RETURN_IF_ERROR(AdaptiveDCSmoothing(

        memory_manager, dec_state_->shared->quantizer.MulDC(),

        &dec_state_->shared_storage.dc_storage, pool_));

  }

  finalized_dc_ = true;

  return true;

}

Status FrameDecoder::AllocateOutput() {

  if (allocated_) return true;

  modular_frame_decoder_.MaybeDropFullImage();

  decoded_->origin = frame_header_.frame_origin;

  JXL_RETURN_IF_ERROR(

      dec_state_->InitForAC(frame_header_.passes.num_passes, nullptr));

  allocated_ = true;

  return true;

}

Status FrameDecoder::ProcessACGlobal(BitReader* br) {

  JXL_ENSURE(finalized_dc_);

  JxlMemoryManager* memory_manager = dec_state_->memory_manager();

  // Decode AC group.

  if (frame_header_.encoding == FrameEncoding::kVarDCT) {

    JXL_RETURN_IF_ERROR(dec_state_->shared_storage.matrices.Decode(

        memory_manager, br, &modular_frame_decoder_));

    JXL_RETURN_IF_ERROR(dec_state_->shared_storage.matrices.EnsureComputed(

        memory_manager, dec_state_->used_acs));

    size_t num_histo_bits =

        CeilLog2Nonzero(dec_state_->shared->frame_dim.num_groups);

    dec_state_->shared_storage.num_histograms =

        1 + br->ReadBits(num_histo_bits);

    JXL_DEBUG_V(3,

                "Processing AC global with %d passes and %" PRIuS

                " sets of histograms",

                frame_header_.passes.num_passes,

                dec_state_->shared_storage.num_histograms);

    dec_state_->code.resize(kMaxNumPasses);

    dec_state_->context_map.resize(kMaxNumPasses);

    // Read coefficient orders and histograms.

    size_t max_num_bits_ac = 0;

    for (size_t i = 0; i < frame_header_.passes.num_passes; i++) {

      uint16_t used_orders = U32Coder::Read(kOrderEnc, br);

      JXL_RETURN_IF_ERROR(DecodeCoeffOrders(

          memory_manager, used_orders, dec_state_->used_acs,

          &dec_state_->shared_storage

               .coeff_orders[i * dec_state_->shared_storage.coeff_order_size],

          br));

      size_t num_contexts =

          dec_state_->shared->num_histograms *

          dec_state_->shared_storage.block_ctx_map.NumACContexts();

      JXL_RETURN_IF_ERROR(DecodeHistograms(memory_manager, br, num_contexts,

                                           &dec_state_->code[i],

                                           &dec_state_->context_map[i]));

      // Add extra values to enable the cheat in hot loop of DecodeACVarBlock.

      dec_state_->context_map[i].resize(

          num_contexts + kZeroDensityContextLimit - kZeroDensityContextCount);

      max_num_bits_ac =

          std::max(max_num_bits_ac, dec_state_->code[i].max_num_bits);

    }

    max_num_bits_ac += CeilLog2Nonzero(frame_header_.passes.num_passes);

    // 16-bit buffer for decoding to JPEG are not implemented.

    // TODO(veluca): figure out the exact limit - 16 should still work with

    // 16-bit buffers, but we are excluding it for safety.

    bool use_16_bit = max_num_bits_ac < 16 && !decoded_->IsJPEG();

    bool store = frame_header_.passes.num_passes > 1;

    size_t xs = store ? kGroupDim * kGroupDim : 0;

    size_t ys = store ? frame_dim_.num_groups : 0;

    if (use_16_bit) {

      JXL_ASSIGN_OR_RETURN(dec_state_->coefficients,

                           ACImageT<int16_t>::Make(memory_manager, xs, ys));

    } else {

      JXL_ASSIGN_OR_RETURN(dec_state_->coefficients,

                           ACImageT<int32_t>::Make(memory_manager, xs, ys));

    }

    if (store) {

      dec_state_->coefficients->ZeroFill();

    }

  }

  // Set JPEG decoding data.

  if (decoded_->IsJPEG()) {

    decoded_->color_transform = frame_header_.color_transform;

    decoded_->chroma_subsampling = frame_header_.chroma_subsampling;

    const std::vector<QuantEncoding>& qe =

        dec_state_->shared_storage.matrices.encodings();

    if (qe.empty() || qe[0].mode != QuantEncoding::Mode::kQuantModeRAW ||

        std::abs(qe[0].qraw.qtable_den - 1.f / (8 * 255)) > 1e-8f) {

      return JXL_FAILURE(

          "Quantization table is not a JPEG quantization table.");

    }

    jpeg::JPEGData* jpeg_data = decoded_->jpeg_data.get();

    size_t num_components = jpeg_data->components.size();

    bool is_gray = (num_components == 1);

    JXL_ENSURE(frame_header_.color_transform != ColorTransform::kXYB);

    auto jpeg_c_map = JpegOrder(frame_header_.color_transform, is_gray);

    size_t qt_set = 0;

    JXL_ENSURE(num_components <= 3);

    JXL_ENSURE(qe[0].qraw.qtable->size() == 3 * 8 * 8);

    int* qtable = qe[0].qraw.qtable->data();

    for (size_t c = 0; c < num_components; c++) {

      // TODO(eustas): why 1-st quant table for gray?

      size_t quant_c = is_gray ? 1 : c;

      size_t qpos = jpeg_data->components[jpeg_c_map[c]].quant_idx;

      JXL_ENSURE(qpos != jpeg_data->quant.size());

      qt_set |= 1 << qpos;

      for (size_t x = 0; x < 8; x++) {

        for (size_t y = 0; y < 8; y++) {

          jpeg_data->quant[qpos].values[x * 8 + y] =

              qtable[quant_c * 64 + y * 8 + x];

        }

      }

    }

    for (size_t i = 0; i < jpeg_data->quant.size(); i++) {

      if (qt_set & (1 << i)) continue;

      if (i == 0) return JXL_FAILURE("First quant table unused.");

      // Unused quant table is set to copy of previous quant table

      for (size_t j = 0; j < 64; j++) {

        jpeg_data->quant[i].values[j] = jpeg_data->quant[i - 1].values[j];

      }

    }

  }

  decoded_ac_global_ = true;

  return true;

}

Status FrameDecoder::ProcessACGroup(size_t ac_group_id, PassesReaders& br,

                                    size_t num_passes, size_t thread,

                                    bool force_draw, bool dc_only) {

  size_t group_dim = frame_dim_.group_dim;

  const size_t gx = ac_group_id % frame_dim_.xsize_groups;

  const size_t gy = ac_group_id / frame_dim_.xsize_groups;

  const size_t x = gx * group_dim;

  const size_t y = gy * group_dim;

  JxlMemoryManager* memory_manager = dec_state_->memory_manager();

  JXL_DEBUG_V(3,

              "Processing AC group %" PRIuS "(%" PRIuS ",%" PRIuS

              ") group_dim: %" PRIuS " decoded passes: %u new passes: %" PRIuS,

              ac_group_id, gx, gy, group_dim,

              decoded_passes_per_ac_group_[ac_group_id], num_passes);

  RenderPipelineInput render_pipeline_input =

      dec_state_->render_pipeline->GetInputBuffers(ac_group_id, thread);

  bool should_run_pipeline = true;

  if (frame_header_.encoding == FrameEncoding::kVarDCT) {

    JXL_RETURN_IF_ERROR(group_dec_caches_[thread].InitOnce(

        memory_manager, frame_header_.passes.num_passes, dec_state_->used_acs));

    JXL_RETURN_IF_ERROR(DecodeGroup(

        frame_header_, br.data(), num_passes, ac_group_id, dec_state_,

        &group_dec_caches_[thread], thread, render_pipeline_input,

        decoded_->jpeg_data.get(), decoded_passes_per_ac_group_[ac_group_id],

        force_draw, dc_only, &should_run_pipeline));

  }

  // don't limit to image dimensions here (is done in DecodeGroup)

  const Rect mrect(x, y, group_dim, group_dim);

  bool modular_ready = false;

  size_t pass0 = decoded_passes_per_ac_group_[ac_group_id];

  size_t pass1 =

      force_draw ? frame_header_.passes.num_passes : pass0 + num_passes;

  for (size_t i = pass0; i < pass1; ++i) {

    int minShift;

    int maxShift;

    frame_header_.passes.GetDownsamplingBracket(i, minShift, maxShift);

    bool modular_pass_ready = true;

    JXL_DEBUG_V(2, "Decoding modular in group %d pass %d",

                static_cast<int>(ac_group_id), static_cast<int>(i));

    if (i < pass0 + num_passes) {  // i.e. i - pass0 < num_passes

      BitReader* r = br[i - pass0];

      JXL_ENSURE(r);

      JXL_DEBUG_V(2, "Bit reader position: %" PRIuS " / %" PRIuS,

                  r->TotalBitsConsumed(), r->TotalBytes() * kBitsPerByte);

      JXL_RETURN_IF_ERROR(modular_frame_decoder_.DecodeGroup(

          frame_header_, mrect, r, minShift, maxShift,

          ModularStreamId::ModularAC(ac_group_id, i),

          /*zerofill=*/false, dec_state_, &render_pipeline_input,

          /*allow_truncated=*/false, &modular_pass_ready));

    } else {

      JXL_RETURN_IF_ERROR(modular_frame_decoder_.DecodeGroup(

          frame_header_, mrect, nullptr, minShift, maxShift,

          ModularStreamId::ModularAC(ac_group_id, i), /*zerofill=*/true,

          dec_state_, &render_pipeline_input,

          /*allow_truncated=*/false, &modular_pass_ready));

    }

    if (modular_pass_ready) modular_ready = true;

  }

  decoded_passes_per_ac_group_[ac_group_id] += num_passes;

  if ((frame_header_.flags & FrameHeader::kNoise) != 0) {

    PrepareNoiseInput(*dec_state_, frame_dim_, frame_header_, ac_group_id,

                      thread);

  }

  if (!modular_frame_decoder_.UsesFullImage() && !decoded_->IsJPEG()) {

    if (should_run_pipeline && modular_ready) {

      JXL_RETURN_IF_ERROR(render_pipeline_input.Done());

    } else if (force_draw) {

      return JXL_FAILURE("Modular group decoding failed.");

    }

  }

  return true;

}

void FrameDecoder::MarkSections(const SectionInfo* sections, size_t num,

                                const SectionStatus* section_status) {

  num_sections_done_ += num;

  for (size_t i = 0; i < num; i++) {

    if (section_status[i] != SectionStatus::kDone) {

      processed_section_[sections[i].id] = JXL_FALSE;

      num_sections_done_--;

    }

  }

}

Status FrameDecoder::ProcessSections(const SectionInfo* sections, size_t num,

                                     SectionStatus* section_status) {

  if (num == 0) return true;  // Nothing to process

  std::fill(section_status, section_status + num, SectionStatus::kSkipped);

  size_t dc_global_sec = num;

  size_t ac_global_sec = num;

  std::vector<size_t> dc_group_sec(frame_dim_.num_dc_groups, num);

  std::vector<std::vector<size_t>> ac_group_sec(

      frame_dim_.num_groups,

      std::vector<size_t>(frame_header_.passes.num_passes, num));

  // This keeps track of the number of ac passes we want to process during this

  // call of ProcessSections.

  std::vector<size_t> desired_num_ac_passes(frame_dim_.num_groups);

  bool single_section =

      frame_dim_.num_groups == 1 && frame_header_.passes.num_passes == 1;

  if (single_section) {

    JXL_ENSURE(num == 1);

    JXL_ENSURE(sections[0].id == 0);

    if (processed_section_[0] == JXL_FALSE) {

      processed_section_[0] = JXL_TRUE;

      ac_group_sec[0].resize(1);

      dc_global_sec = ac_global_sec = dc_group_sec[0] = ac_group_sec[0][0] = 0;

      desired_num_ac_passes[0] = 1;

    } else {

      section_status[0] = SectionStatus::kDuplicate;

    }

  } else {

    size_t ac_global_index = frame_dim_.num_dc_groups + 1;

    for (size_t i = 0; i < num; i++) {

      JXL_ENSURE(sections[i].id < processed_section_.size());

      if (processed_section_[sections[i].id]) {

        section_status[i] = SectionStatus::kDuplicate;

        continue;

      }

      if (sections[i].id == 0) {

        dc_global_sec = i;

      } else if (sections[i].id < ac_global_index) {

        dc_group_sec[sections[i].id - 1] = i;

      } else if (sections[i].id == ac_global_index) {

        ac_global_sec = i;

      } else {

        size_t ac_idx = sections[i].id - ac_global_index - 1;

        size_t acg = ac_idx % frame_dim_.num_groups;

        size_t acp = ac_idx / frame_dim_.num_groups;

        if (acp >= frame_header_.passes.num_passes) {

          return JXL_FAILURE("Invalid section ID");

        }

        ac_group_sec[acg][acp] = i;

      }

      processed_section_[sections[i].id] = JXL_TRUE;

    }

    // Count number of new passes per group.

    for (size_t g = 0; g < ac_group_sec.size(); g++) {

      size_t j = 0;

      for (; j + decoded_passes_per_ac_group_[g] <

             frame_header_.passes.num_passes;

           j++) {

        if (ac_group_sec[g][j + decoded_passes_per_ac_group_[g]] == num) {

          break;

        }

      }

      desired_num_ac_passes[g] = j;

    }

  }

  if (dc_global_sec != num) {

    Status dc_global_status = ProcessDCGlobal(sections[dc_global_sec].br);

    if (dc_global_status.IsFatalError()) return dc_global_status;

    if (dc_global_status) {

      section_status[dc_global_sec] = SectionStatus::kDone;

    } else {

      section_status[dc_global_sec] = SectionStatus::kPartial;

    }

  }

  if (decoded_dc_global_) {

    const auto process_section = [this, &dc_group_sec, &num, &sections,

                                  &section_status](size_t i,

                                                   size_t thread) -> Status {

      if (dc_group_sec[i] != num) {

        JXL_RETURN_IF_ERROR(ProcessDCGroup(i, sections[dc_group_sec[i]].br));

        section_status[dc_group_sec[i]] = SectionStatus::kDone;

      }

      return true;

    };

    JXL_RETURN_IF_ERROR(RunOnPool(pool_, 0, dc_group_sec.size(),

                                  ThreadPool::NoInit, process_section,

                                  "DecodeDCGroup"));

  }

  if (!HasDcGroupToDecode() && !finalized_dc_) {

    PassesDecoderState::PipelineOptions pipeline_options;

    pipeline_options.use_slow_render_pipeline = use_slow_rendering_pipeline_;

    pipeline_options.coalescing = coalescing_;

    pipeline_options.render_spotcolors = render_spotcolors_;

    pipeline_options.render_noise = true;

    JXL_RETURN_IF_ERROR(dec_state_->PreparePipeline(

        frame_header_, &frame_header_.nonserialized_metadata->m, decoded_,

        pipeline_options));

    JXL_RETURN_IF_ERROR(FinalizeDC());

    JXL_RETURN_IF_ERROR(AllocateOutput());

    if (progressive_detail_ >= JxlProgressiveDetail::kDC) {

      MarkSections(sections, num, section_status);

      return true;

    }

  }

  if (finalized_dc_ && ac_global_sec != num && !decoded_ac_global_) {

    JXL_RETURN_IF_ERROR(ProcessACGlobal(sections[ac_global_sec].br));

    section_status[ac_global_sec] = SectionStatus::kDone;

  }

  if (progressive_detail_ >= JxlProgressiveDetail::kLastPasses) {

    // Mark that we only want the next progression pass.

    size_t target_complete_passes = NextNumPassesToPause();

    for (size_t i = 0; i < ac_group_sec.size(); i++) {

      desired_num_ac_passes[i] =

          std::min(desired_num_ac_passes[i],

                   target_complete_passes - decoded_passes_per_ac_group_[i]);

    }

  }

  if (decoded_ac_global_) {

    // Mark all the AC groups that we received as not complete yet.

    for (size_t i = 0; i < ac_group_sec.size(); i++) {

      if (desired_num_ac_passes[i] != 0) {

        dec_state_->render_pipeline->ClearDone(i);

      }

    }

    const auto prepare_storage = [this](size_t num_threads) -> Status {

      JXL_RETURN_IF_ERROR(

          PrepareStorage(num_threads, decoded_passes_per_ac_group_.size()));

      return true;

    };

    const auto process_group = [this, &ac_group_sec, &desired_num_ac_passes,

                                &num, &sections, &section_status](

                                   size_t g, size_t thread) -> Status {

      if (desired_num_ac_passes[g] == 0) {

        // no new AC pass, nothing to do

        return true;

      }

      size_t first_pass = decoded_passes_per_ac_group_[g];

      PassesReaders readers = {};

      for (size_t i = 0; i < desired_num_ac_passes[g]; i++) {

        JXL_ENSURE(ac_group_sec[g][first_pass + i] != num);

        readers[i] = sections[ac_group_sec[g][first_pass + i]].br;

      }

      JXL_RETURN_IF_ERROR(ProcessACGroup(

          g, readers, desired_num_ac_passes[g], GetStorageLocation(thread, g),

          /*force_draw=*/false, /*dc_only=*/false));

      for (size_t i = 0; i < desired_num_ac_passes[g]; i++) {

        section_status[ac_group_sec[g][first_pass + i]] = SectionStatus::kDone;

      }

      return true;

    };

    JXL_RETURN_IF_ERROR(RunOnPool(pool_, 0, ac_group_sec.size(),

                                  prepare_storage, process_group,

                                  "DecodeGroup"));

  }

  MarkSections(sections, num, section_status);

  return true;

}

Status FrameDecoder::Flush() {

  bool has_blending = frame_header_.blending_info.mode != BlendMode::kReplace ||

                      frame_header_.custom_size_or_origin;

  for (const auto& blending_info_ec :

       frame_header_.extra_channel_blending_info) {

    if (blending_info_ec.mode != BlendMode::kReplace) has_blending = true;

  }

  // No early Flush() if blending is enabled.

  if (has_blending && !is_finalized_) {

    return false;

  }

  // No early Flush() - nothing to do - if the frame is a kSkipProgressive

  // frame.

  if (frame_header_.frame_type == FrameType::kSkipProgressive &&

      !is_finalized_) {

    return true;

  }

  if (decoded_->IsJPEG()) {

    // Nothing to do.

    return true;

  }

  JXL_RETURN_IF_ERROR(AllocateOutput());

  uint32_t completely_decoded_ac_pass = *std::min_element(

      decoded_passes_per_ac_group_.begin(), decoded_passes_per_ac_group_.end());

  if (completely_decoded_ac_pass < frame_header_.passes.num_passes) {

    // We don't have all AC yet: force a draw of all the missing areas.

    // Mark all sections as not complete.

    for (size_t i = 0; i < decoded_passes_per_ac_group_.size(); i++) {

      if (decoded_passes_per_ac_group_[i] < frame_header_.passes.num_passes) {

        dec_state_->render_pipeline->ClearDone(i);

      }

    }

    const auto prepare_storage = [this](const size_t num_threads) -> Status {

      JXL_RETURN_IF_ERROR(

          PrepareStorage(num_threads, decoded_passes_per_ac_group_.size()));

      return true;

    };

    const auto process_group = [this](const uint32_t g,

                                      size_t thread) -> Status {

      if (decoded_passes_per_ac_group_[g] == frame_header_.passes.num_passes) {

        // This group was drawn already, nothing to do.

        return true;

      }

      PassesReaders readers = {};

      JXL_RETURN_IF_ERROR(ProcessACGroup(

          g, readers, /*num_passes=*/0, GetStorageLocation(thread, g),

          /*force_draw=*/true, /*dc_only=*/!decoded_ac_global_));

      return true;

    };

    JXL_RETURN_IF_ERROR(RunOnPool(pool_, 0, decoded_passes_per_ac_group_.size(),

                                  prepare_storage, process_group,

                                  "ForceDrawGroup"));

  }

  // undo global modular transforms and copy int pixel buffers to float ones

  JXL_RETURN_IF_ERROR(modular_frame_decoder_.FinalizeDecoding(

      frame_header_, dec_state_, pool_, is_finalized_));

  return true;

}

int FrameDecoder::SavedAs(const FrameHeader& header) {

  if (header.frame_type == FrameType::kDCFrame) {

    // bits 16, 32, 64, 128 for DC level

    return 16 << (header.dc_level - 1);

  } else if (header.CanBeReferenced()) {

    // bits 1, 2, 4 and 8 for the references

    return 1 << header.save_as_reference;

  }

  return 0;

}

bool FrameDecoder::HasEverything() const {

  if (!decoded_dc_global_) return false;

  if (!decoded_ac_global_) return false;

  if (HasDcGroupToDecode()) return false;

  for (const auto& nb_passes : decoded_passes_per_ac_group_) {

    if (nb_passes < frame_header_.passes.num_passes) return false;

  }

  return true;

}

int FrameDecoder::References() const {

  if (is_finalized_) {

    return 0;

  }

  if (!HasEverything()) return 0;

  int result = 0;

  // Blending

  if (frame_header_.frame_type == FrameType::kRegularFrame ||

      frame_header_.frame_type == FrameType::kSkipProgressive) {

    bool cropped = frame_header_.custom_size_or_origin;

    if (cropped || frame_header_.blending_info.mode != BlendMode::kReplace) {

      result |= (1 << frame_header_.blending_info.source);

    }

    const auto& extra = frame_header_.extra_channel_blending_info;

    for (const auto& ecbi : extra) {

      if (cropped || ecbi.mode != BlendMode::kReplace) {

        result |= (1 << ecbi.source);

      }

    }

  }

  // Patches

  if (frame_header_.flags & FrameHeader::kPatches) {

    result |= dec_state_->shared->image_features.patches.GetReferences();

  }

  // DC Level

  if (frame_header_.flags & FrameHeader::kUseDcFrame) {

    // Reads from the next dc level

    int dc_level = frame_header_.dc_level + 1;

    // bits 16, 32, 64, 128 for DC level

    result |= (16 << (dc_level - 1));

  }

  return result;

}

Status FrameDecoder::FinalizeFrame() {

  if (is_finalized_) {

    return JXL_FAILURE("FinalizeFrame called multiple times");

  }

  is_finalized_ = true;

  if (decoded_->IsJPEG()) {

    // Nothing to do.

    return true;

  }

  // undo global modular transforms and copy int pixel buffers to float ones

  JXL_RETURN_IF_ERROR(

      modular_frame_decoder_.FinalizeDecoding(frame_header_, dec_state_, pool_,

                                              /*inplace=*/true));

  if (frame_header_.CanBeReferenced()) {

    auto& info = dec_state_->shared_storage

                     .reference_frames[frame_header_.save_as_reference];

    *info.frame = std::move(dec_state_->frame_storage_for_referencing);

    info.ib_is_in_xyb = frame_header_.save_before_color_transform;

  }

  return true;

}

}  // namespace jxl