[Vxl-users] vnl_sparse_lu least squares

[David D. <do...@rp...>]

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 non-square
vnl_sparse_matrix? And should we add an error if it is passed a
non-square matrix?

Thanks,

David