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             : #include <utility>
       8             : 
       9             : #include "DataStructures/DataVector.hpp"
      10             : #include "DataStructures/Variables.hpp"
      11             : #include "Domain/Structure/Direction.hpp"
      12             : #include "Domain/Structure/Side.hpp"
      13             : #include "NumericalAlgorithms/Spectral/BoundaryInterpolationTerm.hpp"
      14             : #include "NumericalAlgorithms/Spectral/Mesh.hpp"
      15             : #include "NumericalAlgorithms/Spectral/Quadrature.hpp"
      16             : #include "Utilities/Algorithm.hpp"
      17             : #include "Utilities/Gsl.hpp"
      18             : 
      19             : namespace dg {
      20             : namespace detail {
      21             : template <size_t Dim>
      22             : void interpolate_dt_terms_gauss_points_impl(
      23             :     gsl::not_null<double*> volume_dt_vars, size_t num_independent_components,
      24             :     const Mesh<Dim>& volume_mesh, size_t dimension, size_t num_boundary_pts,
      25             :     const gsl::span<const double>& dt_corrections,
      26             :     const DataVector& boundary_interpolation_term);
      27             : }  // namespace detail
      28             : 
      29             : /*!
      30             :  * \brief Interpolate the Bjorhus/time derivative corrections to the volume time
      31             :  * derivatives in the specified direction.
      32             :  *
      33             :  * The general interpolation term (for the \f$+\xi\f$-dimension) is:
      34             :  *
      35             :  * \f{align*}{
      36             :  *   \partial_t u_{\alpha\breve{\imath}\breve{\jmath}\breve{k}}=\cdots
      37             :  *   +\ell^{\mathrm{Gauss-Lobatto}}_{N}
      38             :  *    \left(\xi_{\breve{\imath}}^{\mathrm{Gauss}}\right)
      39             :  *    \partial_t u^{\mathrm{BC}}_{\alpha\breve{\jmath}\breve{k}},
      40             :  * \f}
      41             :  *
      42             :  * where \f$\breve{\imath}\f$, \f$\breve{\jmath}\f$, and \f$\breve{k}\f$ are
      43             :  * indices in the logical \f$\xi\f$, \f$\eta\f$, and \f$\zeta\f$ dimensions.
      44             :  * \f$\partial_t u^{\mathrm{BC}}\f$ is the time derivative correction, and
      45             :  * the  function Spectral::boundary_interpolation_term() is used to compute and
      46             :  * cache the terms from the lifting terms.
      47             :  */
      48             : template <size_t Dim, typename DtTagsList>
      49           1 : void interpolate_dt_terms_gauss_points(
      50             :     const gsl::not_null<Variables<DtTagsList>*> dt_vars,
      51             :     const Mesh<Dim>& volume_mesh, const Direction<Dim>& direction,
      52             :     const Variables<DtTagsList>& dt_corrections) {
      53             :   ASSERT(alg::all_of(volume_mesh.quadrature(),
      54             :                      [](const Spectral::Quadrature quadrature) {
      55             :                        return quadrature == Spectral::Quadrature::Gauss;
      56             :                      }),
      57             :          "Must use Gauss points in all directions but got the mesh: "
      58             :              << volume_mesh);
      59             :   const Mesh<Dim - 1> boundary_mesh =
      60             :       volume_mesh.slice_away(direction.dimension());
      61             :   const Mesh<1> volume_stripe_mesh =
      62             :       volume_mesh.slice_through(direction.dimension());
      63             :   const size_t num_boundary_grid_points = boundary_mesh.number_of_grid_points();
      64             :   detail::interpolate_dt_terms_gauss_points_impl(
      65             :       make_not_null(dt_vars->data()), dt_vars->number_of_independent_components,
      66             :       volume_mesh, direction.dimension(), num_boundary_grid_points,
      67             :       gsl::make_span(dt_corrections.data(), dt_corrections.size()),
      68             :       direction.side() == Side::Upper
      69             :           ? Spectral::boundary_interpolation_term(volume_stripe_mesh).second
      70             :           : Spectral::boundary_interpolation_term(volume_stripe_mesh).first);
      71             : }
      72             : }  // namespace dg