// SPDX-FileCopyrightText: Copyright (c) Ken Martin, Will Schroeder, Bill Lorensen
// SPDX-License-Identifier: BSD-3-Clause

#include "vtkDataObjectMeshCache.h"

#include "vtkCellData.h"
#include "vtkCompositeDataSet.h"
#include "vtkDataArrayRange.h"
#include "vtkDataObjectTree.h"
#include "vtkDataObjectTreeRange.h"
#include "vtkDataSet.h"
#include "vtkPointData.h"

VTK_ABI_NAMESPACE_BEGIN
vtkStandardNewMacro(vtkDataObjectMeshCache);

/**
 * Interface to dispatch work over every contained vtkDataSet.
 * If input is a vtkDataSet subclass, forward it directly to
 * ComputeDataSet.
 * If input is a vtkDataObjectTree subclass, iterate over
 * inner non empty vtkDataSet leaves.
 */
struct GenericDataObjectWorker
{
  virtual ~GenericDataObjectWorker() = default;

  /**
   * Entry point. In the end, call ComputeDataSet for every
   * contained vtkDataSet.
   */
  void Compute(vtkDataObject* dataobject)
  {
    auto dataset = vtkDataSet::SafeDownCast(dataobject);
    if (dataset)
    {
      this->ComputeDataSet(dataset);
      return;
    }
    auto composite = vtkDataObjectTree::SafeDownCast(dataobject);
    if (composite)
    {
      this->ComputeComposite(composite);
      return;
    }

    this->SkippedData = true;
  }

  /**
   * Iterate over inner vtkDataSet to call ComputeDataSet
   */
  void ComputeComposite(vtkDataObjectTree* composite)
  {
    auto options = vtk::DataObjectTreeOptions::TraverseSubTree |
      vtk::DataObjectTreeOptions::SkipEmptyNodes | vtk::DataObjectTreeOptions::VisitOnlyLeaves;
    for (auto dataLeaf : vtk::Range(composite, options))
    {
      auto dataset = vtkDataSet::SafeDownCast(dataLeaf);
      if (dataset)
      {
        this->ComputeDataSet(dataset);
      }
      else
      {
        this->SkippedData = true;
      }
    }
  }

  /**
   * To be reimplemented to do the actual work.
   * Will be called multiple times for composite.
   */
  virtual void ComputeDataSet(vtkDataSet* dataset) = 0;

  bool SkippedData = false;
};

/**
 * Worker to compute mesh mtime.
 * For composite, return the max value.
 */
struct MeshMTimeWorker : public GenericDataObjectWorker
{
  ~MeshMTimeWorker() override = default;

  void ComputeDataSet(vtkDataSet* dataset) override
  {
    this->MeshTime = std::max(this->MeshTime, dataset->GetMeshMTime());
  }

  vtkMTimeType MeshTime = 0;
};

/**
 * Worker to verify if data is supported.
 * If input is not a dataset, Supported() will return false.
 * If any inner dataset is unsupported, Supported() will return false.
 * Otherwise, Supported() return true;
 */
struct SupportedDataWorker : public GenericDataObjectWorker
{
  ~SupportedDataWorker() override = default;

  bool Supported() { return !this->SkippedData; }

  void ComputeDataSet(vtkDataSet* vtkNotUsed(dataset)) override {}
};

/**
 * Worker to verify that data has requested arrays.
 * OriginalIdsName is the map of requested array names per attribute type.
 * HasRequestedIds is set to false if a requested array is not found.
 */
struct RequestedIdsWorker : public GenericDataObjectWorker
{
  ~RequestedIdsWorker() override = default;

  void ComputeDataSet(vtkDataSet* dataset) override
  {
    for (const auto& attribute : this->OriginalIdsName)
    {
      vtkDataSetAttributes* field = dataset->GetAttributes(attribute.first);
      if (!field)
      {
        this->HasRequestedIds = false;
        continue;
      }

      if (!field->GetArray(attribute.second.c_str()))
      {
        this->HasRequestedIds = false;
      }
    }
  }

