5dd346007b
clang-format BinPack settings have been disabled to make sure that the VTK-m style guideline is obeyed.
252 lines
12 KiB
C++
252 lines
12 KiB
C++
//============================================================================
|
|
// 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.
|
|
//
|
|
// Copyright 2016 Sandia Corporation.
|
|
// Copyright 2016 UT-Battelle, LLC.
|
|
// Copyright 2016 Los Alamos National Security.
|
|
//
|
|
// Under the terms of Contract DE-AC04-94AL85000 with Sandia Corporation,
|
|
// the U.S. Government retains certain rights in this software.
|
|
//
|
|
// Under the terms of Contract DE-AC52-06NA25396 with Los Alamos National
|
|
// Laboratory (LANL), the U.S. Government retains certain rights in
|
|
// this software.
|
|
//============================================================================
|
|
#ifndef vtk_m_exec_CellFace_h
|
|
#define vtk_m_exec_CellFace_h
|
|
|
|
#include <vtkm/Assert.h>
|
|
#include <vtkm/CellShape.h>
|
|
#include <vtkm/Types.h>
|
|
#include <vtkm/exec/FunctorBase.h>
|
|
#include <vtkm/internal/Assume.h>
|
|
|
|
namespace vtkm
|
|
{
|
|
namespace exec
|
|
{
|
|
|
|
namespace detail
|
|
{
|
|
|
|
static const vtkm::IdComponent MAX_FACE_SIZE = 4;
|
|
static const vtkm::IdComponent MAX_NUM_FACES = 6;
|
|
|
|
VTKM_EXEC_CONSTANT
|
|
static const vtkm::IdComponent NumFaces[vtkm::NUMBER_OF_CELL_SHAPES] = {
|
|
0, // 0: CELL_SHAPE_EMPTY
|
|
0, // 1: CELL_SHAPE_VERTEX
|
|
0, // 2: Unused
|
|
0, // 3: CELL_SHAPE_LINE
|
|
0, // 4: Unused
|
|
0, // 5: CELL_SHAPE_TRIANGLE
|
|
0, // 6: Unused
|
|
0, // 7: CELL_SHAPE_POLYGON
|
|
0, // 8: Unused
|
|
0, // 9: CELL_SHAPE_QUAD
|
|
4, // 10: CELL_SHAPE_TETRA
|
|
0, // 11: Unused
|
|
6, // 12: CELL_SHAPE_HEXAHEDRON
|
|
5, // 13: CELL_SHAPE_WEDGE
|
|
5 // 14: CELL_SHAPE_PYRAMID
|
|
};
|
|
|
|
VTKM_EXEC_CONSTANT
|
|
static const vtkm::IdComponent NumPointsInFace[vtkm::NUMBER_OF_CELL_SHAPES][MAX_NUM_FACES] = {
|
|
{ -1, -1, -1, -1, -1, -1 }, // 0: CELL_SHAPE_EMPTY
|
|
{ -1, -1, -1, -1, -1, -1 }, // 1: CELL_SHAPE_VERTEX
|
|
{ -1, -1, -1, -1, -1, -1 }, // 2: Unused
|
|
{ -1, -1, -1, -1, -1, -1 }, // 3: CELL_SHAPE_LINE
|
|
{ -1, -1, -1, -1, -1, -1 }, // 4: Unused
|
|
{ -1, -1, -1, -1, -1, -1 }, // 5: CELL_SHAPE_TRIANGLE
|
|
{ -1, -1, -1, -1, -1, -1 }, // 6: Unused
|
|
{ -1, -1, -1, -1, -1, -1 }, // 7: CELL_SHAPE_POLYGON
|
|
{ -1, -1, -1, -1, -1, -1 }, // 8: Unused
|
|
{ -1, -1, -1, -1, -1, -1 }, // 9: CELL_SHAPE_QUAD
|
|
{ 3, 3, 3, 3, -1, -1 }, // 10: CELL_SHAPE_TETRA
|
|
{ -1, -1, -1, -1, -1, -1 }, // 11: Unused
|
|
{ 4, 4, 4, 4, 4, 4 }, // 12: CELL_SHAPE_HEXAHEDRON
|
|
{ 3, 3, 4, 4, 4, -1 }, // 13: CELL_SHAPE_WEDGE
|
|
{ 4, 3, 3, 3, 3, -1 } // 14: CELL_SHAPE_PYRAMID
|
|
};
|
|
|
|
VTKM_EXEC_CONSTANT
|
|
static const vtkm::IdComponent PointsInFace[vtkm::NUMBER_OF_CELL_SHAPES][MAX_NUM_FACES]
|
|
[MAX_FACE_SIZE] = {
|
|
// 0: CELL_SHAPE_EMPTY
|
|
{ { -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 } },
|
|
// 1: CELL_SHAPE_VERTEX
|
|
{ { -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 } },
|
|
// 2: Unused
|
|
{ { -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 } },
|
|
// 3: CELL_SHAPE_LINE
|
|
{ { -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 } },
|
|
// 4: Unused
|
|
{ { -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 } },
|
|
// 5: CELL_SHAPE_TRIANGLE
|
|
{ { -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 } },
|
|
// 6: Unused
|
|
{ { -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 } },
|
|
// 7: CELL_SHAPE_POLYGON
|
|
{ { -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 } },
|
|
// 8: Unused
|
|
{ { -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 } },
|
|
// 9: CELL_SHAPE_QUAD
|
|
{ { -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 } },
|
|
// 10: CELL_SHAPE_TETRA
|
|
{ { 0, 1, 3, -1 },
|
|
{ 1, 2, 3, -1 },
|
|
{ 2, 0, 3, -1 },
|
|
{ 0, 2, 1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 } },
|
|
// 11: Unused
|
|
{ { -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 },
|
|
{ -1, -1, -1, -1 } },
|
|
// 12: CELL_SHAPE_HEXAHEDRON
|
|
{ { 0, 4, 7, 3 },
|
|
{ 1, 2, 6, 5 },
|
|
{ 0, 1, 5, 4 },
|
|
{ 3, 7, 6, 2 },
|
|
{ 0, 3, 2, 1 },
|
|
{ 4, 5, 6, 7 } },
|
|
// 13: CELL_SHAPE_WEDGE
|
|
{ { 0, 1, 2, -1 },
|
|
{ 3, 5, 4, -1 },
|
|
{ 0, 3, 4, 1 },
|
|
{ 1, 4, 5, 2 },
|
|
{ 2, 5, 3, 0 },
|
|
{ -1, -1, -1, -1 } },
|
|
// 14: CELL_SHAPE_PYRAMID
|
|
{ { 0, 3, 2, 1 },
|
|
{ 0, 1, 4, -1 },
|
|
{ 1, 2, 4, -1 },
|
|
{ 2, 3, 4, -1 },
|
|
{ 3, 0, 4, -1 },
|
|
{ -1, -1, -1, -1 } }
|
|
};
|
|
|
|
} // namespace detail
|
|
|
|
template <typename CellShapeTag>
|
|
static inline VTKM_EXEC vtkm::IdComponent CellFaceNumberOfFaces(CellShapeTag shape,
|
|
const vtkm::exec::FunctorBase&)
|
|
{
|
|
return detail::NumFaces[shape.Id];
|
|
}
|
|
|
|
template <typename CellShapeTag>
|
|
static inline VTKM_EXEC vtkm::IdComponent CellFaceNumberOfPoints(
|
|
vtkm::IdComponent faceIndex,
|
|
CellShapeTag shape,
|
|
const vtkm::exec::FunctorBase& worklet)
|
|
{
|
|
VTKM_ASSUME(faceIndex >= 0);
|
|
VTKM_ASSUME(faceIndex < detail::MAX_NUM_FACES);
|
|
if (faceIndex >= vtkm::exec::CellFaceNumberOfFaces(shape, worklet))
|
|
{
|
|
worklet.RaiseError("Invalid face number.");
|
|
return 0;
|
|
}
|
|
|
|
return detail::NumPointsInFace[shape.Id][faceIndex];
|
|
}
|
|
|
|
template <typename CellShapeTag>
|
|
static inline VTKM_EXEC vtkm::UInt8 CellFaceShape(vtkm::IdComponent faceIndex,
|
|
CellShapeTag shape,
|
|
const vtkm::exec::FunctorBase& worklet)
|
|
{
|
|
VTKM_ASSUME(faceIndex >= 0);
|
|
VTKM_ASSUME(faceIndex < detail::MAX_NUM_FACES);
|
|
switch (CellFaceNumberOfPoints(faceIndex, shape, worklet))
|
|
{
|
|
case 3:
|
|
return vtkm::CELL_SHAPE_TRIANGLE;
|
|
case 4:
|
|
return vtkm::CELL_SHAPE_QUAD;
|
|
default:
|
|
return vtkm::CELL_SHAPE_POLYGON;
|
|
}
|
|
}
|
|
|
|
template <typename CellShapeTag>
|
|
static inline VTKM_EXEC vtkm::VecCConst<vtkm::IdComponent> CellFaceLocalIndices(
|
|
vtkm::IdComponent faceIndex,
|
|
CellShapeTag shape,
|
|
const vtkm::exec::FunctorBase& worklet)
|
|
{
|
|
vtkm::IdComponent numPointsInFace = vtkm::exec::CellFaceNumberOfPoints(faceIndex, shape, worklet);
|
|
if (numPointsInFace < 1)
|
|
{
|
|
// An invalid face. We should already have gotten an error from
|
|
// CellFaceNumberOfPoints.
|
|
return vtkm::VecCConst<vtkm::IdComponent>();
|
|
}
|
|
|
|
return vtkm::make_VecC(detail::PointsInFace[shape.Id][faceIndex], numPointsInFace);
|
|
}
|
|
}
|
|
} // namespace vtkm::exec
|
|
|
|
#endif //vtk_m_exec_CellFace_h
|