x, y, z = (x * scale_x + offset_x, y * scale_y + offset_y, z * scale_z +
    offset_z)
distances = np.sqrt((x - x0) ** 2 + (y - y0) ** 2 + (z - z0) ** 2)
Fastest way to fill numpy array with distances from a point