Method RequestData

IO/Chemistry/vtkGaussianCubeReader.cxx:54–237

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

//------------------------------------------------------------------------------
int vtkGaussianCubeReader::RequestData(vtkInformation* vtkNotUsed(request),
  vtkInformationVector** vtkNotUsed(inputVector), vtkInformationVector* outputVector)
{
  vtkInformation* outInfo = outputVector->GetInformationObject(0);
  vtkPolyData* output = vtkPolyData::SafeDownCast(outInfo->Get(vtkDataObject::DATA_OBJECT()));

  FILE* fp;
  char title[256];
  char data_name[256];
  double elements[16];
  int JN1, N1N2, n1, n2, n3, i, j, k;
  bool orbitalCubeFile = false;
  int numberOfOrbitals;

  // Output 0 (the default is the polydata)
  // Output 1 will be the gridded Image data

  vtkImageData* grid = this->GetGridOutput();

  if (!this->FileName)
  {
    return 0;
  }

  if ((fp = vtksys::SystemTools::Fopen(this->FileName, "r")) == nullptr)
  {
    vtkErrorMacro(<< "File " << this->FileName << " not found");
    return 0;
  }

  if (!fgets(title, 256, fp))
  {
    vtkErrorMacro("GaussianCubeReader error reading file: "
      << this->FileName << " Premature EOF while reading title.");
    fclose(fp);
    return 0;
  }

  // TODO: SystemTools::Split should be replaced by a SystemTools::SplitN call
  // which only splits up to N times as soon as it exists
  std::vector<std::string> tokens;
  vtksys::SystemTools::Split(title, tokens, ':');
  if (tokens.size() > 2)
  {
    for (std::size_t token = 3; token < tokens.size(); ++token)
    {
      tokens[2] += ":" + tokens[token];
    }
    strcpy(data_name, tokens[2].c_str());
    vtk::print(stderr, "label = {:s}\n", data_name);
  }

  if (!fgets(title, 256, fp))
  {
    vtkErrorMacro("GaussianCubeReader error reading file: "
      << this->FileName << " Premature EOF while reading title.");
    fclose(fp);
    return 0;