Method operator()

Parallel/DIY/vtkDIYUtilities.cxx:75–131

  template <typename ArrayT>
  void operator()(ArrayT* array)
  {
    using ValueType = vtk::GetAPIType<ArrayT>;

    int numberOfComponents;
    vtkIdType numberOfTuples;
    std::string name;
    diy::load(this->BB, numberOfComponents);
    diy::load(this->BB, numberOfTuples);
    diy::load(this->BB, name);

    array->SetNumberOfComponents(numberOfComponents);
    array->SetNumberOfTuples(numberOfTuples);

    if (name != "_vtkArrayNoName")
    {
      array->SetName(name.c_str());
    }

    vtkNew<vtkAOSDataArrayTemplate<ValueType>> aosArray;
    if (array->HasStandardMemoryLayout())
    {
      aosArray->ShallowCopy(array);
    }
    else
    {
      // create a temporary array for loading
      aosArray->SetNumberOfComponents(numberOfComponents);
      aosArray->SetNumberOfTuples(numberOfTuples);
    }
    ValueType* data = aosArray->GetPointer(0);

    if (array->GetNumberOfValues())
    {
      diy::load(this->BB, data, array->GetNumberOfValues());
    }

    if (!array->HasStandardMemoryLayout())
    {
      // read the data into the non-standard array.
      auto range = vtk::DataArrayTupleRange(array);

      using TupleRef = typename decltype(range)::TupleReferenceType;
      using ComponentRef = typename decltype(range)::ComponentReferenceType;

      vtkIdType i(0);
      for (TupleRef tpl : range)
      {
        for (ComponentRef comp : tpl)
        {
          comp = data[i++];
        }
      }

      assert(i == numberOfComponents * numberOfTuples);
    }
  }