Procs
proc gelsd[T: SomeFloat](a, b: Tensor[T]; solution, residuals: var Tensor[T]; singular_values: var Tensor[T]; matrix_rank: var int; rcond = -1.T)
-
Wrapper around LAPACK's gelsd taking tensors and preparing the required work space.
rcond is the condition for singular values to be considered zero, s(i) <= rcond * s(i) are treated as zero.
If rcond = -1 is used, it determines the size automatically (to the machine precision).
Source Edit