本征增量列一

Question

如何将动态矩阵的列增加一，作为就地操作（不创建 copies/intermediates）？

尝试：

#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;
}

Answer 1

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

这类似于

#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;
}

本征增量列一

Eigen Increment Column by One

c++

eigen

eigen3