ThreadIndicesTopologyMap.h

Go to the documentation of this file.

//============================================================================
//  The contents of this file are covered by the Viskores license. See
//  LICENSE.txt for details.
//
//  By contributing to this file, all contributors agree to the Developer
//  Certificate of Origin Version 1.1 (DCO 1.1) as stated in DCO.txt.
//============================================================================
 
//============================================================================
//  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 viskores_exec_arg_ThreadIndicesTopologyMap_h
#define viskores_exec_arg_ThreadIndicesTopologyMap_h
 
#include <viskores/exec/arg/ThreadIndicesBasic.h>
 
#include <viskores/exec/ConnectivityPermuted.h>
#include <viskores/exec/ConnectivityStructured.h>
 
namespace viskores
{
namespace exec
{
namespace arg
{
 
namespace detail
{
 
static inline VISKORES_EXEC viskores::Id3 InflateTo3D(viskores::Id3 index)
{
  return index;
}
 
static inline VISKORES_EXEC viskores::Id3 InflateTo3D(viskores::Id2 index)
{
  return viskores::Id3(index[0], index[1], 0);
}
 
static inline VISKORES_EXEC viskores::Id3 InflateTo3D(viskores::Vec<viskores::Id, 1> index)
{
  return viskores::Id3(index[0], 0, 0);
}
 
static inline VISKORES_EXEC viskores::Id3 InflateTo3D(viskores::Id index)
{
  return viskores::Id3(index, 0, 0);
}
 
static inline VISKORES_EXEC viskores::Id3 Deflate(const viskores::Id3& index, viskores::Id3)
{
  return index;
}
 
static inline VISKORES_EXEC viskores::Id2 Deflate(const viskores::Id3& index, viskores::Id2)
{
  return viskores::Id2(index[0], index[1]);
}
 
} // namespace detail
 
struct DefaultScatterAndMaskTag
{
};
 
struct CustomScatterOrMaskTag
{
};
 
template <typename ConnectivityType, typename ScatterAndMaskMode>
class ThreadIndicesTopologyMap : public viskores::exec::arg::ThreadIndicesBasic
{
  using Superclass = viskores::exec::arg::ThreadIndicesBasic;
 
public:
  using IndicesIncidentType = typename ConnectivityType::IndicesType;
  using CellShapeTag = typename ConnectivityType::CellShapeTag;
  using Connectivity = ConnectivityType;
 
  VISKORES_SUPPRESS_EXEC_WARNINGS
  VISKORES_EXEC ThreadIndicesTopologyMap(viskores::Id threadIndex,
                                         viskores::Id inputIndex,
                                         viskores::IdComponent visitIndex,
                                         viskores::Id outputIndex,
                                         const ConnectivityType& connectivity)
    : Superclass(threadIndex, inputIndex, visitIndex, outputIndex)
    // The connectivity is stored in the invocation parameter at the given
    // input domain index. If this class is being used correctly, the type
    // of the domain will match the connectivity type used here. If there is
    // a compile error here about a type mismatch, chances are a worklet has
    // set its input domain incorrectly.
    , IndicesIncident(connectivity.GetIndices(inputIndex))
    , CellShape(connectivity.GetCellShape(inputIndex))
  {
  }
 
  VISKORES_EXEC
  const IndicesIncidentType& GetIndicesIncident() const { return this->IndicesIncident; }
 
  VISKORES_EXEC
  const IndicesIncidentType* GetIndicesIncidentPointer() const { return &this->IndicesIncident; }
 
