Matrix Rank Calculation and GSL

[Lopaka(Rob) Lee]

I'm a student learning numerical methods with the GSL.

I'd like to be able to calculate matrix rank - however, no GSL routines 
directly support this.  I realize that the non-zero singular values from 
SVD will give me matrix rank, but the GSL SVD does not support cases where M<N!

The most useful algorithm for obtaining matrix rank I've used is 
the QR-decomposition sub 'dqrdc2' (from LINPACK I believe).

Any suggestions?

-L