Method RequestData

Filters/ParallelStatistics/vtkGenerateStatistics.cxx:604–670

}

int vtkGenerateStatistics::RequestData(
  vtkInformation* vtkNotUsed(request), vtkInformationVector** input, vtkInformationVector* output)
{
  vtkDataObject* dataObjIn = vtkDataObject::GetData(input[0], 0);
  if (!dataObjIn)
  {
    // Silently ignore missing data.
    return 1;
  }

  int numArrays = this->GetNumberOfInputArraySpecifications();
  if (numArrays <= 0)
  {
    // If we have variables specified by the EnableAttributeArray API,
    // use it to populate "SetInputArraysToProcess".
    if (!this->P->Buffer.empty())
    {
      this->P->AddBufferToRequests();
    }
    int aa = 0;
    for (const auto& request : this->P->Requests)
    {
      for (const auto& name : request)
      {
        this->SetInputArrayToProcess(
          aa, 0, 0, vtkDataObject::FIELD_ASSOCIATION_POINTS_THEN_CELLS, name.c_str());
        ++aa;
      }
    }
    // Silently ignore empty requests.
    numArrays = this->GetNumberOfInputArraySpecifications();
    if (numArrays <= 0)
    {
      return 1;
    }
  }

  // Get output model data and sci-viz data.
  auto* modelObjOu = vtkPartitionedDataSetCollection::GetData(output, 0);
  if (!modelObjOu)
  {
    // Silently ignore missing data.
    return 1;
  }

  // TODO: Perform a "pre-run" stage to compute total number of samples across
  //       all blocks in all partitions on all ranks? If so, we should then
  //       have a target sample size on a per-leaf basis. This stage is only
  //       needed if users are allowed to specify a fixed sample size rather
  //       than a training *fraction*.

  int stat = this->RequestLocalDataDispatch(dataObjIn, modelObjOu);

  // The RequestLocalDataDispatch() above requires no communication and does
  // all the local model aggregation possible. The remaining global model(s)
  // then need to be collectively aggregated (in pairs of ranks with ⌈log₂(N)⌉
  // merges for N ranks) resulting in the final model aggregated on rank 0.
  // The final model(s) should then be broadcast from rank 0 to all ranks so