I'm having quite a headache with the equality constraints lagrange multipliers coming from the quadprog solver. The example to reproduce the issue is as follows: import numpy as np from quadprog ...