Re: New keyword name for linalg.lu

Ilhan Polat <[email protected]> Mon, 22 May 2023 09:58:17 +0200
Newsgroups gmane.comp.python.scientific.devel
Message-ID <CAEBuzr-vn-KeB-xJ=6aJB4QE8WMZ+6v06oMojrhtx0nwf7Miwg@mail.gmail.com>
Yes indeed that's something I checked but in sparse case there is no
possibility of full array return due to the size constraints hence, say in
SuperLU, they only have "perm_r" and "perm_c" integer 1D arrays anyways for
row and column permutation indices. Also permutation index is a thing in
combinatorics so I have

p_as_vector = False/True
return_indices = False/True

so far that doesn't sound confusing and also not that ugly.


On Mon, May 22, 2023 at 8:37 AM Todd Bailey <[email protected]>
wrote:

> Conceptually this seems related to sparse matrices so I wonder if there is
> some helpful terminology there.  Your 1D option would return “column
> indices” for the non-zero entries in each row of the 2D option.
>
> Best,
> Todd
>
> On 21 May 2023, at 16:08, Ilhan Polat <[email protected]> wrote:
>
> [snip]
>
> 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.
> [snip]
> Could you please offer some alternatives even just for inspiration?
>
> 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]