- we have found that the 1e-7 comparison can cause some problems. sometimes it is True, sometimes it is False. like 2 > 2 - 1e-7 is true, 4 > 4 - 1e-7 is false..?
- it seems that the first element of indptr will not be sampled...?
for (EdgeIDType k = e_search; k > std::max(s_search - 1, e_search - neighs); k--)