Add timer to interpolation

Author: nilsvu

inputs/z4c/spectre_bbh/z4c_spectre_bbh.athinput

@@ -6,7 +6,7 @@ basename  = bbh      # problem ID: basename of output filenames
 
 <problem>
 pgen_name = z4c_spectre_bbh
-id_filename_glob = BbhVolume*.h5
+id_filename_glob = /home/nilsvu/2024-09-bbh-athenak/ID/BbhVolume*.h5
 id_subfile_name = VolumeData
 id_observation_step = -1
 

src/pgen/z4c_spectre_bbh.cpp

@@ -10,6 +10,7 @@
 #include <algorithm>
 #include <cmath>
 #include <iostream>
+#include <Kokkos_Timer.hpp>
 #include <sstream>
 #include <string>
 #include <vector>
@@ -135,6 +136,7 @@ void LoadSpectreInitialData(MeshBlockPack *pmbp, const std::string &filename_glo
   x[2].resize(sz);
 
   int nmb = pmbp->nmb_thispack;
+  Kokkos::Timer timer{};
   for (int m = 0; m < nmb; ++m) {
     // Get the coordinates for this meshblock
     Real &x1min = size.h_view(m).x1min;
@@ -155,8 +157,9 @@ void LoadSpectreInitialData(MeshBlockPack *pmbp, const std::string &filename_glo
     }
 
     // Interpolate data to the coordinates
+    timer.reset();
     std::cout << "Interpolating initial data for meshblock " << m << "/"
-              << nmb - 1 << " with " << sz << " points " << std::endl;
+              << nmb - 1 << " with " << sz << " points..." << std::endl;
     const auto data = spectre::Exporter::interpolate_to_points<3>(
         filename_glob, subfile_name, spectre::Exporter::ObservationStep{observation_step},
         {"SpatialMetric_xx", "SpatialMetric_yx", "SpatialMetric_yy",
@@ -165,6 +168,8 @@ void LoadSpectreInitialData(MeshBlockPack *pmbp, const std::string &filename_glo
          "ExtrinsicCurvature_zx", "ExtrinsicCurvature_zy", "ExtrinsicCurvature_zz"},
         /* target_points */ x,
         /* extrapolate_into_excisions */ true);
+    std::cout << "  done in " << timer.seconds() << " seconds." << std::endl;
+
     const auto &gxx = data[0];
     const auto &gyx = data[1];
     const auto &gyy = data[2];