1234567891011121314151617181920212223242526272829303132333435 |
- #include "median.h"
- #include "matrix_to_list.h"
- #include <vector>
- #include <algorithm>
- IGL_INLINE bool igl::median(const Eigen::VectorXd & V, double & m)
- {
- using namespace std;
- if(V.size() == 0)
- {
- return false;
- }
- vector<double> vV;
- matrix_to_list(V,vV);
-
- size_t n = vV.size()/2;
- nth_element(vV.begin(),vV.begin()+n,vV.end());
- if(vV.size()%2==0)
- {
- nth_element(vV.begin(),vV.begin()+n-1,vV.end());
- m = 0.5*(vV[n]+vV[n-1]);
- }else
- {
- m = vV[n];
- }
- return true;
- }
|