Start of BinaryWithGravitationalWaves class for XCTS

- class to implement background data
- right now it implements superposition of two boosted black holes
- no radiative content is added for now

Additional BoundaryCondition class to set on inner boundaries: SuperposedBoostedBinary

- Dirichlet is implemented
- Neumann is not implemented (set to zero)
- Neumann will be implemented later

docs/References.bib

@@ -457,6 +457,20 @@ @article{Buchman:2012dw
     year = "2012"
 }
 
+@article{Buonanno2005xu,
+    author = "Buonanno, Alessandra and Chen, Yanbei and Damour, Thibault",
+    title = "{Transition from inspiral to plunge in precessing binaries of
+              spinning black holes}",
+    eprint = "gr-qc/0508067",
+    archivePrefix = "arXiv",
+    doi = "10.1103/PhysRevD.74.104005",
+    journal = "Phys. Rev. D",
+    volume = "74",
+    pages = "104005",
+    year = "2006",
+    url = "https://doi.org/10.1103/PhysRevD.74.104005"
+}
+
 @article{Casoni2012,
   author  = {Casoni, E. and Peraire, J. and Huerta, A.},
   title   = {One-dimensional shock-capturing for high-order discontinuous
@@ -1345,6 +1359,21 @@ @Article{Hou2007
   url     = "https://doi.org/10.1007/s10915-006-9105-9"
 }
 
+@article{Jaranowski1997ky,
+    author = "Jaranowski, Piotr and Schaefer, Gerhard",
+    title = "{Third postNewtonian higher order ADM Hamilton dynamics for
+             two-body point mass systems}",
+    eprint = "gr-qc/9712075",
+    archivePrefix = "arXiv",
+    doi = "10.1103/PhysRevD.57.7274",
+    url = "https://doi.org/10.1103/PhysRevD.57.7274",
+    journal = "Phys. Rev. D",
+    volume = "57",
+    pages = "7274--7291",
+    year = "1998",
+    note = "[Erratum: Phys.Rev.D 63, 029902 (2001)]"
+}
+
 @article{Jiang1996,
   author  = {Guang-Shan Jiang and Chi-Wang Shu},
   doi     = {https://doi.org/10.1006/jcph.1996.0130},
@@ -1375,6 +1404,22 @@ @article{Kastaun2020uxr
   year =         2021
 }
 
+@article{Kelly2007uc,
+    author = "Kelly, Bernard J. and Tichy, Wolfgang and Campanelli,
+             Manuela and Whiting, Bernard F.",
+    title = "{Black hole puncture initial data with realistic gravitational
+             wave content}",
+    eprint = "0704.0628",
+    archivePrefix = "arXiv",
+    primaryClass = "gr-qc",
+    doi = "10.1103/PhysRevD.76.024008",
+    url = "https://doi.org/10.1103/PhysRevD.76.024008",
+    journal = "Phys. Rev. D",
+    volume = "76",
+    pages = "024008",
+    year = "2007"
+}
+
 @article{Kennedy2003,
   author = "Kennedy, Christopher A. and Carpenter, Mark H.",
   title = "Additive {Runge-Kutta} schemes for convection-diffusion-reaction
@@ -2006,6 +2051,22 @@ @article{Muhlberger2014pja
   year =         2014
 }
 
+@article{Mundim2010hu,
+    author = "Mundim, Bruno C. and Kelly, Bernard J. and Zlochower,
+              Yosef and Nakano, Hiroyuki and Campanelli, Manuela",
+    editor = "Lehner, Luis and Pfeiffer, Harald P. and Poisson, E.",
+    title = "{Hybrid black-hole binary initial data}",
+    eprint = "1012.0886",
+    archivePrefix = "arXiv",
+    primaryClass = "gr-qc",
+    doi = "10.1088/0264-9381/28/13/134003",
+    url = "https://doi.org/10.1088/0264-9381/28/13/134003",
+    journal = "Class. Quant. Grav.",
+    volume = "28",
+    pages = "134003",
+    year = "2011"
+}
+
 @article{Nee2024bur,
   author        = {Nee, Peter James and Lara, Guillermo and Pfeiffer, Harald
                     P. and Vu, Nils L.},
@@ -2620,6 +2681,21 @@ @book{ThorneBlandford2017
   publisher  = "Princeton University Press",
 }
 
+@article{Tichy2002ec,
+    author = "Tichy, Wolfgang and Bruegmann, Bernd and Campanelli, Manuela
+             and Diener, Peter",
+    title = "{Binary black hole initial data for numerical general relativity
+             based on postNewtonian data}",
+    eprint = "gr-qc/0207011",
+    archivePrefix = "arXiv",
+    reportNumber = "UTBRG-2002-06, AEI-2002-048",
+    doi = "10.1103/PhysRevD.67.064008",
+    journal = "Phys. Rev. D",
+    volume = "67",
+    pages = "064008",
+    year = "2003"
+}
+
 @article{Toro1994,
   author         = "Toro, E. F. and Spruce, M. and Speares, W.",
   title          = "{Restoration of the Contact Surface in the HLL–Riemann

src/Elliptic/Executables/Xcts/SolveXcts.hpp

@@ -43,6 +43,7 @@
 #include "ParallelAlgorithms/LinearSolver/Multigrid/ElementsAllocator.hpp"
 #include "ParallelAlgorithms/LinearSolver/Multigrid/Tags.hpp"
 #include "PointwiseFunctions/AnalyticData/Xcts/Binary.hpp"
+#include "PointwiseFunctions/AnalyticData/Xcts/BinaryWithGravitationalWaves.hpp"
 #include "PointwiseFunctions/AnalyticSolutions/Xcts/Factory.hpp"
 #include "PointwiseFunctions/Hydro/LowerSpatialFourVelocity.hpp"
 #include "PointwiseFunctions/Hydro/Tags.hpp"
@@ -118,6 +119,9 @@ struct Metavariables {
         Xcts::Solutions::all_analytic_solutions,
         Xcts::AnalyticData::Binary<elliptic::analytic_data::AnalyticSolution,
                                    Xcts::Solutions::all_analytic_solutions>,
+        Xcts::AnalyticData::BinaryWithGravitationalWaves<
+            elliptic::analytic_data::AnalyticSolution,
+            Xcts::Solutions::all_analytic_solutions>,
         elliptic::analytic_data::NumericData>;
 
     using factory_classes = tmpl::map<

src/Elliptic/Systems/Xcts/BoundaryConditions/CMakeLists.txt

@@ -10,6 +10,7 @@ spectre_target_sources(
   PRIVATE
   ApparentHorizon.cpp
   Flatness.cpp
+  SuperposedBoostedBinary.cpp
   Robin.cpp
   )
 
@@ -20,6 +21,7 @@ spectre_target_headers(
   ApparentHorizon.hpp
   Factory.hpp
   Flatness.hpp
+  SuperposedBoostedBinary.hpp
   Robin.hpp
   )
 

src/Elliptic/Systems/Xcts/BoundaryConditions/Factory.hpp

@@ -7,13 +7,17 @@
 #include "Elliptic/Systems/Xcts/BoundaryConditions/ApparentHorizon.hpp"
 #include "Elliptic/Systems/Xcts/BoundaryConditions/Flatness.hpp"
 #include "Elliptic/Systems/Xcts/BoundaryConditions/Robin.hpp"
+#include "Elliptic/Systems/Xcts/BoundaryConditions/SuperposedBoostedBinary.hpp"
+#include "PointwiseFunctions/AnalyticSolutions/Xcts/Factory.hpp"
+#include "PointwiseFunctions/InitialDataUtilities/AnalyticSolution.hpp"
 #include "Utilities/TMPL.hpp"
 
 namespace Xcts::BoundaryConditions {
 template <typename System>
-using standard_boundary_conditions =
-    tmpl::list<elliptic::BoundaryConditions::AnalyticSolution<System>,
-               Flatness<System::enabled_equations>,
-               Robin<System::enabled_equations>,
-               ApparentHorizon<System::conformal_geometry>>;
+using standard_boundary_conditions = tmpl::list<
+    elliptic::BoundaryConditions::AnalyticSolution<System>,
+    Flatness<System::enabled_equations>, Robin<System::enabled_equations>,
+    ApparentHorizon<System::conformal_geometry>,
+    SuperposedBoostedBinary<elliptic::analytic_data::AnalyticSolution,
+                            Xcts::Solutions::all_analytic_solutions>>;
 }