[SciPy-User] Q: Provide .H for numpy arrays?

Pauli Virtanen pav@iki...
Thu Oct 14 14:50:37 CDT 2010


Thu, 14 Oct 2010 21:18:56 +0200, Nico Schlömer wrote:
>> I think the proposal needs further details. ".T" is NOT completely
>> equivalent because no copying takes place.
> 
> That's right, I was more thinking about semantic, mathematical
> equivalence. Looking at the code, this is certainly a different matter
> about which I, as a Scipy user, don't know much.
> 
> As mentioned in the previous mail, at this moment there appears to be no
> possibility to efficiently block-dot-product two complex-valued
> blockvectors as vdot() doesn't work for blocks, and
> 
>     np.dot( z.T.conjugate(), z )
> 
> unnecessarily creates and intermediate vector. Correct? The suggested
> view option would certainly help out here.

You can probably write a couple of lines function calling ZGEMM from 
scipy.lib.blas. But yes, nothing out of the box.

-- 
Pauli Virtanen



More information about the SciPy-User mailing list