https://fairytalefashionheadband.blogspot.com/2023/10/avoiding-double-for-loops-in-numpy.html