Given a matrix M, let M' be the matrix created by M.rref() Let E be a matrix made from composing the elementary row operations which took M to M.rref(). In particular, E satisfies M.rref() = EM I now wish to compute E (actually, EA, for any matrix A of the right dimension). Does Sagemath offer anything which allows me to bypass manually computing E?