Gibt es eine Möglichkeit, nur den eigentlichen Teil eines Numpy -Matmuls zu berechnen? ⇐ Python

[Guest]

Nehmen wir an, ich habe zwei Arrays A und b beide mit dType np.complex128 und ich möchte c = np.matmul (a, b) .real . < /p>
Das heißt, ich kümmere mich nicht um den imaginären Teil, nur den realen Teil. Gibt es einen besseren Weg, dies zu tun? Es scheint, als würde eine solche Methode eine anständige Beschleunigung liefern, wenn man bedenkt, dass tatsächlich nur etwa die Hälfte der schwimmenden Punktmultiplikationen erforderlich ist. Realer Teil des Einsum-NP-INeinsum ("AB, BC, CD, DA-> A", ...). Real mit ziemlich großen Operanden, also hoffentlich sind Methoden für die grundlegende Frage verallgemeinerbar auf Einsum.