Line data    Source code

       1            0 : // Distributed under the MIT License.
       2              : // See LICENSE.txt for details.
       3              : 
       4              : #pragma once
       5              : 
       6              : #include <ostream>
       7              : 
       8              : #include "DataStructures/DataBox/Access.hpp"
       9              : #include "DataStructures/Variables.hpp"
      10              : #include "NumericalAlgorithms/DiscontinuousGalerkin/Formulation.hpp"
      11              : #include "Utilities/ErrorHandling/Assert.hpp"
      12              : #include "Utilities/Gsl.hpp"
      13              : #include "Utilities/TMPL.hpp"
      14              : 
      15              : namespace evolution::dg::subcell {
      16              : template <typename... CorrectionTags, typename BoundaryCorrection,
      17              :           typename... PackageFieldTags, typename... BoundaryTermsVolumeTags>
      18            0 : void compute_boundary_terms(
      19              :     const gsl::not_null<Variables<tmpl::list<CorrectionTags...>>*>
      20              :         boundary_corrections_on_face,
      21              :     const BoundaryCorrection& boundary_correction,
      22              :     const Variables<tmpl::list<PackageFieldTags...>>& upper_packaged_data,
      23              :     const Variables<tmpl::list<PackageFieldTags...>>& lower_packaged_data,
      24              :     const db::Access& box, tmpl::list<BoundaryTermsVolumeTags...> /*meta*/) {
      25              :   ASSERT(
      26              :       upper_packaged_data.number_of_grid_points() ==
      27              :           lower_packaged_data.number_of_grid_points(),
      28              :       "The number of grid points must be the same for the upper packaged data ("
      29              :           << upper_packaged_data.number_of_grid_points()
      30              :           << ") and the lower packaged data ("
      31              :           << lower_packaged_data.number_of_grid_points() << ')');
      32              :   ASSERT(upper_packaged_data.number_of_grid_points() ==
      33              :              boundary_corrections_on_face->number_of_grid_points(),
      34              :          "The number of grid points must be the same for the packaged data ("
      35              :              << upper_packaged_data.number_of_grid_points()
      36              :              << ") and the boundary corrections on the faces ("
      37              :              << boundary_corrections_on_face->number_of_grid_points() << ')');
      38              :   boundary_correction.dg_boundary_terms(
      39              :       make_not_null(&get<CorrectionTags>(*boundary_corrections_on_face))...,
      40              :       get<PackageFieldTags>(upper_packaged_data)...,
      41              :       get<PackageFieldTags>(lower_packaged_data)...,
      42              :       // FD schemes are basically weak form FV scheme
      43              :       ::dg::Formulation::WeakInertial,
      44              :       db::get<BoundaryTermsVolumeTags>(box)...);
      45              : }
      46              : }  // namespace evolution::dg::subcell