Add pivoting

[rreusser]

QR with pivoting is a thing that exists, but it's not a thing I'm doing. Is this just a matter of defining and using a permutation vector? I expect it's a minimally to moderately invasive change.