ttkIcospheresFromPoints.cpp

Go to the documentation of this file.

#include <ttkIcospheresFromPoints.h>
 
#include <vtkInformation.h>
#include <vtkObjectFactory.h>
 
#include <vtkPointData.h>
#include <vtkPointSet.h>
 
#include <ttkUtils.h>
 
vtkStandardNewMacro(ttkIcospheresFromPoints);
 
ttkIcospheresFromPoints::ttkIcospheresFromPoints() : ttkIcosphere() {
  this->setDebugMsgPrefix("IcospheresFromPoints");
  this->SetNumberOfInputPorts(1);
}
ttkIcospheresFromPoints::~ttkIcospheresFromPoints() = default;
 
int ttkIcospheresFromPoints::FillInputPortInformation(int port,
                                                      vtkInformation *info) {
  if(port == 0)
    info->Set(vtkAlgorithm::INPUT_REQUIRED_DATA_TYPE(), "vtkPointSet");
  else
    return 0;
  return 1;
}
 
template <typename VTK_TT>
int copyArrayData(vtkDataArray *oldArray,
                  vtkDataArray *newArray,
                  const size_t &nSpheres,
                  const size_t &nVerticesPerSphere,
                  const size_t &nComponents) {
  auto oldData = ttkUtils::GetPointer<VTK_TT>(oldArray);
  auto newData = ttkUtils::GetPointer<VTK_TT>(newArray);
 
  for(size_t i = 0; i < nSpheres; i++) {
    size_t sphereIndex = i * nVerticesPerSphere * nComponents;
    for(size_t j = 0; j < nComponents; j++) {
      const auto &value = oldData[i * nComponents + j];
      size_t newIndex = sphereIndex + j;
      for(size_t k = 0; k < nVerticesPerSphere; k++) {
        newData[newIndex] = value;
        newIndex += nComponents;
      }
    }
  }
  return 1;
}
 
int ttkIcospheresFromPoints::RequestData(vtkInformation *request,
                                         vtkInformationVector **inputVector,
                                         vtkInformationVector *outputVector) {
  auto input = vtkPointSet::GetData(inputVector[0], 0);
  if(!input)
    return 1;
  size_t nPoints = input->GetNumberOfPoints();
  if(nPoints < 1)
    return 1;
 
  this->SetCenters(input->GetPoints()->GetData());
 
  // compute spheres
  int status
    = this->ttkIcosphere::RequestData(request, inputVector, outputVector);
  if(!status)
    return 0;
 
  size_t nVertices = 0;
  size_t nTriangles = 0;
  this->computeNumberOfVerticesAndTriangles(
    nVertices, nTriangles, this->GetNumberOfSubdivisions());
 
  auto output = vtkDataSet::GetData(outputVector);
  auto outputPD = output->GetPointData();
 
  // copy point data
  if(this->CopyPointData) {
 
    auto inputPD = input->GetPointData();
    for(size_t i = 0, n = inputPD->GetNumberOfArrays(); i < n; i++) {
      auto oldArray = vtkDataArray::SafeDownCast(inputPD->GetAbstractArray(i));
      if(!oldArray || oldArray->GetName() == nullptr) {
        continue;
      }
      std::string oldArrayName(oldArray->GetName());
      if(this->GetComputeNormals() && oldArrayName.compare("Normals") == 0) {
        continue;
      }
      auto newArray
        = vtkSmartPointer<vtkDataArray>::Take(oldArray->NewInstance());
      size_t nComponents = oldArray->GetNumberOfComponents();
      newArray->SetName(oldArray->GetName());
      newArray->SetNumberOfComponents(nComponents);
      newArray->SetNumberOfTuples(nVertices * nPoints);
 
      switch(newArray->GetDataType()) {
        vtkTemplateMacro((copyArrayData<VTK_TT>(
          oldArray, newArray, nPoints, nVertices, nComponents)));
      }
      outputPD->AddArray(newArray);
    }
  }
 
  return 1;
}