Function TestDataArrayIterators

Common/Core/Testing/Cxx/TestDataArrayIterators.cxx:17–120

#define SILENT

int TestDataArrayIterators(int, char*[])
{
  vtkIdType numComps = 4;
  vtkIdType numValues = 100000000; // 10 million
  assert(numValues % numComps == 0);
  vtkIdType numTuples = numValues / numComps;

  vtkNew<vtkFloatArray> arrayContainer;
  vtkFloatArray* array = arrayContainer;
  array->SetNumberOfComponents(numComps);
  array->SetNumberOfTuples(numTuples);
  for (vtkIdType i = 0; i < numValues; ++i)
  {
    // Just fill with consistent data
    array->SetValue(i, i % 97);
  }

  // Create the vtkAOSDataArrayTemplate testing implementation:
  vtkNew<vtkAOSDataArrayTemplate<float>> tdaContainer;
  vtkAOSDataArrayTemplate<float>* tda = tdaContainer;
  tda->SetNumberOfComponents(numComps);
  tda->SetArray(arrayContainer->GetPointer(0), numValues, 1);

  // should be vtkAOSDataArrayTemplate<float>::Iterator (float*):
  vtkFloatArray::Iterator datBegin = array->Begin();
  vtkFloatArray::Iterator datIter = array->Begin();
  if (typeid(datBegin) != typeid(float*))
  {
    std::cerr << "Error: vtkFloatArray::Iterator is not a float*.";
    return EXIT_FAILURE;
  }

  // should be vtkTypedDataArrayIterator<float>:
  vtkAOSDataArrayTemplate<float>::Iterator tdaBegin =
    vtkAOSDataArrayTemplate<float>::FastDownCast(tda)->Begin();
  vtkAOSDataArrayTemplate<float>::Iterator tdaIter =
    vtkAOSDataArrayTemplate<float>::FastDownCast(tda)->Begin();
  if (typeid(tdaBegin) != typeid(float*))
  {
    std::cerr << "Error: vtkAOSDataArrayTemplate<float>::Iterator is not a "
                 "float*.";
    return EXIT_FAILURE;
  }

  // Validate that the iterators return the same values from operator[] and
  // operator* as GetValue;
  for (vtkIdType i = 0; i < numValues; ++i)
  {
    float lookup = array->GetValue(i);
    if (lookup != datBegin[i] || lookup != tdaBegin[i] || lookup != *datIter || lookup != *tdaIter)
    {
      std::cerr << "Mismatch at " << i << ":"
                << " GetValue(i)=" << lookup << " datBegin[i]=" << datBegin[i]
                << " tdaBegin[i]=" << tdaBegin[i] << " *datIter=" << *datIter
                << " *tdaIter=" << *tdaIter << std::endl;
      return EXIT_FAILURE;
    }
    ++datIter;