Null space of a sparse matrix

calculates null space and range of a rectangular sparse matrix

finds null space and range of a sparse matrix using lu decomposition.

This routine is devoted to find null space of possibly large matrix fast.
Results may be inaccurate if zero and nonzero singular values are not well separated.