Python/Numpy - Kreuzprodukt von übereinstimmenden Zeilen in zwei Arrays

Question

Was ist der beste Weg, das Kreuzprodukt jeder entsprechenden Zeile zwischen zwei Arrays zu nehmen? Zum Beispiel:

a = 20x3 array 
b = 20x3 array 
c = 20x3 array = some_cross_function(a, b) where: 
c[0] = np.cross(a[0], b[0]) 
c[1] = np.cross(a[1], b[1]) 
c[2] = np.cross(a[2], b[2]) 
...etc... 

Ich weiß, dass dies mit einer einfachen Python-Schleife durchgeführt werden kann oder numpy des apply_along_axis verwenden, aber ich frage mich, ob es eine gute Möglichkeit ist, dies innerhalb des zugrunde liegenden C-Code von numpy völlig zu tun. Ich benutze derzeit eine einfache Schleife, aber das ist bei weitem der langsamste Teil meines Codes (meine tatsächlichen Arrays sind Zehntausende von Zeilen lang).

Answer 1

Ich werde wahrscheinlich diese Antwort in ein paar Minuten löschen müssen, wenn ich meinen Fehler erkenne, aber funktioniert das Offensichtliche nicht?

>>> a = np.random.random((20,3)) 
>>> b = np.random.random((20,3)) 
>>> c = np.cross(a,b) 
>>> c[0], np.cross(a[0], b[0]) 
(array([-0.02469147, 0.52341148, -0.65514102]), array([-0.02469147, 0.52341148, -0.65514102])) 
>>> c[1], np.cross(a[1], b[1]) 
(array([-0.0733347 , -0.32691093, 0.40987079]), array([-0.0733347 , -0.32691093, 0.40987079])) 
>>> all((c[i] == np.cross(a[i], b[i])).all() for i in range(len(c))) 
True