CD_PolyhedralEBUtils.H

Go to the documentation of this file.

/*
 * SPDX-FileCopyrightText: 2021-2026 SINTEF Energy Research
 *
 * SPDX-License-Identifier: GPL-3.0-or-later
 */
 
#ifndef CD_POLYHEDRALEBUTILS_H
#define CD_POLYHEDRALEBUTILS_H
 
// Chombo includes
#include <RealVect.H>
 
// Our includes
#include <CD_CutCellSurface.H>
#include <CD_NamespaceHeader.H>
 
namespace PolyhedralEB {
 
inline bool
isFluid(const Real a_value) noexcept;
 
namespace detail {
 
#if CH_SPACEDIM == 3
constexpr int s_transverse[3][2] = {{1, 2}, {0, 2}, {0, 1}};
#endif
 
inline int
edgeDirection(const int a_edge) noexcept;
 
inline void
edgeOrigin(const int a_edge, int a_offset[SpaceDim]) noexcept;
 
inline void
edgeCorners(const int a_edge, int& a_lo, int& a_hi) noexcept;
 
inline RealVect
cornerPosition(const int a_corner) noexcept;
 
inline void
faceCorners(const int a_dir, const int a_side, int a_corner[1 << (SpaceDim - 1)]) noexcept;
 
inline int
edgeIndex(const int a_dir, const int a_offset[SpaceDim]) noexcept;
 
inline RealVect
crossingPosition(const CutCellSurface& a_surface, const int a_edge, const Real a_tolerance) noexcept;
 
#if CH_SPACEDIM == 3
inline void
faceEdges(const int a_dir, const int a_side, int a_edge[4]) noexcept;
 
inline void
polygonMoments(const RealVect* a_vertex,
               const int       a_num,
               Real&           a_area,
               RealVect&       a_vector,
               RealVect&       a_centroid) noexcept;
 
inline int
facePairs(const int a_dir, const int a_side, const CutCellSurface& a_surface, int a_pair[2][2]) noexcept;
 
inline int
crossingLoops(const CutCellSurface& a_surface,
              int                   a_loop[CutCellSurface::s_numEdges],
              int                   a_start[CutCellSurface::s_numEdges + 1]) noexcept;
#endif
constexpr Real s_weldTolerance = 1.0E-10;
 
inline bool
sameVertex(const RealVect& a_a, const RealVect& a_b) noexcept;
 
} // namespace detail
} // namespace PolyhedralEB
 
#include <CD_NamespaceFooter.H>
 
#include <CD_PolyhedralEBUtilsImplem.H>
 
#endif