A robust implementation for solving the $S$-unit equation and several applications