DGELSx family of functions: overdetermined matrix eqns

[RRLapackuser]

I'm failry new to LAPACK and have a quick question on the DGELSx family of functions for solving an overdetermined set of equations (Ax=b). It's a linear least squares (LLS) problem.

For my case, matrix A has more rows than columns.
b and x are column vectors

I'd like to solve them only if the rank of A = #columns in A
I also want to know the error (residual) in the solution (for full rank)

DGELS doesn't give me info about rank although error norm info is available (overwritten on B)

DGELSY gives me the rank. But is there a way to find out the error in the solution?

DGELSS seems like it might work for me. It gives rank and error info (overwritten on B) but it actually gives me a solution that is very off. Besides it computes SVD which I don't need.

Any suggestions for which function I should use?