Update interpolator to 2nd order and fix bug

GRTLCollaboration · Dec 3, 2024 · b6806f6 · b6806f6
1 parent 58e4d3d
commit b6806f6
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git a/Examples/ScalarFieldBH/CustomExtraction.hpp b/Examples/ScalarFieldBH/CustomExtraction.hpp
@@ -58,7 +58,7 @@ class CustomExtraction
     ~CustomExtraction() {}
 
     //! Execute the query
-    void execute_query(AMRInterpolator<Lagrange<4>> *a_interpolator,
+    void execute_query(AMRInterpolator<Lagrange<2>> *a_interpolator,
                        std::string a_file_prefix) const
     {
         CH_TIME("CustomExtraction::execute_query");

diff --git a/Examples/ScalarFieldBH/ScalarFieldLevel.cpp b/Examples/ScalarFieldBH/ScalarFieldLevel.cpp
@@ -83,7 +83,7 @@ void ScalarFieldLevel::postRestart()
 
     // Use AMR Interpolator and do lineout data extraction
     // pass the boundary params so that we can use symmetries
-    AMRInterpolator<Lagrange<4>> interpolator(
+    AMRInterpolator<Lagrange<2>> interpolator(
         m_gr_amr, m_p.origin, m_p.dx, m_p.boundary_params, m_p.verbosity);
 
     // this should fill all ghosts including the boundary ones according
@@ -115,7 +115,7 @@ void ScalarFieldLevel::postRestart()
                                         m_dt, m_time);
         ham_extraction.execute_query(&interpolator,
                                      m_p.data_path + "Ham_Lineout");
-        CustomExtraction mom_extraction(c_Ham, num_points, m_p.L, m_p.center,
+        CustomExtraction mom_extraction(c_Mom, num_points, m_p.L, m_p.center,
                                         m_dt, m_time);
         mom_extraction.execute_query(&interpolator,
                                      m_p.data_path + "Mom_Lineout");