solve_nivp: A Python toolkit for integrating nonsmooth dynamical systems (Journal of Open Source Software, 2026)