Which way is better?
Show older comments
Note: This is not about preallocating variable in a loop. This is for general case.
For instance:
a(10,10) = 10; % the rest are filled with zeros
%or
a = zeros(10,10);
a(10,10) = 10;
4 Comments
Walter Roberson
on 15 Aug 2019
There was a Question in which people did some timing. The answer is not always obvious.
The larger question at the time also asked about the timing of assigning a constant to the matrix. For example is
5*ones(10,10)
faster than
5+zeros(10, 10l
Or
repmat(5, 10, 10)
madhan ravi
on 15 Aug 2019
Walter Roberson
on 15 Aug 2019
madhan ravi
on 15 Aug 2019
Answers (1)
In terms of speed, the 2nd option is 3.0 times faster on average (p<0.0001, Wilcoxon Signed Rank; 1 million individually timed iterations of the following 2 options).
a = []; % to reset a upon each loop
a(10,10) = 10;
% versus
a = []; % to reset a upon each loop
a = zeros(10,10);
a(10,10) = 10;
Declaring 'a' as zeros and then filling in the index (10,10) could also avoid other errors that option 1 may impose if 'a' is already a variable or function that exists.
Lastly, I think the 2nd option is more readable. Upon first glace of the first option, I'm wondering where 'a' came from.
3 Comments
madhan ravi
on 15 Aug 2019
Rik
on 15 Aug 2019
There is nothing inherently wrong with the first option, but because it can cause confusion, you should add a comment there. It is also important to make sure that a doesn't exist before this executes.
madhan ravi
on 15 Aug 2019
Categories
Find more on Matrix Indexing in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!