Real ErrorVector::median() { const unsigned int n = this->size(); if (n == 0) return 0.; // Build a StatisticsVector<ErrorVectorReal> containing // only our active entries and take its mean StatisticsVector<ErrorVectorReal> sv; sv.reserve (n); for (unsigned int i=0; i<n; i++) if(this->is_active_elem(i)) sv.push_back((*this)[i]); return sv.median(); }
Real StatisticsVector<T>::median() const { StatisticsVector<T> sv = (*this); return sv.median(); }