def get_all_delaunay_2d(vertices, N_s, k, w, r):
tangents = find_tangents(vertices, k)
curve = find_curve(vertices, tangents, k, w)
edges, v_s, sparse = delaunay_complex_2d(vertices, N_s)
t_s = tangents[sparse]
c_s = curve[sparse]
return v_s, t_s, c_s, edges
def get_all_witness_4d(vertices, N_s, k, w, r):
tangents = find_tangents(vertices, k)
curve = find_curve(vertices, tangents, k, w)
edges, v_s, sparse = witness_complex_4d(vertices, tangents, w=w, N_s=N_s)
t_s = tangents[sparse]
c_s = curve[sparse]
return v_s, t_s, c_s, edges
def get_all_rips_2d(vertices, N_s, k, w, r):
tangents = find_tangents(vertices, k)
curve = find_curve(vertices, tangents, k, w)
sparse = sparse_sample(vertices, N_s)
v_s = vertices[sparse, :]
t_s = tangents[sparse]
c_s = curve[sparse]
edges = rips_complex_2d(v_s, r=r)
return v_s, t_s, c_s, edges