def METHOD_NAME(self):
    # One line is fully covered by another
    p_1 = np.ones(3)
    p_2 = 2 * p_1
    p_3 = 0 * p_1
    p_4 = 3 * p_1
    p_int = pp.intersections.segments_3d(p_1, p_2, p_3, p_4)
    p_known_1 = p_1.reshape((-1, 1))
    p_known_2 = p_2.reshape((-1, 1))
    self.assertTrue(np.min(np.sum(np.abs(p_int - p_known_1), axis=0)) < 1e-8)
    self.assertTrue(np.min(np.sum(np.abs(p_int - p_known_2), axis=0)) < 1e-8)