Question

我是Eigen的新手，我想创建10个mxn矩阵。由于某些原因，我使用以下方法：

Matrix<double, m, n*10>

这似乎内存分配类似于以下内容：

 _______________________________________________________________
|M1(1,1)|M2(1,1)|...|M10(1,1)|.....|M1(1,n)|M2(1,n)|...|M10(1,n)|
|  . .                                                          |
|  .                                                            |

现在如何才能为这10个矩阵中的每一个创建一个参考矩阵（通过引用表示而不复制数据）？

Answer 1

我建议使用动态分配的矩阵，因为m和n可能很大。此外，当默认为列主要时，您似乎假设矩阵存储器是行主要的。在下面的例子中，我明确地将它们排成了主要部分。

您可以像这样使用系列Eigen::Map<MatrixXd>：

#include <Eigen/Core>
#include <iostream>

using namespace Eigen;

int main(void)
{
    int m = 3;
    int n = 4;
    int x = 6;
    typedef Matrix < double, Dynamic, Dynamic, RowMajor > ourMat;
    ourMat  M1(m, n * x);
    M1.setConstant(9.9);

    for (int i = 0; i < x; i++)
    {
        Eigen::Map<ourMat, 0, Stride<Dynamic, Dynamic>> m_i(M1.data() + i,
                                                            m, n,
                                   Stride<Dynamic, Dynamic>(n*x,x));
        m_i.setConstant(double(i));
        std::cout << m_i << std::endl;
        std::cout << M1 << "\n" << std::endl;

    }

    Eigen::Map<VectorXd> m_i(M1.data(), m * n * x);
    std::cout << m_i.transpose() << std::endl;

    return 0;
}

子特征在特征库中大步前进

1 个答案: