Add scipy.linalg.lu() decomposition support

[abagusetty]

This PR adds dpnp.scipy.linalg.lu() with support for all three output modes: default (P, L, U), permute_l=True (PL, U), and p_indices=True (p, L, U), including batched inputs.

Fixes: #2786

• Have you provided a meaningful PR description?
• Have you added a test, reproducer or referred to an issue with a reproducer?
• Have you tested your changes locally for CPU and GPU devices?
• Have you made sure that new changes do not introduce compiler warnings?
• Have you checked performance impact of proposed changes?
• Have you added documentation for your changes, if necessary?
• Have you added your changes to the changelog?

[intel-python-devops]

Can one of the admins verify this patch?