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 comentarios

Stephen23
Stephen23 el 8 de Oct. de 2018
Editada: Stephen23 el 8 de Oct. de 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?
Iori Yagami Kusanagi
Iori Yagami Kusanagi el 8 de Oct. de 2018
When I change the value of c, implementing linear search, takes a lot of time this procedure.
Bruno Luong
Bruno Luong el 8 de Oct. de 2018
Editada: Bruno Luong el 8 de Oct. de 2018
Are you awared that you creates arrays of few hundred of Tb?
Iori Yagami Kusanagi
Iori Yagami Kusanagi el 8 de Oct. de 2018
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 el 9 de Oct. de 2018
Editada: Bruno Luong el 9 de Oct. de 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?

Iniciar sesión para comentar.

Respuestas (1)

Jan
Jan el 9 de Oct. de 2018

0 votos

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.

Categorías

Más información sobre Creating and Concatenating Matrices en Centro de ayuda y File Exchange.

Productos

Etiquetas

Preguntada:

el 8 de Oct. de 2018

Editada:

el 9 de Oct. de 2018

Community Treasure Hunt

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

Start Hunting!

Translated by