How vectorize this operation

2018 10 8
1 个回答
更新时间：2018 10 9
1 次查看（30 天）

0 个投票

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 个评论
显示 3更早的评论 隐藏 3更早的评论

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?
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?
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 Exchange 中查找有关 Creating and Concatenating Matrices 的更多信息

产品

标签

Community Treasure Hunt

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

Start Hunting!

Translated by