Metodo efficiente per la ricerca di elementi nella matrice MATLAB

Question

Vorrei sapere come si può trattare il collo di bottiglia nel pezzo di codice indicato.

%% Points is an Nx3 matrix having the coordinates of N points where N ~ 10^6 
Z = points(:,3) 
listZ = (Z >= a & Z < b); % Bottleneck 
np = sum(listZ); % For later usage 
slice = points(listZ,:); 

Attualmente per N ~ 10^6, np ~ 1000 e number of calls to this part of code = 1000, la dichiarazione collo di bottiglia sta prendendo circa 10 secondi in totale, che è un grande pezzo di tempo rispetto al resto del mio codice.

Alcuni ulteriori screenshot per un codice di esempio per solo la dichiarazione di indicizzazione come richiesto dal @EitanT

Answer 1

Se l'uguaglianza da un lato non è importante che si può riformulare in un confronto unilaterale e diventa un ordine di grandezza più veloce:

Z = rand(1e6,3); 
a=0.5; b=0.6; 
c=(a+b)/2; 
d=abs(a-b)/2; 
tic 
for k=1:100, 
    listZ1 = (Z >= a & Z < b); % Bottleneck 
end 
toc 

tic 
for k=1:100, 
    listZ2 = (abs(Z-c)<d); 
end 
toc 

isequal(listZ1, listZ2) 

rendimenti

Elapsed time is 5.567460 seconds. 
Elapsed time is 0.625646 seconds. 

ans = 

    1 

Answer 2

provare a fare qualcosa di simile:

for i = 1:1000 
    x = (a >= 0.5); 
    x = (x < 0.6); 
end 

Ho trovato che fosse più veloce di:

for i = 1:1000 
    x = (a >= 0.5 & a < 0.6); 
end 

da circa 4 secondi:

Elapsed time is 0.985001 seconds. (first one) 
Elapsed time is 4.888243 seconds. (second one) 

Penso che il motivo per il rallentamento è l'elemento saggio & operazione.

Answer 3

Supponendo che il caso peggiore:

• elemento-saggio & non è in corto circuito internamente
• i confronti sono single-threaded

Si sta facendo 2*1e6*1e3 = 2e9 confronti in ~ 10 secondi . Si tratta di ~ 200 milioni di confronti al secondo (~ 200 MFLOPS).

Considerando che è possibile fare un po 'di 1.7 GFLops on a single core, questo sembra piuttosto basso.

Stai utilizzando Windows 7? In tal caso, hai controllato le impostazioni di alimentazione? Sei su un processore mobile, quindi mi aspetto che per impostazione predefinita, ci sarà uno schema di consumo a basso consumo. Ciò consente alle finestre di ridimensionare la velocità di elaborazione, quindi ... verificarlo.

A parte questo .... non ne ho idea.