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 + 1元素必须在i + 2元素初始化之前初始化？还是有更微妙的事情发生？

谢谢

Answer 1

简单的拼写错误。

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

应该是

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

对于循环的每次迭代，在纸上进行处理都会揭示这一点。

Answer 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 ; 
    }   
}

在上一次循环执行期间采取您的示例：

4x2大小：i-1 = 3，_rows = 4且j-1 = 1，因此您可以访问大小为8的向量的第13个元素。
4000x4000大小，您可以访问超过16000000的（4000-1）* 4000 +（4000-1）= 15999999元素，因此不会超出界限。

最后一个例子也是如此。

std :: vector的奇怪行为

2 个答案: