From: David Doria <doriad@rp...>  20110306 23:21:29

It seems like vnl_sparse_lu expects a square matrix  is this correct? I suspect this because in spOutput.c I see things like Size = Matrix>Size; and printf("Size of matrix = %1d x %1d.\n", Size, Size); indicating that it doesn't have a concept of two different dimensions. Is this just an error in the output function? Or is it really only supposed to work with square matrices? If so, is there a way to find the least squares solution to a system Ax = b where A is a nonsquare vnl_sparse_matrix? And should we add an error if it is passed a nonsquare matrix? Thanks, David 