ttkLDistance.cpp

Go to the documentation of this file.

#include <vtkDataArray.h>
#include <vtkDataSet.h>
#include <vtkDoubleArray.h>
#include <vtkInformation.h>
#include <vtkObjectFactory.h>
#include <vtkPointData.h>
#include <vtkSmartPointer.h>
 
#include "ttkLDistance.h"
#include <ttkUtils.h>
 
vtkStandardNewMacro(ttkLDistance);
 
ttkLDistance::ttkLDistance() {
  this->SetNumberOfInputPorts(1);
  this->SetNumberOfOutputPorts(1);
}
 
int ttkLDistance::FillInputPortInformation(int port, vtkInformation *info) {
  if(port == 0) {
    info->Set(vtkAlgorithm::INPUT_REQUIRED_DATA_TYPE(), "vtkDataSet");
    return 1;
  }
  return 0;
}
 
int ttkLDistance::FillOutputPortInformation(int port, vtkInformation *info) {
  if(port == 0) {
    info->Set(ttkAlgorithm::SAME_DATA_TYPE_AS_INPUT_PORT(), 0);
    return 1;
  }
  return 0;
}
 
int ttkLDistance::RequestData(vtkInformation *ttkNotUsed(request),
                              vtkInformationVector **inputVector,
                              vtkInformationVector *outputVector) {
 
  const auto input = vtkDataSet::GetData(inputVector[0]);
  auto output = vtkDataSet::GetData(outputVector);
 
  // Test validity of datasets (must present the same number of points).
#ifndef TTK_ENABLE_KAMIKAZE
  if(!input) {
    this->printErr("Input pointer is NULL.");
    return 0;
  }
  if(!input->GetNumberOfPoints()) {
    this->printErr("Input has no point.");
    return 0;
  }
  if(!input->GetPointData()) {
    this->printErr("Input has no point data.");
    return 0;
  }
  if(!output) {
    this->printErr("Output pointer is NULL.");
    return 0;
  }
#endif
 
  // Use a pointer-base copy for the input.
  output->ShallowCopy(input);
 
  const auto inputScalarField1 = this->GetInputArrayToProcess(0, input);
  const auto inputScalarField2 = this->GetInputArrayToProcess(1, input);
 
#ifndef TTK_ENABLE_KAMIKAZE
  if(!inputScalarField1 || !inputScalarField2
     || inputScalarField1->GetDataType() != inputScalarField2->GetDataType()) {
    this->printErr("Input scalar fields are NULL or have different types.");
    return 0;
  }
#endif
 
  // Allocate memory for the output scalar field, based on the first input.
  vtkSmartPointer<vtkDataArray> const outputScalarField{
    inputScalarField1->NewInstance()};
 
  ttk::SimplexId const numberOfPoints = input->GetNumberOfPoints();
 
  outputScalarField->SetNumberOfTuples(numberOfPoints);
  outputScalarField->SetName(DistanceFieldName.data());
  output->GetPointData()->AddArray(outputScalarField);
 
  // Calling the executing package.
  switch(inputScalarField1->GetDataType()) {
    vtkTemplateMacro(this->execute<VTK_TT>(
      static_cast<VTK_TT *>(ttkUtils::GetVoidPointer(inputScalarField1)),
      static_cast<VTK_TT *>(ttkUtils::GetVoidPointer(inputScalarField2)),
      static_cast<VTK_TT *>(ttkUtils::GetVoidPointer(outputScalarField)),
      DistanceType, numberOfPoints));
  }
 
  vtkNew<vtkDoubleArray> meanDistanceArray{};
  const std::string arrayName = "L" + DistanceType + "-distance";
  meanDistanceArray->SetName(arrayName.c_str());
  meanDistanceArray->SetNumberOfTuples(1);
  meanDistanceArray->SetTuple1(0, this->getResult());
  output->GetFieldData()->AddArray(meanDistanceArray);
 
  return 1;
}