import numpy as np
from scipy.spatial.distance import cdist, pdist, squareform

np.set_printoptions(precision=3, suppress=True)

points = np.array([
    [0.0, 0.0],
    [3.0, 4.0],
    [6.0, 8.0],
])

condensed = pdist(points, metric="euclidean")

print("pdist euclidean")
print(condensed)

print("\nsquareform rows")
for row in squareform(condensed):
    print(row)

queries = np.array([
    [0.0, 4.0],
    [9.0, 12.0],
])

print("\ncdist euclidean rows")
for row in cdist(queries, points, metric="euclidean"):
    print(row)

print("\npdist cityblock")
print(pdist(points, metric="cityblock"))