Изчисляването на евклидови разстояния с Python работи твърде бавно

Чета набори от данни от файл в numpy масиви като този:

def read_data(filename):
   data = np.empty(shape=[0, 65], dtype=int)
   with open(filename) as f:
       for line in f:
           data = np.vstack((data, np.array(list(map(int, line.split(','))), dtype=int)))
   return data

Използвам numpy за изчисляване на евклидовото разстояние между два списъка:

def euclidean_distance(x, z):
   return np.linalg.norm(x-z)

След това изчислявам евклидовите разстояния по следния начин:

for data in testing_data:
   for data2 in training_data:
       dist = euclidean_distance(data, data2)

Проблемът ми е, че този код работи много бавно, отнема около ~10 минути, за да завърши. Как мога да подобря това, какво пропускам?
Трябва да използвам разстоянията в друг алгоритъм, така че скоростта е много важна.