The other option was using sqlsq from scipy, the only method that accepts equality constraints. However, it was so painfully slow that it became unusable. Why don't we pass it without the option and then work on a PR to add it?