def METHOD_NAME(kernel_size=21, sigma=3):
    """
    Returns a 2D Gaussian kernel
    :param kernel_size: should not be even
    :param sigma: sigma of the gaussian
    :return: array containing the gaussian values
    """
    X = np.arange(kernel_size)
    return Gauss2D(X, X, 1.0, kernel_size//2, kernel_size//2, sigma)