Changeset 1039 for trunk/yat/statistics

Timestamp:

Feb 6, 2008, 1:20:31 AM (16 years ago)

Author:

Peter

Message:

fixes #276

File:

• trunk/yat/statistics/utility.h (modified) (1 diff)

trunk/yat/statistics/utility.h

@@ -181 +181 @@
     std::copy(first, last, std::back_inserter(v_copy));
     double j = p/100 * (v_copy.size()-1);
-    int i = static_cast<int>(j);
-    std::partial_sort(v_copy.begin(),v_copy.begin()+i+2 , v_copy.end());
+    size_t i = static_cast<size_t>(j);
+    if (i+2 < v_copy.size()) {
+      utility::yat_assert<std::out_of_range>(i+2 < v_copy.size(), 
+                                             "in utility::percentile");
+      std::partial_sort(v_copy.begin(), v_copy.begin()+i+2, v_copy.end());
+    }
+    else
+      std::sort(v_copy.begin(), v_copy.end());
+      
     return percentile(v_copy.begin(), v_copy.end(), p, true);
   }