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