Allocating elements in a large matrix takes a VERY long time. Why?
7 ビュー (過去 30 日間)
古いコメントを表示
MATLAB is doing something I have never seen before. Can anyone explain to me the following behaviour:
I have some code:
function myFuncMain()
myMat = NaN(1E5, 5E2);
for t = 1: T
%carry out some tests
%allocate elements to the matrix, eg
myMat(i, j) = 10;
myMat(z,y) = now;
%etc
%Now go into an external function
myMat = myFunc(myMat, variables);
end
end
function myMat = myFunc(myMat, variables)
%carry out some more tests
%allocate some more elements to the matrix, eg
myMat(ii, jj) = 10;
myMat(zz,yy) = now;
%etc
end
When I look at this code in the profiler, I see that >99.99% of the time is spent on the first allocation to the matrix in myFunc.m
Everytime, this function is called it takes a HUGE amount of time, according to the profiler, for the first element to be assigned to the matrix. Then according to the profiler, all the other assignations are instantaneous as expected.
whats going on? Is this something to do with the size of the matrix? Is it be copied somewhere behind the scenes?
I use a 64bit, 16GB machine.
0 件のコメント
採用された回答
Matt J
2013 年 2 月 14 日
編集済み: Matt J
2013 年 2 月 14 日
Because you are changing myMat inside myFunc, MATLAB assumes that a new copy of the entire myMat matrix is required. This happens the instant myFunc makes the first change, resulting in the large memory allocation time that you see.
If MATLAB knew that you were planning to overwrite your original myMat in the calling workspace with the one generated inside myFunc, it would not need to create a new copy, but its parsing mechanism isn't smart enough to know that.
4 件のコメント
Matt J
2013 年 2 月 14 日
編集済み: Matt J
2013 年 2 月 14 日
Just to elaborate a little on my previous Comment, using the Reference class, your code would get modified as follows
function myFuncMain()
obj=Reference;
myMat = NaN(1E5, 5E2);
for t = 1: T
myMat(i, j) = 10;
myMat(z,y) = now;
%etc
%Now go into an external function
obj.X=myMat; clear myMat
myMat = myFunc(obj, variables);
end
end
function myMat = myFunc(obj, variables)
%carry out some more tests
%allocate some more elements to the matrix, eg
myMat=obj.X; obj.X=[];
myMat(ii, jj) = 10;
myMat(zz,yy) = now;
%etc
end
その他の回答 (1 件)
James Tursa
2013 年 2 月 14 日
You might also consider using this fast matrix allocation function UNINIT from the FEX:
The UNINT function allocates arrays of uninitialized values and can be useful in cases where you know the values will be overwritten downstream before their use.
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!