Sorry for reviving this thread, but indeed "p_indices" sounding uglier the more I look at it. Just to summarize what it does to save you from the thread above;

when

P, L, U = scipy.linalg.lu(A)

is run, currently, P is returning a full 2D array. If A is a tall array say, (25, 5) then P is necessarily (25, 25). And it is just a permutaiton matrix, a row shuffled np.eye(25). Instead, you can ask with this new keyword to return that shuffle pattern. as a 1D array and hence P becomes (25, ) array. And there are more than rare cases that only L and U is needed for other purposes.

So we would like to make that matrix vs array switch with a keyword. I thought "p_indices = True/False" was a good idea but it became a bit disturbing. Other matlab following software offers a string selector with "outputForm='matrix'/'vector'" which forces users to remember the options hence I am voting for a TrueFalse flag. However I can't come up with the right keyword.

Could you please offer some alternatives even just for inspiration?

Best,

ilhan