The orthogonal projection matrix \(P\) has unit determinant and
transforms an n by d matrix by taking \(x * P\).
The components of the first column vector of the resulting matrix are strictly positive.
proj_hs(beta, inv = FALSE)a d by d orthogonal projection matrix
vector defining the half-space
logical; if TRUE, return the inverse matrix