Conjugate transpose l'opérateur ".H" en numpy

Il est très pratique dans numpy pour utiliser le .T attribut pour obtenir une version transposée d'une ndarray. Cependant, il n'existe aucune autre façon d'obtenir le conjugué transposer. Numpy est la matrice de la classe a l' .H opérateur, mais pas ndarray. Parce que j'aime un code lisible, et parce que je suis trop paresseux pour écrire toujours .conj().Tje voudrais le .H propriété d'être toujours disponible pour moi. Comment puis-je ajouter cette fonctionnalité? Est-il possible de l'ajouter, de sorte qu'il est brainlessly disponible à chaque fois que numpy est importé?

(Même question pourrait être posée à propos de la .I inverse de l'opérateur.)

source d'informationauteur benpro

  1. 8

    Vous pouvez sous-classe de la ndarray objet comme:

    from numpy import ndarray

class myarray(ndarray):    
    @property
    def H(self):
        return self.conj().T

    tels que:

    a = np.random.random((3, 3)).view(myarray)
a.H

    vous donnera le comportement désiré...

  2. 7

    En général, la difficulté de ce problème est que Numpy est une C-extension, qui ne peut pas être singe patché...ou peut-il? Le forbiddenfruit module permet de le faire, bien qu'il se sent un peu comme jouer avec des couteaux.

    Voici donc ce que j'ai fait:

    1. Installer le très simple forbiddenfruit paquet

    2. Déterminer l'utilisateur le répertoire de personnalisation:

      import site
print site.getusersitepackages()

    3. Dans ce répertoire, modifier usercustomize.py à inclure les éléments suivants:

      from forbiddenfruit import curse
from numpy import ndarray
from numpy.linalg import inv
curse(ndarray,'H',property(fget=lambda A: A.conj().T))
curse(ndarray,'I',property(fget=lambda A: inv(A)))

    4. Tester:

      python -c python -c "import numpy as np; A = np.array([[1,1j]]);  print A; print A.H"

      Résultats dans:

      [[ 1.+0.j  0.+1.j]]
[[ 1.-0.j]
 [ 0.-1.j]]