Function computeL1Norm

Common/Core/vtkArrayComponents.cxx:120–141

};

vtkSmartPointer<vtkDataArray> computeL1Norm(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 sum = 0;
        for (int jj = 0; jj < numComponents; ++jj)
        {
          sum += absolute(tuple[jj]);
        }
        norm->SetTuple1(ii, sum);
      }
    });
  return norm;
}

vtkSmartPointer<vtkDataArray> computeL2Norm(vtkDataArray* array)
{