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