Function MedianAbsoluteDeviation

src/colmap/math/math.h:242–249

template <typename T>
std::pair<double, double> MedianAbsoluteDeviation(std::vector<T>& elems) {
  const double median = Median(elems);
  std::vector<double> abs_deviations(elems.size());
  for (size_t i = 0; i < elems.size(); i++) {
    abs_deviations[i] = std::abs(static_cast<double>(elems[i]) - median);
  }
  return {median, Median(abs_deviations)};
}

template <typename T>
std::pair<double, double> MedianAbsoluteDeviation(std::vector<T>&& elems) {