проверить эвклидовое расстояние через периодические границы

Question

Каков наилучший способ проверить эвклидовое расстояние между двумя точками, когда граничные условия являются периодическими? Я попытался взять минимум расстояния и (дальность - расстояние), но я думаю, из-за того, как я структурировал программу, это дает странный результат. Я надеюсь, что есть еще один хороший способ приблизиться к этому, что я могу принять, а не переоценивать остальную часть функции.

Это в Python кстати. В настоящее время поиск эвклидовой дистанции с numpy.linalg.norm, хотя есть подпрограммы piist SciPy, которые делают то же самое, что я могу использовать, я думаю.

Answer 1

Вы должны отдельно проверять размер каждого измерения на то, что меньше между dx и range-dx.

def distance(p1, p2): 
    total = 0 
    for i, (a, b) in enumerate(zip(p1, p2)): 
     delta = abs(b - a) 
     if delta > dimension[i] - delta: 
      delta = dimension[i] - delta 
     total += delta ** 2 
    return total ** 0.5