Method ProcessNode

Filters/HyperTree/vtkHyperTreeGridEvaluateCoarse.cxx:142–187

------------------------------------------------------------------------------

//------------------------------------------------------------------------------
void vtkHyperTreeGridEvaluateCoarse::ProcessNode(vtkHyperTreeGridNonOrientedCursor* outCursor)
{
  if (this->CheckAbort())
  {
    return;
  }

  vtkIdType currentId = outCursor->GetGlobalNodeIndex();

  // Leaf/Masked cell: data does not change
  if (outCursor->IsLeaf() || outCursor->IsMasked())
  {
    this->OutData->CopyData(this->InData, currentId, currentId);
    return;
  }

  // Coarse cell: recurse and retrieve values
  int nbArray = this->InData->GetNumberOfArrays();
  std::vector<std::vector<std::vector<double>>> childrenValues(outCursor->GetNumberOfChildren());

  for (int ichild = 0; ichild < outCursor->GetNumberOfChildren(); ++ichild)
  {
    this->ProcessChild(outCursor, ichild, childrenValues[ichild]);
  }

  // Reduction operation over the resulting array
  for (int arrayId = 0; arrayId < nbArray; ++arrayId)
  {
    if (this->CheckAbort())
    {
      break;
    }
    vtkDataArray* arr = this->OutData->GetArray(arrayId);
    int nbComponents = arr->GetNumberOfComponents();
    for (int componentID = 0; componentID < nbComponents; ++componentID)
    {
      // Copy child values from every child array
      std::vector<double> childVals(outCursor->GetNumberOfChildren(), 0.0);
      for (unsigned char ichild = 0; ichild < outCursor->GetNumberOfChildren(); ichild++)
      {
        childVals[ichild] = childrenValues[ichild][arrayId][componentID];
      }
      arr->SetComponent(currentId, componentID, EvalCoarse(childVals));
    }
  }
}

//------------------------------------------------------------------------------
void vtkHyperTreeGridEvaluateCoarse::ProcessChild(vtkHyperTreeGridNonOrientedCursor* outCursor,