docs: def lu_solve(LU, b):
      r"""Solve for LU decomposition.
  
      Solve the linear equations :math:`\mathrm A \mathbf x = \mathbf b`,
      given the LU factorization of :math:`\mathrm A`.
  
      Args:
          LU (array_like): LU decomposition.
          b (array_like): Right-hand side.
  
      Returns:
          :class:`numpy.ndarray`: The solution to the system
          :math:`\mathrm A \mathbf x = \mathbf b`.
  
      See Also
      --------
      scipy.linalg.lu_factor : LU decomposition.
      scipy.linalg.lu_solve : Solve linear equations given LU factorization.
      """
      from scipy.linalg import lu_solve as sp_lu_solve
  
      LU = (asarray(LU[0], float), asarray(LU[1], float))
      b = asarray(b, float)
      return sp_lu_solve(LU, b, check_finite=False)