如何处理本征(C++)中的复数？

Question

我正在成功地与艾根合作，我试图用复数来理解一些细节。

1. 如何将矩阵或向量乘以复常数？乘成复矩阵很容易；同样地，与复向量和矩阵的内积一样。
2. 如何将复值赋给矩阵元素？我试过：

这很好()

kx.real()(0, 0) = 1.0;

这会引发编译器错误。

kz_r.imag()(0, ii) =1.0

我得到的错误是：

严重程度代码描述项目文件行抑制状态错误C2440‘返回’：不能从'double‘转换为'double &’\eigen\src\mathfunctions.h.h 919

Answer 1

只需使用std::complex<double> (或float)：

std::complex<double> c(1,1);
Eigen::MatrixXd R1; R1.setRandom(2,2);
Eigen::MatrixXcd C1 = c*R1; // multiply complex*real
Eigen::MatrixXcd C2 = c*C1; // complex scalar times complex matrix
C1(0,0) = c; // assign complex value.