Revert "timing test file [skip ci]"

This reverts commit 46a764c.

Author: lums658

src/include/scoring.h

@@ -27,6 +27,8 @@
  *
  * @section DESCRIPTION
  *
+ * Gemm-based scoring.
+ *
  */
 
 #ifndef TDB_SCORING_H
@@ -73,12 +75,103 @@ class without_ids {};
 // Distance functions
 // ----------------------------------------------------------------------------
 
+/**
+ * @brief Compute sum of squares distance between two vectors.
+ * @tparam V
+ * @tparam U
+ * @param a
+ * @param b
+ * @return
+ */
+#if 0
+template <class V, class U>
+inline auto sum_of_squares(V const& a, U const& b) {
+  float sum{0.0};
+  size_t size_a = size(a);
+
+  if constexpr (std::is_same_v<decltype(a[0]),decltype(b[0])>) {
+    for (size_t i = 0; i < size_a; ++i) {
+      float diff = a[i]- b[i];
+      sum += diff * diff;
+    }
+  } else {
+    for (size_t i = 0; i < size_a; ++i) {
+      float diff = ((float)a[i]) - ((float)b[i]);
+      sum += diff * diff;
+    }
+  }
+  return sum;
+}
+#else
+template <class V, class U>
+inline auto sum_of_squares(V const& a, U const& b) {
+  float sum{0.0};
+  size_t size_a = size(a);
+
+  if constexpr (
+      std::unsigned_integral<std::remove_reference_t<decltype(a[0])>> ||
+      std::unsigned_integral<std::remove_reference_t<decltype(b[0])>>) {
+    for (size_t i = 0; i < size_a; ++i) {
+      // float diff = (float)a[i] - (float)b[i];  // converting to float is slow
+      float diff = (float)a[i] - (float)b[i];
+      sum += diff * diff;
+    }
+  } else {
+    for (size_t i = 0; i < size_a; ++i) {
+      // float diff = (float)a[i] - (float)b[i];  // converting to float is slow
+      float diff = a[i] - b[i];
+      sum += diff * diff;
+    }
+  }
+  return sum;
+}
 
-// sum_of_squares(a, b)
-// sum_of_squares(a)
-// sub_sum_of_squares(a, b, start, stop)
+template <class V>
+inline auto sum_of_squares(V const& a) {
+  float sum{0.0};
+  size_t size_a = size(a);
+
+  if constexpr (std::unsigned_integral<
+                    std::remove_reference_t<decltype(a[0])>>) {
+    for (size_t i = 0; i < size_a; ++i) {
+      // float diff = (float)a[i] - (float)b[i];  // converting to float is slow
+      float diff = (float)a[i];
+      sum += diff * diff;
+    }
+  } else {
+    for (size_t i = 0; i < size_a; ++i) {
+      // float diff = (float)a[i] - (float)b[i];  // converting to float is slow
+      float diff = a[i];
+      sum += diff * diff;
+    }
+  }
+  return sum;
+}
 
+template <class V, class U>
+inline auto sub_sum_of_squares(
+    V const& a, U const& b, size_t start, size_t end) {
+  float sum{0.0};
+
+  if constexpr (
+      std::unsigned_integral<std::remove_reference_t<decltype(a[0])>> ||
+      std::unsigned_integral<std::remove_reference_t<decltype(b[0])>>) {
+    for (size_t i = start; i < end; ++i) {
+      // float diff = (float)a[i] - (float)b[i];  // converting to float is slow
+      float diff = (float)a[i] - (float)b[i];
+      sum += diff * diff;
+    }
+  } else {
+    for (size_t i = start; i < end; ++i) {
+      // float diff = (float)a[i] - (float)b[i];  // converting to float is slow
+      float diff = a[i] - b[i];
+      sum += diff * diff;
+    }
+  }
+  return sum;
+}
 
+#endif
 
 /**
  * @brief Compute L2 distance between two vectors.

src/include/test/CMakeLists.txt

@@ -53,8 +53,6 @@ endmacro (kmeans_add_test)
 # kmeans_add_test(time_l2)
 # kmeans_add_test(time_open)
 
-kmeans_add_test(time_scoring)
-
 kmeans_add_test(unit_vamana_index)
 
 kmeans_add_test(unit_vamana_group)

src/include/test/test_utils.h

@@ -58,11 +58,6 @@ void randomize(R& r, std::tuple<int, int> range = {0, 128}) {
   }
 }
 
-template <std::ranges::range R>
-void randomize(R&& r, std::tuple<int, int> range = {0, 128}) {
-  return randomize(r, range);
-}
-
 // Fill a matrix with sequentially increasing values. Will delete data from the
 // URI if it exists.
 template <class Matrix>