From 349cfd5090f586b7ec189b72c00786522199fe34 Mon Sep 17 00:00:00 2001
From: Douglas Rumbaugh <dbr4@psu.edu>
Date: Wed, 1 May 2024 16:03:19 -0400
Subject: TS BSM Adjustments

---
 benchmarks/include/standard_benchmarks.h | 39 +++++++++++++++++++++++++++++++-
 benchmarks/include/triespline_bsm.h      | 31 ++++++++++++++++++++-----
 2 files changed, 63 insertions(+), 7 deletions(-)

(limited to 'benchmarks/include')

diff --git a/benchmarks/include/standard_benchmarks.h b/benchmarks/include/standard_benchmarks.h
index 76b5f7c..1261a4c 100644
--- a/benchmarks/include/standard_benchmarks.h
+++ b/benchmarks/include/standard_benchmarks.h
@@ -22,6 +22,19 @@
 static size_t g_deleted_records = 0;
 static double delete_proportion = 0.05;
 
+static volatile size_t total = 0;
+
+template<typename DE, typename QP, typename R>
+static void run_queries(DE *extension, DE *ghost, std::vector<QP> &queries) {
+    for (size_t i=0; i<queries.size(); i++) {
+        std::vector<R> res = extension->query(&queries[i]);
+        std::vector<R> negres = ghost->query(&queries[i]);
+        auto result = res[0].first - negres[0].first;
+        total = result;
+    }
+}
+
+
 template<typename DE, typename QP, bool BSM=false>
 static void run_queries(DE *extension, std::vector<QP> &queries) {
     for (size_t i=0; i<queries.size(); i++) {
@@ -63,6 +76,7 @@ static void run_queries(DE *extension, std::vector<QP> &queries) {
                 fflush(stdout);
                 #endif
             } else {
+                total = res.size();
                 #ifdef BENCH_PRINT_RESULTS
                     fprintf(stdout, "\n\n");
                     for (int i=res.size()-1; i>=0; i--) {
@@ -123,7 +137,6 @@ static void insert_records(psudb::bsm::BentleySaxe<R, DS, MDSP> *extension,
                            size_t start, size_t stop, std::vector<R> &records) {
 
     psudb::progress_update(0, "Insert Progress");
-    size_t reccnt = 0;
     for (size_t i=start; i<stop; i++) {
         extension->insert(records[i]);
     }
@@ -132,6 +145,30 @@ static void insert_records(psudb::bsm::BentleySaxe<R, DS, MDSP> *extension,
 }
 
 
+template<typename DS, typename R, bool MDSP=false>
+static void insert_records(psudb::bsm::BentleySaxe<R, DS, MDSP> *extension, 
+                           psudb::bsm::BentleySaxe<R, DS, MDSP> *ghost, 
+                           size_t start, size_t stop, std::vector<R> &records,
+                           std::vector<size_t> &to_delete, size_t &delete_idx, 
+                           gsl_rng *rng) {
+
+    psudb::progress_update(0, "Insert Progress");
+    size_t reccnt = 0;
+    for (size_t i=start; i<stop; i++) {
+
+        extension->insert(records[i]);
+
+        if (gsl_rng_uniform(rng) <= delete_proportion && to_delete[delete_idx] <= i) {
+            ghost->insert(records[to_delete[delete_idx]]);
+            delete_idx++;
+            g_deleted_records++;
+        }
+
+    }
+
+}
+
+
 template<typename DE, typename R>
 static void insert_records(DE *structure, size_t start, size_t stop, 
                            std::vector<R> &records, std::vector<size_t> &to_delete, 
diff --git a/benchmarks/include/triespline_bsm.h b/benchmarks/include/triespline_bsm.h
index dfddc43..eaf1079 100644
--- a/benchmarks/include/triespline_bsm.h
+++ b/benchmarks/include/triespline_bsm.h
@@ -18,11 +18,19 @@ public:
 
 public:
     static BSMTrieSpline *build(std::vector<R> &records) {
+        if (records.size() == 0) {
+            return nullptr;
+        }
+
         std::sort(records.begin(), records.end());
         return new BSMTrieSpline(records);
     }
 
     static BSMTrieSpline *build_presorted(std::vector<R> &records) {
+        if (records.size() == 0) {
+            return nullptr;
+        }
+
         return new BSMTrieSpline(records);
     }
 
@@ -44,8 +52,9 @@ public:
         size_t idx = lower_bound(parms->lower_bound);
 
         size_t cnt = 0;
-        while (idx < m_data.size() && m_data[idx].first < parms->upper_bound) {
+        while (idx < m_data.size() && m_data[idx].first <= parms->upper_bound) {
             cnt++;
+            idx++;
         }
 
         rs.push_back({cnt, 0});
@@ -54,6 +63,10 @@ public:
     }
 
     std::vector<R> query_merge(std::vector<R> &rsa, std::vector<R> &rsb, void *parms) {
+        /* initialize rsa on the first merge */
+        if (rsa.size() == 0) {
+            rsa.push_back({0, 0});
+        }
         rsa[0].first += rsb[0].first;
         return std::move(rsa);
     }
@@ -75,6 +88,10 @@ private:
         m_min_key = m_data[0].first;
         m_max_key = m_data[m_data.size() - 1].first;
 
+        if (m_data.size() < 50) {
+            return;
+        }
+
         auto bldr = ts::Builder<K>(m_min_key, m_max_key, E);
         for (size_t i=0; i<m_data.size(); i++) {
             bldr.AddKey(m_data[i].first);
@@ -88,12 +105,14 @@ private:
     size_t lower_bound(K key) {
         if (m_data.size() == 0) {
             return 1;
-        } else if (m_data.size() == 1) {
-            if (m_data[0].first < key) {
-                return 1;
-            } else {
-                return 0;
+        } else if (m_data.size() < 50) {
+            for (size_t i=0; i<m_data.size(); i++) {
+                if (m_data[i].first >= key) {
+                    return i;
+                }
             }
+
+            return m_data.size();
         }
 
         auto bound = m_ts.GetSearchBound(key);
-- 
cgit v1.2.3