GSL currently does not support a singular value decomposition when M<N. Any alternative, and GSL-supported, methods to find the entire solution space of a system where M<N? Brian?