Matrix-Zeilen/-Spalten effizient in MATLAB löschen

Question

Tipps zum effizienten und schnellen Löschen einer gegebenen Zeile/Spalte aus einer Matrix?

Ich hatte anfangs geglaubt, dass das Löschen der letzten Spalte einer gegebenen Matrix effizienter wäre als die erste Spalte, und alle Spaltenoperationen effizienter wären als Zeilenoperationen (MATLABs spaltenbasierter Speicher), was ich konnte Bestätigung durch Tests Allerdings war die Leistung, die ich bekam, ziemlich unglücklich.

someB = rand(4,50000); 
someC = someB.'; 

tic 
while size(someB,2) > 2 
    someB(:,size(someB,2)) = []; 
end 
toc 

tic 
while size(someC,1) > 2 
    someC(size(someC,1),:) = []; 
end 
toc 

%Elapsed time is 13.869280 seconds. 
%Elapsed time is 10.198270 seconds. 

Ich habe eine schnelle Suche und in this MATLAB newsgroup discussion fand ich die Hoffnung, dass durch externe C-MEX-Funktionen gibt es in der Tat eine Art und Weise effizient die letzte Spalte einer Matrix schnell zu löschen sein. Der Code ist unten angehängt.

#include "mex.h" 

// You may need to uncomment the next line 
//#define mwSize int 

void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]) 
{ 
    mwSize n; 
    if(n = mxGetN(prhs[0])) 
     mxSetN(prhs[0], n - 1); 
} 

Allerdings war ich nicht in der Lage an den Code zu bekommen, mich laufen. Wenn Sie sich die Ergebnisse, die der Autor gefunden hat, kurz ansehen, werden Sie bemerkenswerte Leistungen feststellen. Ich bin selbst nicht so gut in MEX; Würde jemand wissen, wie man obigen Code repariert, so dass er läuft, oder alternativ einen gleichwertigen/fast gleichwertigen MEX-Code/MATLAB-Code leistungsmäßig haben?

Danke!

Answer 1

Nun, zum einen, die MEX-Lösung, die Sie gebucht haben, ist eigentlich nicht löschen die Spalte, so ist es nicht wirklich ein fairer Vergleich. Die Lösung, "mxSetN" zu verwenden, setzt nur den internen Header für das mxArray, um zu denken, dass es N Spalten hat. Aus der Dokumentation:

Normalerweise verwenden Sie mxSetN, um die Form eines vorhandenen mxArray zu ändern. Die Funktion "mxSetN" weist keinen Speicherplatz für die Arrays "pr", "pi", "ir" oder "jc" zu.

In Bezug auf die Timing-Ergebnisse, die Sie sehen, was passiert, wenn Sie diese mehrmals ausführen? Als ich das erste Mal lief, hatte ich ähnliche Ergebnisse wie bei Ihnen, aber beim zweiten Durchlauf waren die Timings für beide Methoden ähnlich.

Beachten Sie, dass Sie in jedem Fall eine große Speichermenge neu zuordnen müssen, da Sie nur 4 Elemente gleichzeitig löschen. Dies bedeutet, dass die früheren Schleifeniterationen mehr kosten als die späteren Iterationen.

Die große Frage hier ist, was Sie letztlich versuchen zu tun. Vielleicht können Sie die Löschung vermeiden und nur den Teil der Matrix verwenden, den Sie brauchen, oder vielleicht gibt es eine Möglichkeit, zu vermeiden, dass Sie bei jeder Iteration eine einzelne Spalte löschen müssen (und stattdessen mehrere Spalten verwenden).