TransportTagTopologyFieldIn.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_cont_arg_TransportTagTopologyFieldIn_h
#define vtk_m_cont_arg_TransportTagTopologyFieldIn_h
 
#include <vtkm/TopologyElementTag.h>
#include <vtkm/Types.h>
 
#include <vtkm/cont/ArrayHandle.h>
#include <vtkm/cont/CellSet.h>
 
#include <vtkm/cont/arg/Transport.h>
 
namespace vtkm
{
namespace cont
{
namespace arg
{
 
template <typename TopologyElementTag>
struct TransportTagTopologyFieldIn
{
};
 
namespace detail
{
 
VTKM_CONT
inline static vtkm::Id TopologyDomainSize(const vtkm::cont::CellSet& cellSet,
                                          vtkm::TopologyElementTagPoint)
{
  return cellSet.GetNumberOfPoints();
}
 
VTKM_CONT
inline static vtkm::Id TopologyDomainSize(const vtkm::cont::CellSet& cellSet,
                                          vtkm::TopologyElementTagCell)
{
  return cellSet.GetNumberOfCells();
}
 
VTKM_CONT
inline static vtkm::Id TopologyDomainSize(const vtkm::cont::CellSet& cellSet,
                                          vtkm::TopologyElementTagFace)
{
  return cellSet.GetNumberOfFaces();
}
 
VTKM_CONT
inline static vtkm::Id TopologyDomainSize(const vtkm::cont::CellSet& cellSet,
                                          vtkm::TopologyElementTagEdge)
{
  return cellSet.GetNumberOfEdges();
}
 
} // namespace detail
 
template <typename TopologyElementTag, typename ContObjectType, typename Device>
struct Transport<vtkm::cont::arg::TransportTagTopologyFieldIn<TopologyElementTag>,
                 ContObjectType,
                 Device>
{
  VTKM_IS_ARRAY_HANDLE(ContObjectType);
 
 
  using ExecObjectType = decltype(
    std::declval<ContObjectType>().PrepareForInput(Device(), std::declval<vtkm::cont::Token&>()));
 
  VTKM_CONT
  ExecObjectType operator()(const ContObjectType& object,
                            const vtkm::cont::CellSet& inputDomain,
                            vtkm::Id,
                            vtkm::Id,
                            vtkm::cont::Token& token) const
  {
    if (object.GetNumberOfValues() != detail::TopologyDomainSize(inputDomain, TopologyElementTag()))
    {
      throw vtkm::cont::ErrorBadValue("Input array to worklet invocation the wrong size.");
    }
 
    return object.PrepareForInput(Device(), token);
  }
};
}
}
} // namespace vtkm::cont::arg
 
#endif //vtk_m_cont_arg_TransportTagTopologyFieldIn_h