Function computeL2Norm

Common/Core/vtkArrayComponents.cxx:143–164

}

vtkSmartPointer<vtkDataArray> computeL2Norm(vtkDataArray* array)
{
  auto norm = vtkSmartPointer<vtkDoubleArray>::New();
  norm->SetNumberOfTuples(array->GetNumberOfTuples());
  vtkSMPTools::For(0, array->GetNumberOfTuples(),
    [&norm, array](vtkIdType begin, vtkIdType end)
    {
      int numComponents = array->GetNumberOfComponents();
      std::vector<double> tuple(numComponents);
      for (vtkIdType ii = begin; ii < end; ++ii)
      {
        array->GetTuple(ii, tuple.data());
        double norm2 = 0;
        for (int jj = 0; jj < numComponents; ++jj)
        {
          norm2 += tuple[jj] * tuple[jj];
        }
        norm->SetValue(ii, std::sqrt(norm2));
      }
    });
  return norm;
}

vtkSmartPointer<vtkDataArray> computeLInfNorm(vtkDataArray* array)
{