eigen 矩阵块操作

eigen 矩阵块操作

在Eigen中最基本的快操作运算是用.block()完成的。提取的子矩阵同样分为动态大小和固定大小。

块操作	构建动态大小子矩阵
提取块大小为(p,q),起始于(i,j)	matrix.block(i,j,p,q)

同样需要注意的是在Eigen中，索引是从0开始。所有的操作方法都可以适用于Array.同样使用固定大小的操作方式在小型矩阵运算时更加的快，但要求在编译时就要知道矩阵的大小。

        Matrix16x16 E = Matrix16x16::Identity();
        Matrix16x16 Q = Matrix16x16::Identity();
        
        for (int i = 0; i < A.Rows; i++)
        {
            for (int j = 0; j < A.Rows; i++)
            {
                if (i==j)
                {
                    Q.block(i*4, j*4, 4, 4) = A(i, j)*E - B.matrix.transpose();
                }
                else
                {
                    Q.block(i*4, j*4, 4, 4) = A(i, j)*E;
                }
                
            }
        }

.block() 可以为左值，也可以为右值。

 

参考：

(5条消息) Eigen子矩阵操作(block)_不二青衣的博客-CSDN博客