将向量的值传递给特征库格式[重复]

Question

我正在尝试使用 Eigen 对 A 作为二维平方向量的能力来求解线性方程 Ax=b。我将 A 和 b 分别作为基于 C++ 的 2D 向量和 1D 向量。但是，我找不到将它们的值传递给特征格式矩阵和向量的方法。你能告诉我如何以 Eigen 格式复制变量吗？ 此外，一开始应该包括什么才能使用 Map 类作为可能的溶剂？！ 代码如下：

#include <iostream>
#include <vector>

#include "Eigen/Dense"

using namespace std;
using namespace Eigen;
int main()
{
    // Let's consider that the A and b are following CPP based vectors:
    vector<vector<double>> mainA= { { 10.,11.,12. },{ 13.,14.,15. },{ 16.,17.,18. } };  
    vector<double> mainB = { 2.,5.,8. };

    // ??? Here I need to do something to pass the values to the following Eigen 
    //format matrix and vector

    MatrixXf A;
    VectorXf b;

    cout << "Here is the matrix A:\n" << A << endl;
    cout << "Here is the vector b:\n" << b << endl;
    Vector3f x = A.colPivHouseholderQr().solve(b);
    cout << "The solution is:\n" << x << endl;

}

Answer 1

正如 cmets 中提到的，Eigen::Map 应该可以解决问题。 通常不使用Unaligned 就可以逃脱，但为了正确/稳定，最好使用它：

auto A = Eigen::Map<Eigen::MatrixXd, Eigen::Unaligned>(mainA.data(), mainA.size(), mainA[0].size())
auto b = Eigen::Map<Eigen::VectorXd, Eigen::Unaligned>(mainB.data(), mainB.size());

Vector3d x = A.colPivHouseholderQr().solve(b);

要回答以下关于健壮性的问题：这可以使用辅助函数来完成：

template <typename T, int Align>
Eigen::Map<Eigen::Matrix<T, -1, -1>, Align> CreateMatrix(const std::vector<std::vector<T>>& x)
{
    int R = x.size();
    assert(!x.empty());
    int C = x[0].size();
    #ifndef NDEBUG
    for(int r=1; r < R; r++)
        assert(x[r].size() == x[0].size());
    #endif
    return auto A = Eigen::Map<Eigen::Matrix<T,-1,-1>, Align>(x.data(), R, C);
}

虽然它看起来很冗长，但它会进行大量的健全性检查，然后您可以在进行单元测试的所有代码中依赖这些检查。