I was looking for a port of from cupyx.scipy.linalg import lu_factor, lu_solve for LU factorization & eq. solve via LU factorization

Apparently `lu_factor does exist luckily (thanks!) but not exposed since it is a private function: https://github.com/IntelPython/dpnp/blob/master/dpnp/linalg/dpnp_utils_linalg.py#L874

I would like to request a feature to expose this method and also implement a lu_solve feature possibly via getrs: Ref: https://github.com/cupy/cupy/blob/v13.4.1/cupyx/scipy/linalg/_decomp_lu.py#L259