Functionality related to discontinuous Galerkin schemes. More...
Namespaces | |
| protocols | |
| Protocols related to Discontinuous Galerkin functionality | |
Classes | |
| struct | SimpleBoundaryData |
| Distinguishes between field data, which can be projected to a mortar, and extra data, which will not be projected. More... | |
| class | SimpleMortarData |
| Storage of boundary data on two sides of a mortar. More... | |
Typedefs | |
| template<size_t VolumeDim> | |
| using | MortarId = std::pair<::Direction< VolumeDim >, ElementId< VolumeDim > > |
| template<size_t MortarDim> | |
| using | MortarSize = std::array< Spectral::MortarSize, MortarDim > |
| template<size_t VolumeDim, typename ValueType > | |
| using | MortarMap = std::unordered_map< MortarId< VolumeDim >, ValueType, boost::hash< MortarId< VolumeDim > >> |
Enumerations | |
| enum | Formulation { StrongInertial, WeakInertial } |
| The DG formulation to use. More... | |
Functions | |
| std::ostream & | operator<< (std::ostream &os, Formulation t) noexcept |
| template<typename InboxTag , size_t Dim, typename TemporalIdType , typename... InboxTags> | |
| bool | has_received_from_all_mortars (const TemporalIdType &temporal_id, const Element< Dim > &element, const tuples::TaggedTuple< InboxTags... > &inboxes) noexcept |
Determines if data on all mortars has been received for the InboxTag at time temporal_id. More... | |
| template<size_t Dim> | |
| void | metric_identity_det_jac_times_inv_jac (gsl::not_null< InverseJacobian< DataVector, Dim, Frame::Logical, Frame::Inertial > * > det_jac_times_inverse_jacobian, const Mesh< Dim > &mesh, const tnsr::I< DataVector, Dim, Frame::Inertial > &inertial_coords, const Jacobian< DataVector, Dim, Frame::Logical, Frame::Inertial > &jacobian) noexcept |
| Compute the Jacobian determinant times the inverse Jacobian so that the result is divergence-free. More... | |
| template<size_t Dim> | |
| void | metric_identity_jacobian_quantities (gsl::not_null< InverseJacobian< DataVector, Dim, Frame::Logical, Frame::Inertial > * > det_jac_times_inverse_jacobian, gsl::not_null< InverseJacobian< DataVector, Dim, Frame::Logical, Frame::Inertial > * > inverse_jacobian, gsl::not_null< Jacobian< DataVector, Dim, Frame::Logical, Frame::Inertial > * > jacobian, gsl::not_null< Scalar< DataVector > * > det_jacobian, const Mesh< Dim > &mesh, const tnsr::I< DataVector, Dim, Frame::Inertial > &inertial_coords) noexcept |
| Compute the Jacobian, inverse Jacobian, and determinant of the Jacobian so that they satisfy the metric identities. More... | |
| template<size_t Dim> | |
| Mesh< Dim > | mortar_mesh (const Mesh< Dim > &face_mesh1, const Mesh< Dim > &face_mesh2) noexcept |
| template<size_t Dim> | |
| MortarSize< Dim - 1 > | mortar_size (const ElementId< Dim > &self, const ElementId< Dim > &neighbor, size_t dimension, const OrientationMap< Dim > &orientation) noexcept |
| template<typename Tags , size_t Dim> | |
| Variables< Tags > | project_to_mortar (const Variables< Tags > &vars, const Mesh< Dim > &face_mesh, const Mesh< Dim > &mortar_mesh, const MortarSize< Dim > &mortar_size) noexcept |
| template<typename Tags , size_t Dim> | |
| Variables< Tags > | project_from_mortar (const Variables< Tags > &vars, const Mesh< Dim > &face_mesh, const Mesh< Dim > &mortar_mesh, const MortarSize< Dim > &mortar_size) noexcept |
Functionality related to discontinuous Galerkin schemes.
|
noexcept |
Apply the DG mass matrix to the data, in the diagonal mass-matrix approximation ("mass-lumping")
The DG mass matrix is:
\begin{equation} M_{pq} = \int_{\Omega_k} \psi_p(\xi) \psi_q(\xi) \mathrm{d}V \end{equation}
where \(\psi_p(\xi)\) are the basis functions on the element \(\Omega_k\). In the diagonal mass-matrix approximation ("mass-lumping") we evaluate the integral directly on the collocation points, i.e. with a Gauss or Gauss-Lobatto quadrature determined by the element mesh. Then it reduces to:
\begin{equation} M_{pq} \approx \delta_{pq} \prod_{i=1}^d w_{p_i} \end{equation}
where \(d\) is the spatial dimension and \(w_{p_i}\) are the quadrature weights in dimension \(i\). To apply the mass matrix in coordinates different than logical, or to account for a curved background metric, the data can be pre-multiplied with the Jacobian determinant and/or a metric determinant.
|
noexcept |
Apply the DG mass matrix to the data, in the diagonal mass-matrix approximation ("mass-lumping")
The DG mass matrix is:
\begin{equation} M_{pq} = \int_{\Omega_k} \psi_p(\xi) \psi_q(\xi) \mathrm{d}V \end{equation}
where \(\psi_p(\xi)\) are the basis functions on the element \(\Omega_k\). In the diagonal mass-matrix approximation ("mass-lumping") we evaluate the integral directly on the collocation points, i.e. with a Gauss or Gauss-Lobatto quadrature determined by the element mesh. Then it reduces to:
\begin{equation} M_{pq} \approx \delta_{pq} \prod_{i=1}^d w_{p_i} \end{equation}
where \(d\) is the spatial dimension and \(w_{p_i}\) are the quadrature weights in dimension \(i\). To apply the mass matrix in coordinates different than logical, or to account for a curved background metric, the data can be pre-multiplied with the Jacobian determinant and/or a metric determinant.