https://keroseneheaterreview.blogspot.com/2024/01/optimize-this-function-with-numpy-or.html