Lowell Johnson writes: > I think there is an unnecessary limitation in the QR solution code > (linalg/qr.c). The functions gsl_linalg_QR_svx() and gsl_linalg_QR_solve() > only support square matrices. Shouldn't the input matrix be allowed to be M x > N, M >= N ? > Good point -- I will add that. Thanks.