From 9269777672f336e02aeb7da2197abb00624ea2bd Mon Sep 17 00:00:00 2001
From: Stan Jenkins <mc.stanislaw@gmail.com>
Date: Sun, 29 Sep 2024 20:34:19 -0400
Subject: [PATCH 1/3] [Jenkins] auto-formatting by clang-format version
 10.0.0-4ubuntu1

---
 src/stan/mcmc/chainset.hpp | 4 +---
 1 file changed, 1 insertion(+), 3 deletions(-)

diff --git a/src/stan/mcmc/chainset.hpp b/src/stan/mcmc/chainset.hpp
index 44dc06d053..d9a019685c 100644
--- a/src/stan/mcmc/chainset.hpp
+++ b/src/stan/mcmc/chainset.hpp
@@ -249,9 +249,7 @@ class chainset {
    * @param index parameter index
    * @return median
    */
-  double median(const int index) const {
-    return quantile(index, 0.5);
-  }
+  double median(const int index) const { return quantile(index, 0.5); }
 
   /**
    * Compute median value of specified parameter across all chains.

From 197542ff7632e5a12b99aada8719be0785a348f6 Mon Sep 17 00:00:00 2001
From: Mitzi Morris <mitzi212@gmail.com>
Date: Mon, 30 Sep 2024 13:38:18 -0400
Subject: [PATCH 2/3] using nth_element

---
 src/stan/mcmc/chainset.hpp | 29 +++++++++++++----------------
 1 file changed, 13 insertions(+), 16 deletions(-)

diff --git a/src/stan/mcmc/chainset.hpp b/src/stan/mcmc/chainset.hpp
index d9a019685c..d7d4a463e7 100644
--- a/src/stan/mcmc/chainset.hpp
+++ b/src/stan/mcmc/chainset.hpp
@@ -9,8 +9,7 @@
 #include <boost/accumulators/accumulators.hpp>
 #include <boost/accumulators/statistics/stats.hpp>
 #include <boost/accumulators/statistics/mean.hpp>
-#include <boost/accumulators/statistics/tail_quantile.hpp>
-#include <boost/accumulators/statistics/p_square_quantile.hpp>
+#include <boost/accumulators/statistics/median.hpp>
 #include <boost/accumulators/statistics/variance.hpp>
 #include <boost/accumulators/statistics/covariance.hpp>
 #include <boost/accumulators/statistics/variates/covariate.hpp>
@@ -249,8 +248,10 @@ class chainset {
    * @param index parameter index
    * @return median
    */
-  double median(const int index) const { return quantile(index, 0.5); }
-
+  double median(const int index) const {
+    return (quantile(index, 0.5));
+  }
+  
   /**
    * Compute median value of specified parameter across all chains.
    *
@@ -275,9 +276,9 @@ class chainset {
     auto center = median(index);
     Eigen::MatrixXd abs_dev = (draws.array() - center).abs();
     size_t idx = static_cast<size_t>(0.5 * (abs_dev.size() - 1));
-    std::vector<double> sorted(abs_dev.data(), abs_dev.data() + abs_dev.size());
-    std::nth_element(sorted.begin(), sorted.begin() + idx, sorted.end());
-    return 1.4826 * sorted[idx];
+    std::vector<double> copy(abs_dev.data(), abs_dev.data() + abs_dev.size());
+    std::nth_element(copy.begin(), copy.begin() + idx, copy.end());
+    return 1.4826 * copy[idx];
   }
 
   /**
@@ -311,10 +312,10 @@ class chainset {
       throw std::out_of_range("Probability must be between 0 and 1.");
     }
     Eigen::MatrixXd draws = samples(index);
-    size_t idx = static_cast<size_t>(prob * (draws.size() - 1));
-    std::vector<double> sorted(draws.data(), draws.data() + draws.size());
-    std::nth_element(sorted.begin(), sorted.begin() + idx, sorted.end());
-    return sorted[idx];
+    std::vector<double> copy(draws.data(), draws.data() + draws.size());
+    int idx = static_cast<int>(prob * (copy.size() - 1));
+    std::nth_element(copy.begin(), copy.begin() + idx, copy.end());
+    return copy[idx];
   }
 
   /**
@@ -349,12 +350,8 @@ class chainset {
     if (probs.minCoeff() < 0.0 || probs.maxCoeff() > 1.0) {
       throw std::out_of_range("Probabilities must be between 0 and 1.");
     }
-    Eigen::MatrixXd draws = samples(index);
-    std::vector<double> sorted(draws.data(), draws.data() + draws.size());
-    std::sort(sorted.begin(), sorted.end());
     for (size_t i = 0; i < probs.size(); ++i) {
-      size_t idx = static_cast<size_t>(probs[i] * (sorted.size() - 1));
-      quantiles[i] = sorted[idx];
+      quantiles[i] = quantile(index, probs[i]);
     }
     return quantiles;
   }

From 021ac040c0d18dd97c5311b033b5e58c21502731 Mon Sep 17 00:00:00 2001
From: Stan Jenkins <mc.stanislaw@gmail.com>
Date: Mon, 30 Sep 2024 13:39:03 -0400
Subject: [PATCH 3/3] [Jenkins] auto-formatting by clang-format version
 10.0.0-4ubuntu1

---
 src/stan/mcmc/chainset.hpp | 6 ++----
 1 file changed, 2 insertions(+), 4 deletions(-)

diff --git a/src/stan/mcmc/chainset.hpp b/src/stan/mcmc/chainset.hpp
index d7d4a463e7..1d17f1b395 100644
--- a/src/stan/mcmc/chainset.hpp
+++ b/src/stan/mcmc/chainset.hpp
@@ -248,10 +248,8 @@ class chainset {
    * @param index parameter index
    * @return median
    */
-  double median(const int index) const {
-    return (quantile(index, 0.5));
-  }
-  
+  double median(const int index) const { return (quantile(index, 0.5)); }
+
   /**
    * Compute median value of specified parameter across all chains.
    *