特征增量列一个

时间:2016-10-03 15:40:54

标签: c++ eigen eigen3

如何将动态矩阵的列增加1,作为就地操作(不创建副本/中间体)?

尝试:

#include <Eigen/Dense>
#include <iostream>
#include <stdint.h>
int main(void){
    Eigen::MatrixXf A;
    A = Eigen::MatrixXf::Random(3, 5);
    std::cout << A << std::endl << std::endl;
    A.col(1) = A.col(1)*2; //this works.
    A.col(1) = A.col(1) + 1; //this doesn't work.
    std::cout << A << std::endl;
}

1 个答案:

答案 0 :(得分:2)

我找到了一种方法来做到这一点。但我不知道这项行动是否到位。

这类似于eigen: Subtracting a scalar from a vector

#include <Eigen/Dense>
#include <iostream>
int main(void){
    Eigen::MatrixXf A;
    A = Eigen::MatrixXf::Random(3, 5);
    std::cout << A << std::endl << std::endl;

    A.col(1) = A.col(1)*2;
    A.col(1) = A.col(1) + Eigen::VectorXf::Ones(3);
    std::cout << A << std::endl;
}

另一种方法是使用数组操作。这种方式似乎更好(我猜)。

https://eigen.tuxfamily.org/dox/group__TutorialArrayClass.html

#include <Eigen/Dense>
#include <iostream>
int main(void){
    Eigen::MatrixXf A;
    A = Eigen::MatrixXf::Random(3, 5);
    std::cout << A << std::endl << std::endl;

    A.array() += 1;
    A.col(1).array() += 100;

    std::cout << A << std::endl;
}