title: How to compute laplacian of a field?
body: <p>I'm trying to compute the laplacian of a 2d field <strong>A</strong> using <em>scipy.ndimage.convolve</em>.</p>
  
  <pre><code>stencil = numpy.array([[0, 1, 0],[1, -4, 1], [0, 1, 0]])
  scipy.ndimage.convolve(A, stencil, mode='wrap')
  </code></pre>
  
  <p>This doesn't seem to give me the right answer though. Any ideas where I'm going wrong, or are there better ways of computing the laplacian in numpy?</p>