  std::map<int, std::string> OriginalIdsName;
  bool HasRequestedIds = true;
};

/**
 * Worker to clear dataset attributes from data.
 */
struct ClearAttributesWorker : public GenericDataObjectWorker
{
  ~ClearAttributesWorker() override = default;

  void ComputeDataSet(vtkDataSet* dataset) override
  {
    for (int attribute = vtkDataObject::POINT; attribute < vtkDataObject::NUMBER_OF_ATTRIBUTE_TYPES;
         attribute++)
    {
      vtkFieldData* field = dataset->GetAttributesAsFieldData(attribute);
      if (field)
      {
        field->Initialize();
      }
    }
  }
};

/**
 * Worker to count number of datasets.
 */
struct NumberOfDataSetWorker : public GenericDataObjectWorker
{
  ~NumberOfDataSetWorker() override = default;

  void ComputeDataSet(vtkDataSet* vtkNotUsed(dataset)) override { this->NumberOfDataSets++; }
  vtkIdType NumberOfDataSets = 0;
};

//------------------------------------------------------------------------------
void vtkDataObjectMeshCache::PrintSelf(ostream& os, vtkIndent indent)
{
  this->Superclass::PrintSelf(os, indent);

  os << indent << "Cache:";
  if (this->Cache)
  {
    os << endl;
    this->Cache->PrintSelf(os, indent.GetNextIndent());
  }
  else
  {
    os << "(none)\n";
  }

  os << indent << "CachedOriginalMeshTime: " << this->CachedOriginalMeshTime << "\n";
  os << indent << "CachedConsumerTime: " << this->CachedConsumerTime << "\n";

  os << indent << "OriginalIdsName:\n";
  for (const auto& attribute : this->OriginalIdsName)
  {
    os << indent.GetNextIndent() << vtkDataObject::GetAssociationTypeAsString(attribute.first)
       << " " << attribute.second << "\n";
  }

  Status status = this->GetStatus();
  os << indent << "Status\n:";
  status.PrintSelf(os, indent.GetNextIndent());
}

//------------------------------------------------------------------------------
bool vtkDataObjectMeshCache::IsSupportedData(vtkDataObject* dataobject) const
{
  SupportedDataWorker supportWorker;
  supportWorker.Compute(dataobject);

  vtkDebugMacro(" return IsSupportedData: " << supportWorker.Supported());
  return supportWorker.Supported();
}

//------------------------------------------------------------------------------
void vtkDataObjectMeshCache::SetOriginalDataObject(vtkDataObject* input)
{
  if (!input)
  {
    vtkWarningMacro("Invalid original dataobject: nullptr");
    return;
  }

  if (this->IsSupportedData(input))
  {
    this->OriginalDataSet = vtkDataSet::SafeDownCast(input);
    this->OriginalCompositeDataSet = vtkCompositeDataSet::SafeDownCast(input);
    vtkDebugMacro(" set OriginalDataObject: " << input);
    this->Modified();
    return;
  }

  // Clear existing dataset ptrs
  this->OriginalCompositeDataSet = nullptr;
  this->OriginalDataSet = nullptr;

  if (vtkCompositeDataSet::SafeDownCast(input))
  {
    vtkWarningMacro("Composite " << input->GetClassName() << " has unsupported block(s).");
  }
  else
  {
    vtkWarningMacro("Unsupported input type: " << input->GetClassName());
  }
}

//------------------------------------------------------------------------------
void vtkDataObjectMeshCache::ClearOriginalIds()
{
  this->OriginalIdsName.clear();
  vtkDebugMacro(" clear OriginalIdsName");
  this->Modified();
}

//------------------------------------------------------------------------------
void vtkDataObjectMeshCache::AddOriginalIds(int attribute, const std::string& name)
{
  if (attribute < 0 || attribute >= vtkDataObject::NUMBER_OF_ATTRIBUTE_TYPES)
  {
    vtkWarningMacro("Invalid attribute type: " << attribute);
    return;
  }

  this->OriginalIdsName[attribute] = name;
  vtkDebugMacro(" set OriginalIds: " << attribute << " array name to " << name.c_str());
  this->Modified();
}

