Pre allocation do not work ...
2 ビュー (過去 30 日間)
古いコメントを表示
Hi everyone,
I am a bit lost,
I ran a code this morning and the pre allocation do not work, it changes the size of my matrices with any reason,
So I did a small try like this :
n=3
x = zeros(1,n);
for ii=1:10
x( ii ) = ( ii );
end
%
Before I got an error because the size just exceeds... now it works, it changes the size of x, how to avoid that ?
Best regards,
1 件のコメント
Adam
2020 年 2 月 3 日
Pre-size it correctly! You pre-size it to 3 then put it in a for loop up to 10 assignging values to it. Where do you expect the next 7 elements to come from? You didn't pre-allocate those so the array grows in the loop.
採用された回答
Raj
2020 年 2 月 3 日
Question is not clear. What exactly are you trying to do here?The pre allocation works fine. You can check it by putting breakpoints in your code. The loop simply overwrites the pre allocated matrix. If you want to stop the loop when the execution reaches the preallocated matrix size then use something like this:
n=3
x = zeros(1,n);
temp=length(x);
for ii=1:10
x( ii ) = ( ii );
if ii==temp
break
end
end
%
4 件のコメント
Rik
2020 年 2 月 3 日
For 1xN it doesn't matter, but as soon as your input is not a vector you will get into trouble. As you can see in the 'answer' below, x is actually an array. If you want to do something with all elements, use numel, otherwise you should use size with the second argument. If the user supplies a column vector and your code expects a row vector, the result could be a memory error. How should your user understand that?
bigN=100000;
A=rand(bigN,1);
B1=rand(bigN,1);
B2=B1';%user didn't read the doc and supplies the wrong shape
no_crash=A.*B1;
crash=A.*B2;
その他の回答 (1 件)
Gauthier Briere
2020 年 2 月 3 日
3 件のコメント
Raj
2020 年 2 月 3 日
編集済み: Raj
2020 年 2 月 3 日
I still don't get it!! Why are you expecting an error? The preallocation simply helps in speeding up the loop. See this answer for details:
Other than that the for loop will always overwrite the preallocated loop. If you MUST HAVE an error you can insert it manually like this:
n=3
x = zeros(1,n);
temp=length(x);
for ii=1:10
if ii>temp
error('Index exceeds preallocated matrix size')
break
end
x( ii ) = ( ii );
end
%
参考
カテゴリ
Help Center および File Exchange で Matrix Indexing についてさらに検索
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!