LSP

Author: MSABuschmann

source/sledgehamr.cpp

@@ -97,6 +97,7 @@ void Sledgehamr::Evolve() {
         amrex::Print() << "Full step took " << utils::DurationSeconds(timer)
                        << "s.\n"
                        << std::endl;
+        last_full_time = utils::DurationSeconds(timer);
         io_module->Write();
 
 #ifdef AMREX_MEM_PROFILING

source/sledgehamr.h

@@ -245,7 +245,7 @@ class Sledgehamr : public amrex::AmrCore {
     /** @brief Initialize project specific details. To be overriden by each
      *         project if required.
      */
-    virtual void Init(){};
+    virtual void Init() {};
 
     /** @brief Function to check whether a given level is allowed to be created
      *         by a given time. Will always allow a level to be created by
@@ -269,7 +269,7 @@ class Sledgehamr : public amrex::AmrCore {
     /** @brief Any work that should be done before a coarse level time step.
      *         To be overriden by the project class.
      */
-    virtual void BeforeTimestep(const double time){};
+    virtual void BeforeTimestep(const double time) {};
 
     /** @brief Stopping condition of the simulation. Can be overriden by the
      *         project class.
@@ -284,7 +284,7 @@ class Sledgehamr : public amrex::AmrCore {
      * @param   lev     Current level.
      */
     virtual void SetParamsRhs(std::vector<double> &params, const double time,
-                              const int lev){};
+                              const int lev) {};
 
     /** @brief Function to be overriden by the project class to set user-defined
      *         parameters that will be passed to the Rhs computation of
@@ -295,7 +295,7 @@ class Sledgehamr : public amrex::AmrCore {
      */
     virtual void SetParamsGravitationalWaveRhs(std::vector<double> &params,
                                                const double time,
-                                               const int lev){};
+                                               const int lev) {};
 
     /** @brief Function to be overriden by the project class to set user-defined
      *         parameters that will be passed to the tagging procedure.
@@ -305,7 +305,7 @@ class Sledgehamr : public amrex::AmrCore {
      */
     virtual void SetParamsTagCellForRefinement(std::vector<double> &params,
                                                const double time,
-                                               const int lev){};
+                                               const int lev) {};
 
     /** @brief Function to be overriden by the project class to set user-defined
      *         parameters that will be available when modifying the truncation
@@ -316,7 +316,7 @@ class Sledgehamr : public amrex::AmrCore {
      */
     virtual void SetParamsTruncationModifier(std::vector<double> &params,
                                              const double time,
-                                             const int lev){};
+                                             const int lev) {};
 
     /** @brief Function to be overriden by the project class to set user-defined
      *         parameters that will be passed to the spectrum computation.
@@ -325,7 +325,7 @@ class Sledgehamr : public amrex::AmrCore {
      * @param   lev     Current level.
      */
     virtual void SetParamsSpectra(std::vector<double> &params,
-                                  const double time){};
+                                  const double time) {};
 
     /** @brief Function to be overriden by the project class to set user-defined
      *         parameters that will be passed to the computation of projections.
@@ -334,7 +334,7 @@ class Sledgehamr : public amrex::AmrCore {
      * @param   lev     Current level.
      */
     virtual void SetParamsProjections(std::vector<double> &params,
-                                      const double time){};
+                                      const double time) {};
 
     void CreateShadowLevel();
 
@@ -443,6 +443,8 @@ class Sledgehamr : public amrex::AmrCore {
      */
     bool restart_sim = false;
 
+    double last_full_time = 0;
+
   private:
     void DoErrorEstCpu(int lev, amrex::TagBoxArray &tags, double time);
     void DoErrorEstGpu(int lev, amrex::TagBoxArray &tags, double time);

source/timer.cpp

@@ -8,7 +8,7 @@ void Timer::Start() {
     if (is_running)
         return;
 
-    //amrex::ParallelDescriptor::Barrier();
+    // amrex::ParallelDescriptor::Barrier();
     start_time = std::chrono::steady_clock::now();
     is_running = true;
 }
@@ -19,9 +19,9 @@ void Timer::Stop() {
     if (!is_running)
         return;
 
-    //amrex::ParallelDescriptor::Barrier();
+    // amrex::ParallelDescriptor::Barrier();
     CheckClock();
-    total_micro +=  static_cast<double>(last_duration_micro.count());
+    total_micro += static_cast<double>(last_duration_micro.count());
     is_running = false;
 }
 
@@ -36,27 +36,26 @@ double Timer::GetTotalTimeSeconds() {
         extra_micro = static_cast<double>(last_duration_micro.count());
     }
 
-    return (total_micro + extra_micro)/1e6;
+    return (total_micro + extra_micro) / 1e6;
 };
 
-/** @brief Same as GetTotalTimeSeconds() but only time since the timer was 
+/** @brief Same as GetTotalTimeSeconds() but only time since the timer was
  *         last started.
  */
 double Timer::GetLastDurationSeconds() {
     if (is_running)
         CheckClock();
 
-    return static_cast<double>(last_duration_micro.count())/1e6;
+    return static_cast<double>(last_duration_micro.count()) / 1e6;
 };
 
 /** @brief Computes time passed between now and the last starting time point
            internally.
  */
 void Timer::CheckClock() {
     stop_time = std::chrono::steady_clock::now();
-    last_duration_micro =
-            std::chrono::duration_cast<std::chrono::microseconds>(
-                stop_time - start_time);
+    last_duration_micro = std::chrono::duration_cast<std::chrono::microseconds>(
+        stop_time - start_time);
 }
 
 }; // namespace sledgehamr