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 Comments

Stephen23
Stephen23 on 8 Oct 2018
Edited: Stephen23 on 8 Oct 2018
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 on 8 Oct 2018
Edited: Bruno Luong on 8 Oct 2018
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 on 9 Oct 2018
Edited: Bruno Luong on 9 Oct 2018
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?

Sign in to comment.

Answers (1)

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.

Categories

Products

Edited:

on 9 Oct 2018

Community Treasure Hunt

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

Start Hunting!