Вы можете сделать это без for для l oop, используя функцию numpy.linalg.norm
x = np.array([[5, 4, 1],
[3, 2, 1]])
y = np.array([[5, 4, 2],
[3, 3, 2]])
np.linalg.norm(x - y, axis=1)
>> array([1. , 1.41421356])
Чтобы сделать то же самое, но по столбцам, используйте axis=0
np.linalg.norm(x - y, axis=0)
>>array([0. , 1. , 1.41421356])