//------------------------------------------------------------------------------
void vtkDataObjectMeshCache::RemoveOriginalIds(int attribute)
{
  if (attribute < 0 || attribute >= vtkDataObject::NUMBER_OF_ATTRIBUTE_TYPES)
  {
    vtkWarningMacro("Invalid attribute type: " << attribute);
    return;
  }

  this->OriginalIdsName.erase(attribute);
  vtkDebugMacro(" remove OriginalIdsName: " << attribute);
  this->Modified();
}

//------------------------------------------------------------------------------
void vtkDataObjectMeshCache::UpdateCache(vtkDataObject* output)
{
  if (!output)
  {
    vtkWarningMacro("Cannot update from nullptr");
    return;
  }

  bool isSupported = this->IsSupportedData(output);
  if (!isSupported)
  {
    vtkWarningMacro("Cannot update from unsupported data type: " << output->GetClassName());
    return;
  }

  this->Cache.TakeReference(output->NewInstance());
  this->Cache->ShallowCopy(output);
  this->CachedOriginalMeshTime = this->GetOriginalMeshTime();
  this->CachedConsumerTime = this->Consumer->GetMTime();

  vtkDebugMacro(" update Cache: " << this->Cache.GetPointer());
  this->Modified();
}

//------------------------------------------------------------------------------
void vtkDataObjectMeshCache::InvalidateCache()
{
  this->Cache = nullptr;
  this->CachedOriginalMeshTime = 0;
  this->CachedConsumerTime = 0;
  vtkDebugMacro(" invalidate Cache");
  this->Modified();
}

//------------------------------------------------------------------------------
vtkIdType vtkDataObjectMeshCache::GetNumberOfDataSets(vtkDataObject* dataobject) const
{
  NumberOfDataSetWorker countWorker;
  countWorker.Compute(dataobject);
  return countWorker.NumberOfDataSets;
}

//------------------------------------------------------------------------------
vtkMTimeType vtkDataObjectMeshCache::GetOriginalMeshTime() const
{
  MeshMTimeWorker meshtime;
  meshtime.Compute(this->GetOriginalDataObject());
  return meshtime.MeshTime;
}

//------------------------------------------------------------------------------
vtkDataObject* vtkDataObjectMeshCache::GetOriginalDataObject() const
{
  if (this->OriginalDataSet)
  {
    return this->OriginalDataSet;
  }

  if (this->OriginalCompositeDataSet)
  {
    return this->OriginalCompositeDataSet;
  }

  return nullptr;
}

//------------------------------------------------------------------------------
void vtkDataObjectMeshCache::Status::PrintSelf(ostream& os, vtkIndent indent)
{
  os << indent << "OriginalDataDefined: " << this->OriginalDataDefined << "\n";
  os << indent << "ConsumerDefined: " << this->ConsumerDefined << "\n";
  os << indent << "CacheDefined: " << this->CacheDefined << "\n";
  os << indent << "OriginalMeshUnmodified: " << this->OriginalMeshUnmodified << "\n";
  os << indent << "ConsumerUnmodified: " << this->ConsumerUnmodified << "\n";
  os << indent << "AttributesIdsExists: " << this->AttributesIdsExists << "\n";
}

