Here is one approach that reuses the array in every call of rng.standard_normal:
rng.standard_normal
import numpy as np rng = np.random.default_rng() size = 1000 arr = np.empty(size) rng.standard_normal(size, out=arr)