def test_find_pore_to_pore_distance(self): net = op.network.Cubic(shape=[3, 3, 3], connectivity=6) dm = topotools.find_pore_to_pore_distance(network=net, pores1=net.pores('left'), pores2=net.pores('right')) a = np.unique(dm) b = np.array([2., 2.23606798, 2.44948974, 2.82842712, 3., 3.46410162]) assert np.allclose(a, b)