https://keroseneheaterreview.blogspot.com/2022/07/removing-loops-in-numpy-for-simple.html