I could not find any mention of how to compute the least norm solution of an underdetermined linear system. The column-pivoted QR decomposition is a building block, but does not take the additional steps of computing a complete orthogonal decomposition. What is the best way to go about implementing this?