Python: Fast alternatives to GEKKO for solving Algebraic differential equations

Question

To simulate power flows in grids, I use GEKKO to solve my algebraic differential equation system.

For smaller simulations, it works pretty well, but running it several ti