Logo Questions Linux Laravel Mysql Ubuntu Git Menu
 

Loopless function calls on vector/matrix members in Matlab/Octave

I came into matrix world from loop world (C, etc)

I would like to call a function on each individual member of a vector/matrix, and return the resulting vector/matrix.

This is how I currently do it:

function retval = gauss(v, a, b, c)
  for i = 1:length(v)
    retval(i) = a*(e^(-(v(i)-b)*(v(i)-b)/(2*c*c)));
  endfor
endfunction

Example usage:

octave:47> d=[1:1000];
octave:48> mycurve=gauss(d, 1, 500, 100);

Now, all advice on MATLAB/Octave says: STOP whenever you catch yourself using loops and think of a better way of doing it.

Thus, my question: Can one call a function on each member of a vector/matrix and return the result in a new vector/matrix all at once without using explicit loops?

That is I am looking for something like this:

 function retval = newfun(v)
    retval = 42*(v^23); 
endfunction

Perhaps, it is just syntactic sugar, that is all, but still would be useful to know.

like image 679
Sint Avatar asked Mar 18 '10 15:03

Sint


2 Answers

The function should look like this:

function retval = gauss(v, a, b, c)
  retval = a*exp(-(v-b).^2/(2*c^2));

I would recommend you to read MATLAB documentation on how to vectorize the code and avoid loops:

Code Vectorization Guide

Techniques for Improving Performance

Also remember that sometime code with loops can be more clear that vectorized one, and with recent introduction of JIT compiler MATLAB deals with loops pretty well.

like image 189
yuk Avatar answered Sep 28 '22 04:09

yuk


In matlab the '.' prefix on operators is element-wise operation.

Try something like this:

function r = newfun(v)
 r = a.*exp(-(v-b).^2./(2*c^2))
end
like image 35
dkantowitz Avatar answered Sep 28 '22 03:09

dkantowitz