访问 boost稀疏矩阵的元素似乎使程序停止运行

Question

我有一个奇怪的错误，我希望更有经验的程序员能够对此有所了解。我正在使用 boost ublas 稀疏矩阵，特别是mapped_matrix，并且最终会出现间歇性错误，但不会在程序的初始阶段出现。这是一个很大的程序，所以我不能发布所有代码，但核心思想是我调用属于特定类的函数：

bool MyClass::get_cell(unsigned int i, unsigned int j) const
{
    return c(i,j);
}

变量 c 被定义为该类的成员

boost::numeric::ublas::mapped_matrix<bool> c;

当错误发生时，程序似乎停止（但不会崩溃）。使用 Eclipse 进行调试，我可以看到程序进入了 boost Matrix 代码，并继续向下进入 std::map、std::_Rb_tree 和 std::less 几个级别。此外，程序有时会追溯到 std::map、std::_Rb_tree 和 std::_Select1st。当代码正在执行并且 _Rb_tree 中内存中的活动行发生变化时，执行似乎永远不会返回到 std::map 级别。程序卡在 std::map 中的行是以下函数的返回。

const_iterator
find(const key_type& __x) const
{ return _M_t.find(__x); }

在我看来，程序正在寻找 c 矩阵中的某些元素，但不知何故底层存储机制“放错了它”。但是，我不确定为什么或如何解决它。这也可能是完全错误的。

您能提供的任何帮助将不胜感激。如果我在这个问题中没有包含正确的信息，请告诉我我遗漏了什么。谢谢。

Answer 1

尝试调试代码的一些事情（不一定是永久性更改）：

• 将 c 的矩阵类型中的 bool 更改为 int，以查看矩阵是否需要数字类型。
• 将矩阵类型更改为具有类似接口的另一种矩阵类型，可能是普通的旧矩阵。
• Valgrind 该应用程序（如果您使用的是 Linux）来检查您没有损坏内存。

如果失败，您可以尝试在每次修改矩阵时调用 get_cell 以查看可能导致问题的原因。

如果做不到这一点，您可能必须尝试将问题减少到更小的代码子集，您可以将其发布在此处。

如果您告诉我们您正在使用什么编译器和操作系统，这可能会有所帮助。

Answer 2

这是多线程程序的一部分吗？

我之所以这么问，是因为通常当我在 STL 中看到问题时，最终都会发现是访问不同步的问题。