VM
Я думал так ```
n = 1```
for i = -n:n
for j = -n:n
for k = -n:n
Vt = a1 + [i, j, k]
end
end
end
Size: a a a
2020 January 20
VM
Но так не получается
АО
Куча вариантов есть.
Например так
Например так
translations(x, n = 1) = reshape(map(t -> x .+ t, Iterators.product(-n:n, -n:n, -n:n)), :)
translations((0, 0, 0)) # Vector{Tuple{Int64,Int64,Int64}} with 27 elements
АО
1000? Это тяжело... Получается 2001^3 трансляций = 8*10^9
VM
Блин, прям трансляция есть... Это не язык, а находка!
АО
Возможно стоит оптимизировать алгоритм как-то.
VM
1000? Это тяжело... Получается 2001^3 трансляций = 8*10^9
Возможно стоит, но до этого я еще доберусь
VM
1000? Это тяжело... Получается 2001^3 трансляций = 8*10^9
при 50 Scilab ложится и дергается в конвульсиях
АО
50 -> 10^6 трансляций
x = translations((0, 0, 0), 50) # Vector{Tuple{Int64,Int64,Int64}} with 1030301 elements
Base.summarysize(x)/1024/1024 # 23.6 MbАО
То есть трансляции одной точки занимают 23.6 мегабайта.
АО
При достаточно большом наборе точек очень лего OOM словить.
KT
вы уверены, что вам нужно их генерировать все?
KT
Может, вы просто отождествите те точки, которые сдвигами друг в друга переходят
KT
проверить это гораздо проще, чем в явном виде это множество генерить
VM
проверить это гораздо проще, чем в явном виде это множество генерить
Попробую. Спасибо!
KT
т.е. написать тип Point{T} свой и переопределить для него isequal
VM
т.е. написать тип Point{T} свой и переопределить для него isequal
До этого доберусь попозже. Пока мне нужно n =1
VM
Позже построенную ячейку буду размножать, вот тогда воспользуюсь советом. Спасибо!
VM
