Can anyone tell me a general value for using in DGELSY and DGELSD routines? I suppose that the value depends on the type of problem and real data, but I suppose too that for most problems a common value can be used. For example, in GNU Octave, for solving a least squares problem via DGELSD, a RCOND value of -1 is set (http://hg.savannah.gnu.org/hgweb/octave ... x.cc#l2430). This imply that via DLALSD the actual RCOND value will be the machine epsilon via DLAMCH('Epsilon') (http://www.netlib.org/lapack/double/dlalsd.f). Could be epsilon a good value for common problems?