//------------------------------------------------------------------------------
vtkDataObjectMeshCache::Status vtkDataObjectMeshCache::GetStatus() const
{
  Status status;
  status.OriginalDataDefined =
    this->OriginalDataSet != nullptr || this->OriginalCompositeDataSet != nullptr;
  if (!status.OriginalDataDefined)
  {
    vtkDebugMacro("OriginalDataObject is not set.");
  }

  status.ConsumerDefined = this->Consumer != nullptr;
  if (!status.ConsumerDefined)
  {
    vtkDebugMacro("Consumer is nullptr.");
  }

  status.CacheDefined = this->Cache != nullptr;
  if (!status.CacheDefined)
  {
    vtkDebugMacro("Cache is uninitialized.");
    return status;
  }

  status.ConsumerUnmodified = this->Consumer->GetMTime() <= this->CachedConsumerTime;
  if (!status.ConsumerUnmodified)
  {
    vtkDebugMacro("Consumer modification time has changed.");
  }

  // Be sure that original data and cache are of same class.
  if (this->OriginalCompositeDataSet)
  {
    status.OriginalMeshUnmodified =
      strcmp(this->OriginalCompositeDataSet->GetClassName(), this->Cache->GetClassName()) == 0;
  }
  else if (this->OriginalDataSet)
  {
    status.OriginalMeshUnmodified =
      strcmp(this->OriginalDataSet->GetClassName(), this->Cache->GetClassName()) == 0;
  }

  status.OriginalMeshUnmodified &= this->GetNumberOfDataSets(this->Cache) ==
    this->GetNumberOfDataSets(this->GetOriginalDataObject());

  if (!status.OriginalMeshUnmodified)
  {
    vtkDebugMacro("Input structure has changed.");
  }

  auto originalMeshMTime = this->GetOriginalMeshTime();
  status.OriginalMeshUnmodified &= originalMeshMTime > 0;
  if (!status.OriginalMeshUnmodified)
  {
    vtkDebugMacro(
      "Invalid input mesh time. Input may be of unsupported type or has no valid mesh.");
  }

  status.OriginalMeshUnmodified &= (originalMeshMTime == this->CachedOriginalMeshTime);
  if (!status.OriginalMeshUnmodified)
  {
    vtkDebugMacro("Input mesh time has changed.");
  }

  status.AttributesIdsExists = this->CacheHasRequestedIds();
  if (!status.AttributesIdsExists)
  {
    vtkDebugMacro("Cache does not have requested ids");
  }

  vtkDebugMacro(" returning status");
  return status;
}

//------------------------------------------------------------------------------
void vtkDataObjectMeshCache::CopyCacheToDataObject(vtkDataObject* output)
{
  if (!output)
  {
    vtkWarningMacro("Cannot copy to nullptr");
    return;
  }
  if (!this->Cache)
  {
    vtkWarningMacro("Cannot copy from nullptr");
    return;
  }

  if (!this->IsSupportedData(output))
  {
    vtkWarningMacro("Cannot copy to unsupported data type: " << output->GetClassName());
    return;
  }
  vtkSmartPointer<vtkDataObject> input = nullptr;
  if (this->HasConsumerNoInputPort())
  {
    input = vtkSmartPointer<vtkDataObject>::Take(output->NewInstance());
    input->ShallowCopy(output);
  }
  else if (this->OriginalDataSet)
  {
    input = this->OriginalDataSet.Get();
  }
  else if (this->OriginalCompositeDataSet)
  {
    input = this->OriginalCompositeDataSet.Get();
  }

  vtkDebugMacro(" copy Cache to data object");
  output->ShallowCopy(this->Cache);
  this->ClearAttributes(output);

  auto outputDataSet = vtkDataSet::SafeDownCast(output);
  auto outputComposite = vtkCompositeDataSet::SafeDownCast(output);
  if (outputDataSet)
  {
    auto cacheDataSet = vtkDataSet::SafeDownCast(this->Cache);
    auto inputDataSet = vtkDataSet::SafeDownCast(input);
    this->ForwardAttributesToDataSet(inputDataSet, cacheDataSet, outputDataSet);
  }
  else if (outputComposite)
  {
    auto inputComposite = vtkCompositeDataSet::SafeDownCast(input);
    this->ForwardAttributesToComposite(inputComposite, outputComposite);
  }
}

//------------------------------------------------------------------------------
void vtkDataObjectMeshCache::ForwardAttributesToDataSet(
  vtkDataSet* input, vtkDataSet* cache, vtkDataSet* outputDataSet)
{
  if (!input || !cache)
  {
    return;
  }

  for (const auto& attribute : this->OriginalIdsName)
  {
    this->ForwardAttributes(input, cache, outputDataSet, attribute.first, attribute.second);
  }

  outputDataSet->GetFieldData()->PassData(input->GetFieldData());
}

