Added fix to array lifetimes

CHANGES.md

@@ -5,6 +5,7 @@
 #### v0.12.1 (2024-12-02)
 
 Fixes:
+* Added fix so that arrays returned by `WrapCySolverResult` class are not destroyed when class is. Instead they are managed by python's garbage collector once they are held by a python variable.
 
 Tests:
 * Added test to check that arrays stay alive after underlying `WrapCySolverResult` class is destroyed.

CyRK/cy/cysolver_api.pyx

@@ -124,11 +124,11 @@ cdef class WrapCySolverResult:
 
     @property
     def t(self):
-        return np.asarray(self.time_view, dtype=np.float64, order='C')
+        return np.copy(np.asarray(self.time_view, dtype=np.float64, order='C'))
 
     @property
     def y(self):
-        return np.asarray(self.y_view, dtype=np.float64, order='C').reshape((self.size, self.num_dy)).T
+        return np.copy(np.asarray(self.y_view, dtype=np.float64, order='C')).reshape((self.size, self.num_dy)).T
 
     @property
     def size(self):

pyproject.toml

@@ -1,6 +1,6 @@
 [project]
 name='CyRK'
-version = '0.12.0'
+version = '0.12.1'
 description='Runge-Kutta ODE Integrator Implemented in Cython and Numba.'
 authors= [
     {name = 'Joe P. Renaud', email = '[email protected]'}