wip
return pivot array
Merge branch 'lu' into solve
returns permutation matrix