Update MovingPunctureGauge adaptive parameters in BinaryBHLevel.cpp

GRTLCollaboration · Apr 2, 2021 · a803fe5 · a803fe5
1 parent 6483081
commit a803fe5
Showing 1 changed file with 6 additions and 0 deletions.
diff --git a/Examples/BinaryBH/BinaryBHLevel.cpp b/Examples/BinaryBH/BinaryBHLevel.cpp
@@ -49,6 +49,10 @@ void BinaryBHLevel::initialData()
     // Can't use simd with this initial data
     BoxLoops::loop(two_punctures_initial_data, m_state_new, m_state_new,
                    INCLUDE_GHOST_CELLS, disable_simd());
+    m_p.ccz4_params.puncture_masses[0] = m_tp_amr.m_two_punctures.mm;
+    m_p.ccz4_params.puncture_masses[1] = m_tp_amr.m_two_punctures.mp;
+    m_p.ccz4_params.puncture_coords =
+        m_bh_amr.m_puncture_tracker.get_puncture_coords();
 #else
     // Set up the compute class for the BinaryBH initial data
     BinaryBH binary(m_p.bh1_params, m_p.bh2_params, m_dx);
@@ -199,6 +203,8 @@ void BinaryBHLevel::specificPostTimeStep()
         bool write_punctures = at_level_timestep_multiple(coarsest_level);
         m_bh_amr.m_puncture_tracker.execute_tracking(m_time, m_restart_time,
                                                      m_dt, write_punctures);
+        m_p.ccz4_params.puncture_coords =
+            m_bh_amr.m_puncture_tracker.get_puncture_coords();
     }
 }