'Helper function of function diagonal_solve,
without the error checks, to be used privately.'

def _diagonal_solve(self, rhs):

    return self._new(rhs.rows, rhs.cols, (lambda i, j: (rhs[(i, j)] / self[(i, i)])))
