12345678910111213141516171819202122232425262728293031323334353637383940414243444546 |
- #include "median.h"
- #include "matrix_to_list.h"
- #include <vector>
- #include <algorithm>
- template <typename DerivedV, typename mType>
- IGL_INLINE bool igl::median(
- const Eigen::MatrixBase<DerivedV> & V, mType & m)
- {
- using namespace std;
- if(V.size() == 0)
- {
- return false;
- }
- vector<typename DerivedV::Scalar> 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;
- }
- #ifdef IGL_STATIC_LIBRARY
- template bool igl::median<Eigen::Block<Eigen::Matrix<float, -1, -1, 0, -1, -1>, -1, 1, true>, float>(Eigen::MatrixBase<Eigen::Block<Eigen::Matrix<float, -1, -1, 0, -1, -1>, -1, 1, true> > const&, float&);
- template bool igl::median<Eigen::Block<Eigen::Matrix<double, -1, -1, 0, -1, -1>, -1, 1, true>, double>(Eigen::MatrixBase<Eigen::Block<Eigen::Matrix<double, -1, -1, 0, -1, -1>, -1, 1, true> > const&, double&);
- template bool igl::median<Eigen::Matrix<double, -1, 1, 0, -1, 1>, double>(Eigen::MatrixBase<Eigen::Matrix<double, -1, 1, 0, -1, 1> > const&, double&);
- #endif
|