Line data    Source code

       1            0 : // Distributed under the MIT License.
       2              : // See LICENSE.txt for details.
       3              : 
       4              : #pragma once
       5              : 
       6              : #include <cstddef>
       7              : 
       8              : #include "DataStructures/DataVector.hpp"
       9              : 
      10              : namespace elliptic::dg {
      11              : 
      12              : /*!
      13              :  * \brief The penalty factor in internal penalty fluxes
      14              :  *
      15              :  * The penalty factor is computed as
      16              :  *
      17              :  * \f{equation}
      18              :  * \sigma = C \frac{N_\text{points}^2}{h}
      19              :  * \f}
      20              :  *
      21              :  * where \f$N_\text{points} = 1 + N_p\f$ is the number of points (or one plus
      22              :  * the polynomial degree) and \f$h\f$ is a measure of the element size. Both
      23              :  * quantities are taken perpendicular to the face of the DG element that the
      24              :  * penalty is being computed on. \f$C\f$ is the "penalty parameter". For details
      25              :  * see section 7.2 in \cite HesthavenWarburton.
      26              :  */
      27            1 : DataVector penalty(const DataVector& element_size, size_t num_points,
      28              :                    double penalty_parameter);
      29              : 
      30              : }  // namespace elliptic::dg