关于高斯核函数（rbf）径向基函数的代码

在scdn上找了好久，想找一个比较的好用的python代码就是没有。索性自己写一个吧；我把高斯核函数的公式放在这里
这个公式有多简单，就三个输入，x1,x2,sigma。为什么给我搞这么复杂？为什么。就不能简单一点，各位看看这个好吗：

def kernelfun(xtrain,xt,gamma):
    xtrain = xtrain.T
    xt = xt.T
    nb_data = xtrain.shape[0]
    omega = -2*xtrain*xt.T
    xtrain = np.dot(np.sum(xtrain**2,axis=1),np.ones((1,xt.shape[0])))
    xt = np.dot(np.sum(xt**2,axis=1),np.ones((1,nb_data)))
    omega = omega+xtrain+xt.T
    omega = np.exp(-omega/(2*gamma**2))
    return omega

用法也很简单，怎么用看这里

k = kernelfun(a,b,2)