//------------------------------------------------------------------------------
void vtkDataObjectMeshCache::ForwardAttributesToComposite(
  vtkCompositeDataSet* input, vtkCompositeDataSet* output)
{
  auto inputDataTree = vtkDataObjectTree::SafeDownCast(input);
  auto outputDataTree = vtkDataObjectTree::SafeDownCast(output);
  auto cacheDataTree = vtkDataObjectTree::SafeDownCast(this->Cache);

  if (!inputDataTree || !outputDataTree || !cacheDataTree)
  {
    vtkWarningMacro("Only vtkDataObjectTree are supported for now");
    return;
  }

  auto options = vtk::DataObjectTreeOptions::TraverseSubTree |
    vtk::DataObjectTreeOptions::SkipEmptyNodes | vtk::DataObjectTreeOptions::VisitOnlyLeaves;
  auto inputDataRange = vtk::Range(inputDataTree, options);
  auto outputDataRange = vtk::Range(outputDataTree, options);
  auto cacheDataRange = vtk::Range(cacheDataTree, options);

  auto inputBlock = inputDataRange.begin();
  auto outputBlock = outputDataRange.begin();
  auto cacheBlock = cacheDataRange.begin();
  while (outputBlock != outputDataRange.end() && inputBlock != inputDataRange.end() &&
    cacheBlock != cacheDataRange.end())
  {
    auto inputDataSet = vtkDataSet::SafeDownCast(*inputBlock);
    auto outputDataSet = vtkDataSet::SafeDownCast(*outputBlock);
    auto cacheDataSet = vtkDataSet::SafeDownCast(*cacheBlock);
    if (outputDataSet && cacheDataSet)
    {
      this->ForwardAttributesToDataSet(inputDataSet, cacheDataSet, outputDataSet);
    }

    cacheBlock++;
    outputBlock++;
    inputBlock++;
  }

  outputDataTree->GetFieldData()->PassData(inputDataTree->GetFieldData());
}

//------------------------------------------------------------------------------
void vtkDataObjectMeshCache::ForwardAttributes(
  vtkDataSet* input, vtkDataSet* cache, vtkDataSet* output, int attribute, const std::string& name)
{
  vtkDebugMacro("Forward attribute " << vtkDataObject::GetAssociationTypeAsString(attribute));

  auto inAttribute = input->GetAttributes(attribute);
  auto outAttribute = output->GetAttributes(attribute);
  auto cacheAttribute = cache->GetAttributes(attribute);

  auto originalIds = cacheAttribute->GetArray(name.c_str());
  if (!originalIds)
  {
    vtkDebugMacro(
      "Global Ids not found for " << vtkDataObject::GetAssociationTypeAsString(attribute));
    return;
  }

  outAttribute->CopyAllOn();
  outAttribute->CopyAllocate(inAttribute);

  // NOTE potential optimization:
  // this copy may be replaced by an (optional ?) use of the implicit vtkIndexedArray
  auto ptsIdsRange = vtk::DataArrayValueRange(originalIds);
  vtkIdType outId = 0;
  for (auto originalId : ptsIdsRange)
  {
    outAttribute->CopyData(inAttribute, originalId, outId);
    outId++;
  }
}

//------------------------------------------------------------------------------
bool vtkDataObjectMeshCache::CacheHasRequestedIds() const
{
  if (this->OriginalIdsName.empty())
  {
    return true;
  }

  RequestedIdsWorker idsWorker;
  idsWorker.OriginalIdsName = this->OriginalIdsName;
  idsWorker.Compute(this->Cache);
  return idsWorker.HasRequestedIds;
}

//------------------------------------------------------------------------------
void vtkDataObjectMeshCache::ClearAttributes(vtkDataObject* dataobject)
{
  ClearAttributesWorker clearWorker;
  clearWorker.Compute(dataobject);
}

//------------------------------------------------------------------------------
bool vtkDataObjectMeshCache::HasConsumerNoInputPort() const
{
  return this->Consumer->GetNumberOfInputPorts() == 0;
}

VTK_ABI_NAMESPACE_END
