Re: New keyword name for linalg.lu

Ilhan Polat <[email protected]> Mon, 22 May 2023 21:06:43 +0200
Newsgroups gmane.comp.python.scientific.devel
Message-ID <CAEBuzr9xv+8pnF2sZYwOBdTVQRP-S6U8uFLRL_RJSvRE2kGk+w@mail.gmail.com>
In the end I went with "matrix_p = True" to ring like "permute_l" which was
I think what you were suggesting anyways. Hopefully some day in the future
we can change this to False and be done with it.

On Sun, May 21, 2023 at 11:31 PM Andras Deak <[email protected]> wrote:

> On Sun, May 21, 2023, at 17:08, Ilhan Polat wrote:
> > 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?
>
> How about reversing the meaning of the flag and naming it something like
> `p_full` or `full_p` with a `True` default? I also considered "squareform"
> instead of "full" as in `spatial.distance.squareform`, because there's a
> similar 1d <-> 2d square thing here, but the semantics are different enough
> that this name would probably be more confusing than helpful.
>
> András
>
> >
> > Best,
> > ilhan
> >
> >
> >
> > On Tue, Apr 25, 2023 at 9:34 AM Jake Bowhay <[email protected]>
> wrote:
> >> It would be nice to add a quick note to the docs explaining when/which
> you should use. Currently both state "Compute pivoted LU decomposition of a
> matrix." which while true isn't very helpful for a user trying to decide
> which function to pick!
> >> _______________________________________________
> >> SciPy-Dev mailing list -- [email protected]
> >> To unsubscribe send an email to [email protected]
> >> https://mail.python.org/mailman3/lists/scipy-dev.python.org/
> >> Member address: [email protected]
> > _______________________________________________
> > SciPy-Dev mailing list -- [email protected]
> > To unsubscribe send an email to [email protected]
> > https://mail.python.org/mailman3/lists/scipy-dev.python.org/
> > Member address: [email protected]
> _______________________________________________
> SciPy-Dev mailing list -- [email protected]
> To unsubscribe send an email to [email protected]
> https://mail.python.org/mailman3/lists/scipy-dev.python.org/
> Member address: [email protected]
>

_______________________________________________
SciPy-Dev mailing list -- [email protected]
To unsubscribe send an email to [email protected]
https://mail.python.org/mailman3/lists/scipy-dev.python.org/
Member address: [email protected]