Я недавно начал использовать numpy и пока что не являюсь опытным пользователем, нужны ваши советы) Есть две матрицы А и В, каждая строчка в которых - это координаты определённого вектора (по сути, есть два типа атомов и наборы их координат). Количество атомов в каждой матрице(то есть количество строк каждой матрицы) в районе 10000. Но в дальнейшем, возможно, придётся делать расчеты для нескольких сотен тысяч атомов. Каким способом (в плане быстродействия и оптимального использования расчетных ресурсов) лучше всего найти все расстояния между каждым вектором матрицы А и каждым вектором матрицы В?
Пока что написал такую функцию, которая вычисляет расстояния и возвращает количество расстояний, меньших заранее заданной длины (defaultLength). Но мне кажется, что это не самый быстрый способ расчёта, возможно ли как-то ускорить код?
def Bond(a, b, defaultLength): #a и b - матрицы с координатами атомов A и B
def Length(x):
return np.sqrt(x[0]**2 + x[2]**2 + x[2]**2)
i = 0
for row in b:
c = a - row
d = np.apply_along_axis(Length, axis=1, arr=c)
for item in d:
if item < defaultLength: i = i + 1
return i