Method copyContinuumData

IO/Parallel/vtkNek5000Reader.cxx:1643–1739

}

void vtkNek5000Reader::copyContinuumData(vtkUnstructuredGrid* pv_ugrid)
{
#ifndef NDEBUG
  int my_rank;
  vtkMultiProcessController* ctrl = vtkMultiProcessController::GetGlobalController();
  if (ctrl != nullptr)
  {
    my_rank = ctrl->GetLocalProcessId();
  }
  else
  {
    my_rank = 0;
  }
#endif
  int index = 0;
  int num_verts = this->myNumBlocks * this->totalBlockSize;

  // for each variable
  for (auto v_index = 0; v_index < this->num_vars; v_index++)
  {
    if (this->GetPointArrayStatus(v_index))
    {
      // if this is a scalar
      if (this->var_length[v_index] == 1)
      {
        index = 0;
        vtkFloatArray* scalars = vtkFloatArray::New();
        // scalars->SetNumberOfComponents(1);
        // scalars->SetNumberOfTuples(num_verts);
        scalars->SetName(this->var_names[v_index]);
        scalars->SetArray(
          this->dataArray[v_index], num_verts, 0, vtkDataArray::VTK_DATA_ARRAY_DELETE);
        /*
        for (int b_index = 0; b_index < this->myNumBlocks; ++b_index)
        {
          for (int p_index = 0; p_index < this->totalBlockSize; ++p_index)
          {
            scalars->SetValue(index, this->dataArray[v_index][index]);
            index++;
          }
        }*/
        this->UGrid->GetPointData()->AddArray(scalars);
        scalars->Delete();
      }
      // if this is a vector
      else if (this->var_length[v_index] > 1)
      {
        index = 0;
        vtkFloatArray* vectors = vtkFloatArray::New();
        vectors->SetNumberOfComponents(3);
        vectors->SetNumberOfTuples(num_verts);
        vectors->SetName(this->var_names[v_index]);

        // for each  element/block in the continuum mesh
        for (int b_index = 0; b_index < this->myNumBlocks; ++b_index)
        {
          // for every point in this element/block
#ifndef NDEBUG