The recursive version

recursive.py

def get_descendants(node, A, descendants=None):
    if not descendants:
        descendants = []

    descendants.append(node)

    children = A.getrow(node).indices
    for child in children:
        descendants += get_descendants(child, A)

    return descendants