We introduce solve_nivp, a Python library for robust time integration of nonsmooth ODE/DAE systems with impacts, switching, and inequality constraints. Such systems are widespread: frictional contact in mechanics, earthquakes, piecewise and switching behaviour in circuits, sliding–mode control, discontinuous rules in finance and energy markets…
By embedding nonsmooth rules directly into an implicit time-stepping scheme, the approach avoids regularisation and restrictive step sizes required by classical solvers. The package is supported by documentation, examples, and tests for reproducible use.
For more details: https://joss.theoj.org/papers/10.21105/joss.09775 & https://github.com/ERC-INJECT/solve_nivp & https://pypi.org/project/solve-nivp/

#Python #ScientificComputing #NumericalMethods #ComputationalMechanics #Earthquakes #ContactMechanics #Circuits #SlidingModeControl #Markets #ODE/DAE @DavidRiley @European Research Council (ERC), @Hi! PARIS Center – AI for Science, Business & Society, @ENSTA