Я пытаюсь найти медиану векторов (x, y) точек, используя nth_element
cv::Point2f medOffset;
vector<float> tempOffsetsX = offsetsX;
int medLoc = tempOffsetsX.size()/2;
nth_element(tempOffsetsX.begin(), tempOffsetsX.begin()+medLoc, tempOffsetsX.end());
// sort(tempOffsetsX.begin(), tempOffsetsX.end());
medOffset.x = tempOffsetsX[medLoc];
vector<float> tempOffsetsY = offsetsY;
****** debug out line 1 *********
nth_element(tempOffsetsY.begin(), tempOffsetsY.begin()+medLoc, tempOffsetsY.end());
// sort(tempOffsetsY.begin(), tempOffsetsY.end());
medOffset.y = tempOffsetsY[medLoc];
****** debug out line 2 *********
tempOffsetsX работает просто отлично, но иногда tempOffsetsY дает очень странные результаты после nth_element. Вот пример вывода на отмеченные строки отладки
tempOffsetsY1: 5.184135 -1.564125 3.751759 0.221855 -0.742348 1.737648
tempOffsetsY2: -0.742348 -1.564125 -8885092352.000000 -8850636800.000000 0.000000 0.000000
Результаты повторяются до тех пор, пока я не перекомпилирую, после чего специфика изменится, и общая проблема останется. Ясно, что вектор как-то искажается, но я не могу придумать, как это сделать.
Кроме того, если я использую сортировку вместо nth_element, это работает без проблем. Для отладки я попытался выполнить сортировку, а затем nth_element, которая работала просто отлично. Так или иначе переупорядочение, которое происходит внутри nth_element, испортилось, но я не могу придумать, как это сделать.
Есть идеи, как это происходит?
редактировать — Больше информации о моем окружении. Я использую Arch Linux. Я только что сделал обновление системы. Я должен отметить, что этот же код работал без проблем до обновления, и я впервые запускаю его после обновления. Но это разрыв в несколько дней, и я всегда не решаюсь указывать на системные библиотеки, что обычно является моей собственной проблемой.
[]$ uname -r
3.11.6-1-ARCH
[]$ gcc -v
Using built-in specs.
COLLECT_GCC=gcc
COLLECT_LTO_WRAPPER=/usr/lib/gcc/x86_64-unknown-linux-gnu/4.8.2/lto-wrapper
Target: x86_64-unknown-linux-gnu
Configured with: /build/gcc/src/gcc-4.8.2/configure --prefix=/usr --libdir=/usr/lib --libexecdir=/usr/lib --mandir=/usr/share/man --infodir=/usr/share/info --with-bugurl=https://bugs.archlinux.org/ --enable-languages=c,c++,ada,fortran,go,lto,objc,obj-c++ --enable-shared --enable-threads=posix --with-system-zlib --enable-__cxa_atexit --disable-libunwind-exceptions --enable-clocale=gnu --disable-libstdcxx-pch --enable-gnu-unique-object --enable-linker-build-id --enable-cloog-backend=isl --disable-cloog-version-check --enable-lto --enable-gold --enable-ld=default --enable-plugin --with-plugin-ld=ld.gold --with-linker-hash-style=gnu --disable-install-libiberty --disable-multilib --disable-libssp --disable-werror --enable-checking=release
Thread model: posix
gcc version 4.8.2 (GCC)
[]$ pacman -Qi glibc
Name : glibc
Version : 2.18-8
....
Вполне вероятно, что вы столкнулись с недавней ошибкой в libstdc ++, которая сломала функцию nth_element. С тех пор это было исправлено, но некоторые выпуски Linux поставлялись с неработающей версией (например, Ubuntu 13.10)
Патч и обсуждение ошибки можно найти на трекере GCC: http://gcc.gnu.org/bugzilla/show_bug.cgi?id=58800
Других решений пока нет …