return pivot array

`factorize_lu` returns the pivot array alongwith the LU decomposition if `:pivot` is passed as an argument.