Cholesky decomposition: Difference between revisions

Content deleted Content added

Inline

@@ Line 398: / Line 398: @@
 <lang clojure>(defn cholesky
   [matrix]
-  (let [cols (count matrix)
+  (let [n (count matrix)
         A (to-array-2d matrix)
-        L (make-array Double/TYPE cols cols)]
+        L (make-array Double/TYPE n n)]
-    (doseq [i (range cols) j (range (inc i))]
+    (doseq [i (range n) j (range (inc i))]
       (let [s (reduce + (for [k (range j)] (* (aget L i k) (aget L j k))))]
         (aset L i j (if (= i j)