Джулия: распределение случайных полиномиальных чисел

Мне нужно сгенерировать вектор случайного числа, взятый из полиномиального распределения. Поскольку мне нужно сделать это внутри цикла, я хотел бы уменьшить количество выделений. Это мой код:

pop = zeros(100)
p = rand(100)
p /= sum(p)

#imagine the next line inside a for loop
@time pop = rand(Multinomial(100, p)) #test1
@time pop .= rand(Multinomial(100, p)) #test2

Почему test1 — это 2 распределения, а test2 — 4? Есть ли способ сделать это с 0 распределениями?

Заранее спасибо.