[SciPy-User] Q: Provide .H for numpy arrays?
Pauli Virtanen
pav at iki.fi
Thu Oct 14 15:50:37 EDT 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