Function computeLInfNorm

Common/Core/vtkArrayComponents.cxx:166–192

}

vtkSmartPointer<vtkDataArray> computeLInfNorm(vtkDataArray* array)
{
  vtkSmartPointer<vtkDataArray> norm;
  norm.TakeReference(array->NewInstance());
  norm->SetNumberOfTuples(array->GetNumberOfTuples());
  if (!vtkArrayDispatch::Dispatch2::Execute(array, norm, InfNorm{}))
  {
    vtkSMPTools::For(0, array->GetNumberOfTuples(),
      [&norm, array](vtkIdType begin, vtkIdType end)
      {
        int numComponents = array->GetNumberOfComponents();
        assert(numComponents > 1);
        std::vector<double> tuple(numComponents);
        for (vtkIdType ii = begin; ii < end; ++ii)
        {
          array->GetTuple(ii, tuple.data());
          double max = std::abs(tuple[0]);
          for (int jj = 1; jj < numComponents; ++jj)
          {
            max = std::max(max, std::abs(tuple[jj]));
          }
          norm->SetTuple1(ii, max);
        }
      });
  }
  return norm;
}

} // anonymous namespace