If you're pulling in LAPACK anyway, wouldn't just calling dgels to compute the least-squares solution directly be easier?