Accumarray application for rectangle submatices
조회 수: 1 (최근 30일)
이전 댓글 표시
Is there an elegant way to exclude the for loop from the following function using the accumarray function, for example?
function K = accum(K_size, ...
Pi0, Pi1, Pj0, Pj1, ...
Ki0, Ki1, Kj0, Kj1, ...
psf_arr, psf_ind)
K = zeros(K_size);
for s = 1 : numel(psf_ind)
K( Ki0(s) : Ki1(s), Kj0(s) : Kj1(s) ) = ...
K( Ki0(s) : Ki1(s), Kj0(s) : Kj1(s) ) + ...
psf_arr( Pi0(s) : Pi1(s), Pj0(s) : Pj1(s), psf_ind(s) );
end
end
댓글 수: 6
Jan
2021년 8월 6일
편집: Jan
2021년 8월 6일
@Grigorii Nefedov: Please provide some input data with usual sizes. It matters, if this hould be optimized for millions of loops over small blocks or some loops over huge blocks.
It is tedious to guess, how usual inputs looks in oyur code, but you should be able to create some "meaningful" data with some calls of rand().
Do you want to get an elegenat code, or a fast one? I do not see how avoiding a loop is useful in general.
Matt J
2021년 8월 6일
If all the rectangles were the same size, I can see how you might make it faster. Otherwise, I think the for-loop is your best bet.
답변 (0개)
참고 항목
카테고리
Help Center 및 File Exchange에서 Loops and Conditional Statements에 대해 자세히 알아보기
제품
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!