https://keroseneheaterreview.blogspot.com/2023/10/vectorizing-triple-for-loop-in.html