<br><div class="gmail_quote">On Wed, Mar 25, 2009 at 3:56 PM, David Baddeley <span dir="ltr">&lt;<a href="mailto:david_baddeley@yahoo.com.au">david_baddeley@yahoo.com.au</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">

<br>
There&#39;s always:<br>
<br>
X, Y = numpy.mgrid[-128:128, -128:128]<br>
R2 = X**2 + Y**2<br>
p = numpy.exp(-0.01*R2)<br>
<br>
which uses a bit more memory, but avoids loops.<br>
<br>
David<br>
</blockquote><div><br>If you use ogrid instead of mgrid, you&#39;ll avoid the increase in memory, since the arrays are created in a form ready for broadcasting.<br><br>Ryan <br></div></div><br>-- <br>Ryan May<br>Graduate Research Assistant<br>

School of Meteorology<br>University of Oklahoma<br>