fix(approx_fns): use exact percentile when no compression

Author: aryan-212

datafusion/functions-aggregate-common/src/tdigest.rs

@@ -440,6 +440,10 @@ impl TDigest {
             return 0.0;
         }
 
+        // No compression happened since each centroid = one data point, so use exact percentile instead of interpolation
+        if self.count as usize == self.centroids.len() {
+            return self.exact_quantile(q);
+        }
         let rank = q * self.count as f64;
 
         let mut pos: usize;
@@ -515,6 +519,20 @@ impl TDigest {
         Self::clamp(value, min, max)
     }
 
+    fn exact_quantile(&self, q: f64) -> f64 {
+        if q <= 0.0 {
+            return self.min();
+        }
+        if q >= 1.0 {
+            return self.max();
+        }
+
+        let n = self.centroids.len();
+        let idx = (q * n as f64).ceil() as usize;
+        let idx = idx.saturating_sub(1).min(n - 1);
+        self.centroids[idx].mean()
+    }
+
     /// This method decomposes the [`TDigest`] and its [`Centroid`] instances
     /// into a series of primitive scalar values.
     ///