code: def _dm_rref_GJ_sparse(M):
      """Compute RREF using sparse Gauss-Jordan elimination with division."""
      M_rref_d, pivots, _ = sdm_irref(M.rep)
      M_rref_sdm = SDM(M_rref_d, M.shape, M.domain)
      pivots = tuple(pivots)
      return M.from_rep(M_rref_sdm), pivots
  
file_path: sympy/sympy/polys/matrices/rref.py
imports: from sympy.polys.domains import ZZ
  from sympy.polys.matrices.sdm import SDM, sdm_irref, sdm_rref_den
  from sympy.polys.matrices.ddm import DDM
  from sympy.polys.matrices.dense import ddm_irref, ddm_irref_den
label_desc: Use image node_id 6 for calling a global function with example usage: _dm_rref_GJ_sparse(M) and returns: M,  pivots