MATLAB 행렬에서 원소를 찾는 효율적인 방법

Question

주어진 코드에서 어떻게 병목 현상을 처리 할 수 ​​있는지 알고 싶습니다.

%% 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,:); 

현재 N ~ 10^6, np ~ 1000 및 number of calls to this part of code = 1000를 들어, 병목 문은 내 코드의 나머지 부분에 비해 시간의 큰 덩어리 인, 총 약 10 초하고있다.

만 인덱싱 문에 대한 샘플 코드를 좀 더 스크린 샷

@EitanT의 요청에 따라하는 것은

Answer 1

:

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) 

반환

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

ans = 

    1 

Answer 2

같은 일을보십시오 :

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

나는 그것이보다 빠르다는 것을 발견했다 :

01 약 4 초에 의해 23,516,

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

는 :

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

나는 당신의 둔화에 대한 이유는 요소 현명한 & 작업입니다 생각합니다. 최악의 경우 가정

Answer 3

:

• 요소가 많다는 & 단락하지 않습니다를 내부적으로
• 비교 식 단일 스레드 당신이 2*1e6*1e3 = 2e9 비교를하고있는

를 ~ 10 초입니다 . 이는 초당 ~ 2 억 회의 비교 (~ 200 MFLOPS)입니다.

일부 1.7 GFLops on a single core을 수행 할 수 있다고 생각하면 실제로는 다소 저조한 것으로 보입니다.

Windows 7을 사용하고 있습니까? 그렇다면 전원 설정을 확인 했습니까? 모바일 프로세서를 사용하고 있으므로 기본적으로 저전력 소모 방식이 효과적 일 것으로 기대됩니다. 이것은 창문이 처리 속도를 축소 할 수있게 해주므로 ... 확인하십시오.

그 외 ... 나는 실마리가 없습니다. 한쪽 평등 당신이 일방적 인 비교로 재구성 할 수 있습니다 중요하지 않습니다 그것은 빠른 진도의 한 순서를 얻는 경우