Eigen :: VectorXdにアクセスする際に2番目のインデックスとしてゼロを使用するのは安全ですか？

Question

Eigen::VectorXdはScalar operator()(Index i)を持ち、ベクトルのインデックスiの係数を返します。しかし、Eigen::VectorXdはEigen::Matrixの特殊タイプ、すなわちタイプEigen::Matrix<Scalar, Eigen::Dynamic, 1>;であるので、Scalar operator()(Index i, Index j)もある。

質問：

は私がゼロにjを設定した場合、第2バージョンを使用する（すなわちなし未定義の動作）安全であることを前提とすることはできますか？言い換えれば、以下のコードはOKですか？

Eigen::VectorXd v(4); 
v << 1, 2, 3, 4; 
std::cout << v(2, 0); // displays 3 

それはOKだように見え、そこにすべての警告をデバッグモードでコンパイルされたときに何も失敗したアサーションや警告はありませんが、私は100％を確認していません。

Answer 1

はv限り安全である例えば、両方の列と行ベクトルに対してv(i)作品を使用し、一方、列ベクトルである：

template<typename T> 
void foo(const T &v) { 
    v(2); // OK 
    v(2,0); // -> out of bounds runtime assertion 
} 
MatrixXd mat(10,10); 
foo(mat.row(5)); 

Answer 2

Iは@ggaels answer時に言えます。あなたはDenseCoeffsBase.hでoperator()定義（私は3.2.10を引用しています）を見れば、あなたは彼らが両方PlainObjectBase.hでcoeffRefの定義を見てみるとcoeff（またはcoeffRef）

EIGEN_STRONG_INLINE CoeffReturnType operator()(Index row, Index col) const 
{ 
    eigen_assert(row >= 0 && row < rows() 
     && col >= 0 && col < cols()); 
    return derived().coeff(row, col); 
} 

EIGEN_STRONG_INLINE CoeffReturnType 
operator()(Index index) const 
{ 
    eigen_assert(index >= 0 && index < size()); 
    return derived().coeff(index); 
} 

を呼び出すことがわかります

EIGEN_STRONG_INLINE Scalar& coeffRef(Index rowId, Index colId) 
{ 
    if(Flags & RowMajorBit) 
    return m_storage.data()[colId + rowId * m_storage.cols()]; 
    else // column-major 
    return m_storage.data()[rowId + colId * m_storage.rows()]; 
} 

EIGEN_STRONG_INLINE Scalar& coeffRef(Index index) 
{ 
    return m_storage.data()[index]; 
} 

ので、行ベクトルの場合には、あなたが境界エラーのうちの可能なアサーションの失敗を/避けるためにv(0,2)を記述する必要があります：私たちは、オフセットを簡単に計算されていることがわかります。