  VISKORES_EXEC
  CellShapeTag GetCellShape() const { return this->CellShape; }
 
private:
  IndicesIncidentType IndicesIncident;
  CellShapeTag CellShape;
};
 
template <typename VisitTopology, typename IncidentTopology, viskores::IdComponent Dimension>
class ThreadIndicesTopologyMap<
  viskores::exec::ConnectivityStructured<VisitTopology, IncidentTopology, Dimension>,
  CustomScatterOrMaskTag>
{
  using ConnectivityType =
    viskores::exec::ConnectivityStructured<VisitTopology, IncidentTopology, Dimension>;
 
public:
  using IndicesIncidentType = typename ConnectivityType::IndicesType;
  using CellShapeTag = typename ConnectivityType::CellShapeTag;
  using LogicalIndexType = typename ConnectivityType::SchedulingRangeType;
  using Connectivity =
    viskores::exec::ConnectivityStructured<VisitTopology, IncidentTopology, Dimension>;
 
  VISKORES_EXEC ThreadIndicesTopologyMap(viskores::Id threadIndex,
                                         viskores::Id inputIndex,
                                         viskores::IdComponent visitIndex,
                                         viskores::Id outputIndex,
                                         const ConnectivityType& connectivity)
  {
    this->ThreadIndex = threadIndex;
    this->InputIndex = inputIndex;
    this->VisitIndex = visitIndex;
    this->OutputIndex = outputIndex;
    this->LogicalIndex = connectivity.FlatToLogicalVisitIndex(this->InputIndex);
    this->IndicesIncident = connectivity.GetIndices(this->LogicalIndex);
    this->CellShape = connectivity.GetCellShape(this->InputIndex);
  }
 
  VISKORES_EXEC ThreadIndicesTopologyMap(const viskores::Id3& threadIndex3D,
                                         viskores::Id threadIndex1D,
                                         const ConnectivityType& connectivity)
  {
    // This constructor handles multidimensional indices on one-to-one input-to-output
    auto logicalIndex = detail::Deflate(threadIndex3D, LogicalIndexType());
    this->ThreadIndex = threadIndex1D;
    this->LogicalIndex = logicalIndex;
    this->IndicesIncident = connectivity.GetIndices(logicalIndex);
    this->CellShape = connectivity.GetCellShape(threadIndex1D);
  }
 
  VISKORES_EXEC ThreadIndicesTopologyMap(const viskores::Id3& threadIndex3D,
                                         viskores::Id threadIndex1D,
                                         viskores::Id inIndex,
                                         viskores::IdComponent visitIndex,
                                         viskores::Id outIndex,
                                         const ConnectivityType& connectivity)
  {
    // This constructor handles multidimensional indices on many-to-many input-to-output
    auto logicalIndex = detail::Deflate(threadIndex3D, LogicalIndexType());
    this->ThreadIndex = threadIndex1D;
    this->InputIndex = inIndex;
    this->VisitIndex = visitIndex;
    this->OutputIndex = outIndex;
    this->LogicalIndex = logicalIndex;
    this->IndicesIncident = connectivity.GetIndices(logicalIndex);
    this->CellShape = connectivity.GetCellShape(threadIndex1D);
  }
 
 
  VISKORES_EXEC
  viskores::Id GetThreadIndex() const { return this->ThreadIndex; }
 
  VISKORES_EXEC
  LogicalIndexType GetIndexLogical() const { return this->LogicalIndex; }
 
  VISKORES_EXEC
  viskores::Id GetInputIndex() const { return this->InputIndex; }
 
  VISKORES_EXEC
  viskores::Id3 GetInputIndex3D() const { return detail::InflateTo3D(this->GetIndexLogical()); }
 
  VISKORES_EXEC
  viskores::Id GetOutputIndex() const { return this->OutputIndex; }
 
  VISKORES_EXEC
  viskores::IdComponent GetVisitIndex() const { return this->VisitIndex; }
 
  VISKORES_EXEC
  const IndicesIncidentType& GetIndicesIncident() const { return this->IndicesIncident; }
 
  VISKORES_EXEC
  const IndicesIncidentType* GetIndicesIncidentPointer() const { return &this->IndicesIncident; }
 
  VISKORES_EXEC
  CellShapeTag GetCellShape() const { return this->CellShape; }
 
private:
  viskores::Id ThreadIndex;
  viskores::IdComponent VisitIndex;
  LogicalIndexType LogicalIndex;
  IndicesIncidentType IndicesIncident;
  CellShapeTag CellShape;
  viskores::Id InputIndex;
  viskores::Id OutputIndex;
};
 
template <typename VisitTopology, typename IncidentTopology, viskores::IdComponent Dimension>
class ThreadIndicesTopologyMap<
  viskores::exec::ConnectivityStructured<VisitTopology, IncidentTopology, Dimension>,
  DefaultScatterAndMaskTag>
{
  using ConnectivityType =
    viskores::exec::ConnectivityStructured<VisitTopology, IncidentTopology, Dimension>;
 
public:
  using IndicesIncidentType = typename ConnectivityType::IndicesType;
  using CellShapeTag = typename ConnectivityType::CellShapeTag;
  using LogicalIndexType = typename ConnectivityType::SchedulingRangeType;
  using Connectivity = ConnectivityType;
 
  VISKORES_EXEC ThreadIndicesTopologyMap(viskores::Id threadIndex,
                                         viskores::Id inputIndex,
                                         viskores::IdComponent viskoresNotUsed(visitIndex),
                                         viskores::Id viskoresNotUsed(outputIndex),
                                         const ConnectivityType& connectivity)
  {
    this->ThreadIndex = threadIndex;
    this->LogicalIndex = connectivity.FlatToLogicalVisitIndex(inputIndex);
    this->IndicesIncident = connectivity.GetIndices(this->LogicalIndex);
    this->CellShape = connectivity.GetCellShape(inputIndex);
  }
 
 
 
 
  VISKORES_EXEC ThreadIndicesTopologyMap(const viskores::Id3& threadIndex3D,
                                         viskores::Id threadIndex1D,
                                         const ConnectivityType& connectivity)
  {
    // This constructor handles multidimensional indices on one-to-one input-to-output
    auto logicalIndex = detail::Deflate(threadIndex3D, LogicalIndexType());
    this->ThreadIndex = threadIndex1D;
    this->LogicalIndex = logicalIndex;
    this->IndicesIncident = connectivity.GetIndices(logicalIndex);
    this->CellShape = connectivity.GetCellShape(threadIndex1D);
  }
 
  VISKORES_EXEC ThreadIndicesTopologyMap(const viskores::Id3& threadIndex3D,
                                         viskores::Id threadIndex1D,
                                         viskores::Id viskoresNotUsed(inIndex),
                                         viskores::IdComponent viskoresNotUsed(visitIndex),
                                         viskores::Id viskoresNotUsed(outIndex),
                                         const ConnectivityType& connectivity)
  {
    // This constructor handles multidimensional indices on many-to-many input-to-output
    auto logicalIndex = detail::Deflate(threadIndex3D, LogicalIndexType());
    this->ThreadIndex = threadIndex1D;
    this->LogicalIndex = logicalIndex;
    this->IndicesIncident = connectivity.GetIndices(logicalIndex);
    this->CellShape = connectivity.GetCellShape(threadIndex1D);
  }
 
  VISKORES_EXEC
  viskores::Id GetThreadIndex() const { return this->ThreadIndex; }
 
  VISKORES_EXEC
  LogicalIndexType GetIndexLogical() const { return this->LogicalIndex; }
 
  VISKORES_EXEC
  viskores::Id GetInputIndex() const { return this->ThreadIndex; }
 
  VISKORES_EXEC
  viskores::Id3 GetInputIndex3D() const { return detail::InflateTo3D(this->GetIndexLogical()); }
 
  VISKORES_EXEC
  viskores::Id GetOutputIndex() const { return this->ThreadIndex; }
 
  VISKORES_EXEC
  viskores::IdComponent GetVisitIndex() const { return 0; }
 
  VISKORES_EXEC
  const IndicesIncidentType& GetIndicesIncident() const { return this->IndicesIncident; }
 
  VISKORES_EXEC
  const IndicesIncidentType* GetIndicesIncidentPointer() const { return &this->IndicesIncident; }
 
  VISKORES_EXEC
  CellShapeTag GetCellShape() const { return this->CellShape; }
 
private:
  viskores::Id ThreadIndex;
  LogicalIndexType LogicalIndex;
  IndicesIncidentType IndicesIncident;
  CellShapeTag CellShape;
};
 
 
template <typename PermutationPortal, viskores::IdComponent Dimension>
class ThreadIndicesTopologyMap<
  viskores::exec::ConnectivityPermutedVisitCellsWithPoints<
    PermutationPortal,
    viskores::exec::ConnectivityStructured<viskores::TopologyElementTagCell,
                                           viskores::TopologyElementTagPoint,
                                           Dimension>>,
  CustomScatterOrMaskTag>
{
  using PermutedConnectivityType = viskores::exec::ConnectivityPermutedVisitCellsWithPoints<
    PermutationPortal,
    viskores::exec::ConnectivityStructured<viskores::TopologyElementTagCell,
                                           viskores::TopologyElementTagPoint,
                                           Dimension>>;
  using ConnectivityType = viskores::exec::ConnectivityStructured<viskores::TopologyElementTagCell,
                                                                  viskores::TopologyElementTagPoint,
                                                                  Dimension>;
 
public:
  using IndicesIncidentType = typename ConnectivityType::IndicesType;
  using CellShapeTag = typename ConnectivityType::CellShapeTag;
  using LogicalIndexType = typename ConnectivityType::SchedulingRangeType;
  using Connectivity = viskores::exec::ConnectivityPermutedVisitCellsWithPoints<
    PermutationPortal,
    viskores::exec::ConnectivityStructured<viskores::TopologyElementTagCell,
                                           viskores::TopologyElementTagPoint,
                                           Dimension>>;
 
  VISKORES_EXEC ThreadIndicesTopologyMap(viskores::Id threadIndex,
                                         viskores::Id inputIndex,
                                         viskores::IdComponent visitIndex,
                                         viskores::Id outputIndex,
                                         const PermutedConnectivityType& permutation)
  {
    this->ThreadIndex = threadIndex;
    this->InputIndex = inputIndex;
    this->VisitIndex = visitIndex;
    this->OutputIndex = outputIndex;
 
    const viskores::Id permutedIndex = permutation.Portal.Get(this->InputIndex);
    this->LogicalIndex = permutation.Connectivity.FlatToLogicalVisitIndex(permutedIndex);
    this->IndicesIncident = permutation.Connectivity.GetIndices(this->LogicalIndex);
    this->CellShape = permutation.Connectivity.GetCellShape(permutedIndex);
  }
 
  VISKORES_EXEC
  viskores::Id GetThreadIndex() const { return this->ThreadIndex; }
 
  VISKORES_EXEC
  LogicalIndexType GetIndexLogical() const { return this->LogicalIndex; }
 
  VISKORES_EXEC
  viskores::Id GetInputIndex() const { return this->InputIndex; }
 
  VISKORES_EXEC
  viskores::Id3 GetInputIndex3D() const { return detail::InflateTo3D(this->GetIndexLogical()); }
 
  VISKORES_EXEC
  viskores::Id GetOutputIndex() const { return this->OutputIndex; }
 
  VISKORES_EXEC
  viskores::IdComponent GetVisitIndex() const { return this->VisitIndex; }
 
  VISKORES_EXEC
  const IndicesIncidentType& GetIndicesIncident() const { return this->IndicesIncident; }
 
  VISKORES_EXEC
  const IndicesIncidentType* GetIndicesIncidentPointer() const { return &this->IndicesIncident; }
 
  VISKORES_EXEC
  CellShapeTag GetCellShape() const { return this->CellShape; }
 
private:
  viskores::Id ThreadIndex;
  viskores::Id InputIndex;
  viskores::IdComponent VisitIndex;
  viskores::Id OutputIndex;
  LogicalIndexType LogicalIndex;
  IndicesIncidentType IndicesIncident;
  CellShapeTag CellShape;
};
}
}
} // namespace viskores::exec::arg
 
#endif //viskores_exec_arg_ThreadIndicesTopologyMap_h