https://keroseneheaterreview.blogspot.com/2024/02/removing-loops-in-numpy-for-simple.html