I need to calculate several products of matrices where at least one of them is always hermitian. The function zhemm (in blas, level 3) seems to directly do that in an efficient manner.
However ... how can i access that function and dirctly apply it on numpy arrays?
If you know alternatives that are equivalent or even faster, please let me know.
Any help is highly appreciated. Q