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              : #include "DataStructures/DataVector.hpp"
       9              : #include "DataStructures/Tensor/Tensor.hpp"
      10              : #include "Elliptic/Systems/Punctures/Tags.hpp"
      11              : #include "Utilities/Gsl.hpp"
      12              : #include "Utilities/TMPL.hpp"
      13              : 
      14              : namespace Punctures {
      15              : 
      16              : /// @{
      17              : /*!
      18              :  * \brief The volume integrand for the ADM mass $M_\mathrm{ADM}$
      19              :  *
      20              :  * The ADM mass for Punctures is (Eq. (12.56) in \cite BaumgarteShapiro)
      21              :  *
      22              :  * \begin{equation}
      23              :  * M_\mathrm{ADM} = \sum_I M_I + \frac{1}{2\pi} \int
      24              :  *   \beta \left(\alpha \left(1 + u\right) + 1\right)^{-7}
      25              :  * \end{equation}
      26              :  *
      27              :  * \see Punctures
      28              :  */
      29            1 : void adm_mass_integrand(const gsl::not_null<Scalar<DataVector>*> result,
      30              :                         const Scalar<DataVector>& field,
      31              :                         const Scalar<DataVector>& alpha,
      32              :                         const Scalar<DataVector>& beta);
      33              : 
      34            1 : Scalar<DataVector> adm_mass_integrand(const Scalar<DataVector>& field,
      35              :                                       const Scalar<DataVector>& alpha,
      36              :                                       const Scalar<DataVector>& beta);
      37              : /// @}
      38              : 
      39              : namespace Tags {
      40              : 
      41              : /// @{
      42            0 : struct AdmMassIntegrand : db::SimpleTag {
      43            0 :   using type = Scalar<DataVector>;
      44              : };
      45              : 
      46              : /*!
      47              :  * \brief The volume integrand for the ADM mass $M_\mathrm{ADM}$
      48              :  *
      49              :  * \see adm_mass_integrand
      50              :  */
      51            1 : struct AdmMassIntegrandCompute : AdmMassIntegrand, db::ComputeTag {
      52            0 :   using base = AdmMassIntegrand;
      53            0 :   using argument_tags = tmpl::list<Field, Alpha, Beta>;
      54            0 :   using return_type = Scalar<DataVector>;
      55            0 :   static constexpr auto function =
      56              :       static_cast<void (*)(gsl::not_null<Scalar<DataVector>*>,
      57              :                            const Scalar<DataVector>&, const Scalar<DataVector>&,
      58              :                            const Scalar<DataVector>&)>(&adm_mass_integrand);
      59              : };
      60              : /// @}
      61              : 
      62              : }  // namespace Tags
      63              : }  // namespace Punctures