[Numpy-discussion] Efficient Matrix-matrix product of hermitian matrices, zhemm (blas) and numpy

[qubax at gmx.at]

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.