与std :: vector的奇怪行为

Question

我想实现一个矩阵模板类使用std :: vector。

代码：

// matrix.h 
template <class T> 
class matrix 
{ 
public: 
    ~matrix(void); 
    matrix(int rows, int cols):_rows(rows),_cols(cols){ _size = _rows*_cols;init();} // 

private: 
    matrix(void); 
    void init(); // sets up _matirx 

    //DATA 
    size_t _rows; 
    size_t _cols; 
    std::vector<T> _matrix; 
} 

// continued in matrix.tpp file 
template <class T> 
void matrix<T>::init(){ 
    _matrix = std::vector<T>(_rows*_cols); 
    for(size_t i = 1; i <= _rows ;i++){ 
     for(size_t j = 1 ; j <= _cols ; j++){ 
      _matrix[(i - 1)*_rows + (j - 1)] = 0 ; 
     } 
    } 
} 


template <class T> 
matrix<T>::matrix(const matrix<T>& rhs) 
{ 
    _matrix = rhs._matrix; 
    _rows = rhs._rows; 
    _cols = rhs._cols; 
} 



//in Source.cpp 

matrix<int> ABC = matrix<int>(4,2) ; 
// Gives Debug Assertion Failed , subscript error in VS 

matrix<int> ABC = matrix<int>(4000,4000) ;// Works , No Error 

matrix<int> ABC = matrix<int>(2,4) ; // Works No Error 

我知道如何使用的push_back，我会重新使用它实现类，但我想知道，为什么它在过去两年的情况下，在第一种情况下不呢？我的直觉是，在第一种情况下，一些元素没有被初始化。 std :: vector中是否存在对索引i的限制，在i + 2元素初始化之前，i + 1个元素是否必须初始化？ 还是有更微妙的事情呢？

• 谢谢

Answer 1

在这里，你可能有超出越界访问：

• 与4×2尺寸：

  _matrix = std::vector<T>(_rows*_cols); 
  for(size_t i = 1; i <= _rows ;i++){ 
      for(size_t j = 1 ; j <= _cols ; j++){ 
       _matrix[(i - 1)*_rows + (j - 1)] = 0 ; 
      } 
  } 
  

  最后的循环执行过程中把你的例子：i-1 = 3，_rows = 4和j-1 = 1，因此您正在访问大小为8的向量的第13个元素。

• 与4000x4000大小，您访问超过16000000（4000-1）* 4000 +（4000-1）= 15999999 th元素，因此没有出界限制访问。

最后一个例子也是如此。

Answer 2

简单的错字。

_matrix[(i - 1)*_rows + (j - 1)] = 0 ; 

应该

_matrix[(i - 1)*_cols + (j - 1)] = 0; 

通过工作在纸上的循环的每次迭代会发现这一点。