def absorbing_gaussian(x,sigma):
'''
Applies a gaussian convolution to 2d array `x` with absorbing
boundary conditions.
Parameters
----------
Returns
-------
'''
support = 1+sigma*6
normalization = np.zeros(x.shape,'double')
result = np.zeros(x.shape,'double')
kernel = gaussian_2D_kernel(sigma)
return convolve2d(x, kernel, mode='same', boundary='symm')
评论列表
文章目录