1

我正在尝试使用 nth_element 找到 (x,y) 点向量的中位数

        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

在我重新编译之前，结果是相当可重复的，此时细节发生了变化，但一般问题仍然存在。显然，向量以某种方式被破坏了，但我想不出如何。

此外，如果我使用 sort 而不是 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
....
4

1 回答 1

1

您很可能遇到了 libstdc++ 中最近的一个错误，它破坏了 nth_element 函数。这已经得到修复，但一些 Linux 版本附带了损坏的版本（例如 Ubuntu 13.10）

可以在 GCC 跟踪器上找到该错误的补丁和讨论：http ://gcc.gnu.org/bugzilla/show_bug.cgi?id=58800

于 2014-02-04T15:15:09.403 回答