def __call__(self, field):
<EOL>tile = peri.util.Tile(field.shape)<EOL>rx, ry = tile.kvectors(norm=<NUM_LIT:1.0>/tile.shape)<EOL>sx, sy = self.values<EOL>psf = np.exp(-((rx/sx)**<NUM_LIT:2> + (ry/sy)**<NUM_LIT:2>)/<NUM_LIT:2>)<EOL>psf = psf / psf.sum()<EOL>self.psf = psf<EOL>out = fft.fftn(fft.ifftn(field)*fft.ifftn(psf))<EOL>return fftnorm(np.real(out))<EOL>
Accept a field, apply the point-spread-function, and return
the resulting image of a blurred field