added BoostedBHComplexProca example

GRTLCollaboration · Mar 20, 2024 · 02b659c · 02b659c
1 parent 6bf5db5
commit 02b659c
Show file tree

Hide file tree

Showing 8 changed files with 25 additions and 25 deletions.
diff --git a/Examples/BoostedBHComplexProca/BoostedBHComplexProcaLevel.cpp b/Examples/BoostedBHComplexProca/BoostedBHComplexProcaLevel.cpp
@@ -20,14 +20,14 @@
 #include "FixedGridsTaggingCriterion.hpp"
 
 // Problem specific includes
+#include "ComplexProcaField.hpp"
 #include "EnergyConservation.hpp"
-#include "ExcisionDiagnostics.hpp"
 #include "ExcisionComplexProcaEvolution.hpp"
+#include "ExcisionDiagnostics.hpp"
 #include "FluxExtraction.hpp"
 #include "InitialComplexProcaData.hpp"
 #include "LinearMomConservation.hpp"
 #include "ProcaConstraint.hpp"
-#include "ComplexProcaField.hpp"
 
 // Diagnostics
 #include "ProcaDiagnostic.hpp"
@@ -155,8 +155,8 @@ void BoostedBHComplexProcaLevel::specificPostTimeStep()
 
 // Things to do in RHS update, at each RK4 step
 void BoostedBHComplexProcaLevel::specificEvalRHS(GRLevelData &a_soln,
-                                          GRLevelData &a_rhs,
-                                          const double a_time)
+                                                 GRLevelData &a_rhs,
+                                                 const double a_time)
 {
     // Calculate right hand side with matter_t = ComplexProcaField
     // and background_t = BoostedBH

diff --git a/Examples/BoostedBHComplexProca/DiagnosticVariables.hpp b/Examples/BoostedBHComplexProca/DiagnosticVariables.hpp
@@ -16,7 +16,7 @@ enum
     c_fluxLinMom,
     c_sourceLinMom,
     c_Zvec_Re_out, // Auxiliary Z field. Used for tracking (minus) the Gauss
-                // constraint
+                   // constraint
     c_Zvec_Im_out,
 
     NUM_DIAGNOSTIC_VARS

diff --git a/Examples/BoostedBHComplexProca/InitialComplexProcaData.hpp b/Examples/BoostedBHComplexProca/InitialComplexProcaData.hpp
@@ -9,8 +9,8 @@
 #include "ADMFixedBGVars.hpp"
 #include "BoostedBH.hpp"
 #include "Cell.hpp"
-#include "Coordinates.hpp"
 #include "ComplexProcaField.hpp"
+#include "Coordinates.hpp"
 #include "Tensor.hpp"
 #include "TensorAlgebra.hpp"
 #include "UserVariables.hpp" //This files needs NUM_VARS - total no. components
@@ -37,13 +37,12 @@ class InitialProcaData
   public:
     //! The constructor for the class
     InitialProcaData(
-        const double a_amplitude,
-        const double a_mass,
+        const double a_amplitude, const double a_mass,
         const std::array<double, CH_SPACEDIM> a_center,
         const BoostedBH::params_t a_bg_params, const double a_dx,
         const std::string a_proca_initial_data_profile = "uniform-x")
-        : m_dx(a_dx), m_amplitude(a_amplitude), m_mass{a_mass}, m_center(a_center),
-          m_bg_params(a_bg_params),
+        : m_dx(a_dx), m_amplitude(a_amplitude), m_mass{a_mass},
+          m_center(a_center), m_bg_params(a_bg_params),
           m_proca_initial_data_profile(a_proca_initial_data_profile)
     {
     }
@@ -84,15 +83,16 @@ class InitialProcaData
             vars.Avec_Re[2] = m_amplitude;
         }
 
-        //Since Evec is minus the conjugate momentum to Avec, we set E = - dAvec/dt
-        // Assuming Avec is a phasor, i.e. Avec[i] ~ e^(-i*w*t)
+        // Since Evec is minus the conjugate momentum to Avec, we set E = -
+        // dAvec/dt
+        //  Assuming Avec is a phasor, i.e. Avec[i] ~ e^(-i*w*t)
         FOR1(i)
         {
             vars.Evec_Im[i] = 0;
-            
+
             FOR1(j)
             {
-                vars.Evec_Im[i] += - gamma_UU[i][j] * m_mass * vars.Avec_Re[j];
+                vars.Evec_Im[i] += -gamma_UU[i][j] * m_mass * vars.Avec_Re[j];
             }
         }
 

