我需要帮助解决一个奇怪的运行时错误。这是给出它的代码:
JacobianCol &diag_index_column = J[diag_index];
JacobianColData::iterator &diagonal_element = diag_index_column.find(diag_index);
Jacobian J2 = J; //added to reveal the problem
J[diag_index].divide_by(diagonal_element);
我想要的是。我想保存迭代器对角元素并将其传递给divide_by函数。但是当我调用J变量时,迭代器会下降。指向内存的指针仍然存在(我已在调试器中检查过),但迭代器的内容已损坏(未引用的变量)。
我究竟做错了什么?
更多代码:
雅可比 J:
class Jacobian
{
private:
std::vector<JacobianCol> _J;
...
public:
...
JacobianCol& operator [](int i); //return _J[i];
};
雅可比列:
typedef std::map<int, Submatrix> JacobianColData;
class JacobianCol
{
private:
...
JacobianColData _col_data;
public:
...
JacobianColData::iterator &find(int key, bool need_append = false);
};
找到实现:
JacobianColData::iterator &JacobianCol::find(int key, bool need_append)
{
if(need_append)
this->insert(key);
JacobianColData::iterator &res = this->_col_data.find(key);
return res;
}