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 -- The king who needs to remind his people of his rank, is no king. A beggar's mistake harms no one but the beggar. A king's mistake, however, harms everyone but the king. Too often, the measure of power lies not in the number who obey your will, but in the number who suffer your stupidity.