log a warning when an incoming dt gets restricted

reneSchm · reneSchm · commit b93c0e6294c3 · 2024-03-08T15:25:13.000+01:00
diff --git a/cpp/memilio/math/adapt_rk.cpp b/cpp/memilio/math/adapt_rk.cpp
@@ -76,6 +76,11 @@ bool RKIntegratorCore::step(const DerivFunction& f, Eigen::Ref<const Eigen::Vect
 {
     assert(0 <= m_dt_min);
     assert(m_dt_min <= m_dt_max);
+
+    if (dt < m_dt_min || dt > m_dt_max) {
+        mio::log_warning("IntegratorCore: Restricting given step size dt = {} to [{}, {}].", dt, m_dt_min, m_dt_max);
+    }
+
     dt = std::min(dt, m_dt_max);
 
     double t_eval; // shifted time for evaluating yt
diff --git a/cpp/memilio/math/stepper_wrapper.h b/cpp/memilio/math/stepper_wrapper.h
@@ -22,6 +22,7 @@
 
 #include "memilio/utils/compiler_diagnostics.h"
 #include "memilio/math/integrator.h"
+#include "memilio/utils/logging.h"
 
 GCC_CLANG_DIAGNOSTIC(push)
 GCC_CLANG_DIAGNOSTIC(ignored "-Wshadow")
@@ -56,7 +57,6 @@ class ControlledStepperWrapper : public mio::IntegratorCore
     static_assert(!is_fsal_stepper,
                   "FSAL steppers cannot be used until https://github.com/boostorg/odeint/issues/72 is resolved.");
 
-
 public:
     /**
      * @brief Set up the integrator
@@ -90,6 +90,11 @@ class ControlledStepperWrapper : public mio::IntegratorCore
         using boost::numeric::odeint::fail;
         assert(0 <= m_dt_min);
         assert(m_dt_min <= m_dt_max);
+
+        if (dt < m_dt_min || dt > m_dt_max) {
+            mio::log_warning("IntegratorCore: Restricting given step size dt = {} to [{}, {}].", dt, m_dt_min,
+                             m_dt_max);
+        }
         // set initial values for exit conditions
         auto step_result = fail;
         bool is_dt_valid = true;
