OK, I will give it a try.
One delicate thing (At least when using direct solvers), The LS approach will always result in a solution with the equality forced. The solution with \boldsymbol{E}, in case no solution in the range of the model matrix will result in LS solution. So the equality will be approximated.
I will create a short script and share results.