Come calcolare solo la diagonale di un prodotto matrice in Octave?

Question

Esiste un modo in Octave per calcolare e memorizzare solo la diagonale di un prodotto matrice?

Fondamentalmente come fare: vector = diag(A*B);

Non mi importa niente dei valori di A*B ad eccezione di quelli sulla diagonale. Le dimensioni della matrice sono intorno a 80k x 12 e 12 x 80k, quindi anche se non mi interessa la velocità/memoria extra semplicemente non si adatta alla RAM.

Strano, dal momento che Octave è un pacchetto per enormi set di dati e diagonali sono molto importanti, quindi dovrebbe essere possibile.

Answer 1

Il primo elemento nella diagonale è il prodotto scalare della prima riga di A con la prima colonna di B. Il secondo elemento nella diagonale è il prodotto scalare della seconda riga di A con la seconda colonna di B.

in altre parole:

vector = sum(A.*B',2); 

Answer 2

questo è come si potrebbe fare in MATLAB (probabilmente simile a sintassi Octave):

vector = sum(A.*B',2); 

consente di calcolare solo la diagonale risultante dell'operazione A*B come vettore colonna vector.

Answer 3

realtà credo sia il prodotto scalare della prima riga di A con il primo colonna di B ... secondo elemento diagonale è il prodotto scalare della seconda riga e la seconda colonna ... ecc