Hi, <br>This is a really quick and stupid question. Let&#39;s say I have a function and I want to calculate its value for a set of points. What&#39;s the best way of quickly doing this avoiding loops?<br><br>def MyFunction ( p1, p2, p3, t, y)<br>
  return *some_value*<br><br>t = numpy.array ( ..... )<br>y = numpy.array ( ....)<br>a1,a2,a3 = numpy.mgrid [ 0:100, 0:1:100j, 0:100 ]<br><br>evaluated_function = MyFunction (a1, a2, a3, t, y) <br><br>evaluated_function.shape would be (100,100,100), and evaluated_function [m, n, k] should be equal to<br>
MyFunction ( a1[m], a2[n], a3[k], t, y )<br><br>I think mgrid (ogrid) are the right tools, but I don&#39;t really understand how I can get it to do what I want. I have read &lt;<a href="http://docs.scipy.org/doc/scipy/reference/tutorial/basic.html#id5">http://docs.scipy.org/doc/scipy/reference/tutorial/basic.html#id5</a>&gt;, but that doesn&#39;t clear things up. I guess that what I want to do is just mentioned at the end of this Section of the tutorial, but can anyone maybe give an example?<br>
<br>many thanks!<br>J<br>