performance - Clojure 性能，如何将提示键入 r/map

Question

下面，我有 2 个函数计算它们参数的平方和。第一个很好且实用，但比第二个慢 20 倍。我认为 r/map 没有利用 aget 从双数组中检索元素，而我在函数 2 中明确地这样做了。

有什么方法可以进一步输入提示或帮助 r/map r/fold 更快地执行？

(defn sum-of-squares
  "Given a vector v, compute the sum of the squares of elements."
  ^double [^doubles v]
  (r/fold + (r/map #(* % %) v)))

(defn sum-of-squares2
  "This is much faster than above.  Post to stack-overflow to see."
  ^double [^doubles v]
  (loop [val 0.0
         i (dec (alength v))]
    (if (neg? i)
      val
      (let [x (aget v i)]
        (recur (+ val (* x x)) (dec i))))))

(def a (double-array (range 10)))
(quick-bench (sum-of-squares a))

800 纳秒

(quick-bench (sum-of-squares2 a))

40 纳秒

Answer 1

在实验之前，我在 project.clj 中添加了下一行:

:jvm-opts ^:replace [] ; Makes measurements more accurate

基本测量:

(def a (double-array (range 1000000))) ; 10 is too small for performance measurements
(quick-bench (sum-of-squares a)) ; ... Execution time mean : 27.617748 ms ...
(quick-bench (sum-of-squares2 a)) ; ... Execution time mean : 1.259175 ms ...

这或多或少与问题中的时差一致。让我们尽量不使用 Java 数组(这对于 Clojure 来说并不是真正惯用的):

(def b (mapv (partial * 1.0) (range 1000000))) ; Persistent vector
(quick-bench (sum-of-squares b)) ; ... Execution time mean : 14.808644 ms ...

快了将近 2 倍。现在让我们删除类型提示:

(defn sum-of-squares3
"Given a vector v, compute the sum of the squares of elements."
[v]
(r/fold + (r/map #(* % %) v)))

(quick-bench (sum-of-squares3 a)) ; Execution time mean : 30.392206 ms
(quick-bench (sum-of-squares3 b)) ; Execution time mean : 15.583379 ms

与带有类型提示的版本相比，执行时间仅略微增加。顺便说一下，版本是 transducers具有非常相似的性能并且更干净:

(defn sum-of-squares3 [v]
  (transduce (map #(* % %)) + v))

现在关于附加类型提示。我们确实可以先优化 sum-of-squares执行:

(defn square ^double [^double x] (* x x))

(defn sum-of-squares4
  "Given a vector v, compute the sum of the squares of elements."
  [v]
  (r/fold + (r/map square v)))

(quick-bench (sum-of-squares4 b)) ; ... Execution time mean : 12.891831 ms ...

(defn pl
  (^double [] 0.0)
  (^double [^double x] (+ x))
  (^double [^double x ^double y] (+ x y)))

(defn sum-of-squares5
  "Given a vector v, compute the sum of the squares of elements."
  [v]
  (r/fold pl (r/map square v)))

(quick-bench (sum-of-squares5 b)) ; ... Execution time mean : 9.441748 ms ...

注意#1:关于 sum-of-squares4 的参数和返回值的类型提示和 sum-of-squares5没有额外的性能优势。

注释 #2 :从 optimizations 开始通常是不好的做法.直白版 (apply + (map square v))在大多数情况下都会有足够好的性能。 sum-of-squares2与惯用语相去甚远，并且几乎不使用任何 Clojure 概念。如果这真的是性能关键代码 - 最好在 Java 中实现它并使用互操作。尽管有 2 种语言，但代码会更清晰。或者甚至在非托管代码(C、C++)中实现它并使用 JNI(不是真正可维护的，但如果正确实现，可以提供最佳性能)。