CD_MirrorDeposition.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_MIRRORDEPOSITION_H
#define CD_MIRRORDEPOSITION_H
 
// Chombo includes
#include <REAL.H>
#include <RealVect.H>
#include <SPACE.H>
 
// Our includes
#include <CD_NamespaceHeader.H>
 
namespace MirrorDeposition {
 
constexpr int numTangents = SpaceDim - 1;
 
constexpr int numShapeComp = (SpaceDim * (SpaceDim + 1)) / 2;
 
constexpr int compStatus = 0;
 
constexpr int compCentroid = 1;
 
constexpr int compNormal = compCentroid + SpaceDim;
 
constexpr int compShape = compNormal + SpaceDim;
 
constexpr int numComp = compShape + numShapeComp;
 
enum class Status
{
  None = 0,
 
  Fitted = 1,
 
  Planar = 2
};
 
enum class Refusal
{
  None,
 
  SmallDenominator,
 
  LargeJacobian,
 
  NonPositiveJacobian
};
 
inline void
liftShapeOperator(const RealVect (&a_tangents)[numTangents],
                  const Real (&a_S)[numTangents][numTangents],
                  Real (&a_shape)[numShapeComp]) noexcept;
 
inline int
shapeIndex(const int a_i, const int a_j) noexcept;
 
inline Real
shapeEntry(const Real* a_shape, const int a_i, const int a_j) noexcept;
 
inline RealVect
applyShapeOperator(const Real* a_shape, const RealVect& a_vec) noexcept;
 
inline Real
meanCurvatureTimesTwo(const Real* a_shape) noexcept;
 
inline Real
gaussianCurvature(const Real* a_shape) noexcept;
 
inline Real
jacobian(const Real a_twoH, const Real a_K, const Real a_d, Real& a_denominator) noexcept;
 
inline bool
reflect(const RealVect& a_pos,
        const Real*     a_surfaceData,
        const Real      a_maxJacobian,
        const Real      a_minDenominator,
        RealVect&       a_image,
        Real&           a_jacobian,
        Real&           a_signedDistance,
        Refusal&        a_refusal) noexcept;
 
} // namespace MirrorDeposition
 
#include <CD_NamespaceFooter.H>
 
#include <CD_MirrorDepositionImplem.H>
 
#endif