Sorry for the late reply on my behalf. this is a new forum for me. I really do appreciate all your replies.
Ultimately, I feel like numpy is well versed for floating point operations, but falls short for integer operations when overflow and memory consumption is a concern (a very common concern for image processing).
Honestly, I didn’t know that numpy was able to use gemm with syntax like:
npMat3 = npMat4 = npMat5 = npTmp + npTmp*1j
%timeit npMat3.T @ npMat4 + npMat5
I’ll have to get back to you all with real benchmarks it seems.