proj_hs function

Orthogonal projection matrix onto the half-space