MATLAB Answers

Robin L.
0

Do I need pre-allocation ?

Robin L.
さんによって質問されました 2019 年 3 月 19 日
最新アクティビティ Robin L.
さんによって コメントされました 2019 年 3 月 20 日
Hello everybody !
You helped me a lot these last days, thanks !
I come up today with another question, concerning preallocation.
Let's consider 2 versions of code :
% a and b some matrix 1D
tab = zeros(length(a), length(b)); % PREALLOCATION
for x_idx = 1:length(a)
for y_idx = 1:length(b)
tab(x_idx, y_idx) = a(x_idx) * b(y_idx);
end
end
And :
% x and y some matrix 1D
tab = zeros(length(x), length(y)); % PREALLOCATION
tab = a' .* b;
I am wondering if I really need preallocation in the second version, because I don't access to every key of an "growing up" array, but I use Matlab vectorization..
So I think I can write (without loosing in terms of perf, and event do a bargain) :
% x and y some matrix 1D
tab = a' .* b;
Am I right ?
Robin.

  6 件のコメント

Stephen Cobeldick
2019 年 3 月 19 日
"so the matrix has been created and has been "pre-allocated" with zero values... ?"
Sure, you create an array of zeros... that is not the problem! That line could be array preallocation, if it were followed by some code that accessed that array using some indexing.
But what you actually do on your very next line of code is to simply discard that entire array and replace it entirely with a totally new array. All you are doing is just basic allocatiion of an array to a variable name. And then allocation of another array to that same variable name.
But none of this has anything to do with array preallocation before a loop.
The explanation at the top of that documentation explains what array preallocation is and when it is needed: preallocate before LOOP/S and then inside the loop/s access the array elements using INDEXING. Without the indexing there is no point, you simply discard everything.
Because you simply discard the first array then its class and size are totally irrelevant to the second line. Try it with a cell array, and you will see that MATLAB simply replaces your first (cell) array with the second (numeric) array: no problem, no errors, no preallocation:
>> tab = cell(3,4); % NOT PREALLOCATION
>> tab = [1,2,3] .* [4,5,6]
tab =
4 10 18
Steven Lord
2019 年 3 月 19 日
In general the zeros function is one tool you can use to preallocate.
In your specific example your call to the zeros function is not preallocation.
Preallocating an array means allocating/initializing an array so you can fill it in later on, as opposed to reallocating the array to be a new size every time you add an element to it.
You're allocating an array, then throwing that allocated array away and reusing that name for a whole new array that has no relation with the originally allocated array.
Robin L.
2019 年 3 月 20 日
I thank you Stephen Cobeldick and Steven Lord ! Now I understand what means pre-allocation 😃

サインイン to comment.

製品


リリース

R2018a

2 件の回答

Guillaume
回答者: Guillaume
2019 年 3 月 19 日
編集済み: Guillaume
2019 年 3 月 19 日
 採用された回答

Robin, in the 2nd case, yes you were attempting to preallocate the array. However, the array you preallocated is not the same array that you will end up with. In fact, the array you preallocated is immediately destroyed and replaced by a completely different array on the next line. Hence why Stephen says it's not preallocation.
In fact, more than unnecessary, your preallocation attempt is counter productive. Matlab waste time preallocating an array that is never going to be used and will be destroyed immediately.
There is a big difference between
tab(r, c) = something %indexing
and
tab = something %no indexing
In the first case, using indexing, you're assigning to one or more elements of the array. If the array is not big enough to start with, then matlab waste time making room for the new element(s). Hence preallocation is important.
In the second case, where there's no indexing, you're copying an entire array. Whatever was in the variable before that gets discarded, so preallocation doesn't work.

  1 件のコメント

Robin L.
2019 年 3 月 20 日
Thank you Guillaume ! Now I understand what means pre-allocation 😃

サインイン to comment.


madhan ravi
回答者: madhan ravi
2019 年 3 月 19 日

Preallocation is not necessary in your second version.

  3 件のコメント

Stephen Cobeldick
2019 年 3 月 19 日
The second version has no preallocation anyway. It just allocates an array to a variable. And then allocates another array to that variable. Which has nothing to do with preallocation.
madhan ravi
2019 年 3 月 19 日
I said Preallocation is not needed for the second case , where have I ever said it was preallocation?, one flaw was I never explained the reason for it but Guillaume did so +1.
Robin L.
2019 年 3 月 20 日
Thank you madhan ravi ! Now I understand what means pre-allocation 😃

サインイン to comment.



Translated by