proc gesv[T: SomeFloat](a, b: var Tensor[T]; pivot_indices: var seq[int32])
Wrapper for LAPACK *gesv routines Solve AX = B for general matrix
In-place version, this will overwrite a and b
