Vectorising cell array operations

1 visualización (últimos 30 días)
Simon
Simon el 25 de Jul. de 2013
I have currently a series of polygons in my code, stored as m*2 arrays, where each 1*2 section is a point. these polygons are stored in a cell array called polys.
I need to flip all of these vertically, so the code i am currently using is
for i = 1:numel(polys)
polys{i}(:, 1) = m-polys{i}(:, 1)
end
Where m is the height of the region (and the y coords are stored in the first position (ie a pioint is (y, x), not (x, y))
I want to execute this whole loop in one line (probably using cellfun?), but im not sure how to do it??
  2 comentarios
Daniel Shub
Daniel Shub el 25 de Jul. de 2013
Why do you want to do that? The for loop seems nice and clean and self documenting to me.
Simon
Simon el 26 de Jul. de 2013
Various reasons that aren't worth going into. One if which being I'm doing everything I possibly can to remove all loops, for speed, but there are others

Iniciar sesión para comentar.

Respuestas (1)

Daniel Shub
Daniel Shub el 26 de Jul. de 2013
Editada: Daniel Shub el 26 de Jul. de 2013
Blindly removing loops for speed is an antiquated view of MATLAB based on the performance of versions over 10 years old prior to the introduction of the JIT accelerator. To fully optimize your code you will want to make use of the tools like the profile (but be aware that it disables the JIT accelerator).
Given all you want to do is remove the loop you could look into either doing the task recursively or using CELLFUN. Neither are likely to give you a substantial speed up (and are likely to actually slow things down) and will be less readable than your simple loop.
For example, with cellfun
polys{1} = randn(5);
polys{2} = randn(6);
polys{3} = randn(7);
m = 10;
polys = cellfun(@(x)[m-x(:, 1), x(:, 2:end)], polys, 'UniformOutput', false);
  2 comentarios
Simon
Simon el 26 de Jul. de 2013
I did mention that I wanted to use cellfun, I can't get it to work... I had an anonymous function declared but it wasn't working at all
Daniel Shub
Daniel Shub el 26 de Jul. de 2013
See my edit, but again I see no advantage to doing this. It is unlikely to be faster.

Iniciar sesión para comentar.

Categorías

Más información sobre Loops and Conditional Statements en Help Center y File Exchange.

Community Treasure Hunt

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

Start Hunting!

Translated by