Line data    Source code

       1           0 : // Distributed under the MIT License.
       2             : // See LICENSE.txt for details.
       3             : 
       4             : #pragma once
       5             : 
       6             : #include <array>
       7             : #include <cstddef>
       8             : 
       9             : #include "DataStructures/Tensor/TypeAliases.hpp"
      10             : #include "DataStructures/Variables.hpp"
      11             : #include "Evolution/Systems/GeneralizedHarmonic/ApplyTensorYlmFilter.hpp"
      12             : #include "Evolution/Systems/GeneralizedHarmonic/SphericalShellPowerMonitor.hpp"
      13             : #include "NumericalAlgorithms/Spectral/Mesh.hpp"
      14             : #include "Utilities/Gsl.hpp"
      15             : 
      16             : /// \cond
      17             : class DataVector;
      18             : /// \endcond
      19             : 
      20           0 : namespace gh::power_monitor {
      21             : 
      22             : /// Holds separate GH filled-sphere power monitors for each evolved GH
      23             : /// variable. Index 0 is the radial monitor and index 1 is the angular monitor.
      24           1 : struct GhFilledSpherePowerMonitors {
      25           0 :   std::array<DataVector, 2> spacetime_metric{};
      26           0 :   std::array<DataVector, 2> pi{};
      27           0 :   std::array<DataVector, 2> phi{};
      28             : };
      29             : 
      30             : /*!
      31             :  * \brief Compute radial and angular power monitors for the evolved GH
      32             :  * variables on a filled-sphere element.
      33             :  *
      34             :  * Each monitor is an array with the radial monitor at index 0 and the angular
      35             :  * monitor at index 1. Both monitors are computed after transforming the GH
      36             :  * variables to TensorYlm coefficients in the grid frame using
      37             :  * `jac_inertial_to_grid`.
      38             :  *
      39             :  * **Radial monitor** (index 0): groups squared Jacobi spectral coefficients by
      40             :  * the radial mode index \f$\text{mode} = (n+1)/2\f$ where
      41             :  * \f$n = \ell + 2k\f$ is the total ZernikeB3 degree, \f$\ell\f$ is the
      42             :  * angular degree, and \f$k\f$ is the Jacobi radial index. The result has
      43             :  * \f$n_r\f$ entries.
      44             :  *
      45             :  * **Angular monitor** (index 1): groups squared Jacobi spectral coefficients
      46             :  * by angular degree \f$\ell\f$. Spin-weighted spherical harmonic modes with
      47             :  * \f$\ell < |s|\f$ (where \f$s\f$ is the spin weight of the TensorYlm
      48             :  * component) contribute zero to both the sum and the count, following the same
      49             :  * convention as `gh_shell_power_monitors`. The result has \f$\ell_\max + 1\f$
      50             :  * entries.
      51             :  *
      52             :  * The Cartesian-to-TensorYlm matrices in `cart_to_sphere_matrices` are filled
      53             :  * on first use and reused on subsequent calls.
      54             :  *
      55             :  * \param cart_to_sphere_matrices Cache for Cartesian-to-TensorYlm sparse
      56             :  *   matrices.
      57             :  * \param gh_vars Generalized Harmonic variables at collocation points on the
      58             :  *   B3 ball mesh.
      59             :  * \param mesh The ZernikeB3 mesh.
      60             :  * \param jac_inertial_to_grid Inverse Jacobian mapping inertial-frame spatial
      61             :  *   tensor components to the grid frame.
      62             :  */
      63           1 : GhFilledSpherePowerMonitors gh_filled_sphere_power_monitors(
      64             :     gsl::not_null<CartToSphereMatrices*> cart_to_sphere_matrices,
      65             :     const Variables<ylm::TensorYlm::filter_detail::gh_spacetime_vars_list>&
      66             :         gh_vars,
      67             :     const Mesh<3>& mesh,
      68             :     const InverseJacobian<DataVector, 3, Frame::Inertial, Frame::Grid>&
      69             :         jac_inertial_to_grid);
      70             : 
      71             : }  // namespace gh::power_monitor