Even more print statements

.github/workflows/tests.yml

@@ -24,7 +24,7 @@ jobs:
     strategy:
       fail-fast: false
       matrix:
-        os: [ubuntu-latest]
+        os: []
         # Test other versions of Python in special cases to avoid exploding the
         # matrix size; make sure to test all supported versions in some form.
         python-version: ["3.11"]

qutip/solver/integrator/scipy_integrator.py

@@ -305,6 +305,7 @@ def set_state(self, t, state0):
             state0.to_array().ravel().view(np.float64),
             t
         )
+        print(self._ode_solver._y)
         print(self._ode_solver.y)
         print(' --- --- ')
 
@@ -447,7 +448,11 @@ def _one_step(self, t):
         if t > self._front and t_ode >= self._front:
             # The state is at self._front, do a step
             self._back = self.get_state()
+            print(' --- THERE --- ')
+            print(self._ode_solver._y)
+            print(self._ode_solver.y)
             self._ode_solver.integrate(min(self._front + safe_delta, t))
+            print(' --- --- ')
             self._front = self._ode_solver._integrator.rwork[12]
             # We asked for a fraction of a step, now complete it.
             self._ode_solver.integrate(min(self._front, t))

qutip/tests/solver/test_mcsolve.py

@@ -516,7 +516,9 @@ def test_MCSolver_stepping():
     H = qutip.num(size)
     solver = MCSolver(H, a)
     solver.start(qutip.basis(size, size-1), 0, seed=0)
+    print('a', solver.integrator._ode_solver._y)
     solver.options = {'method': 'lsoda'}
+    print('b', solver.integrator._ode_solver._y)
     state = solver.step(1)
     assert qutip.expect(qutip.qeye(size), state) == pytest.approx(1)
     assert qutip.expect(qutip.num(size), state) == pytest.approx(size - 1)