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 <cstdint>
       8             : 
       9             : /// \cond
      10             : class Matrix;
      11             : namespace Spectral {
      12             : enum class Basis : uint8_t;
      13             : enum class Quadrature : uint8_t;
      14             : }  // namespace Spectral
      15             : /// \endcond
      16             : 
      17             : namespace ader::dg {
      18             : /*!
      19             :  * \brief Computes the matrix applied to the spacetime predictor volume
      20             :  * contributions to give the next initial guess in the iteration.
      21             :  *
      22             :  * Specifically, the returned matrix is:
      23             :  *
      24             :  * \f{align}{
      25             :  * S_{a_0 b_0}=\left(\ell_{a_0}(1)\ell_{c_0}(1) - w_{c_0}
      26             :  * \left(\partial_\tau\ell_{a_0}(\tau)\right)\rvert_{c_0}
      27             :  * \right)^{-1}M_{c_{0}b_{0}}
      28             :  * \f}
      29             :  *
      30             :  * where \f$M_{a_0 b_0}\f$ is the temporal mass matrix.
      31             :  */
      32             : template <Spectral::Basis BasisType, Spectral::Quadrature QuadratureType>
      33           1 : const Matrix& predictor_inverse_temporal_matrix(size_t num_points);
      34             : }  // namespace ader::dg