使用 Eigen 时出错：无效使用不完整类型 'const class Eigen

Question

我有以下有效的代码：

Matrix <float, ny+1, nx> eXX;
eXX.setZero();

Eigen::Matrix< double, (ny+1), (ny)> u; 
    u.setZero(); 
    for(int i = 0; i< nx; i++){
        for(int j = 0; j< ny+1; j++){
            eXX(j + (ny+1)*i) = (i)*2*EIGEN_PI/nx; 
            u(j + (ny+1)*i) = cos(eXX(j + (ny+1)*i));
        }
    }

但是当我编写以下内容时，它不起作用：

Matrix <float, ny+1, nx> eXX;
eXX.setZero();

Eigen::Matrix< double, (ny+1), (ny)> u; 
    u.setZero(); 
    for(int i = 0; i< nx; i++){
        for(int j = 0; j< ny+1; j++){
            eXX(j + (ny+1)*i) = (i)*2*EIGEN_PI/nx; 
        }
    }
u = eXX.matrix().cos();// -or- std::cos(eXX.array());
    std::cout << u << "\n"; //error

完整的错误信息：

Test.cpp:418:23: error: invalid use of incomplete type ‘const class Eigen::MatrixFunctionReturnValue<Eigen::Matrix<float, 11, 10> >’
  418 |  u = eXX.matrix().cos();
      |                       ^
In file included from /mnt/c/Users/eigen-3.4.0/eigen-3.4.0/Eigen/Core:163,
                 from /mnt/c/Users/eigen-3.4.0/eigen-3.4.0/Eigen/Dense:1,
                 from Test.cpp:21:
/mnt/c/Users/eigen-3.4.0/eigen-3.4.0/Eigen/src/Core/util/ForwardDeclarations.h:305:34: note: declaration of ‘class Eigen::MatrixFunctionReturnValue<Eigen::Matrix<float, 11, 10> >’
  305 | template<typename Derived> class MatrixFunctionReturnValue;

我想我可以尝试在不使用 for 循环的情况下重写eXX并通过它，但这也行不通。 另外，我读到有人建议添加类似#include <MatrixFunctionReturnValue>的内容，这实际上使事情变得更糟。 谢谢。

我在这里也添加了我的包含：

#define _USE_MATH_DEFINES 
#include <cmath>
#include<math.h>
#include<stdio.h>
#include "fftw3.h"
#include <cstring>
#include <sstream>
#include <string>
#include <sys/resource.h>
#include <iostream> 
#include <vector>
#include <fstream> 
#include <iomanip>
#include <numeric>
#include <assert.h>
#include <Eigen/Dense>
#include <unsupported/Eigen/FFT>
#include <Eigen/SparseCore> 
#include <Eigen/Sparse>

Answer 1

Matrix类是为线性代数构建的。 当您想要对矩阵的元素进行操作时，您需要使用Array类。 请参阅 有关 Array 的 Eigen 文档。 另一种方法是使用unaryExpr将矩阵的每个元素作为输入。

以下是这两种方法：

#include <iostream>
#include <Eigen/Dense>
....
Eigen::Matrix<double, 3, 3> vals;
vals.setZero();
std::cout << vals << '\n';
std::cout << "******\n";
std::cout << vals.array().cos() << '\n';
std::cout << vals << '\n';
std::cout << "******\n";
Eigen::Matrix<double, 3, 3> res = vals.unaryExpr([](double a) { return cos(a); });
std::cout << res << '\n';
std::cout << vals << '\n';

请注意vals如何随各种操作发生变化（并且不会发生变化）。