Added tracker for steps taken during integration

CHANGES.md

@@ -2,6 +2,11 @@
 
 ## 2024
 
+#### v0.11.4 (2024-11-27)
+
+New:
+* Added a `steps_taken` tracking variable to the C++ class `CySolverResult` and the Cython wrapped `WrapCySolverResult` so that users can see how many steps were taken during integration even when using `t_eval`.
+
 #### v0.11.4 (2024-11-25)
 
 C++ Changes

CyRK/cy/cysolution.cpp

@@ -133,8 +133,9 @@ void CySolverResult::reset()
     this->storage_capacity = this->original_expected_size;
 
     // Reset trackers
-    this->size = 0;
+    this->size             = 0;
     this->num_interpolates = 0;
+    this->steps_taken      = 0;
 
     // Reserve the memory for the vectors
     try

CyRK/cy/cysolution.hpp

@@ -60,8 +60,9 @@ class CySolverResult : public std::enable_shared_from_this<CySolverResult>{
 
     // More status information
     char* message_ptr = &message[0];
-    size_t size = 0;
+    size_t size             = 0;
     size_t num_interpolates = 0;
+    size_t steps_taken      = 0;
 
     // Pointer to storage arrays
     std::vector<double> time_domain_vec        = std::vector<double>();

CyRK/cy/cysolver.cpp

@@ -324,6 +324,7 @@ void CySolverBase::take_step()
 
             this->p_step_implementation();
             this->len_t++;
+            this->storage_sptr->steps_taken++;
 
             // Take care of dense output and t_eval
             int dense_built = 0;

CyRK/cy/cysolver_api.pxd

@@ -85,6 +85,7 @@ cdef extern from "cysolution.cpp" nogil:
             char* message_ptr
             size_t size
             size_t num_interpolates
+            size_t steps_taken
             vector[double] time_domain_vec
             vector[double] time_domain_vec_sorted
             vector[double] solution

CyRK/cy/cysolver_api.pyx

@@ -72,6 +72,10 @@ cdef class WrapCySolverResult:
     def size(self):
         return self.cyresult_ptr.size
 
+    @property
+    def steps_taken(self):
+        return self.cyresult_ptr.steps_taken
+
     @property
     def num_y(self):
         return self.cyresult_ptr.num_y

Tests/D_PySolver_Tests/test_a_pysolve_ivp.py

@@ -199,6 +199,7 @@ def test_pysolve_ivp(use_scipy_style, use_args, use_njit_always,
     assert result.success
     assert result.error_code == 1
     assert result.size > 1
+    assert result.steps_taken > 1
     assert result.message == "Integration completed without issue."
     # Check that the ndarrays make sense
     assert type(result.t) == np.ndarray

pyproject.toml

@@ -1,6 +1,6 @@
 [project]
 name='CyRK'
-version = '0.11.4'
+version = '0.11.5'
 description='Runge-Kutta ODE Integrator Implemented in Cython and Numba.'
 authors= [
     {name = 'Joe P. Renaud', email = 'joe.p.renaud@gmail.com'}