ScatterCounting.h

Go to the documentation of this file.

//============================================================================
//  Copyright (c) Kitware, Inc.
//  All rights reserved.
//  See LICENSE.txt for details.
//
//  This software is distributed WITHOUT ANY WARRANTY; without even
//  the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
//  PURPOSE.  See the above copyright notice for more information.
//============================================================================
#ifndef vtk_m_worklet_ScatterCounting_h
#define vtk_m_worklet_ScatterCounting_h
 
#include <vtkm/worklet/internal/ScatterBase.h>
#include <vtkm/worklet/vtkm_worklet_export.h>
 
#include <vtkm/cont/UnknownArrayHandle.h>
 
#include <sstream>
 
namespace vtkm
{
namespace worklet
{
 
namespace detail
{
 
struct ScatterCountingBuilder;
 
} // namespace detail
 
struct VTKM_WORKLET_EXPORT ScatterCounting : internal::ScatterBase
{
  using CountTypes = vtkm::List<vtkm::Int64,
                                vtkm::Int32,
                                vtkm::Int16,
                                vtkm::Int8,
                                vtkm::UInt64,
                                vtkm::UInt32,
                                vtkm::UInt16,
                                vtkm::UInt8>;
 
  VTKM_CONT ScatterCounting(const vtkm::cont::UnknownArrayHandle& countArray,
                            vtkm::cont::DeviceAdapterId device = vtkm::cont::DeviceAdapterTagAny(),
                            bool saveInputToOutputMap = false)
  {
    this->BuildArrays(countArray, device, saveInputToOutputMap);
  }
  VTKM_CONT ScatterCounting(const vtkm::cont::UnknownArrayHandle& countArray,
                            bool saveInputToOutputMap)
  {
    this->BuildArrays(countArray, vtkm::cont::DeviceAdapterTagAny(), saveInputToOutputMap);
  }
 
  using OutputToInputMapType = vtkm::cont::ArrayHandle<vtkm::Id>;
 
  template <typename RangeType>
  VTKM_CONT OutputToInputMapType GetOutputToInputMap(RangeType) const
  {
    return this->OutputToInputMap;
  }
 
  using VisitArrayType = vtkm::cont::ArrayHandle<vtkm::IdComponent>;
  template <typename RangeType>
  VTKM_CONT VisitArrayType GetVisitArray(RangeType) const
  {
    return this->VisitArray;
  }
 
  VTKM_CONT
  vtkm::Id GetOutputRange(vtkm::Id inputRange) const
  {
    if (inputRange != this->InputRange)
    {
      std::stringstream msg;
      msg << "ScatterCounting initialized with input domain of size " << this->InputRange
          << " but used with a worklet invoke of size " << inputRange << std::endl;
      throw vtkm::cont::ErrorBadValue(msg.str());
    }
    return this->VisitArray.GetNumberOfValues();
  }
  VTKM_CONT
  vtkm::Id GetOutputRange(vtkm::Id3 inputRange) const
  {
    return this->GetOutputRange(inputRange[0] * inputRange[1] * inputRange[2]);
  }
 
  VTKM_CONT
  OutputToInputMapType GetOutputToInputMap() const { return this->OutputToInputMap; }
 
  VTKM_CONT
  vtkm::cont::ArrayHandle<vtkm::Id> GetInputToOutputMap() const { return this->InputToOutputMap; }
 
private:
  vtkm::Id InputRange;
  vtkm::cont::ArrayHandle<vtkm::Id> InputToOutputMap;
  OutputToInputMapType OutputToInputMap;
  VisitArrayType VisitArray;
 
  friend struct detail::ScatterCountingBuilder;
 
  VTKM_CONT void BuildArrays(const vtkm::cont::UnknownArrayHandle& countArray,
                             vtkm::cont::DeviceAdapterId device,
                             bool saveInputToOutputMap);
};
}
} // namespace vtkm::worklet
 
#endif //vtk_m_worklet_ScatterCounting_h