   "Dist matrix"
   ([m] (dist-matrix m m))
   ([m1 m2]
-   (matrix
+   (matrix/matrix
      (map
        (fn [r1]
          (map

        ; next make sure we're not dealing with any clusters that are identical to eachother
        uniq-clusters (uniqify-clusters clusters)
        ; count uniq data points to figure out how many clusters are possible
        possible-clusters (min k (count (distinct (rows (nm/get-matrix data)))))]
    (loop [clusters uniq-clusters]
      ; Whatever the case here, we want to do one more recentering
      (let [clusters (recenter-clusters data clusters :weights weights)]



(defn xy-clusters-to-nmat2 [clusters]
  (nm/named-matrix
    (map :id clusters) ; row names
    [:x :y] ; column names
    (matrix/matrix (map :center clusters))))


