Add equality constraints of type `Ax + Bz = c` to `admm` routines

[qntwrsm]

Currently the admm implementation can only handle equality constraints of type $x - z = 0$.