def reset(self, reset_events=False):
      """Reset the blocks inputs and outputs and also its internal solver, 
      if the block has a solver instance.
      """
      #reset inputs and outputs
      self.inputs.reset()
      self.outputs.reset()

      #reset engine if block has solver
      if self.engine: self.engine.reset()

      #reset operators if defined
      if self.op_alg: self.op_alg.reset()
      if self.op_dyn: self.op_dyn.reset()
      
      #reset events (re-arm event schedulers with a kwarg for backwards compatibility)
      if reset_events:
          for event in self.events:
              event.reset()

from pathsim import Simulation, Connection
from pathsim.blocks import StepSource, Scope
from pathsim.solvers import SSPRK22

step = StepSource(amplitude=1.0, tau=1.0)
scope = Scope()
sim = Simulation([step, scope], [Connection(step[0], scope[0])], Solver=SSPRK22, dt=0.1)

# First run works
sim.run(duration=5.0)
t1, y1 = scope.read()
print(f"First run: max={y1.max()}")  # Should be 1

# Second run after reset produces zeros
sim.reset(time=0.0)
sim.run(duration=5.0)
t2, y2 = scope.read()
print(f"Second run: max={y2.max()}")  # This is now 0 (BUG)