effiziente Methode zwei nicht schwach besetzte Matrizen des Multiplizierens

Question

Ich habe zwei Matrizen vom Typ float:
A der Abmessungen 7000x100000 und
B der Abmessungen 100000x20.
Wenn ich sie multipliziere, verbraucht mein Code meinen gesamten RAM, obwohl die Ausgabe klein ist.
Gibt es eine Möglichkeit, dies effizienter zu machen?
Ich habe versucht, this from Matlab help page zu tun, aber es hat nicht geholfen.

Answer 1

Sie können versuchen, und teilen sich in Blöcke:

durch Reihen von A:

row_blk = 1000; 
C = zeros(size(A,1), size(B,2), class(A)); 
f = 1; t = row_blk; 
while t <= size(A,1) 
    C(f:t,:) = A(f,t,:)*B; 
    f = t+1; 
    t = min(size(A,1), f+row_blk-1); 
end 

von Spalten von A:

col_blk = 10000; 
C = zeros(size(A,1), size(B,2), class(A)); 
f = 1; t = col_blk; 
while t <= size(A,2) 
    C = C + A(:,f:t)*B(f:t,:); 
    f = t+1; 
    t = min(size(A,2), f+col_blk-1); 
end 

Answer 2

ich diesen Code in GNU laufen Oktave und Speicherverbrauch dieses Codes:

a= rand(7000,100000); 
b = rand(100000,20); 

ist:

Absolute running time: 10.05 sec, cpu time: 9.97 sec, memory peak: 5375 Mb 

und Speicherverbrauch dieser Code:

a= rand(7000,100000); 
b = rand(100000,20); 
c = a * b; 

ist:

Absolute running time: 14.26 sec, cpu time: 14.19 sec, memory peak: 5376 Mb 

so kein signifikanter Unterschied beobachtet!