code.clj

(defn gaussian [m mu gamma]
  (let [m (mp/clone m)
        t (/ 1 (* gamma (Math/sqrt (* 2 Math/PI))))]
    (magic/specialize :double
      (loop-over [m]
        (let [x (aget m-data m-idx)
              exparg (* -0.5 (Math/pow (/ (- x mu) gamma) 2))
              g (* t (Math/exp exparg))]
          (aset m-data m-idx g))))))

code.py

def gen_vector(n):
    return np.array([float(x) / n for x in range(n)])

def gaussian1(m, mu, gamma):
    t = 1. / (gamma * (2 * np.pi) ** 0.5)
    exparg = -0.5 * (((m - mu) / gamma) ** 2)
    return np.exp(exparg) * t