FWIW, I'd suggest to reconsider adding the keyword at all.
Maybe return a permutation vector (ideally, a LAPACK convention) and
add a helper function to convert it to a permutation matrix?
To be clear, `lu()` currently returns the full permutation matrix. The new proposed functionality is to control changing the format of this returned object, not to control whether or not something gets returned.