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