How vectorize this operation

Can someone help me please?
x = 10*ones((10^7),1);
L = ((length(x)*2):-2:2)';
M = x.*L;
p= -M./L;
c= 0.25;
x = x+c*p;
Is possible vectorize x = x+c*p;

5 件のコメント

Stephen23
Stephen23 2018 年 10 月 8 日
編集済み: Stephen23 2018 年 10 月 8 日
What does "Is possible vectorize x = x+c*p" mean? That line of code already operates on vectors (of appropriate sizes) and does not throw an error. What are you trying to achieve?
Iori Yagami Kusanagi
Iori Yagami Kusanagi 2018 年 10 月 8 日
When I change the value of c, implementing linear search, takes a lot of time this procedure.
Bruno Luong
Bruno Luong 2018 年 10 月 8 日
編集済み: Bruno Luong 2018 年 10 月 8 日
Are you awared that you creates arrays of few hundred of Tb?
Iori Yagami Kusanagi
Iori Yagami Kusanagi 2018 年 10 月 8 日
Yes, is a optimization problem that I need to solve, and when run on that part take a lot of time
Bruno Luong
Bruno Luong 2018 年 10 月 9 日
編集済み: Bruno Luong 2018 年 10 月 9 日
You know how much memory has your PC? Usually few or few dozen Gb. How you intend to carry out the storage of such matrix on the computer?

サインインしてコメントする。

回答 (1 件)

Jan
Jan 2018 年 10 月 9 日

0 投票

The code works vectorized already. So the problem is not vectorizing.
The code can be simplified slightly:
M = x .* L;
p = -x; % instead of: -M ./ L, but p is not needed at all
c = 0.25;
x = (1-c) * x; % instead of: x + c * p;
The code is not meaningful yet. Maybe there is another way to improve it. But again, vectorization is applied here already.

カテゴリ

ヘルプ センター および File ExchangeCreating and Concatenating Matrices についてさらに検索

製品

タグ

質問済み:

2018 年 10 月 8 日

編集済み:

2018 年 10 月 9 日

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by