Как улучшить производительность scipy.sparse.csgraph.depth_first_order?

Question

import numpy as np
from scipy.sparse import csgraph, dia_matrix

n = 1000
g = dia_matrix((np.ones(n), 1), shape=(n, n)).tocsr()
for i in range(100):
    for j in range(n):
        csgraph.breadth_first_order(g, j)

Я написал этот фрагмент кода и профилировал его.

Функция deep_first_order провела большую часть времени при проверке графа g Я прошел в.

Как улучшить производительность этой функции?

Есть ли способ заставить Сципи проверять график только один раз?