Update dpnp.linalg.solve() to align NumPy 2.0

[vlad-perevezentsev]

This PR proposes to update dpnp.linalg.solve() implementation to align with NumPy >= 2.0 and also Python array API standard .

In numpy 2.0 the broadcast ambiguity has been removed and now b is treaded as a single vector if and only if it is 1-dimensional; for other cases this signature must be followed (..., m, m), (..., m, n) -> (..., m, n)

This is accomplished by using dpnp.broadcast_to() for input arrays in the current suggestion.

Running array api tests for dpnp.linalg.solve() with NumPy 2.1

$ pytest -v array_api_tests/test_linalg.py::test_solve
========================== 1 passed, 1 warning in 0.35s ===================================

• Have you provided a meaningful PR description?
• Have you added a test, reproducer or referred to 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?
• If this PR is a work in progress, are you filing the PR as a draft?

[github-actions]

View rendered docs @ https://intelpython.github.io/dpnp/index.html

[antonwolfy]

Thank you @vlad-perevezentsev