diff --git a/Examples/BoostedBHComplexProca/Main_BoostedBHComplexProca.cpp b/Examples/BoostedBHComplexProca/Main_BoostedBHComplexProca.cpp
@@ -27,8 +27,8 @@ int runGRChombo(int argc, char *argv[])
     // (To simulate a different problem, define a new child of AMRLevel
     // and an associated LevelFactory)
     GRAMR gr_amr;
-    DefaultLevelFactory<BoostedBHComplexProcaLevel> proca_field_level_fact(gr_amr,
-                                                                    sim_params);
+    DefaultLevelFactory<BoostedBHComplexProcaLevel> proca_field_level_fact(
+        gr_amr, sim_params);
     setupAMRObject(gr_amr, proca_field_level_fact);
 
     // call this after amr object setup so grids known

diff --git a/Examples/BoostedBHComplexProca/ProcaConstraint.hpp b/Examples/BoostedBHComplexProca/ProcaConstraint.hpp
@@ -9,8 +9,8 @@
 #include "ADMFixedBGVars.hpp"
 #include "BoostedBH.hpp"
 #include "Cell.hpp"
-#include "Coordinates.hpp"
 #include "ComplexProcaField.hpp"
+#include "Coordinates.hpp"
 #include "Tensor.hpp"
 #include "UserVariables.hpp" //This files needs NUM_VARS - total no. components
 #include "VarsTools.hpp"
@@ -72,11 +72,10 @@ class ProcaConstraint
             Avec0_Re += -d1.Evec_Re[i][i];
             Avec0_Im += -d1.Evec_Im[i][i];
 
-            FOR1(j) 
-            { 
-              Avec0_Re += -chris_phys.ULL[i][i][j] * vars.Evec_Re[j]; 
-              Avec0_Im += -chris_phys.ULL[i][i][j] * vars.Evec_Im[j];
-
+            FOR1(j)
+            {
+                Avec0_Re += -chris_phys.ULL[i][i][j] * vars.Evec_Re[j];
+                Avec0_Im += -chris_phys.ULL[i][i][j] * vars.Evec_Im[j];
             }
         }
         Avec0_Re = Avec0_Re / m_mu / m_mu;

diff --git a/Examples/BoostedBHComplexProca/ProcaDiagnostic.hpp b/Examples/BoostedBHComplexProca/ProcaDiagnostic.hpp
@@ -4,8 +4,8 @@
 #include "ADMFixedBGVars.hpp"
 #include "BoostedBH.hpp"
 #include "Cell.hpp"
-#include "Coordinates.hpp"
 #include "ComplexProcaField.hpp"
+#include "Coordinates.hpp"
 #include "Tensor.hpp"
 #include "UserVariables.hpp" //This files needs NUM_VARS - total no. components
 #include "VarsTools.hpp"

diff --git a/Examples/BoostedBHComplexProca/SimulationParameters.hpp b/Examples/BoostedBHComplexProca/SimulationParameters.hpp
@@ -78,7 +78,8 @@ class SimulationParameters : public FixedBGSimulationParametersBase
         }
     }
 
-    // Problem specific parameters - the radii for the integrations (Should always equal the radii of extraction)
+    // Problem specific parameters - the radii for the integrations (Should
+    // always equal the radii of extraction)
     double inner_r, outer_r;
     double proca_mass, proca_damping, proca_amplitude;
     // Collection of parameters necessary for the initial conditions

diff --git a/Examples/BoostedBHComplexProca/UserVariables.hpp b/Examples/BoostedBHComplexProca/UserVariables.hpp
@@ -47,7 +47,7 @@ static const std::array<std::string, NUM_VARS> variable_names = {
 
     "Evec1_Re", "Evec1_Im", "Evec2_Re", "Evec2_Im", "Evec3_Re", "Evec3_Im",
 
-    "Zvec_Re", "Zvec_Im"};
+    "Zvec_Re",  "Zvec_Im"};
 
 } // namespace UserVariables