WorkletCellNeighborhood.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_WorkletCellNeighborhood_h
#define vtk_m_worklet_WorkletCellNeighborhood_h
 
 
#include <vtkm/exec/arg/ThreadIndicesCellNeighborhood.h>
#include <vtkm/worklet/DispatcherCellNeighborhood.h>
#include <vtkm/worklet/WorkletNeighborhood.h>
 
namespace vtkm
{
namespace worklet
{
 
template <typename WorkletType>
class DispatcherCellNeighborhood;
 
class WorkletCellNeighborhood : public WorkletNeighborhood
{
public:
  template <typename Worklet>
  using Dispatcher = vtkm::worklet::DispatcherCellNeighborhood<Worklet>;
 
#ifdef VTKM_DOXYGEN_ONLY
  // These redeclarations of superclass features are for documentation purposes only.
 
  struct CellSetIn : vtkm::worklet::WorkletNeighborhood::CellSetIn
  {
  };
 
  struct FieldIn : vtkm::worklet::WorkletNeighborhood::FieldIn
  {
  };
 
  struct FieldInNeighborhood : vtkm::worklet::WorkletNeighborhood::FieldInNeighborhood
  {
  };
 
  struct FieldOut : vtkm::worklet::WorkletNeighborhood::FieldOut
  {
  };
 
  struct FieldInOut : vtkm::worklet::WorkletNeighborhood::FieldInOut
  {
  };
 
  struct WholeArrayIn : vtkm::worklet::internal::WorkletBase::WholeArrayIn
  {
  };
 
  struct WholeArrayOut : vtkm::worklet::internal::WorkletBase::WholeArrayOut
  {
  };
 
  struct WholeArrayInOut : vtkm::worklet::internal::WorkletBase::WholeArrayInOut
  {
  };
 
  struct AtomicArrayInOut : vtkm::worklet::internal::WorkletBase::AtomicArrayInOut
  {
  };
 
  template <typename VisitTopology = Cell, typename IncidentTopology = Point>
  struct WholeCellSetIn
    : vtkm::worklet::internal::WorkletBase::WholeCellSetIn<VisitTopology, IncidentTopology>
  {
  };
 
  struct ExecObject : vtkm::worklet::internal::WorkletBase::ExecObject
  {
  };
#endif // VTKM_DOXYGEN_ONLY
 
#ifdef VTKM_DOXYGEN_ONLY
  // These redeclarations of superclass features are for documentation purposes only.
 
  struct _1 : vtkm::worklet::internal::WorkletBase::_1
  {
  };
 
  struct Boundary : vtkm::worklet::WorkletNeighborhood::Boundary
  {
  };
 
  struct WorkIndex : vtkm::worklet::internal::WorkletBase::WorkIndex
  {
  };
 
  struct VisitIndex : vtkm::worklet::internal::WorkletBase::VisitIndex
  {
  };
 
  struct InputIndex : vtkm::worklet::internal::WorkletBase::InputIndex
  {
  };
 
  struct OutputIndex : vtkm::worklet::internal::WorkletBase::OutputIndex
  {
  };
 
  struct ThreadIndices : vtkm::worklet::internal::WorkletBase::ThreadIndices
  {
  };
 
  struct Device : vtkm::worklet::internal::WorkletBase::Device
  {
  };
#endif // VTKM_DOXYGEN_ONLY
 
  VTKM_SUPPRESS_EXEC_WARNINGS
  template <typename OutToInArrayType,
            typename VisitArrayType,
            typename ThreadToOutArrayType,
            vtkm::IdComponent Dimension>
  VTKM_EXEC vtkm::exec::arg::ThreadIndicesCellNeighborhood GetThreadIndices(
    vtkm::Id threadIndex,
    const OutToInArrayType& outToIn,
    const VisitArrayType& visit,
    const ThreadToOutArrayType& threadToOut,
    const vtkm::exec::ConnectivityStructured<vtkm::TopologyElementTagPoint,
                                             vtkm::TopologyElementTagCell,
                                             Dimension>& inputDomain //this should be explicit
  ) const
  {
    const vtkm::Id outIndex = threadToOut.Get(threadIndex);
    return vtkm::exec::arg::ThreadIndicesCellNeighborhood(
      threadIndex, outToIn.Get(outIndex), visit.Get(outIndex), outIndex, inputDomain);
  }
 
 
private:
  static constexpr bool IsScatterIdentity =
    std::is_same<ScatterType, vtkm::worklet::ScatterIdentity>::value;
  static constexpr bool IsMaskNone = std::is_same<MaskType, vtkm::worklet::MaskNone>::value;
 
public:
  template <bool Cond, typename ReturnType>
  using EnableFnWhen = typename std::enable_if<Cond, ReturnType>::type;
 
  VTKM_SUPPRESS_EXEC_WARNINGS
  template <typename OutToInArrayType,
            typename VisitArrayType,
            typename ThreadToOutArrayType,
            typename InputDomainType,
            bool S = IsScatterIdentity,
            bool M = IsMaskNone>
  VTKM_EXEC EnableFnWhen<S && M, vtkm::exec::arg::ThreadIndicesCellNeighborhood> GetThreadIndices(
    vtkm::Id threadIndex1D,
    const vtkm::Id3& threadIndex3D,
    const OutToInArrayType& vtkmNotUsed(outToIn),
    const VisitArrayType& vtkmNotUsed(visit),
    const ThreadToOutArrayType& vtkmNotUsed(threadToOut),
    const InputDomainType& connectivity) const
  {
    return vtkm::exec::arg::ThreadIndicesCellNeighborhood(
      threadIndex3D, threadIndex1D, connectivity);
  }
 
  VTKM_SUPPRESS_EXEC_WARNINGS
  template <typename OutToInArrayType,
            typename VisitArrayType,
            typename ThreadToOutArrayType,
            typename InputDomainType,
            bool S = IsScatterIdentity,
            bool M = IsMaskNone>
  VTKM_EXEC EnableFnWhen<!(S && M), vtkm::exec::arg::ThreadIndicesCellNeighborhood>
  GetThreadIndices(vtkm::Id threadIndex1D,
                   const vtkm::Id3& threadIndex3D,
                   const OutToInArrayType& outToIn,
                   const VisitArrayType& visit,
                   const ThreadToOutArrayType& threadToOut,
                   const InputDomainType& connectivity) const
  {
    const vtkm::Id outIndex = threadToOut.Get(threadIndex1D);
    return vtkm::exec::arg::ThreadIndicesCellNeighborhood(threadIndex3D,
                                                          threadIndex1D,
                                                          outToIn.Get(outIndex),
                                                          visit.Get(outIndex),
                                                          outIndex,
                                                          connectivity);
  }
};
}
}
 
#endif