Method AddVectorArray

Filters/Hybrid/vtkTemporalFractal.cxx:1021–1081

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

//------------------------------------------------------------------------------
void vtkTemporalFractal::AddVectorArray(vtkOverlappingAMR* output)
{
  double* origin = this->GetTopLevelOrigin();

  int levels = output->GetNumberOfLevels();
  int level = 0;
  while (level < levels)
  {
    int blocks = output->GetNumberOfBlocks(level);
    int block = 0;
    while (block < blocks)
    {
      vtkUniformGrid* grid;
      grid = vtkUniformGrid::SafeDownCast(output->GetDataSetAsCartesianGrid(level, block));
      assert("check: grid_exists" && grid != nullptr);

      vtkDoubleArray* array = vtkDoubleArray::New();
      array->SetNumberOfComponents(3);
      int numCells = grid->GetNumberOfCells();
      array->Allocate(numCells);
      array->SetNumberOfTuples(numCells);
      double* arrayPtr = static_cast<double*>(array->GetPointer(0));
      double spacing[3];
      grid->GetSpacing(spacing);
      int x, y, z;
      int ext[6];
      grid->GetExtent(ext);
      // we need cell extents but we just get point extents
      if (ext[5] > 0)
      {
        --ext[5];
      }
      if (ext[3] > 0)
      {
        --ext[3];
      }
      if (ext[1] > 0)
      {
        --ext[1];
      }

      for (z = ext[4]; z <= ext[5]; ++z)
      {
        for (y = ext[2]; y <= ext[3]; ++y)
        {
          for (x = ext[0]; x <= ext[1]; ++x)
          {
            *arrayPtr++ = origin[0] + spacing[0] * (static_cast<double>(x) + 0.5);
            *arrayPtr++ = origin[1] + spacing[1] * (static_cast<double>(y) + 0.5);
            *arrayPtr++ = origin[2] + spacing[2] * (static_cast<double>(z) + 0.5);
          }
        }
      }
      array->SetName("VectorXYZ");
      grid->GetCellData()->AddArray(array);
      array->Delete();
      ++block;
    }