https://discountsnowblowersthisinstant.blogspot.com/2023/11/avoiding-double-for-loops-in-numpy.html